Non-planar/spherical mesh object construction fails

Noticed by FM:
Constructing a non-R2/S2 mesh object with inla.mesh.create because of an 
oversight in the logic dealing with pre-specified triangulations.  fmesher 
cannot handle refining a non-R2/S2 triangulation, and bails out without 
checking that it wasn't asked to refine.

> pts <- matrix(c(0,3,2,1,0,0,1,2,0,0,1,0), ncol=3)
> mtest <- inla.mesh.create(loc=pts, tv=matrix(as.integer(c(2,3,3,4,1,1)), 
ncol=3))
fmesher.cc(830) Points not in the plane or on a sphere, and --smorg not 
specified.
No output generated.
Error in fmesher.read(prefix, "manifold") : 
  File '/tmp/RtmptdqLa9/fmesher47a55a2a.manifold' does not exist.

A temporary workaround is to construct a mesh object with the desired triangles 
but planar coordinates, and then replace the vertex coordinates with the 
intended ones:

> mtest <- 
inla.mesh.create(loc=cbind(matrix(runif(nrow(pts)*2),nrow(pts),2),0), 
tv=matrix(as.integer(c(2,3,3,4,1,1)), ncol=3))
> mtest$loc=pts

Original issue reported on code.google.com by finn.lin...@gmail.com on 16 Dec 2011 at 1:01

kathryntmorrison / inla

Non-planar/spherical mesh object construction fails #17