create.FEM.basis: Create a FEM basis

Description Usage Arguments Value See Also Examples

View source: R/FEMobjects.R

Description

Sets up a Finite Element basis. It requires a mesh.2D or mesh.2.5D object, as input. The basis' functions are globally continuos functions, that are polynomials once restricted to a triangle in the mesh. The current implementation includes linear finite elements (when order = 1 in the input mesh) and quadratic finite elements (when order = 2 in the input mesh).

Usage

1

Arguments

mesh

A mesh.2D or mesh.2.5D object representing the domain triangulation. See create.mesh.2D, create.mesh.2.5D.

Value

A FEMbasis object. This contains the mesh, along with some additional quantities:

See Also

create.mesh.2D, create.mesh.2.5D

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
## Upload the quasicircle2D data
data(quasicircle2D)
boundary_nodes = quasicircle2D$boundary_nodes
boundary_segments = quasicircle2D$boundary_segments
locations = quasicircle2D$locations
data = quasicircle2D$data

## Create the 2D mesh
mesh = create.mesh.2D(nodes = rbind(boundary_nodes, locations), segments = boundary_segments)
## Plot it
plot(mesh)                   
## Create the basis
FEMbasis = create.FEM.basis(mesh)
## Upload the hub2.5D data
data(hub2.5D)
hub2.5D.nodes = hub2.5D$hub2.5D.nodes
hub2.5D.triangles = hub2.5D$hub2.5D.triangles

## Create the 2.5D mesh
mesh = create.mesh.2.5D(nodes = hub2.5D.nodes, triangles = hub2.5D.triangles)
## Plot it
plot(mesh)                   
## Create the basis
FEMbasis = create.FEM.basis(mesh)

fdaPDE documentation built on July 2, 2020, 2:22 a.m.