Getfem++ User Guide

• Introduction
• How to install
• Build a mesh
• • Add an element to a mesh
  • Remove an element from a mesh
  • Simple structured meshes
  • Mesh regions
  • Methods of the getfem::mesh object
  • Using dal::bit_vector
  • Face numbering
  • Save and load meshes
  • • From getfem file format
    • Import a mesh
• Selecting finite element methods
• • Examples
  • Other methods of the mesh_fem object
  • Obtaining generic mesh_fems
• Selecting integration methods
• • Methods of the mesh_im object
• Mesh refinement
• Linear algebra procedures
• Standard assembly procedures
• • Laplacian (Poisson) problem
  • Linear Elasticity problem
  • Stokes Problem with mixed finite element method
  • Assembling a mass matrix
• Compute arbitrary elementary matrices - generic assembly procedures
• Incorporate new finite element methods in getfem++ 
• Incorporate new approximated integration methods in getfem++ 
• Level-sets, Xfem, fictitious domains
• • Representation of level-sets
  • Mesh cut by level-sets
  • Adapted integration methods
  • Discontinuous field across some level-sets
  • Fictitious domain approach with Xfem
• Support for Xfem methods
• Interpolation of a finite element method on non-matching meshes
• • mixed methods with different meshes
  • mortar methods
• Compute L² and H¹ norms
• Compute derivatives
• Export and view a solution
• • Saving mesh and mesh_fem objects for the Matlab interface
  • Producing mesh slices
  • Exporting mesh, mesh_fem or slices to VTK
  • Exporting mesh, mesh_fem or slices to OpenDX
• Interpolation on different meshes
• The model bricks
• • The model state variable
  • Basic properties of a brick
  • Brick parameters
  • generic elliptic brick
  • Source term brick
  • Constraint brick
  • Dirichlet condition brick
  • Example of a complete Poisson problem
  • Predefined solvers
  • Isotropic linearized elasticity brick
  • Qu term brick
  • Helmholtz brick
  • linear incompressibility (or nearly incompressibility) brick
  • Small displacement plasticity brick
  • Contact and friction conditions brick
  • Linearized plate brick
  • Large strain elasticity brick
  • Large strain incompressibility brick
• Parallelization of getfem++ 
• Catch errors
• Example : Laplacian program
• Appendix A. Finite element method list
• • Dof graphical codification
  • Classical P_K Lagrange elements on simplices
  • Classical Lagrange elements on other geometries
  • Elements with hierarchical basis
  • • Hiercarchical elements with respect to the degree
    • Composite elements
    • Hierarchical composite elements
  • Classical vectorial elements
  • • Raviart-Thomas of lowest order elements
    • Nedelec (or Whitney) edge elements
  • Specific elements in dimension 1
  • • GaussLobatto element
    • Hermite element
    • Lagrange element with an additional bubble function
  • Specific elements in dimension 2
  • • Elements with additional bubble functions
    • Non-conforming P₁ element
    • Hermite element
    • Morley element
    • Argyris element
    • Hsieh-Clough-Tocher element
    • A composite C¹ element on quadrilaterals
  • Specific elements in dimension 3
  • • Elements with additional bubble functions
    • Hermite element
• Appendix B. Cubature method list
• • Exact Integration methods
  • Newton cotes Integration methods
  • Gauss Integration methods on dimension 1
  • Gauss Integration methods on dimension 2
  • Gauss Integration methods on dimension 3
  • Direct product of integration methods
  • Composite integration methods
• References
• Index

Incorporate new approximated integration methods in getfem++ 

A perl script automatically incorporates new cubature methods from a description file. You can see in the directory cubature such description files (with extension .IM) . For instance for IM_TETRAHEDRON(5) the following file describes the method:

NAME = IM_TETRAHEDRON(5)
N = 3
GEOTRANS = GT_PK(3,1)
NBPT = 4
0, 0.25, 0.25, 0.25, 0.008818342151675485
1, 0.31979362782962991, 0.31979362782962991, 0.31979362782962991, 0.011511367871045398 
1, 0.091971078052723033, 0.091971078052723033, 0.091971078052723033, 0.01198951396316977
1, 0.056350832689629156, 0.056350832689629156, 0.44364916731037084, 0.008818342151675485
NBF = 4
IM_TRIANGLE(5)
IM_TRIANGLE(5)
IM_TRIANGLE(5)
IM_TRIANGLE(5)

where NAME is the name of the method in getfem++ (constant integer parameter are allowed), N is the dimension, GEOTRANS describes a valid geometric transformation of getfem++ . This geometric transformation just defines the reference element on which the integration method is described. NBPT is the number of integration node definitions. Integration node definitions include a symmetry definition such that the total number of integration nodes would be greater than NBPT.

Composition of the integration node definition :

• an integer : 0 = no symmetry, 1 = full symmetric (x6 for a triangle, x4 for a quadrangle, x24 for a tetrahedron ...),
• the N coordinates of the integration node,
• the load.

NBF is the number of faces of the reference element (should correspond to GEOTRANS). Then follows an already existing integration method for each face (each on a line). This is necessary to make integrations on boundaries.

The file format is inspired from [3].

---