A multidomain spectral method for solving elliptic equations

From MaRDI portal
Publication:709381

DOI10.1016/S0010-4655(02)00847-0zbMATH Open1196.65179arXivgr-qc/0202096OpenAlexW2089021822MaRDI QIDQ709381FDOQ709381

Harald P. Pfeiffer, Mark A. Scheel, Saul A. Teukolsky, Lawrence E. Kidder

Publication date: 18 October 2010

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three distinct features. First, the combined problem of solving the PDE, satisfying the boundary conditions, and matching between different subdomains is cast into one set of equations readily accessible to standard linear and nonlinear solvers. Second, touching as well as overlapping subdomains are supported; both rectangular blocks with Chebyshev basis functions as well as spherical shells with an expansion in spherical harmonics are implemented. Third, the code is very flexible: The domain decomposition as well as the distribution of collocation points in each domain can be chosen at run time, and the solver is easily adaptable to new PDEs. The code has been used to solve the equations of the initial value problem of general relativity and should be useful in many other problems. We compare the new method to finite difference codes and find it superior in both runtime and accuracy, at least for the smooth problems considered here.


Full work available at URL: https://arxiv.org/abs/gr-qc/0202096




Recommendations



Cites Work


Cited In (51)

Uses Software





This page was built for publication: A multidomain spectral method for solving elliptic equations

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