Related papers: Rigorous numerical integration of algebraic functi…
In this paper, we consider spectral-collocation method base on Legendre-Gauss-Lobatto point. We present a computational method for solving a class of fractional integral equation of the second kind. Then based on Legendre-Gauss-Lobatto…
The Lambert W function has utility for solving various exponential and logarithmic equations arranged in the form of $g(x)e^{g(x)}$. Using the Lambert W function and tetration, a variety of categorized inversion formulas are presented.…
We evaluate the performance of novel numerical methods for solving one-dimensional nonlinear fractional dispersive and dissipative evolution equations. The methods are based on affine combinations of time-splitting integrators and…
Low rank matrix approximation is a popular topic in machine learning. In this paper, we propose a new algorithm for this topic by minimizing the least-squares estimation over the Riemannian manifold of fixed-rank matrices. The algorithm is…
Computer Algebra systems are widely spread because of some of their remarkable features such as their ease of use and performance. Nonetheless, this focus on performance sometimes leads to unwanted consequences: algorithms and computations…
An existing solvability result for relaxed one-sided Lipschitz algebraic inclusions is substantially improved. This enhanced solvability result allows the design of a very robust numerical method for the approximation of a solution of the…
In this paper, in order to improve the spatial accuracy, the exponential integrator Fourier Galerkin method (EIFG) is proposed for solving semilinear parabolic equations in rectangular domains. In this proposed method, the spatial…
We present here algorithms for efficient computation of linear algebra problems over finite fields.
In 1970s, a method was developed for integration of nonlinear equations by means of algebraic geometry. Starting from a Lax representation with spectral parameter, the algebro-geometric method allows to solve the system explicitly in terms…
Motivated by many applications in complex domains with boundaries exposed to large topological changes or deformations, fictitious domain methods regard the actual domain of interest as being embedded in a fixed Cartesian background. This…
We represent the Riemann zeta function in the half-plane $\Re s >1$ via series whose terms admit geometrically decreasing bounds. Due to an underlying recurrence relation, which is used to compute coefficients entering into the terms, the…
Recently, Trefethen (SIAM Review 50 (2008), 67--87) and Xiang and Bornemann (SIAM J. Numer. Anal. 50 (2012), 2581--2587) investigated error bounds for n-point Gauss and Clenshaw-Curtis quadrature for the Legendre weight with integrands…
The aim of this work is to certify lower bounds for real-valued multivariate functions, defined by semialgebraic or transcendental expressions. The certificate must be, eventually, formally provable in a proof system such as Coq. The…
Given a model in algebraic statistics and some data, the likelihood function is a rational function on a projective variety. Algebraic algorithms are presented for computing all critical points of this function, with the aim of identifying…
Algorithms for numerical computation of symmetric elliptic integrals of all three kinds are improved in several ways and extended to complex values of the variables (with some restrictions in the case of the integral of the third kind).…
Finite-part integration is a recently introduced method of evaluating convergent integrals by means of the finite part of divergent integrals [E.A. Galapon, {\it Proc. R. Soc. A 473, 20160567} (2017)]. Current application of the method…
We present an efficient implementation of hypergeometric functions in arbitrary-precision interval arithmetic. The functions ${}_0F_1$, ${}_1F_1$, ${}_2F_1$ and ${}_2F_0$ (or the Kummer $U$-function) are supported for unrestricted complex…
The generalised Gegenbauer functions of fractional degree (GGF-Fs), denoted by ${}^{r\!}G^{(\lambda)}_\nu(x)$ (right GGF-Fs) and ${}^{l}G^{(\lambda)}_\nu(x)$ (left GGF-Fs) with $x\in (-1,1),$ $\lambda>-1/2$ and real $\nu\ge 0,$ are special…
Nonlinear elliptic problems arise in many fields, including plasma physics, astrophysics, and optimal transport. In this article, we propose a novel operator-splitting/finite element method for solving such problems. We begin by introducing…
In this paper, we propose a fast and accurate numerical method based on Fourier transform to solve Kolmogorov forward equations of symmetric scalar L\'evy processes. The method is based on the accurate numerical formulas for Fourier…