Multigrid methods for cubic spline solution of two point (and 2D) boundary value problems

In this paper we propose a scheme based on cubic splines for the solution of the second order two point boundary value problems. The solution of the algebraic system is computed by using optimized multigrid methods. In particular the transformation of the stiffness matrix essentially in a block Toeplitz matrix and its spectral analysis allow to choose smoothers able to reduce error components related to the various frequencies and to obtain an optimal method. The main advantages of our strategy can be listed as follows: (i) a fourth order of accuracy combined with a quadratic conditioning matrix, (ii) a resulting matrix structure whose eigenvalues can be compactly described by a symbol (this information is the key for designing an optimal multigrid method). Finally, some numerics that confirm the predicted behavior of the method are presented and a discussion on the two dimensional case is given, together with few 2D numerical experiments.

[1]  Ming-Jun Lai,et al.  Trivariate spline approximations of 3D Navier-Stokes equations , 2004, Math. Comput..

[2]  J. Stoer,et al.  Introduction to Numerical Analysis , 2002 .

[3]  Stefano Serra Capizzano,et al.  V-cycle Optimal Convergence for Certain (Multilevel) Structured Linear Systems , 2004, SIAM J. Matrix Anal. Appl..

[4]  Marco Donatelli,et al.  A V-cycle Multigrid for multilevel matrix algebras: proof of optimality , 2007, Numerische Mathematik.

[5]  Alfio Quarteroni,et al.  Modellistica numerica per problemi differenziali , 2008 .

[6]  Thomas Huckle,et al.  Multigrid Methods for Block Toeplitz Matrices with Small Size Blocks , 2006 .

[7]  J. Rashidinia,et al.  Cubic spline solution of singularly perturbed boundary value problems with significant first derivatives , 2007, Appl. Math. Comput..

[8]  Seymour V. Parter,et al.  Preconditioning cubic spline collocation discretizations of elliptic equations , 1995 .

[9]  Manoj Kumar.  Higher order method for singular boundary-value problems by using spline function , 2007, Appl. Math. Comput..

[10]  Eisa A. Al-Said.  Cubic spline method for solving two-point boundary-value problems , 1998 .

[11]  R. K. Mohanty,et al.  An O(k2+kh2+h2) Accurate Two-level Implicit Cubic Spline Method for One Space Dimensional Quasi-linear Parabolic Equations , 2011, Am. J. Comput. Math..

[12]  Stefano Serra,et al.  Asymptotic Results on the Spectra of Block Toeplitz Preconditioned Matrices , 1999 .

[13]  S. Serra Capizzano,et al.  Distribution results on the algebra generated by Toeplitz sequences: a finite-dimensional approach , 2001 .

[14]  Ulrich Reif,et al.  Multigrid methods with web-splines , 2002, Numerische Mathematik.

[15]  S. R. K. Iyengar,et al.  Spline finite difference methods for singular two point boundary value problems , 1986 .

[16]  Sabine Fenstermacher,et al.  Numerical Approximation Of Partial Differential Equations , 2016 .

[17]  U. Grenader.  ASYMPTOTIC RESULTS ON THE SPECTRA OF BLOCK TOEPLITZ PRECONDITIONED MATRICES , 1998 .

[18]  S. Serra,et al.  Multi-iterative methods , 1993 .

[19]  Paolo Tilli,et al.  A note on the spectral distribution of toeplitz matrices , 1998 .

[20]  Stefano Serra Capizzano,et al.  Numerische Mathematik Convergence analysis of two-grid methods for elliptic Toeplitz and PDEs Matrix-sequences , 2002 .

[21]  B-spline Method for Two-Point Boundary Value Problems , 2014 .

[22]  M. M. Chawla,et al.  A fourth-order spline method for singular two-point boundary-value problems , 1988 .

[23]  E. E. Tyrtyshnikov.  A unifying approach to some old and new theorems on distribution and clustering , 1996 .

[24]  O. Axelsson.  Iterative solution methods , 1995 .

[25]  Marco Donatelli,et al.  An algebraic generalization of local Fourier analysis for grid transfer operators in multigrid based on Toeplitz matrices , 2010, Numer. Linear Algebra Appl..

[26]  S. Serra Capizzano,et al.  Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations , 2003 .