Basic robust integration weights for mesh/polygon intersections
make_stable_integration_points(mesh, bnd, nsub = NULL)
Mesh on which to integrate
inla.mesh.segment
defining the integration domain
number of subdivision points along each triangle edge, giving
(nsub + 1)^2
proto-integration points used to compute
the vertex weights
(default NULL=9
, giving 100 integration points for each triangle)
list
with elements loc
and weight
with
integration points for the intersection of the mesh and polygon