Sets up a Finite Element basis. It requires a
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
quadratic finite elements (when
order = 2 in the input
If saveTree flag is TRUE, it saves the tree mesh information in advance inside mesh object and can be used later on to save mesh construction time.
create.FEM.basis(mesh, saveTree = FALSE)
a flag to decide to save the tree mesh information in advance (default is FALSE)
FEMbasis object. This contains the
mesh, along with some additional quantities:
orderEither "1" or "2" for the 2D and 2.5D case, and "1" for the 3D case.
Order of the Finite Element basis.
nbasisScalar. The number of basis.
library(fdaPDE) ## Upload the quasicircle2D data data(quasicircle2D) ## Create the 2D mesh mesh = create.mesh.2D(nodes = rbind(quasicircle2D$boundary_nodes, quasicircle2D$locations), segments = quasicircle2D$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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.