Closed finnlindgren closed 6 years ago
Since the weights are constructed based on the existing mesh points, the construction should use the mesh boundary to limit the integration domain to the intersection between the mesh and samplers
.
The mesh/polygon intersection code was non-robust and has been replaced by an more approximate but highly robust construction. Fixed in 98337a2
In
int.polygon
, the integration weight construction may end up constructing small negative integration weights, which are then used as exposure input to INLA:Causes may include degenerate triangles, and/or sampling domains that extend outside the domain mesh.