Index: ../trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp
===================================================================
--- ../trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 13595)
+++ ../trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 13596)
@@ -47,7 +47,7 @@
 	return outmatrix;
 }
 /*}}}*/
-/*NewMat(int M,int N,double sparsity){{{*/
+/*NewMat(int M,int N,double sparsity,COMM comm){{{*/
 Mat NewMat(int M,int N,double sparsity,COMM comm){
 
 	/*output:*/
Index: ../trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h
===================================================================
--- ../trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h	(revision 13595)
+++ ../trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h	(revision 13596)
@@ -23,7 +23,6 @@
 
 int VecTranspose(Vec* tvector,Vec vector);
 int VecToMPISerial(double** pgathered_vector, Vec vector,COMM comm);
-void GetOwnershipBoundariesFromRange(int* plower_row,int* pupper_row,int range,COMM comm);
 void MatFree(Mat* pmat);
 void ISFree(IS* pis);
 void VecFree(Vec* pvec);
Index: ../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp
===================================================================
--- ../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 13595)
+++ ../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 13596)
@@ -36,7 +36,7 @@
 /*FUNCTION PetscMat::PetscMat(int M,int N, IssmDouble sparsity){{{*/
 PetscMat::PetscMat(int M,int N, IssmDouble sparsity){
 
-	this->matrix=NewMat(M,N,sparsity);
+	this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm());
 }
 /*}}}*/
 /*FUNCTION PetscMat(IssmDouble* serial_mat,int M,int N,IssmDouble sparsity){{{*/
@@ -53,7 +53,7 @@
 	for(i=0;i<N;i++)idxn[i]=i;
 
 
-	this->matrix=NewMat(M,N,sparsity);
+	this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm());
 	MatSetValues(this->matrix,M,idxm,N,idxn,serial_mat,INSERT_VALUES);
 	MatAssemblyBegin(this->matrix,MAT_FINAL_ASSEMBLY);
 	MatAssemblyEnd(this->matrix,MAT_FINAL_ASSEMBLY);
@@ -66,7 +66,7 @@
 /*FUNCTION PetscMat::PetscMat(int M,int N, int connectivity, int numberofdofspernode){{{*/
 PetscMat::PetscMat(int M,int N, int connectivity,int numberofdofspernode){
 	
-	this->matrix=NewMat(M,N,connectivity,numberofdofspernode);
+	this->matrix=NewMat(M,N,connectivity,numberofdofspernode,IssmComm::GetComm());
 
 }
 /*}}}*/
Index: ../trunk-jpl/src/c/Makefile.am
===================================================================
--- ../trunk-jpl/src/c/Makefile.am	(revision 13595)
+++ ../trunk-jpl/src/c/Makefile.am	(revision 13596)
@@ -770,7 +770,6 @@
 					./toolkits/petsc/patches/VecDuplicatePatch.cpp\
 					./toolkits/petsc/patches/KSPFree.cpp\
 					./toolkits/petsc/patches/MatFree.cpp\
-					./toolkits/petsc/patches/GetOwnershipBoundariesFromRange.cpp\
 					./toolkits/petsc/patches/VecPartition.cpp\
 					./toolkits/petsc/patches/MatInvert.cpp\
 					./toolkits/petsc/patches/MatMultPatch.cpp\