Using hierarchical matrices in the solution of the time-fractional heat equation by multigrid waveform relaxation release_y4pdytyisje3fo5qalihgfbl3i

by Xiaozhe Hu, Carmen Rodrigo, Francisco J. Gaspar

Released as a article .

2020  

Abstract

This work deals with the efficient numerical solution of the time-fractional heat equation discretized on non-uniform temporal meshes. Non-uniform grids are essential to capture the singularities of "typical" solutions of time-fractional problems. We propose an efficient space-time multigrid method based on the waveform relaxation technique, which accounts for the nonlocal character of the fractional differential operator. To maintain an optimal complexity, which can be obtained for the case of uniform grids, we approximate the coefficient matrix corresponding to the temporal discretization by its hierarchical matrix (H-matrix) representation. In particular, the proposed method has a computational cost of O(k N M log(M)), where M is the number of time steps, N is the number of spatial grid points, and k is a parameter which controls the accuracy of the H-matrix approximation. The efficiency and the good convergence of the algorithm, which can be theoretically justified by a semi-algebraic mode analysis, are demonstrated through numerical experiments in both one- and two-dimensional spaces.
In text/plain format

Archived Files and Locations

application/pdf  1.2 MB
file_2qfnpi5pwrcgdn27vt32lj5sla
arxiv.org (repository)
web.archive.org (webarchive)
Read Archived PDF
Preserved and Accessible
Type  article
Stage   submitted
Date   2020-05-07
Version   v3
Language   en ?
arXiv  1706.07632v3
Work Entity
access all versions, variants, and formats of this works (eg, pre-prints)
Catalog Record
Revision: b3548c71-624c-41d6-b615-b44666380181
API URL: JSON