Index: /issm/trunk-jpl/src/c/matlab/io/FetchMatlabData.cpp
===================================================================
--- /issm/trunk-jpl/src/c/matlab/io/FetchMatlabData.cpp	(revision 13374)
+++ /issm/trunk-jpl/src/c/matlab/io/FetchMatlabData.cpp	(revision 13375)
@@ -72,5 +72,9 @@
 		outmatrix=NULL;
 	}
-	else if (mxIsClass(dataref,"double") ){
+	else if( mxIsClass(dataref,"double") || 
+				mxIsClass(dataref,"single") || 
+				mxIsClass(dataref,"int16") || 
+				mxIsClass(dataref,"int8") || 
+				mxIsClass(dataref,"uint8")){
 
 		/*Check dataref is not pointing to NaN: */
@@ -82,5 +86,7 @@
 		}
 		else{
-
+			if(!mxIsClass(dataref,"double") && !mxIsClass(dataref,"single")){
+				_printLine_("Warning: converting matlab data from '" << mxGetClassName(dataref) << "' to 'double'");
+			}
 			/*Convert matlab n-dim array to double* matrix: */
 			MatlabNArrayToNArray(&outmatrix,&outmatrix_numel,&outmatrix_ndims,&outmatrix_size,dataref);
@@ -114,5 +120,9 @@
 		outmatrix=NULL;
 	}
-	else if (mxIsClass(dataref,"double") ){
+	else if( mxIsClass(dataref,"double") || 
+				mxIsClass(dataref,"single") || 
+				mxIsClass(dataref,"int16") || 
+				mxIsClass(dataref,"int8") || 
+				mxIsClass(dataref,"uint8")){
 
 		/*Check dataref is not pointing to NaN: */
@@ -123,5 +133,7 @@
 		}
 		else{
-
+			if(!mxIsClass(dataref,"double") && !mxIsClass(dataref,"single")){
+				_printLine_("Warning: converting matlab data from '" << mxGetClassName(dataref) << "' to 'double'");
+			}
 			/*Convert matlab matrix to double* matrix: */
 			MatlabMatrixToDoubleMatrix(&doublematrix,&outmatrix_rows,&outmatrix_cols,dataref);
Index: /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13374)
+++ /issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp	(revision 13375)
@@ -25,5 +25,5 @@
 void SolverxSeq(SeqVec<IssmDouble>** puf,SeqMat<IssmDouble>* Kff, SeqVec<IssmDouble>* pf, Parameters* parameters){/*{{{*/
 
-	#ifdef _HAVE_GSL_
+#ifdef _HAVE_GSL_
 	/*Intermediary: */
 	int M,N,N2,s;
@@ -35,5 +35,6 @@
 	if(N!=N2)_error_("Right hand side vector of size " << N2 << ", when matrix is of size " << M << "-" << N << " !");
 	if(M!=N)_error_("Stiffness matrix should be square!");
-        IssmDouble *x  = xNew<IssmDouble>(N);
+	IssmDouble *x  = xNew<IssmDouble>(N);
+
 #ifdef _HAVE_ADOLC_
 	SolverxSeq(x,Kff->matrix,pf->vector,N,parameters);
@@ -41,4 +42,5 @@
 	SolverxSeq(x,Kff->matrix,pf->vector,N);
 #endif
+
 	uf=new SeqVec<IssmDouble>(x,N);
 	xDelete(x);
@@ -47,7 +49,7 @@
 	*puf=uf;
 
-	#else
-		_error_("GSL support not compiled in!");
-	#endif
+#else
+	_error_("GSL support not compiled in!");
+#endif
 
 }/*}}}*/