Index: /issm/trunk-jpl/src/c/Makefile.am
===================================================================
--- /issm/trunk-jpl/src/c/Makefile.am	(revision 13595)
+++ /issm/trunk-jpl/src/c/Makefile.am	(revision 13596)
@@ -771,5 +771,4 @@
 					./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\
Index: /issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 13595)
+++ /issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 13596)
@@ -37,5 +37,5 @@
 PetscMat::PetscMat(int M,int N, IssmDouble sparsity){
 
-	this->matrix=NewMat(M,N,sparsity);
+	this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm());
 }
 /*}}}*/
@@ -54,5 +54,5 @@
 
 
-	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);
@@ -67,5 +67,5 @@
 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: /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 13595)
+++ /issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp	(revision 13596)
@@ -48,5 +48,5 @@
 }
 /*}}}*/
-/*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){
 
Index: /issm/trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h	(revision 13595)
+++ /issm/trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h	(revision 13596)
@@ -24,5 +24,4 @@
 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);