Index: /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h	(revision 14857)
+++ /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h	(revision 14858)
@@ -604,6 +604,13 @@
 		IssmAbsVec<IssmDouble>* Solve(IssmAbsVec<IssmDouble>* pfin, Parameters* parameters){
 
-			printf("IssmMpiDenseMat solve not implemented yet!");
-			exit(1);
+			/*Assume we are getting an IssmMpiVec in input, downcast: */
+			IssmMpiVec<IssmDouble>* pf=(IssmMpiVec<IssmDouble>*)pfin;
+
+			/*Let's try and use the MUMPS solver here: */
+			#ifdef _HAVE_MUMPS_
+
+			#else
+			_error_("IssmMpiDenseMat solver requires MUMPS solver");
+			#endif
 
 		}/*}}}*/