surfpy is a Python package for computing surface integrals over smooth embedded manifolds.

Surfpy rests on curved surface triangulations realised due to kth-order interpolation of the closest point projection, extending initial linear surface approximations. It achieves this by employing a novel technique called square-squeezing, which involves transforming the interpolation tasks of triangulated manifolds to the standard hypercube using a cube-to-simplex transformation that has been recently introduced.