Changeset 4898 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

07/30/10 09:00:31 (15 years ago)

Author:

Mathieu Morlighem

Message:

Use Matrix utils to compute determinant and matrix invert. Simple routines for 2x2 and 3x3 matrices have been added (speeds up)

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (3 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -4014 +4014 @@
 }
 /*}}}*/
-/*FUNCTION Penta::GetMatrixInvert {{{1*/
-void Penta::GetMatrixInvert(double*  Ke_invert, double* Ke){
-        /*Inverse a 3 by 3 matrix only */
-
-        double a,b,c,d,e,f,g,h,i;
-        double det;
-        int calculationdof=3;
-
-        /*Take the matrix components: */
-        a=*(Ke+calculationdof*0+0);
-        b=*(Ke+calculationdof*0+1);
-        c=*(Ke+calculationdof*0+2);
-        d=*(Ke+calculationdof*1+0);
-        e=*(Ke+calculationdof*1+1);
-        f=*(Ke+calculationdof*1+2);
-        g=*(Ke+calculationdof*2+0);
-        h=*(Ke+calculationdof*2+1);
-        i=*(Ke+calculationdof*2+2);
-
-        det=a*(e*i-f*h)-b*(d*i-f*g)+c*(d*h-e*g);
-
-        *(Ke_invert+calculationdof*0+0)=(e*i-f*h)/det;
-        *(Ke_invert+calculationdof*0+1)=(c*h-b*i)/det;
-        *(Ke_invert+calculationdof*0+2)=(b*f-c*e)/det;
-        *(Ke_invert+calculationdof*1+0)=(f*g-d*i)/det;
-        *(Ke_invert+calculationdof*1+1)=(a*i-c*g)/det;
-        *(Ke_invert+calculationdof*1+2)=(c*d-a*f)/det;
-        *(Ke_invert+calculationdof*2+0)=(d*h-e*g)/det;
-        *(Ke_invert+calculationdof*2+1)=(b*g-a*h)/det;
-        *(Ke_invert+calculationdof*2+2)=(a*e-b*d)/det;
-
-}
-/*}}}*/
 /*FUNCTION Penta::GetParameterValue(double* pvalue, double* v_list,double* gauss_coord) {{{1*/
 void Penta::GetParameterValue(double* pvalue, double* v_list,double* gauss_coord){
@@ -5186 +5153 @@
 
         /*Inverse the matrix corresponding to bubble part Kbb */
-        GetMatrixInvert(&Kbbinv[0][0], &Kbb[0][0]);
+        Matrix3x3Invert(&Kbbinv[0][0], &Kbb[0][0]);
 
         /*Multiply matrices to create the reduce matrix Ke_reduced */
@@ -5230 +5197 @@
 
         /*Inverse the matrix corresponding to bubble part Kbb */
-        GetMatrixInvert(&Kbbinv[0][0], &Kbb[0][0]);
+        Matrix3x3Invert(&Kbbinv[0][0], &Kbb[0][0]);
 
         /*Multiply matrices to create the reduce matrix Ke_reduced */