Finite element matrices — inla.mesh.fem • INLA

since 23.08.18. Use fmesher::fm_fem() instead.

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

inla.mesh.fem(mesh, order = 2)

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

See also

inla.mesh.1d.fem()

Author

Finn Lindgren finn.lindgren@gmail.com