source: issm/oecreview/Archive/12678-13393/ISSM-12862-12863.diff

u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/issm/SeqVec.cpp

@@ -41 +41 @@
         this->M=pM;
         this->vector=NULL;
         if(this->M){
-                this->vector=xNewZeroInit<IssmDouble>(pM);
+                this->vector=xNew<IssmDouble>(pM);
                 xMemCpy<IssmDouble>(this->vector,buffer,pM);
         }
 }

u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp

@@ -33 +33 @@
 
 }
 /*}}}*/
+/*FUNCTION PetscVec::PetscVec(Vec petsc_vec){{{*/
+PetscVec::PetscVec(Vec petsc_vec){
+
+        if(petsc_vec==NULL){
+                this->vector=NewVec(0);
+        }
+        else{
+                /*copy vector*/
+                VecDuplicate(petsc_vec,&this->vector);
+                VecCopy(petsc_vec,this->vector);
+        }
+
+}
+/*}}}*/
 /*FUNCTION PetscVec::PetscVec(IssmDouble* serial_vec,int M){{{*/
 PetscVec::PetscVec(IssmDouble* serial_vec,int M){
 
@@ -53 +67 @@
     VecFree(&this->vector);
 }
 /*}}}*/
-/*FUNCTION Vector::Vector(Vec petsc_vec){{{*/
-PetscVec::PetscVec(Vec petsc_vec){
 
-        if(petsc_vec==NULL){
-                this->vector=NewVec(0);
-        }
-        else{
-                /*copy vector*/
-                VecDuplicate(petsc_vec,&this->vector);
-                VecCopy(petsc_vec,this->vector);
-        }
-
-}
-/*}}}*/
-
 /*PetscVec specific routines: */
 /*FUNCTION PetscVec::Echo{{{*/
 void PetscVec::Echo(void){