A robust solver for elliptic PDEs in 3D complex geometries

Matthew J. Morse, Abtin Rahimian, Denis Zorin

Research output: Contribution to journalArticlepeer-review


We develop a boundary integral equation solver for elliptic partial differential equations on complex 3D geometries. Our method is efficient, high-order accurate and robustly handles complex geometries. A key component is our singular and near-singular layer potential evaluation scheme, hedgehog: a simple extrapolation of the solution along a line to the boundary. We present a series of geometry-processing algorithms required for hedgehog to run efficiently with accuracy guarantees on arbitrary geometries and an adaptive upsampling scheme based on a iteration-free heuristic for quadrature error. We validate the accuracy and performance with a series of numerical tests and compare our approach to a competing local evaluation method.

Original languageEnglish (US)
Article number110511
JournalJournal of Computational Physics
StatePublished - Oct 1 2021

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'A robust solver for elliptic PDEs in 3D complex geometries'. Together they form a unique fingerprint.

Cite this