Changeset 13901

Timestamp:

11/07/12 16:23:29 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added printsparsity, very useful to debug memory allocation module

Location:

issm/trunk-jpl/src/c/shared/Elements

Files:

• PrintArrays.cpp (modified) (1 diff)
• elements.h (modified) (1 diff)

issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp

@@ -9 +9 @@
                 _printString_("   [ ");
                 for(int j=0;j<cols;j++) _printString_( " " << setw(11) << setprecision (5) << array[i*cols+j]);
+                _printLine_(" ]");
+        }  
+        _printLine_("");
+}
+void printsparsity(IssmPDouble* array,int lines,int cols){
+        _printLine_("");
+        for(int i=0;i<lines;i++){  
+                _printString_("   [ ");
+                for(int j=0;j<cols;j++){
+                        if(array[i*cols+j]!=0.0)
+                         _printString_( " X");
+                        else
+                         _printString_( " .");
+                }
                 _printLine_(" ]");
         }  

issm/trunk-jpl/src/c/shared/Elements/elements.h

@@ -40 +40 @@
 void printarray(IssmPDouble* array,int lines,int cols=1);
 void printarray(int* array,int lines,int cols=1);
+void printsparsity(IssmPDouble* array,int lines,int cols=1);
 void printbinary(int n);
 #endif //ifndef _SHARED_ELEMENTS_H_