Multigrid methods for cubic spline solution of two point (and 2D) boundary value problems. (English) Zbl 1336.65120
Summary: 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.
MSC:
| 65L10 | Numerical solution of boundary value problems involving ordinary differential equations |
| 65N55 | Multigrid methods; domain decomposition for boundary value problems involving PDEs |
| 35J25 | Boundary value problems for second-order elliptic equations |
| 65L60 | Finite element, Rayleigh-Ritz, Galerkin and collocation methods for ordinary differential equations |
| 65N30 | Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs |
Keywords:
cubic splines; finite elements; multigrid methods; spectral analysis; Toeplitz matrices; symbol; numerical examples; second-order two point boundary value problem; numerical experimentsReferences:
| [1] | Aricò, A.; Donatelli, M., A V-cycle multigrid for multilevel matrix algebras: proof of optimality, Numer. Math., 105, 4, 511-547 (2007) · Zbl 1114.65033 |
| [2] | Aricò, A.; Donatelli, M.; Serra-Capizzano, S., V-cycle optimal convergence for certain (multilevel) structured linear systems, SIAM J. Matrix Anal. Appl., 26, 1, 186-214 (2004) · Zbl 1105.65322 |
| [3] | Awanou, G.; Lai, M., Trivariate spline approximations of 3D Navier-Stokes equations, Math. Comput., 74, 250, 585-601 (2005) · Zbl 1085.76053 |
| [4] | Axelsson, O., Iterative Solution Methods (1994), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0795.65014 |
| [5] | Chawla, M. M.; Subramanian, R., A fourth order spline method for singular two-point boundary value problems, Int. Schriftenreihe Numer. Math., 86, 87-93 (1988) · Zbl 0652.65058 |
| [6] | Donatelli, M., An algebraic generalization of local Fourier analysis for grid transfer operators in multigrid based on Toeplitz matrices, Numer. Linear Algebra Appl., 17, 2-3, 179-197 (2010) · Zbl 1240.65353 |
| [7] | Höllig, K.; Reif, U.; Wipper, J., Multigrid methods with web-splines, Numer. Math., 91, 2, 237-255 (2002) · Zbl 0996.65138 |
| [8] | Huckle, T.; Staudacher, J., Multigrid methods for block Toeplitz matrices with small size blocks, BIT, 46, 61-83 (2006) · Zbl 1103.65035 |
| [9] | Iyengar, S. R.K.; Jain, P., Spline finite difference methods for singular two point boundary value problems, Numer. Math., 50, 3, 363-376 (1986) · Zbl 0642.65062 |
| [10] | Kim, S. D.; Parter, S. V., Preconditioning cubic spline collocation discretizations of elliptic equations, Numer. Math., 72, 1, 39-72 (1995) · Zbl 0844.65086 |
| [11] | Kumar, M., Higher order method for singular boundary-value problems by using spline function, Appl. Math. Comput., 192, 1, 175-179 (2007) · Zbl 1193.65127 |
| [12] | Mohanty, R. K.; Dahiva, V., An \(O(k 2 + k h 2 + h 4)\) accurate two-level implicit cubic spline method for one space dimensional quasi-linear parabolic equations, AJCM, 1, 1, 11-17 (2011) |
| [13] | Quarteroni, A., Modellistica Numerica per Problemi Differenziali (2006), Springer-Verlag: Springer-Verlag Milano, Italia · Zbl 1029.65091 |
| [14] | Quarteroni, A.; Valli, A., Numerical Approximation of Partial Differential Equations (1994), Springer-Verlag: Springer-Verlag Berlin · Zbl 0803.65088 |
| [15] | Rashidinia, J.; Mohammadi, R.; Ghasemi, M., Cubic spline solution of singularly perturbed boundary value problems with significant first derivatives, Appl. Math. Comput., 190, 2, 1762-1766 (2007) · Zbl 1227.65065 |
| [16] | Rashidinia, J.; Mohammadi, R.; Jalilian, R., Cubic spline method for two point boundary value problems, IUST Int. J. Eng. Sci., 19, 5-2, 39-43 (2008) · Zbl 1131.65065 |
| [17] | Serra-Capizzano, S., Multi-iterative methods, Comput. Math. Appl., 26, 4, 65-87 (1993) · Zbl 0790.65025 |
| [18] | Serra-Capizzano, S., Asymptotic results on the spectra of block Toeplitz preconditioned matrices, SIAM J. Matrix Anal. Appl., 20, 1, 31-44 (1999) · Zbl 0932.65037 |
| [19] | Serra-Capizzano, S., Distribution results on the algebra generated by Toeplitz sequences: a finite dimensional approach, Linear Algebra Appl., 328, 1-3, 121-130 (2001) · Zbl 1003.15008 |
| [20] | Serra-Capizzano, S., Convergence analysis of two-grid methods for elliptic Toeplitz and PDEs matrix-sequences, Numer. Math., 92, 3, 433-465 (2002) · Zbl 1013.65026 |
| [21] | Serra-Capizzano, S., Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations, Linear Algebra Appl., 366, 371-402 (2003) · Zbl 1028.65109 |
| [22] | Stoer, J.; Bulirsch, R., Introduction to Numerical Analysis (2002), Springer-Verlag: Springer-Verlag New York · Zbl 1004.65001 |
| [23] | Tilli, P., A note on the spectral distribution of Toeplitz matrices, Linear Multilinear Algebra, 45, 2-3, 147-159 (1998) · Zbl 0951.65033 |
| [24] | Tyrtyshnikov, E. E., A unifying approach to some old and new theorems on distribution and clustering, Linear Algebra Appl., 232, 1-43 (1996) · Zbl 0841.15006 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
