Selecting integration methods¶
The description of an integration method on a whole mesh is done thanks to the
structure getfem::mesh_im
, defined in the file getfem/getfem_mesh_im.h
.
Basically, this structure describes the integration method on each element of the
mesh. One can instantiate a getfem::mesh_im
object as follows:
getfem::mesh_im mim(mymesh);
where mymesh
is an already existing mesh. The structure will be linked to this
mesh and will react when modifications will be done on it (for example when the
mesh is refined, the integration method will be also refined).
It is possible to specify element by element the integration method, so that element of mixed types can be treated, even if the dimensions are different.
To select a particular integration method on a given element, one can use:
mim.set_integration_method(i, ppi);
where i
is the index of the element and ppi
is the descriptor of the
integration method. Alternative forms of this member function are:
void mesh_im::set_integration_method(const dal::bit_vector &cvs,
getfem::pintegration_method ppi);
void mesh_im::set_integration_method(getfem::pintegration_method ppi);
which set the integration method for either the convexes listed in the bit_vector cvs, or all the convexes of the mesh.
The list of all available descriptors of integration methods is in the file
getfem/getfem_integration.h
. Descriptors for integration methods are
available thanks to the following function:
getfem::pintegration_method ppi = getfem::int_method_descriptor("name of method");
where "name of method"
is to be chosen among the existing methods. A name of a
method can be retrieved with:
std::string im_name = getfem::name_of_int_method(ppi);
A non exhaustive list (see Appendix B. Cubature method list or
getfem/getfem_integration.h
for exhaustive lists) of integration methods
is given below.
Examples of exact integration methods:
"IM_NONE()"
: Dummy integration method (new in getfem++-1.7)."IM_EXACT_SIMPLEX(n)"
: Description of the exact integration of polynomials on the simplex of reference of dimensionn
."IM_PRODUCT(a, b)"
: Description of the exact integration on the convex which is the direct product of the convex ina
and inb
."IM_EXACT_PARALLELEPIPED(n)"
: Description of the exact integration of polynomials on the parallelepiped of reference of dimensionn
"IM_EXACT_PRISM(n)"
: Description of the exact integration of polynomials on the prism of reference of dimensionn
Examples of approximated integration methods:
"IM_GAUSS1D(k)"
: Description of the Gauss integration on a segment of orderk
. Available for all odd values ofk <= 99
."IM_NC(n,k)"
: Description of the integration on a simplex of reference of dimensionn
for polynomials of degreek
with the Newton Cotes method (based on Lagrange interpolation)."IM_PRODUCT(a,b)"
: Build a method doing the direct product of methodsa
andb
."IM_TRIANGLE(2)"
: Integration on a triangle of order 2 with 3 points."IM_TRIANGLE(7)"
: Integration on a triangle of order 7 with 13 points."IM_TRIANGLE(19)"
: Integration on a triangle of order 19 with 73 points."IM_QUAD(2)"
: Integration on quadrilaterals of order 2 with 3 points."IM_GAUSS_PARALLELEPIPED(2,3)"
: Integration on quadrilaterals of order 3 with 4 points (shortcut for"IM_PRODUCT(IM_GAUSS1D(3),IM_GAUSS1D(3))"
)."IM_TETRAHEDRON(5)"
: Integration on a tetrahedron of order 5 with 15 points.
Note
Note that "IM_QUAD(3)"
is not able to integrate exactly the base functions
of the "FEM_QK(2,3)"
finite element! Since its base function are tensorial
product of 1D polynomials of degree 3, one would need to use "IM_QUAD(7)"
(6 is not available). Hence "IM_GAUSS_PARALLELEPIPED(2,k)"
should always
be preferred over "IM_QUAD(2*k)"
since it has less integration points.
An alternative way to obtain integration methods:
getfem::pintegration_method ppi =
getfem::classical_exact_im(bgeot::pgeometric_trans pgt);
getfem::pintegration_method ppi =
getfem::classical_approx_im(bgeot::pgeometric_trans pgt, dim_type d);
These functions return an exact (i.e. analytical) integration method, or select an
approximate integration method which is able to integrate exactly polynomials of
degree <= d
(at least) for convexes defined with the specified geometric
transformation.
Methods of the mesh_im object¶
Once an integration method is defined on a mesh, it is possible to obtain information on it with the following methods (the list is not exhaustive).
- mim.convex_index()¶
Set of indexes (a
dal::bit_vector
) on which an integration method is defined.
- mim.linked_mesh()¶
Gives a reference to the linked mesh.
- mim.int_method_of_element(i)¶
Gives a descriptor on the integration method defined on element of index
i
.
- mim.clear()¶
Clear the structure. There are no further integration method defined on the mesh.