Changeset 5885

Timestamp:

09/18/10 13:14:56 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added CheckConsitency method, used only in debugging mode

Location:

issm/trunk/src/c/objects/Numerics

Files:

• ElementMatrix.cpp (modified) (1 diff)
• ElementMatrix.h (modified) (1 diff)

issm/trunk/src/c/objects/Numerics/ElementMatrix.cpp

@@ -360 +360 @@
 }
 /*}}}*/
+/*FUNCTION ElementMatrix::CheckConsistency{{{1*/
+void ElementMatrix::CheckConsistency(void){
+        /*Check element matrix values, only in debugging mode*/
+        #ifdef _ISSM_DEBUG_ 
+        for (int i=0;i<this->nrows;i++){
+                for(int j=0;j<this->ncols;j++){
+                        if (isnan(this->values[i*this->ncols+j])) ISSMERROR("NaN found in Element Matrix");
+                        if (fabs(this->values[i*this->ncols+j])>1.e+50) ISSMERROR("Element Matrix values exceeds 1.e+50");
+                }
+        }
+        #endif
+}
+/*}}}*/
 /*FUNCTION ElementMatrix::Echo{{{1*/
 void ElementMatrix::Echo(void){

issm/trunk/src/c/objects/Numerics/ElementMatrix.h

@@ -60 +60 @@
                 void AddToGlobal(Mat Kgg, Mat Kff, Mat Kfs);
                 void Echo(void);
+                void CheckConsistency(void);
                 void Init(ElementMatrix* Ke);
                 /*}}}*/
 };
 #endif //#ifndef _ELEMENT_MATRIX_H_
-