Plot a FEM object

Share:

Description

Three-dimensional plot of a FEM object, generated by FEM or returned by smooth.FEM.basis, smooth.FEM.PDE.basis or smooth.FEM.PDE.sv.basis.

Usage

1
2
## S3 method for class 'FEM'
plot(x, num_refinements, ...)

Arguments

x

A FEM object.

num_refinements

A natural number specifying how many bisections should by applied to each triangular element for plotting purposes. This functionality is useful where a discretization with 2nd order Finite Element is applied.

...

Arguments representing graphical options to be passed to plot3d.

See Also

image.FEM

Examples

1
2
3
4
5
6
7
8
9
## Upload a triangular mesh and plot it
data("mesh.2D.rectangular")
plot(mesh.2D.rectangular)
## Create a linear Finite Element basis
FEMbasis = create.FEM.basis(mesh.2D.rectangular)
## Define a sinusoidal function as expansion of this basis and plot it
coeff <- sin(mesh.2D.rectangular$nodes[,1])*cos(mesh.2D.rectangular$nodes[,2])
FEM_object<- FEM(coeff, FEMbasis)
plot(FEM_object)