14 Novembre 2005	By LPX

neant

neant

GPState :

- Création d'une classe (SfGpState) intermédiaire à GpStateBase et MechanicalGpState ou Thermal GpState contenant

- la matrice Jacobienne (jacobianMatrix) (ATTENTION : selon le contexte, il s'agira de la matrice (GP1) ou de son inverse (GP0)!!!)

- le determinant de celle ci (detJ)

- Suppression des tableaux J1, J0Inv, J1t, detJ, detJt, .... au profit du gpState

Optimisation CPU:

FourthOrderTensor

- FourthOrderTensor3D::operator():

le test sur les i,j,k,l empechait l'optimisation de la fonction (sous pc) !!!

inline double
FourthOrderTensor3D::operator() (int i, int j, int k, int l) const
{
if(i>=0 && i<=2 && j>=0 && j<=2 && k>=0 && k<=2 && l>=0 && l<=2)
{
return x[i][j][k][l];
}
else
{
FATAL_ERROR("");
return 0.;
}
}

résultat Quantify : (routine appelante)

MechanicalVolumeIntegrationMethods<Dimension2D>::getElementaryStiffness
temps fct + descendants = 385845.12

Remplace par :

inline double
FourthOrderTensor3D::operator() (int i, int j, int k, int l) const
{
assert(i>=0 && i<=2 && j>=0 && j<=2 && k>=0 && k<=2 && l>=0 && l<=2);
return x[i][j][k][l];

}

résultat Quantify : (routine appelante)

MechanicalVolumeIntegrationMethods<Dimension2D>::getElementaryStiffness
temps fct + descendants = 55987.20

 

- FourthOrderTensor2D/3D::Q1SQ2 :

Inlining de ces fonctions (qui sont au "top" des fonctions consommant du CPU)

- FourthOrderTensor3D::setToZero
- FourthOrderTensor3D:: operator *=

Explicitation des fonctions (à la place de la quadruple boucle) : gain important d'apres quantify (à vérifier sur batterie)

EAS :

- Remplacement de qq copies de Matr3 en récuperation de leur reference qd cela est possible

- Inlining de fonction dyadic (la fonction Matr3.dyadic(SymMatr3) est "inlinable" alors que SymMatr3.dyadic(Matr3) ne l'est pas)

Gain total observé :

(sachant que le transfert de jaco/detJ coute un peu)

NewFiles :

\oo_meta\mtMaterials\SfGpState.h/cpp

	created :14 Novembre 2005		modified : 14 Novembre 2005
	contact :L.Papeleux@ulg.ac.be