Loading paper
Enforce the Dirichlet boundary condition by volume constraint in Point Integral method | Tomesphere