14 papers · ranked by Valyu relevance
François Alouges, Giovanni Di Fratta, Alberto Fiorenza, Renato Fiorenza
We present a unified fixed-point construction of the elementary transcendental functions, encompassing the real exponential, the complex exponential (sine and cosine), and the natural logarithm. Each function is characterized as the unique solution of a duplication identity established through the Banach contraction…
Dario A. Bini, Massimiliano Fasi, Bruno Iannazzo
We consider the numerical evaluation of the quantity Af(A − 1B), where A is Hermitian positive definite, B is Hermitian, and f is a function defined on the spectrum of A − 1B. We study the conditioning of the problem, and we introduce several algorithms that combine the Schur decomposition with either the matrix square…
Bowen Gao, Daniel Kressner, Meiyue Shao
Mixed precision algorithms can significantly enhance the performance of linear algebra solvers by leveraging increasingly powerful low precision hardware while recovering working precision accuracy through, for example, iterative refinement. In this paper, we propose a novel mixed precision algorithm for computing…
Flavio Barbosa, Fernando Nogueira
This article proposes new perspectives for developing derivative-based numerical algorithms, supported by the introduction of Generalized Derivative operators. It demonstrates that these operators have the potential to enhance and extend existing derivative-based numerical methods. To this end, two iterative…
Chadi Nour
This paper extends the numerical method introduced by de Pinho et al. [22] and later developed in Nour and Zeidan [33, 37] to the nonautonomous case where the nonsmooth sweeping set depends explicitly on time. The moving character of the sweeping set creates substantial additional difficulties, since several geometric…
Rasmus Vikhamar-Sandberg, Michal Repisky
We present an algorithm for efficient evaluation of Boys functions F0, . . . , Fkmax tailored to modern computing architectures, in particular graphical processing units (GPUs), where maximum throughput is high and data movement is costly. The method combines rational minimax approximations with upward and downward…
Jiechang Fan, Weiping Shen Weiping Shen, Yusong Luo, Enping Lou
In this paper, we investigate numerical solutions for inverse singular value problems (for short, ISVPs) arising in various applications. Inspired by the methodologies employed for inverse eigenvalue problems, we propose a Cayleyfree two-step algorithm for solving the ISVP. Compared to the existing twostep algorithms…
Andrea Gallo Rosso
Local polynomial smoothing is a widespread technique in data analysis, and Savitzky-Golay (SG) filters are one of its most well-known realizations. In real settings, the effectiveness of SG filtering depends critically on proper tuning of its parameters, constrained in turn by repeated polynomial fitting over large…
Pierre Borie, Fabian Bastin, Stéphane Dellacherie
We present an algorithm for solving nonlinear least-squares problems subject to a mix of nonlinear and linear constraints. The nonlinear constraints are handled by reformulating the objective as the augmented Lagrangian function while linear constraints are handled directly. Each iteration consists of approximately…
Umair Zulfiqar
This paper considers large-scale nonsymmetric continuous-time algebraic Riccati equations (NAREs) that admit low-rank solutions. Low-rank alternating direction implicit (ADI) methods have proven to be an efficient approach for solving several matrix equations, including Lyapunov equations, Sylvester equations, and…
Ruixue Gu, Qinian Jin
In this paper we develop a stochastic heavy ball method for solving ill-posed inverse problems. The method updates the iterate using only a randomly selected equation at each iteration step while incorporating a momentum term into the process. To facilitate fast convergence, we propose an adaptive strategy for…
J. Ibáñez, J. Sastre, J. M. Alonso, E. Defez
Computing numerical approximations of matrix functions frequently relies on the efficient evaluation of high-degree matrix polynomials. Although computational bounds are historically governed by the Paterson--Stockmeyer (PS) method, recent theoretical developments have demonstrated the viability of evaluation schemes…
de la Hoz, Francisco, Peru Muniain
In this paper, using a pseudospectral approach, we develop operational matrices based on the shifted Chebyshev polynomials to approximate numerically Caputo fractional derivatives and Riemann-Liouville fractional integrals. In order to make the generation of these matrices stable, we use variable precision arithmetic.…
Jordi Canela, Daniel Pérez-Palau
This paper presents a novel application of Jet Transport, a high-order automatic differentiation technique, to enhance classical numerical methods, with a focus on Newton's method. We prove a central theorem establishing that, under appropriate conditions, applying Jet Transport within a Newton iteration doubles the…