inla.mesh.fem: Finite element matrices

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

Description

Constructs finite element matrices for inla.mesh and inla.mesh.1d objects.

Usage

## 2D and 1D meshes
inla.mesh.fem(mesh, order = 2)

## 1D meshes, order 2 models only
inla.mesh.1d.fem(mesh)

Arguments

mesh	An inla.mesh or inla.mesh.1d object.
order	The model order.

Value

A list of sparse matrices based on basis functions psi_i:

c0	c0[i,j] = < psi_i, 1 >
c1	c1[i,j] = < psi_i, psi_j >
g1	g1[i,j] = < grad psi_i, grad psi_j >
g2	g2 = g1 * c0^-1 * g1
gk	gk = g1 * (c0^-1 * g1)^(k-1), up to and including k=order

Author(s)

Finn Lindgren [e-mail]

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