Preconditioned GMRES for oscillatory integrals
نویسنده
چکیده
None of the existing methods for computing the oscillatory integral ∫ b a f(x)e iωg(x) dx achieve all of the following properties: high asymptotic order, stability, avoiding the computation of the path of steepest descent and insensitivity to oscillations in f . We present a new method that satisfies these properties, based on applying the gmres algorithm to a preconditioned differential operator.
منابع مشابه
Preconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملShifted GMRES for oscillatory integrals
None of the existing methods for computing the oscillatory integral ∫ b a f(x)e iωg(x) dx achieve all of the following properties: high asymptotic order, stability, avoiding deformation into the complex plane and insensitivity to oscillations in f . We present a new method that satisfies these properties, based on applying the gmres algorithm to a shifted linear differential operator.
متن کاملGMRES for the Differentiation Operator
We investigate using the gmres method with the differentiation operator. This operator is unbounded, and thus does not fall into the framework of existing Krylov subspace theory. We establish conditions under which a function can be approximated by its own derivatives in a domain of the complex plane. These conditions are used to determine when gmres converges. This algorithm outperforms tradit...
متن کاملTWO LOW-ORDER METHODS FOR THE NUMERICAL EVALUATION OF CAUCHY PRINCIPAL VSlLUE INTEGRALS OF OSCILLATORY KIND
In this paper, we develop two piecewise polynomial methods for the numerical evaluation of Cauchy Principal Value integrals of oscillatory kind. The two piecewisepolynomial quadratures are compact, easy to implement, and are numerically stable. Two numerical examples are presented to illustrate the two rules developed, The convergence of the two schemes is proved and some error bounds obtai...
متن کاملSolving large systems arising from fractional models by preconditioned methods
This study develops and analyzes preconditioned Krylov subspace methods to solve linear systems arising from discretization of the time-independent space-fractional models. First, we apply shifted Grunwald formulas to obtain a stable finite difference approximation to fractional advection-diffusion equations. Then, we employee two preconditioned iterative methods, namely, the preconditioned gen...
متن کامل