Index: /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp
===================================================================
--- /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp	(revision 13226)
+++ /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp	(revision 13227)
@@ -216,6 +216,6 @@
 }
 /*}}}*/
-/*FUNCTION WriteData(mxArray** pdataref,Matrix* matrix){{{*/
-void WriteData(mxArray** pdataref,Matrix<double>* matrix){
+/*FUNCTION WriteData(mxArray** pdataref,SeqMat<double>* matrix){{{*/
+void WriteData(mxArray** pdataref,SeqMat<double>* matrix){
 		
 	int      i,j;
@@ -246,15 +246,14 @@
 		/*Free ressources:*/
 		xDelete<double>(matrix_ptr);
-
-	}
-	else{
-		dataref = mxCreateDoubleMatrix(0,0,mxREAL);
-	}
-
-	*pdataref=dataref;
-}
-/*}}}*/
-/*FUNCTION WriteData(mxArray** pdataref,Vector<double>* vector){{{*/
-void WriteData(mxArray** pdataref,Vector<double>* vector){
+	}
+	else{
+		dataref = mxCreateDoubleMatrix(0,0,mxREAL);
+	}
+
+	*pdataref=dataref;
+}
+/*}}}*/
+/*FUNCTION WriteData(mxArray** pdataref,SeqVec<double>* vector){{{*/
+void WriteData(mxArray** pdataref,SeqVec<double>* vector){
 	
 	mxArray* dataref=NULL;
Index: /issm/trunk-jpl/src/c/matlab/io/matlabio.h
===================================================================
--- /issm/trunk-jpl/src/c/matlab/io/matlabio.h	(revision 13226)
+++ /issm/trunk-jpl/src/c/matlab/io/matlabio.h	(revision 13227)
@@ -18,8 +18,8 @@
 #include <mex.h>
 
-void WriteData(mxArray** pdataref,Matrix<double>* matrix);
+void WriteData(mxArray** pdataref,SeqMat<double>* matrix);
 void WriteData(mxArray** pdataref,double* matrix, int M,int N);
 void WriteData(mxArray** pdataref,int*    matrix, int M,int N);
-void WriteData(mxArray** pdataref,Vector<double>* vector);
+void WriteData(mxArray** pdataref,SeqVec<double>* vector);
 void WriteData(mxArray** pdataref,double* vector, int M);
 void WriteData(mxArray** pdataref,int integer);
@@ -85,4 +85,3 @@
 #endif
 
-
 #endif	/* _IO_H_ */
Index: /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.cpp	(revision 13226)
+++ /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.cpp	(revision 13227)
@@ -20,5 +20,5 @@
 /*}}}*/
 
-void TriMeshx(Matrix<IssmPDouble>** pindex,Vector<IssmPDouble>** px,Vector<IssmPDouble>** py,Matrix<IssmPDouble>** psegments,Vector<IssmPDouble>** psegmentmarkerlist,DataSet* domain,DataSet* rifts,double area){
+void TriMeshx(SeqMat<IssmPDouble>** pindex,SeqVec<IssmPDouble>** px,SeqVec<IssmPDouble>** py,SeqMat<IssmPDouble>** psegments,SeqVec<IssmPDouble>** psegmentmarkerlist,DataSet* domain,DataSet* rifts,double area){
 
 	/*indexing: */
@@ -26,15 +26,15 @@
 
 	/*output: */
-	IssmPDouble* index=NULL;
-	Matrix<IssmPDouble>* index_matrix=NULL;
-	double* x=NULL;
-	double* y=NULL;
-	double* segments=NULL;
-	Matrix<IssmPDouble>* segments_matrix=NULL;
-	double* segmentmarkerlist=NULL;
+	IssmPDouble         *index             = NULL;
+	SeqMat<IssmPDouble> *index_matrix      = NULL;
+	double              *x                 = NULL;
+	double              *y                 = NULL;
+	double              *segments          = NULL;
+	SeqMat<IssmPDouble> *segments_matrix   = NULL;
+	double              *segmentmarkerlist = NULL;
 
 	/*intermediary: */
-	int      counter,counter2,backcounter;
-	Contour<IssmPDouble>* contour=NULL;
+	int counter,counter2,backcounter;
+	Contour<IssmPDouble> *contour = NULL;
 
 	/* Triangle structures needed to call Triangle library routines: */
@@ -192,12 +192,12 @@
 
 	/*Output : */
-	index_matrix=new Matrix<IssmPDouble>(index,out.numberoftriangles,3,1.0,SeqMatType);
+	index_matrix=new SeqMat<IssmPDouble>(index,out.numberoftriangles,3,1.0);
 	*pindex=index_matrix;
 	
-	segments_matrix=new Matrix<IssmPDouble>(segments,out.numberofsegments,3,1.0,SeqMatType);
+	segments_matrix=new SeqMat<IssmPDouble>(segments,out.numberofsegments,3,1.0);
 	*psegments=segments_matrix;
 
-	*px=new Vector<IssmPDouble>(x,out.numberofpoints,SeqMatType);
-	*py=new Vector<IssmPDouble>(y,out.numberofpoints,SeqMatType);
-	*psegmentmarkerlist=new Vector<IssmPDouble>(segmentmarkerlist,out.numberofsegments,SeqMatType);
+	*px=new SeqVec<IssmPDouble>(x,out.numberofpoints);
+	*py=new SeqVec<IssmPDouble>(y,out.numberofpoints);
+	*psegmentmarkerlist=new SeqVec<IssmPDouble>(segmentmarkerlist,out.numberofsegments);
 }
Index: /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.h
===================================================================
--- /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.h	(revision 13226)
+++ /issm/trunk-jpl/src/c/modules/TriMeshx/TriMeshx.h	(revision 13227)
@@ -11,5 +11,5 @@
 
 /* local prototypes: */
-void TriMeshx(Matrix<IssmPDouble>** pindex,Vector<IssmPDouble>** px,Vector<IssmPDouble>** py,Matrix<IssmPDouble>** psegments,Vector<IssmPDouble>** psegmentmarkerlist,DataSet* domain,DataSet* rifts,double area);
+void TriMeshx(SeqMat<IssmPDouble>** pindex,SeqVec<IssmPDouble>** px,SeqVec<IssmPDouble>** py,SeqMat<IssmPDouble>** psegments,SeqVec<IssmPDouble>** psegmentmarkerlist,DataSet* domain,DataSet* rifts,double area);
 
 #endif  /* _TRIMESHX_H */