Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 12856)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp	(revision 12857)
@@ -30,8 +30,7 @@
 /*FUNCTION PetscMat::PetscMat(int M,int N){{{*/
 PetscMat::PetscMat(int M,int N){
 
-	this->matrix=NULL;
-	if(M*N)this->matrix=NewMat(M,N);
+	this->matrix=NewMat(M,N);
 }
 /*}}}*/
 /*FUNCTION PetscMat::PetscMat(int M,int N, IssmDouble sparsity){{{*/
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp	(revision 12856)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp	(revision 12857)
@@ -28,8 +28,9 @@
 /*}}}*/
 /*FUNCTION PetscVec::PetscVec(int M,bool fromlocalsize){{{*/
 PetscVec::PetscVec(int M,bool fromlocalsize){
-	this->vector=NULL;
-	if(M) this->vector=NewVec(M,fromlocalsize);
+	
+	this->vector=NewVec(M,fromlocalsize);
+
 }
 /*}}}*/
 /*FUNCTION PetscVec::PetscVec(IssmDouble* serial_vec,int M){{{*/