{"id":"arxiv/2604.03097","url":"https://arxiv.org/abs/2604.03097","tags":["math.NA"],"score":0,"title":"A High-Order Fast Direct Solver for Surface PDEs on Triangles","domain":"math.NA","source":"arxiv","authors":["Gentian Zavalani"],"reviews":0,"subject":"math.NA","summary":"We develop a triangular formulation of the hierarchical Poincaré-Steklov (HPS) method for elliptic partial differential equations on surfaces, allowing high-order discretizations on unstructured meshes and complex geometries. Classical HPS formulations rely on high-order quadrilateral meshes and tensor-product spectral discretizations, which enable efficient algorithms but restrict applicability to structured geometries. To overcome this restriction, we introduce a triangle-based hierarchical Poincaré-Steklov scheme (THPS) built on orthogonal Dubiner polynomial bases. As in the classical HPS…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-04-03","reactions":{}}