Python GetFEM interface¶
Introduction¶
GetFEM provides an interface to the Python scripting language. Python is a nice, crossplatform, and free language. With the addition of the numpy package, python provides a subset of Matlab functionalities (i.e. dense arrays). The VTK toolkit may provide visualization tools via its python interface (or via MayaVi), and data files for OpenDX may be exported. In this guide, nevertheless, to visualize the results, we will export to Gmsh postprocessing format. The sparse matrix routines are provided by the getfem interface.
The python interface is available via a python module getfem.py. In order to use the interface you have to load it with:
import getfem
m = getfem.Mesh('cartesian', range(0, 3), range(0,3))
or:
from getfem import *
m = Mesh('cartesian', range(0, 3), range(0,3))
If the getfem.py (and the internal _getfem.so) module is not installed in a
standard location for python, you may have to set the PYTHONPATH
environment variable to its location. For example with:
import sys
sys.path.append('.../getfem/getfem++/interface/src/python/')
Parallel version¶
The python interface is the only one for the moment to interface the mpi based parallel version of Getfem. See MPI Parallelization of GetFEM.
Memory Management¶
A nice advantage over the Matlab interface is that you do not have to
explicitly delete objects that are not used any more, this is done
automagically. You can however inspect the content of the getfem workspace
with the function getfem.memstats()
.
Documentation¶
The getfem module is largely documented. This documentation has been extracted into the API reference. The getfemmatlab user guide may also be used, as 95% of its content translates quite directly into python (with the exception of the plotting functions, which are specific to matlab).
Python GetFEM organization¶
The general organization of the pythoninterface is the following:
Each class from the matlab interface has a corresponding class in the python interface: the gfMesh class becomes the getfem.Mesh class in python, the gfSlice becomes the getfem.Slice etc.
Each get and set method of the matlab interface has been translated into a method of the corresponding class in the python interface. For example:
gf_mesh_get(m, 'outer faces'); gf_mesh_get(m, 'pts');becomes:
m.outer_faces(); m.pts();Some methods have been renamed when there was ambiguity, for example
gf_mesh_set(m, 'pts', P)
ism.set_pts(P)
.The other getfemmatlab function have a very simple mapping to their python equivalent:
gf_compute(mf,U,’foo’,…) getfem.compute_foo(mf,U) or getfem.compute(‘foo’,…) gf_asm(‘foobar’,…) getfem.asm_foobar(…) or getfem.asm(‘foobar’,…) gf_linsolve(‘gmres’,…) getfem.linsolve_gmres(…) or getfem.linsolve(‘gmres’,…)

class
CvStruct
(self, *args)¶ Descriptor for a convex structure objects, stores formal information convex structures (nb. of points, nb. of faces which are themselves convex structures)

class
GeoTrans
(self, *args)¶ Descriptor for geometric transformations objects (defines the shape/position of the convexes).

class
Mesh
(self, *args)¶ Descriptor for mesh structure (nodes, convexes, geometric transformations for each convex).

class
Fem
(self, fem_name)¶ Descriptor for FEM (Finite Element Method) objects (one per convex, can be PK, QK, HERMITE, etc…).

class
Integ
(self, *args)¶ Descriptor for Integration Method objects (exact, quadrature formulaldots). Although not linked directly to GeoTrans, an integration method is usually specific to a given convex structure.

class
MeshFem
(self, *args)¶ Descriptor for object linked to a mesh, where each convex has been assigned an FEM.

class
MeshIm
(self, *args)¶ Descriptor for object linked to a mesh, where each convex has been assigned an integration method.

class
Model
(self, *args)¶ Descriptor for model object, holds the global data, variables and description of a model. Evolution of model state and model brick object for 4.0 version of GetFEM.