Methods for projecting to/from an inla.mesh

Description

since ⁠23.08.18⁠. Use fmesher::fm_evaluate() and fmesher::fm_evaluator() instead.

Calculate a lattice projection to/from an inla.mesh().

The call inla.mesh.project(mesh, loc, field=..., ...), is a shortcut to inla.mesh.project(inla.mesh.projector(mesh, loc), field).

Usage

inla.mesh.project(...)

inla.mesh.projector(...)

Arguments

Value

For inla.mesh.project(mesh, ...), a list with projection information. For inla.mesh.projector(mesh, ...), an inla.mesh.projector object. For inla.mesh.project(projector, field, ...), a field projected from the mesh onto the locations given by the projector object.

Author(s)

Finn Lindgren finn.lindgren@gmail.com

See Also

inla.mesh(), inla.mesh.1d(), inla.mesh.lattice()

Examples

n <- 20
loc <- matrix(runif(n * 2), n, 2)
mesh <- inla.mesh.create(loc, refine = list(max.edge = 0.05))
proj <- inla.mesh.projector(mesh)
field <- cos(mesh$loc[, 1] * 2 * pi * 3) * sin(mesh$loc[, 2] * 2 * pi * 7)
image(proj$x, proj$y, inla.mesh.project(proj, field))

if (require(rgl)) {
  plot(mesh, rgl = TRUE, col = field, draw.edges = FALSE, draw.vertices = FALSE)
}