A Sparse Grid Discretization with Variable Coefficient in High Dimensions
Abstract
We present a Ritz-Galerkin discretization on sparse grids using pre-wavelets, which allows to solve elliptic differential equations with variable coefficients for dimension and higher dimensions . The method applies multilinear finite elements. We introduce an efficient algorithm for matrix vector multiplication using a Ritz-Galerkin discretization and semi-orthogonality. This algorithm is based on standard 1-dimensional restrictions and prolongations, a simple pre-wavelet stencil, and the classical operator dependent stencil for multilinear finite elements. Numerical simulation results are presented for a 3-dimensional problem on a curvilinear bounded domain and for a 6-dimensional problem with variable coefficients. Simulation results show a convergence of the discretization according to the approximation properties of the finite element space. The condition number of the stiffness matrix can be bounded below using a standard diagonal preconditioner.
Cite
@article{arxiv.1603.02906,
title = {A Sparse Grid Discretization with Variable Coefficient in High Dimensions},
author = {Rainer Hartmann and Christoph Pflaum},
journal= {arXiv preprint arXiv:1603.02906},
year = {2016}
}