Search

A class of methods constructed to numerically approximate the solution of two-point singularly perturbed boundary value problems of the form ϵu'' + bu' + cu = f use exponentials to mimic exponential behavior of the solution in the boundary layer(s). We refer to them as exponentially fitted methods. Such methods are usually exact on polynomials of certain degree and some exponential functions. Shortly, they are exact on exponential sums. It is often possible that consistency of the method follows from the convergence of the interpolating function standing behind the method. Because of that, we consider the interpolation error for exponential sums. The main result of the paper is an error bound for interpolation by the exponential sum to the solution of the singularly perturbed problem that does not depend on perturbation parameter ϵ when ϵ is small with the respect to mesh width. The numerical experiment implies that the use of a dense mesh in the boundary layer for small meshwidth results in ϵ-uniform convergence. [ABSTRACT FROM AUTHOR]