Computation of the Complex Error Function using Modified Trapezoidal Rules

10/12/2020
by   Mohammad Al Azah, et al.
0

In this paper we propose a method for computing the Faddeeva function w(z) := e^-z^2erfc(-i z) via truncated modified trapezoidal rule approximations to integrals on the real line. Our starting point is the method due to Matta and Reichel (Math. Comp. 25 (1971), pp. 339-344) and Hunter and Regan (Math. Comp. 26 (1972), pp. 339-541). Addressing shortcomings flagged by Weideman (SIAM. J. Numer. Anal. 31 (1994), pp. 1497-1518), we construct approximations which we prove are exponentially convergent as a function of N+1, the number of quadrature points, obtaining error bounds which show that accuracies of 2× 10^-15 in the computation of w(z) throughout the complex plane are achieved with N = 11, this confirmed by computations. These approximations, moreover, provably achieve small relative errors throughout the upper complex half-plane where w(z) is non-zero. Numerical tests show that this new method is competitive, in accuracy and computation times, with existing methods for computing w(z) for complex z.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
09/24/2021

Computing the Barnes G-function in the entire complex plane

We present an algorithm for generating approximations for the logarithm ...
research
09/22/2015

Modifying iterated Laplace approximations

In this paper, several modifications are introduced to the functional ap...
research
11/26/2017

Algorithms for the Computing Determinants in Commutative Rings

Two known computation methods and one new computation method for matrix ...
research
06/28/2019

Evaluation of Abramowitz functions in the right half of the complex plane

A numerical scheme is developed for the evaluation of Abramowitz functio...
research
04/13/2021

Characterizations of the maximum likelihood estimator of the Cauchy distribution

We consider characterizations of the maximal likelihood estimator (MLE) ...
research
09/29/2020

Numerical approximations of chromatographic models

A numerical scheme based on modified method of characteristics with adju...
research
07/11/2018

On error representation in exact-decisions number types

Accuracy-driven computation is a strategy widely used in exact-decisions...

Please sign up or login with your details

Forgot password? Click here to reset