When solving large linear systems stemming from the approximation of elliptic partial differential equations (PDEs), it is known that the V-cycle multigrid method (MGM) can significantly lower the computational cost. Many convergence estimates already exist for the V-cycle MGM: for example, using the regularity or approximation assumptions of the elliptic PDEs, the results are obtained in Bank and Douglas (SIAM J Numer Anal 22:617–633, 1985), Bramble and Pasciak (Math Comp 49:311–329, 1987); in the case of multilevel matrix algebras (like circulant, tau, Hartely) (Aricò et al. in SIAM J Matrix Anal Appl 26:186–214, 2004; Aricò and Donatelli in Numer Math 105:511–547, 2007), special prolongation operators are provided and the related convergence results are rigorously developed, using a functional approach. In this paper we derive new uniform convergence estimates for the V-cycle MGM applied to symmetric positive definite Toeplitz block tridiagonal matrices, by also discussing few connections with previous results. More concretely, the contributions of this paper are as follows: (1) It tackles the Toeplitz systems directly for the elliptic PDEs. (2) Simple (traditional) restriction operator and prolongation operator are employed in order to handle general Toeplitz systems at each level of the recursion. Such a technique is then applied to systems of algebraic equations generated by the difference scheme of the two-dimensional fractional Feynman–Kac equation, which describes the joint probability density function of non-Brownian motion. In particular, we consider the two coarsening strategies, i.e., doubling the mesh size (geometric MGM) and Galerkin approach (algebraic MGM), which lead to the distinct coarsening stiffness matrices in the general case: however, several numerical experiments show that the two algorithms produce almost the same error behaviour.

Uniform Convergence of V-cycle Multigrid Algorithms for Two-Dimensional Fractional Feynman–Kac Equation

Serra-Capizzano, Stefano
2018-01-01

Abstract

When solving large linear systems stemming from the approximation of elliptic partial differential equations (PDEs), it is known that the V-cycle multigrid method (MGM) can significantly lower the computational cost. Many convergence estimates already exist for the V-cycle MGM: for example, using the regularity or approximation assumptions of the elliptic PDEs, the results are obtained in Bank and Douglas (SIAM J Numer Anal 22:617–633, 1985), Bramble and Pasciak (Math Comp 49:311–329, 1987); in the case of multilevel matrix algebras (like circulant, tau, Hartely) (Aricò et al. in SIAM J Matrix Anal Appl 26:186–214, 2004; Aricò and Donatelli in Numer Math 105:511–547, 2007), special prolongation operators are provided and the related convergence results are rigorously developed, using a functional approach. In this paper we derive new uniform convergence estimates for the V-cycle MGM applied to symmetric positive definite Toeplitz block tridiagonal matrices, by also discussing few connections with previous results. More concretely, the contributions of this paper are as follows: (1) It tackles the Toeplitz systems directly for the elliptic PDEs. (2) Simple (traditional) restriction operator and prolongation operator are employed in order to handle general Toeplitz systems at each level of the recursion. Such a technique is then applied to systems of algebraic equations generated by the difference scheme of the two-dimensional fractional Feynman–Kac equation, which describes the joint probability density function of non-Brownian motion. In particular, we consider the two coarsening strategies, i.e., doubling the mesh size (geometric MGM) and Galerkin approach (algebraic MGM), which lead to the distinct coarsening stiffness matrices in the general case: however, several numerical experiments show that the two algorithms produce almost the same error behaviour.
2018
https://rd.springer.com/journal/10915
Block tridiagonal matrix; Fractional Feynman–Kac equation; V-cycle multigrid method; Software; Theoretical Computer Science; Engineering (all); Computational Theory and Mathematics
Chen, Minghua; Deng, Weihua; Serra-Capizzano, Stefano
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/2073483
 Attenzione

L'Ateneo sottopone a validazione solo i file PDF allegati

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 8
  • ???jsp.display-item.citation.isi??? 7
social impact