Index: ../trunk-jpl/src/c/shared/Elements/PrintArrays.cpp
===================================================================
--- ../trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 13900)
+++ ../trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 13901)
@@ -12,6 +12,20 @@
 	}  
 	_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_(" ]");
+	}  
+	_printLine_("");
+}
 void printarray(int* array,int lines,int cols){
 	_printLine_("");
 	for(int i=0;i<lines;i++){  
Index: ../trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- ../trunk-jpl/src/c/shared/Elements/elements.h	(revision 13900)
+++ ../trunk-jpl/src/c/shared/Elements/elements.h	(revision 13901)
@@ -39,5 +39,6 @@
 /*Print arrays*/
 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_