Index: /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp	(revision 18338)
+++ /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp	(revision 18339)
@@ -339,4 +339,20 @@
 	}
 	#endif
+}
+/*}}}*/
+void ElementMatrix::Lump(void){/*{{{*/
+
+	if(!dofsymmetrical) _error_("not supported yet");
+
+	for(int i=0;i<this->nrows;i++){
+		for(int j=0;j<this->ncols;j++){
+			if(i!=j){
+				this->values[i*this->ncols+i] += this->values[i*this->ncols+j];
+				this->values[i*this->ncols+j]  = 0.;
+			}
+		}
+	}
+
+	return;
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h	(revision 18338)
+++ /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h	(revision 18339)
@@ -63,4 +63,5 @@
 		void CheckConsistency(void);
 		void StaticCondensation(int numindices,int* indices);
+		void Lump(void);
 		void Transpose(void);
 		void Init(ElementMatrix* Ke);