Symmetrization of the sinc-Galerkin method for boundary value problems

Lund, John

doi:10.1090/S0025-5718-1986-0856703-9

Symmetrization of the sinc-Galerkin method for boundary value problems
HTML articles powered by AMS MathViewer

by John Lund PDF

Math. Comp. 47 (1986), 571-588 Request permission

Abstract:

The Sinc-Galerkin method developed in [5], when applied to the second-order selfadjoint boundary value problem, gives rise to a nonsymmetric coefficient matrix. The technique in [5] is based on weighting the Galerkin inner products in such a way that the method will handle boundary value problems with regular singular points. In particular, the method does an accurate job of handling problems with singular solutions (the first or a higher derivative of the solution is unbounded at one or both of the boundary points). Using n function evaluations, the method of [5] converges at the rate $\exp ( - \kappa \sqrt n )$, where k is independent of n. In this paper it is shown that, by changing the weight function used in the Galerkin inner products, the coefficient matrix can be made symmetric. This symmetric method is applicable to a slightly more restrictive set of boundary value problems than the method of [5], The present method, however, still handles a wide class of singular problems and also has the same $\exp ( - \kappa \sqrt n )$ convergence rate.

References

Linda Kaufman and Daniel D. Warner, High-order, fast-direct methods for separable elliptic equations, SIAM J. Numer. Anal. 21 (1984), no. 4, 672–694. MR 749364, DOI 10.1137/0721046
John Lund, Sinc function quadrature rules for the Fourier integral, Math. Comp. 41 (1983), no. 163, 103–113. MR 701627, DOI 10.1090/S0025-5718-1983-0701627-8
John R. Lund and Bruce V. Riley, A sinc-collocation method for the computation of the eigenvalues of the radial Schrödinger equation, IMA J. Numer. Anal. 4 (1984), no. 1, 83–98. MR 740786, DOI 10.1093/imanum/4.1.83
S. W. Schoombie and J. F. Botha, Error estimates for the solution of the radial Schrödinger equation by the Rayleigh-Ritz finite element method, IMA J. Numer. Anal. 1 (1981), no. 1, 47–63. MR 607246, DOI 10.1093/imanum/1.1.47
Frank Stenger, A “sinc-Galerkin” method of solution of boundary value problems, Math. Comp. 33 (1979), no. 145, 85–109. MR 514812, DOI 10.1090/S0025-5718-1979-0514812-4
Frank Stenger, Numerical methods based on Whittaker cardinal, or sinc functions, SIAM Rev. 23 (1981), no. 2, 165–224. MR 618638, DOI 10.1137/1023037
A. Weiser, S. C. Eisenstat, and M. H. Schultz, On solving elliptic equations to moderate accuracy, SIAM J. Numer. Anal. 17 (1980), no. 6, 908–929. MR 595453, DOI 10.1137/0717075