Re: [Metil-lmtpp] gauss integration

[Martin G. <mar...@lm...>]

Un point sur l'intégration:

Si dans formulation.py on trouve integration_totale = True (comportement par défaut), alors le code c++ généré lorsque e.integration(expression) est appelé correspond à l'intégration analytique de expression sur l'élément.

Si dans formulation.py on trouve integration_totale = False, alors le code c++ généré lorsque e.integration(expression) est appelé correspond à l'intégration de Gauss de expression sur l'élément. Par défaut, le degré de expression est calculé automatique en fonction de l'élément (2 pour les éléments linéaires, 4 pour les éléments quadratiques, etc). On peut également spécifier le degré en appelant e.integration(expression, degré), par exemple pour faire de la sous-intégration.

>  - Pourquoi, dans certains fichiers element.py —e.g. element_Tetra, element_Wedge, element_Hexa, ...—, certains jeux de points de Gauss —e.g. gauss_points[2] dans element_Hexa— ont-ils été commentés?

Certains jeux de points de Gauss ont été décommentés (e.g. element_Tetra.py, element_Wedge.py, element_Hexa.py, ...).

>  - Serait-ce envisageable de créer dans ce cas des jeux qui pointent vers leur supérieur —e.g. gauss_points[2] = gauss_points[3]—?

Certains jeux de points de Gauss ont été crées qui pointent vers leur supérieur (e.g. element_Tetra.py, element_Wedge.py, element_Hexa.py, ...).

Martin.