inla.mesh.2d: High-quality triangulations

In andrewzm/INLA: Functions which allow to perform full Bayesian analysis of latent Gaussian models using Integrated Nested Laplace Approximaxion

Description

Create a triangle mesh based on initial point locations, specified or automatic boundaries, and mesh quality parameters.

Usage

inla.mesh.2d(loc = NULL,
             loc.domain = NULL,
             offset = NULL,
             n = NULL,
             boundary = NULL,
             interior = NULL,
             max.edge,
             min.angle = NULL,
             cutoff = 1e-12,
             plot.delay = NULL)

Arguments

loc	Matrix of point locations to be used as initial triangulation nodes.
loc.domain	Matrix of point locations used to determine the domain extent.
offset	The automatic extension distance. One or two values, for an inner and an optional outer extension. If negative, interpreted as a factor relative to the approximate data diameter (default=-0.10???)
n	The number of initial nodes in the automatic extensions (default=16)
boundary	A list of one or two inla.mesh.segment objects describing domain boundaries.
interior	An inla.mesh.segment object describing desired interior edges.
max.edge	The largest allowed triangle edge length. One or two values.
min.angle	The smallest allowed triangle angle. One or two values. (Default=21)
cutoff	The minimum allowed distance between points. Point at most as far apart as this are replaced by a single vertex prior to the mesh refinement step.
plot.delay	On Linux (and Mac if appropriate X11 libraries are installed), specifying a nonnegative numeric value activates a rudimentary plotting system in the underlying fmesher program, showing the triangulation algorithm at work, with waiting time factor plot.delay between each step. On all systems, specifying any negative value activates displaying the result after each step of the multi-step domain extension algorithm.

Value

An inla.mesh object.

Author(s)

Finn Lindgren [e-mail]

See Also

inla.mesh.create, inla.delaunay

Examples

loc = matrix(runif(10*2),10,2)

mesh = inla.mesh.2d(loc,
                    boundary=list(inla.nonconvex.hull(loc, 0.1, 0.15),
                                  inla.nonconvex.hull(loc, 0.2, 0.2)),
                    max.edge=c(0.05, 0.1))
plot(mesh)

andrewzm/INLA documentation built on May 10, 2019, 11:12 a.m.