A scalable algorithm for solving linear parabolic evolution equations

by   Raymond van Venetië, et al.

We present an algorithm for the solution of a simultaneous space-time discretization of linear parabolic evolution equations with a symmetric differential operator in space. Building on earlier work, we recast this discretization into a Schur-complement equation whose solution is a quasi-optimal approximation to the weak solution of the equation at hand. Choosing a tensor-product discretization, we arrive at a remarkably simple linear system. Using wavelets in time and standard finite elements in space, we solve the resulting system in optimal linear complexity on a single processor, and in optimal logarithmic complexity when parallelized in both space and time. We complement these theoretical findings with large-scale parallel computations showing the effectiveness of the method.


page 1

page 2

page 3

page 4


A space-time isogeometric method for the partial differential-algebraic system of Biot's poroelasticity model

Biot's equations of poroelasticity contain a parabolic system for the ev...

Preconditioning Richards Equations: spectral analysis and parallel solution at very large scale

We consider here a cell-centered finite difference approximation of the ...

Minimal residual space-time discretizations of parabolic equations: Asymmetric spatial operators

We consider a minimal residual discretization of a simultaneous space-ti...

Solving differential Riccati equations: A nonlinear space-time method using tensor trains

Differential algebraic Riccati equations are at the heart of many applic...

Approximate Solutions to Second-Order Parabolic Equations: evolution systems and discretization

We study the discretization of a linear evolution partial differential e...

Solving the linear semiclassical Schrödinger equation on the real line

The numerical solution of a linear Schrödinger equation in the semiclass...