fm_fem: Compute finite element matrices

View source: R/fem.R

fm_femR Documentation

Compute finite element matrices

Description

(...)

Usage

fm_fem(mesh, order = 2, ...)

## S3 method for class 'fm_mesh_1d'
fm_fem(mesh, order = 2, ...)

## S3 method for class 'fm_mesh_2d'
fm_fem(mesh, order = 2, aniso = NULL, ...)

## S3 method for class 'inla.mesh.1d'
fm_fem(mesh, order = 2, ...)

## S3 method for class 'inla.mesh'
fm_fem(mesh, order = 2, ...)

Arguments

mesh

inla.mesh object

order

integer

...

Currently unused

aniso

If non-NULL, a list(gamma, v). Calculates anisotropic structure matrices (in addition to the regular) for \gamma and v for an anisotropic operator \nabla\cdot H \nabla, where H=\gamma I + v v^\top. Currently (2023-08-05) the fields need to be given per vertex.

Value

fm_fem.fm_mesh_1d: A list with elements c0, c1, g1, g2. When mesh$degree == 2, also g01, g02, and g12.

fm_fem.fm_mesh_2d: A list with elements c0, c1, g1, va, ta, and more if order > 1. When aniso is non-NULL, also g1aniso matrices, etc.

Examples

str(fm_fem(fmexample$mesh))


fmesher documentation built on Nov. 2, 2023, 5:35 p.m.