Index: /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 13900)
+++ /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 13901)
@@ -9,4 +9,18 @@
 		_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_(" ]");
 	}  
Index: /issm/trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 13900)
+++ /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 13901)
@@ -40,4 +40,5 @@
 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_