  Many geometry processing techniques require the solution of partial
differential equations (PDEs) on surfaces. Such surface PDEs often involve
boundary conditions prescribed on the surface, at points or curves on its
interior or along the geometric (exterior) boundary of an open surface.
However, input surfaces can take many forms (e.g., meshes, parametric surfaces,
point clouds, level sets, neural implicits). One must therefore generate a mesh
to apply finite element-type techniques or derive specialized discretization
procedures for each surface representation.
  We propose instead to address such problems through a novel extension of the
closest point method (CPM) to handle interior boundary conditions specified at
surface points or curves. CPM solves the surface PDE by solving a volumetric
PDE defined over the Cartesian embedding space containing the surface; only a
closest point function is required to represent the surface. As such, CPM
supports surfaces that are open or closed, orientable or not, and of any
codimension or even mixed-codimension. To enable support for interior boundary
conditions, we develop a method to implicitly partition the embedding space
across interior boundaries. CPM's finite difference and interpolation stencils
are adapted to respect this partition while preserving second-order accuracy.
Furthermore, an efficient sparse-grid implementation and numerical solver is
developed that can scale to tens of millions of degrees of freedom, allowing
PDEs to be solved on more complex surfaces. We demonstrate our method's
convergence behaviour on selected model PDEs. Several geometry processing
problems are explored: diffusion curves on surfaces, geodesic distance, tangent
vector field design, and harmonic map construction. Our proposed approach thus
offers a powerful and flexible new tool for a range of geometry processing
tasks on general surface representations.
