In this work, we examine the numerical solution of a 1D distributed-order space-fractional diffusion equation. Discretizing the given problem by means of an implicit finite difference scheme based on the shifted Grünwald-Letnikov formula, the resulting linear systems show a Toeplitz structure. Then, by using well-known spectral tools for Toeplitz sequences, we determine the corresponding symbol describing its asymptotic eigenvalue distribution as the matrix size diverges. The spectral analysis is performed under different assumptions with the aim of estimating the intrinsic asymptotic ill-conditioning of the involved matrices. The obtained results suggest to precondition the involved linear systems with either a Laplacian-like preconditioner or with more general τpreconditioners. Due to the symmetric positive definite nature of the coefficient matrices, we opt for the preconditioned conjugate gradient method, and we compare the performances of our proposal with a Strang circulant alternative given in the literature.

Symbol-based preconditioning for riesz distributed-order space-fractional diffusion equations

Mazza M.
;
Serra Capizzano S.
;
2021-01-01

Abstract

In this work, we examine the numerical solution of a 1D distributed-order space-fractional diffusion equation. Discretizing the given problem by means of an implicit finite difference scheme based on the shifted Grünwald-Letnikov formula, the resulting linear systems show a Toeplitz structure. Then, by using well-known spectral tools for Toeplitz sequences, we determine the corresponding symbol describing its asymptotic eigenvalue distribution as the matrix size diverges. The spectral analysis is performed under different assumptions with the aim of estimating the intrinsic asymptotic ill-conditioning of the involved matrices. The obtained results suggest to precondition the involved linear systems with either a Laplacian-like preconditioner or with more general τpreconditioners. Due to the symmetric positive definite nature of the coefficient matrices, we opt for the preconditioned conjugate gradient method, and we compare the performances of our proposal with a Strang circulant alternative given in the literature.
2021
2021
Fractional diffusion equations; Preconditioning; Spectral distribution; Toeplitz matrices
Mazza, M.; Serra Capizzano, S.; Usman, M.
File in questo prodotto:
File Dimensione Formato  
mazza 2021 etna.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: DRM non definito
Dimensione 801.13 kB
Formato Adobe PDF
801.13 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

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/2119488
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 4
  • ???jsp.display-item.citation.isi??? 3
social impact