Numerical approximation of fractional powers of elliptic operators

From MaRDI portal
(Redirected from Publication:5264121)
Numerical approximation of fractional powers of elliptic operators (scientific article; zbMATH DE number 6460482)



Abstract: We present and study a novel numerical algorithm to approximate the action of where L is a symmetric and positive definite unbounded operator on a Hilbert space H0. The numerical method is based on a representation formula for in terms of Bochner integrals involving (I+t2L)−1 for tin(0,infty). To develop an approximation to , we introduce a finite element approximation Lh to L and base our approximation to on . The direct evaluation of is extremely expensive as it involves expansion in the basis of eigenfunctions for Lh. The above mentioned representation formula holds for and we propose three quadrature approximations denoted generically by . The two results of this paper bound the errors in the H0 inner product of and where pih is the H0 orthogonal projection into the finite element space. We note that the evaluation of involves application of (I+(ti)2Lh)−1 with ti being either a quadrature point or its inverse. Efficient solution algorithms for these problems are available and the problems at different quadrature points can be straightforwardly solved in parallel. Numerical experiments illustrating the theoretical estimates are provided for both the quadrature error and the finite element error .




Cites work


Cited in
(only showing first 100 items - show all)








This page was built for publication: Numerical approximation of fractional powers of elliptic operators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5264121)