Re: [Metil-lmtpp] gauss integration
Status: Pre-Alpha
Brought to you by:
hugo_lec
From: Martin G. <mar...@lm...> - 2009-12-16 17:25:57
|
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. |