Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
skip to main content
research-article

Matrix Decomposition Algorithms for Modified Spline Collocation for Helmholtz Problems

Published: 01 January 2003 Publication History
  • Get Citation Alerts
  • Abstract

    We consider the solution of various boundary value problems for the Helmholtz equation in the unit square using a nodal cubic spline collocation method and modifications of it which produce optimal (fourth-) order approximations. For the solution of the collocation equations, we formulate matrix decomposition algorithms, fast direct methods which employ fast Fourier transforms and require O(N2 log N) operations on an N × N uniform partition of the unit square. A computational study confirms the published analysis for the Dirichlet problem and indicates that similar results hold for Neumann, mixed, and periodic boundary conditions. The numerical results also exhibit superconvergence phenomena not reported in earlier studies.

    References

    [1]
    D. Archer, Some Collocation Methods for Differential Equations, Ph.D. thesis, Rice University, Houston, TX, 1973.
    [2]
    David Archer, An O(h4) cubic spline collocation method for quasilinear parabolic equations, SIAM J. Numer. Anal., 14 (1977), 620–637
    [3]
    Bernard Bialecki, Graeme Fairweather, Matrix decomposition algorithms for separable elliptic boundary value problems in two space dimensions, J. Comput. Appl. Math., 46 (1993), 369–386
    [4]
    B. Bialecki, G. Fairweather, Orthogonal spline collocation methods for partial differential equations, J. Comput. Appl. Math., 128 (2001), 55–82, Numerical analysis 2000, Vol. VII, Partial differential equations
    [5]
    C. de Boor, The Method of Projections as Applied to the Numerical Solution of Two Point Boundary Value Problems Using Cubic Splines, Ph.D. thesis, University of Michigan, Ann Arbor, MI, 1966.
    [6]
    C. Bottcher and M. R. Strayer, The basis spline method and associated techniques, in Computational Atomic and Nuclear Physics, C. Bottcher, M. R. Strayer, and J. B. McGrory, eds., World Scientific, Singapore, 1990, pp. 217–240.
    [7]
    C. Bottcher, M. Strayer, Spline methods for conservationequations, North‐Holland, Amsterdam, 1993, 317–338
    [8]
    R. Carlson, C. Hall, Error bounds for bicubic spline interpolation, J. Approximation Theory, 7 (1973), 41–47
    [9]
    C. C. Christara, Schur complement preconditioned conjugate gradient methods for spline collocation equations, Computer Architecture News, 18 (1990), pp. 108–120.
    [10]
    Christina Christara, Quadratic spline collocation methods for elliptic partial differential equations, BIT, 34 (1994), 33–61
    [11]
    C. C. Christara, Parallel solvers for spline collocation equations, Advances in Engineering Software, 27 (1996), pp. 71–89.
    [12]
    Christina Christara, Barry Smith, Multigrid and multilevel methods for quadratic spline collocation, BIT, 37 (1997), 781–803
    [13]
    J. Daniel, B. Swartz, Extrapolated collocation for two‐point boundary‐value problems using cubic splines, J. Inst. Math. Appl., 16 (1975), 161–174
    [14]
    G. Fairweather, A. Karageorghis, and B. Bialecki, Matrix Decomposition Algorithms for Modified Spline Collocation for Separable Elliptic Boundary Value Problems. I. Helmholtz Problems, Technical report TR/01/2000, Department of Mathematics and Statistics, University of Cyprus, Nicosia, Cyprus, revised, 2001.
    [15]
    D. Fyfe, The use of cubic splines in the solution of two‐point boundary value problems, Comput. J., 12 (1969/1970), 188–192
    [16]
    A. Hadjidimos, E. Houstis, J. Rice, E. Vavalis, Iterative line cubic spline collocation methods for elliptic partial differential equations in several dimensions, SIAM J. Sci. Comput., 14 (1993), 715–734
    [17]
    A. Hadjidimos, E. Houstis, J. Rice, E. Vavalis, Analysis of iterative line spline collocation methods for elliptic partial differential equations, SIAM J. Matrix Anal. Appl., 21 (1999), 508–521
    [18]
    E. Houstis, C. Christara, J. Rice, Quadratic‐spline collocation methods for two‐point boundary value problems, Internat. J. Numer. Methods Engrg., 26 (1988), 935–952
    [19]
    E. N. Houstis, J. R. Rice, and E. A. Vavalis, Parallelization of a new class of cubic spline collocation methods, in Advances in Computer Methods for Partial Differential Equations VI, R. Vichnevetsky and R. S. Stepleman, eds., IMACS, New Brunswick, NJ, 1987, pp. 167–174.
    [20]
    E. Houstis, E. Vavalis, J. Rice, Convergence of O(h4) cubic spline collocation methods for elliptic partial differential equations, SIAM J. Numer. Anal., 25 (1988), 54–74
    [21]
    M. Irodotou‐Ellina, E. Houstis, An O(h6) quintic spline collocation method for fourth order two‐point boundary value problems, BIT, 28 (1988), 288–301
    [22]
    D. Kegley, Jr., V. Oberacker, M. Strayer, A. Umar, J. Wells, Basis spline collocation method for solving the Schrödinger equation in axillary symmetric systems, J. Comput. Phys., 128 (1996), 197–208
    [23]
    P. Prenter, Splines and variational methods, Wiley‐Interscience [John Wiley & Sons], 1975xi+323, Pure and Applied Mathematics
    [24]
    D. Rutherford, Some continuant determinants arising in physics and chemistry. II, Proc. Roy. Soc. Edinburgh. Sect. A., 63 (1952), 232–241
    [25]
    A. A. Samarskii and E. S. Nikolaev, Numerical Methods for Grid Equations, Vol. I, Direct Methods, Birkhäuser Verlag, Basel, Boston, Berlin, 1989.
    [26]
    P. N. Swarztrauber, FFTPACK, NCAR, Boulder, CO, 1985;
    also available online from [email protected].
    [27]
    P. Tsompanopoulou, E. Vavalis, ADI methods for cubic spline collocation discretizations of elliptic PDEs, SIAM J. Sci. Comput., 19 (1998), 341–363
    [28]
    A. S. Umar, Three‐dimensional HF and TDHF calculations with the basis‐spline collocation technique, in Computational Atomic and Nuclear Physics, C. Bottcher, M. R. Strayer, and J. B. McGrory, eds., World Scientific, Singapore, 1990, pp. 377–390.
    [29]
    A. Umar, J. Wu, M. Strayer, C. Bottcher, Basis‐spline collocation method for the lattice solution of boundary value problems, J. Comput. Phys., 93 (1991), 426–448
    [30]
    Charles Van Loan, Computational frameworks for the fast Fourier transform, Frontiers in Applied Mathematics, Vol. 10, Society for Industrial and Applied Mathematics (SIAM), 1992xiv+273

    Cited By

    View all

    Recommendations

    Comments

    Information & Contributors

    Information

    Published In

    cover image SIAM Journal on Scientific Computing
    SIAM Journal on Scientific Computing  Volume 24, Issue 5
    2003
    373 pages

    Publisher

    Society for Industrial and Applied Mathematics

    United States

    Publication History

    Published: 01 January 2003

    Author Tags

    1. 65N35
    2. 65N22

    Author Tags

    1. spline collocation
    2. Helmholtz equation
    3. boundary conditions
    4. tensor product
    5. fast Fourier transforms
    6. convergence rates
    7. superconvergence

    Qualifiers

    • Research-article

    Contributors

    Other Metrics

    Bibliometrics & Citations

    Bibliometrics

    Article Metrics

    • Downloads (Last 12 months)0
    • Downloads (Last 6 weeks)0
    Reflects downloads up to 12 Aug 2024

    Other Metrics

    Citations

    Cited By

    View all
    • (2022)Matrix decomposition algorithms for the C 0-quadratic finite element Galerkin methodBIT10.1007/s10543-009-0233-049:3(509-526)Online publication date: 11-Mar-2022
    • (2022)Optimal superconvergent one step quadratic spline collocation methods BIT10.1007/s10543-008-0188-648:3(449-472)Online publication date: 11-Mar-2022
    • (2019)A quadratic spline collocation method for the Dirichlet biharmonic problemNumerical Algorithms10.1007/s11075-019-00676-z83:1(165-199)Online publication date: 19-Feb-2019
    • (2015)A fast direct solver for a fourth order finite difference scheme for Poisson's equation on the unit disc in polar coordinatesNumerical Algorithms10.1007/s11075-015-9971-z70:4(727-751)Online publication date: 1-Dec-2015
    • (2013)Modified nodal cubic spline collocation for three-dimensional variable coefficient second order partial differential equationsNumerical Algorithms10.1007/s11075-012-9669-464:2(349-383)Online publication date: 1-Oct-2013
    • (2011)Matrix decomposition algorithms for elliptic boundary value problemsNumerical Algorithms10.1007/s11075-010-9384-y56:2(253-295)Online publication date: 1-Feb-2011
    • (2005)Optimal Superconvergent One Step Nodal Cubic Spline Collocation MethodsSIAM Journal on Scientific Computing10.1137/04060979327:2(575-598)Online publication date: 1-Jan-2005

    View Options

    View options

    Get Access

    Login options

    Media

    Figures

    Other

    Tables

    Share

    Share

    Share this Publication link

    Share on social media