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 traditional quadrature schemes for a large class of highly oscillatory integrals, even when the kernel of oscillations is unknown.

[1]  Daan Huybrechs,et al.  On the Evaluation of Highly Oscillatory Integrals by Analytic Continuation , 2006, SIAM J. Numer. Anal..

[2]  N. S. Bakhvalov,et al.  Evaluation of the integrals of oscillating functions by interpolation at nodes of gaussian quadratures , 1968 .

[3]  C. W. Clenshaw,et al.  A method for numerical integration on an automatic computer , 1960 .

[4]  Andreas Frommer,et al.  Restarted GMRES for Shifted Linear Systems , 1998, SIAM J. Sci. Comput..

[5]  Lloyd N. Trefethen,et al.  Is Gauss Quadrature Better than Clenshaw-Curtis? , 2008, SIAM Rev..

[6]  Lloyd N. Trefethen,et al.  An Extension of MATLAB to Continuous Functions and Operators , 2004, SIAM J. Sci. Comput..

[7]  Y. Saad,et al.  Arnoldi methods for large Sylvester-like observer matrix equations, and an associated algorithm for partial spectrum assignment , 1991 .

[8]  F. Olver.  Asymptotics and Special Functions , 1974 .

[9]  A. Iserles,et al.  Efficient quadrature of highly oscillatory integrals using derivatives , 2005, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[10]  Sheehan Olver,et al.  Moment-free numerical approximation of highly oscillatory integrals with stationary points , 2007, European Journal of Applied Mathematics.

[11]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[12]  O. Nevanlinna.  Convergence of Iterations for Linear Equations , 1993 .

[13]  Sheehan Olver,et al.  Moment-free numerical integration of highly oscillatory functions , 2006 .

[14]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[15]  P. Miller.  Applied asymptotic analysis , 2006 .

[16]  Philipp Birken,et al.  Numerical Linear Algebra , 2011, Encyclopedia of Parallel Computing.