Changeset 5371

Timestamp:

08/18/10 12:03:05 (15 years ago)

Author:

seroussi

Message:

added transpose function for matrices

Location:

issm/trunk/src/c/shared/Matrix

Files:

• MatrixUtils.cpp (modified) (1 diff)
• matrix.h (modified) (1 diff)

issm/trunk/src/c/shared/Matrix/MatrixUtils.cpp

@@ -399 +399 @@
 
 }/*}}}*/
+/*FUNCTION MatrixTranspose(double* Adet,double* A) {{{1*/
+void MatrixTranspose(double* tA,double* A, int nrows, int ncols){
+        /*Transpose a n*m matrix*/
+
+        int i,j;
+
+        /*Check the the matrix has been allocated*/
+        ISSMASSERT(tA);
+
+        /*Transpose*/
+        for(i=0;i<ncols;i++){
+                for(j=0;j<nrows;j++){
+                        tA[i*nrows+j]=A[j*ncols+i];
+                }
+        }
+
+}
+/*}}}*/

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

@@ -13 +13 @@
 void Matrix3x3Invert(double* Ainv, double* A);
 void Matrix3x3Determinant(double* Adet,double* A);
+void MatrixTranspose(double* tA,double* A,int nrows, int ncols);
 
 #endif //ifndef _MATRIXUTILS_H_