Abstract In this paper we discuss classical sufficient conditions to be satisfied from the grid transfer operators in order to obtain optimal two-grid and V-cycle multigrid methods utilizing the theory for Toeplitz matrices. We derive relaxed conditions that allow the construction of special grid transfer operators that are computationally less expensive while preserving optimality. This is particularly useful when the generating symbol of the system matrix has a zero of higher order, like in the case of higher order PDEs. These newly derived conditions allow the use of rank deficient grid transfer operators. In this case the use of a pre-relaxation iteration that is lacking the smoothing property is proposed. Combining these pre-relaxations with the new rank deficient grid transfer operators yields a substantial reduction of the convergence rate and of the computational cost at each iteration compared with the classical choice for Toeplitz matrices. The proposed strategy, i.e. a rank deficient grid transfer operator plus a specific pre-relaxation, is applied to linear systems whose system matrix is a Toeplitz matrix where the generating symbol is a high-order polynomial. The necessity of using high-order polynomials as generating symbols for the grid transfer operators usually destroys the Toeplitz structure on the coarser levels. Therefore, we discuss some effective and computational cheap coarsening strategies found in the literature. In particular, we present numerical results showing near-optimal behavior while keeping the Toeplitz structure on the coarser levels.

Generalized grid transfer operators for multigrid methods applied on Toeplitz matrices

DONATELLI, MARCO;
2015-01-01

Abstract

Abstract In this paper we discuss classical sufficient conditions to be satisfied from the grid transfer operators in order to obtain optimal two-grid and V-cycle multigrid methods utilizing the theory for Toeplitz matrices. We derive relaxed conditions that allow the construction of special grid transfer operators that are computationally less expensive while preserving optimality. This is particularly useful when the generating symbol of the system matrix has a zero of higher order, like in the case of higher order PDEs. These newly derived conditions allow the use of rank deficient grid transfer operators. In this case the use of a pre-relaxation iteration that is lacking the smoothing property is proposed. Combining these pre-relaxations with the new rank deficient grid transfer operators yields a substantial reduction of the convergence rate and of the computational cost at each iteration compared with the classical choice for Toeplitz matrices. The proposed strategy, i.e. a rank deficient grid transfer operator plus a specific pre-relaxation, is applied to linear systems whose system matrix is a Toeplitz matrix where the generating symbol is a high-order polynomial. The necessity of using high-order polynomials as generating symbols for the grid transfer operators usually destroys the Toeplitz structure on the coarser levels. Therefore, we discuss some effective and computational cheap coarsening strategies found in the literature. In particular, we present numerical results showing near-optimal behavior while keeping the Toeplitz structure on the coarser levels.
2015
BIT
http://www.kluweronline.com/issn/0006-3835
Grid transfer operators; Multigrid methods; Toeplitz matrices; Computer Networks and Communications; Software; Applied Mathematics; Computational Mathematics
Bolten, Matthias; Donatelli, Marco; Huckle, Thomas; Kravvaritis, Christos
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/2022634
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

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