Changeset 4898 for issm/trunk/src/c/shared/Matrix/matrix.h

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/shared/Matrix/matrix.h (modified) (1 diff)

issm/trunk/src/c/shared/Matrix/matrix.h

@@ -4 +4 @@
 
 #ifndef _MATRIXUTILS_H_
-#define  _MATRIXUTILS_H_
+#define _MATRIXUTILS_H_
 
-int TripleMultiply( double* a, int nrowa, int ncola, int itrna, double* b, int nrowb, int ncolb, int itrnb, double* c, int nrowc, int ncolc, int itrnc, double* d, int iaddd);
-int MatrixMultiply( double* a, int nrowa, int ncola, int itrna, double* b, int nrowb, int ncolb, int itrnb, double* c, int iaddc );
-int MatrixInverse( double* a, int ndim, int nrow, double* b, int nvec, double* pdet );
-        
+int  TripleMultiply( double* a, int nrowa, int ncola, int itrna, double* b, int nrowb, int ncolb, int itrnb, double* c, int nrowc, int ncolc, int itrnc, double* d, int iaddd);
+int  MatrixMultiply( double* a, int nrowa, int ncola, int itrna, double* b, int nrowb, int ncolb, int itrnb, double* c, int iaddc );
+int  MatrixInverse( double* a, int ndim, int nrow, double* b, int nvec, double* pdet );
+void Matrix2x2Invert(double* Ainv, double* A);
+void Matrix2x2Determinant(double* Adet,double* A);
+void Matrix3x3Invert(double* Ainv, double* A);
+void Matrix3x3Determinant(double* Adet,double* A);
 
 #endif //ifndef _MATRIXUTILS_H_