Fast solvers for two-dimensional fractional diffusion equations using rank structured matrices
From MaRDI portal
(Redirected from Publication:5230663)
Matrix equations and identities (15A24) Fractional partial differential equations (35R11) Iterative numerical methods for linear systems (65F10) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs (65M22) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60)
Abstract: We consider the discretization of time-space diffusion equations with fractional derivatives in space and either 1D or 2D spatial domains. The use of implicit Euler scheme in time and finite differences or finite elements in space, leads to a sequence of dense large scale linear systems describing the behavior of the solution over a time interval. We prove that the coefficient matrices arising in the 1D context are rank structured and can be efficiently represented using hierarchical formats (-matrices, HODLR). Quantitative estimates for the rank of the off-diagonal blocks of these matrices are presented. We analyze the use of HODLR arithmetic for solving the 1D case and we compare this strategy with existing methods that exploit the Toeplitz-like structure to precondition the GMRES iteration. The numerical tests demonstrate the convenience of the HODLR format when at least a reasonably low number of time steps is needed. Finally, we explain how these properties can be leveraged to design fast solvers for problems with 2D spatial domains that can be reformulated as matrix equations. The experiments show that the approach based on the use of rank-structured arithmetic is particularly effective and outperforms current state of the art techniques.
Recommendations
- Low-rank solvers for fractional differential equations
- A fast finite difference method for two-dimensional space-fractional diffusion equations
- Fast solvers for finite difference scheme of two-dimensional time-space fractional differential equations
- Fast iterative solvers for linear systems arising from time-dependent space-fractional diffusion equations
- A fast finite difference method for three-dimensional time-dependent space-fractional diffusion equations and its efficient implementation
Cites work
- A Block Lanczos Method for Computing the Singular Values and Corresponding Singular Vectors of a Matrix
- A circulant preconditioner for fractional diffusion equations
- A class of second order difference approximations for solving space fractional diffusion equations
- A direct \(O(N \log ^{2} N)\) finite difference method for fractional diffusion equations
- A finite element formulation preserving symmetric and banded diffusion stiffness matrix characteristics for fractional differential equations
- A New Iterative Method for Solving Large-Scale Lyapunov Matrix Equations
- A sparse matrix arithmetic based on \({\mathfrak H}\)-matrices. I: Introduction to \({\mathfrak H}\)-matrices
- Adaptive finite element method for fractional differential equations using hierarchical matrices
- Adaptive geometrically balanced clustering of \({\mathcal H}\)-matrices
- An error analysis for rational Galerkin projection applied to the Sylvester equation
- An Extension of Chebfun to Two Dimensions
- Boundary Problems for the Fractional and Tempered Fractional Operators
- Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in \(\mathbb R^2\)
- Computational Methods for Linear Matrix Equations
- Convergence analysis of the extended Krylov subspace method for the Lyapunov equation
- Discontinuous Galerkin time stepping method for solving linear space fractional partial differential equations
- Extended Arnoldi methods for large low-rank Sylvester matrix equations
- Fast algorithms for hierarchically semiseparable matrices
- Fast iterative solvers for linear systems arising from time-dependent space-fractional diffusion equations
- Fast preconditioned iterative methods for finite volume discretization of steady-state space-fractional diffusion equations
- Finite difference approximations for fractional advection-dispersion flow equations
- Fractional elliptic equations, Caccioppoli estimates and regularity
- Generalized locally Toeplitz sequences: theory and applications. Volume I
- Generalized Rational Krylov Decompositions with an Application to Rational Approximation
- Hierarchical matrices: algorithms and analysis
- Infinitely Divisible Matrices
- Low-Rank Matrix Approximation Using the Lanczos Bidiagonalization Process with Applications
- Low-rank solvers for fractional differential equations
- Low-Rank Updates and a Divide-And-Conquer Method for Linear Matrix Equations
- Matrix computations and semiseparable matrices. Vol. 1: Linear systems.
- Multigrid method for fractional diffusion equations
- Notes on Hilbert and Cauchy matrices
- On the singular values of matrices with displacement structure
- Preconditioned iterative methods for two-dimensional space-fractional diffusion equations
- Quasi-Toeplitz matrix arithmetic: a MATLAB toolbox
- Separable type representations of matrices and fast algorithms. Volume 2. Eigenvalue method
- Solving rank-structured Sylvester and Lyapunov equations
- Spectral analysis and multigrid preconditioners for two-dimensional space-fractional diffusion equations
- Spectral analysis and structure preserving preconditioners for fractional diffusion equations
- Spectral properties of Toeplitz-plus-Hankel matrices
- Variational formulation for the stationary fractional advection dispersion equation
- Waiting-times and returns in high-frequency financial data: An empirical study
Cited in
(20)- Spectral properties of flipped Toeplitz matrices and related preconditioning
- Space-fractional diffusion with variable order and diffusivity: discretization and direct solution strategies
- All-at-once method for variable-order time fractional diffusion equations
- A fast solver for spectral elements applied to fractional differential equations using hierarchical matrix approximation
- Iterative low-rank approximation solvers for the extension method for fractional diffusion
- Rational Krylov and ADI iteration for infinite size quasi-Toeplitz matrix equations
- A preconditioning technique for all-at-once system from the nonlinear tempered fractional diffusion equation
- Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations
- Rational Krylov for Stieltjes matrix functions: convergence and pole selection
- Low-rank solvers for fractional differential equations
- Nonlocal pagerank
- Divide-and-conquer methods for functions of matrices with banded or hierarchical low-rank structure
- Hm-toolbox: MATLAB software for HODLR and HSS matrices
- Low-rank tensor structure preservation in fractional operators by means of exponential sums
- Numerical study of a fast two-level Strang splitting method for spatial fractional Allen-Cahn equations
- Computing functions of symmetric hierarchically semiseparable matrices
- Approximate parameterized splitting preconditioning for anisotropic space-fractional diffusion equation
- Mixed precision HODLR matrices
- A nested divide-and-conquer method for tensor Sylvester equations with positive definite hierarchically semiseparable coefficients
- Quasi-optimal hierarchically semi-separable matrix approximation
This page was built for publication: Fast solvers for two-dimensional fractional diffusion equations using rank structured matrices
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5230663)