Fourier Methods for Piecewise Hermite Bicubic Or- Thogonal Spline Collocation

| Matrix decomposition algorithms employing fast Fourier transforms were developed recently by the authors to solve the systems of linear algebraic equations that arise when piecewise Hermite bicubic orthogonal spline collocation (OSC) is applied to certain separable elliptic boundary value problems on a rectangle. In this paper, these algorithms are interpreted as Fourier methods in analogy with Fourier methods for nite dierence schemes. This interpretation not only simplies the presentation of the OSC matrix decomposition algorithms but provides a natural unifying framework that closely parallels the Fourier method for the continuous problem. Like their nite dierence counterparts, OSC Fourier methods are based on knowledge of the solution of an eigenvalue problem associated with the corresponding discretization of a two-point boundary value problem. The eigenvalues and the associated matrix representation of the eigenfunctions are presented for various choices of boundary conditions.

[1]  Graeme Fairweather,et al.  Matrix decomposition algorithms for separable elliptic boundary value problems in two space dimensions , 1993 .

[2]  Bernard Bialecki,et al.  Orthogonal spline collocation Laplace-modified and alternating-direction methods for parabolic problems on rectangles , 1993 .

[3]  Gene H. Golub,et al.  On Fourier-Toeplitz methods for separable elliptic problems , 1974 .

[4]  Patrick Keast,et al.  Algorithm 603: COLROW and ARCECO: FORTRAN Packages for Solving Certain Almost Block Diagonal Linear Systems by Modified Alternate Row and Column Elimination , 1983, TOMS.

[5]  Bernard Bialecki,et al.  H 1 -norm error bounds for piecewise Hermite bicubic orthogonal spline collocation schemes for elliptic boundary value problems , 1994 .

[6]  Jim Douglas,et al.  Collocation Methods for Parabolic Equations in a Single Space Variable , 1974 .

[7]  Mary F. Wheeler,et al.  A $C^1 $ Finite Element Collocation Method for Elliptic Equations , 1980 .

[8]  Graeme Fairweather,et al.  Fast direct solvers for piecewise Hermite bicubic orthogonal spline collocation equations , 1992 .

[9]  Morgan Pickering An Introduction to Fast Fourier Transform Methods for Partial Differential Equations, with Applications , 1986 .

[10]  C. Loan Computational Frameworks for the Fast Fourier Transform , 1992 .

[11]  Bernard Bialecki,et al.  Preconditioned Richardson and Minimal Residual Iterative Methods for Piecewise Hermite Bicubic Orthogonal Spline Collocation Equations , 1994, SIAM J. Sci. Comput..

[12]  Roger W. Hockney,et al.  A Fast Direct Solution of Poisson's Equation Using Fourier Analysis , 1965, JACM.

[13]  Bernard Bialecki,et al.  Fourier Matrix Decomposition Methods for the Least Squares Solution of Singular Neumann and Periodic Hermite Bicubic Collocation Problems , 1995, SIAM J. Sci. Comput..

[14]  K. R. Bennett,et al.  Parallel collocation methods for boundary value problems , 1991 .