source: issm/oecreview/Archive/19101-20495/ISSM-19976-19977.diff@ 20498

../trunk-jpl/src/c/toolkits/objects/Vector.h

@@ -280 +280 @@
 
                 }
                 /*}}}*/
+                void Shift(doubletype shift){_assert_(this);/*{{{*/
+
+                        if(type==PetscVecType){
+                                #ifdef _HAVE_PETSC_
+                                this->pvector->Shift(shift);
+                                #endif
+                        }
+                        else this->ivector->Shift(shift);
+                }
+                /*}}}*/
                 void Copy(Vector* to){_assert_(this);/*{{{*/
 
                         if(type==PetscVecType){

../trunk-jpl/src/c/toolkits/issm/IssmAbsVec.h

@@ -41 +41 @@
                 virtual void AXPY(IssmAbsVec* X, doubletype a)=0;
                 virtual void AYPX(IssmAbsVec* X, doubletype a)=0;
                 virtual doubletype* ToMPISerial(void)=0;
+                virtual void Shift(doubletype shift)=0;
                 virtual void Copy(IssmAbsVec* to)=0;
                 virtual doubletype Norm(NormMode mode)=0;
                 virtual void Scale(doubletype scale_factor)=0;

../trunk-jpl/src/c/toolkits/issm/IssmVec.h

@@ -179 +179 @@
                         return vector->ToMPISerial();
                 }
                 /*}}}*/
+                void Shift(doubletype shift){/*{{{*/
+                        vector->Shift(shift);
+                }
+                /*}}}*/
                 void Copy(IssmVec* to){/*{{{*/
                         vector->Copy(to->vector);
                 }

../trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h

@@ -395 +395 @@
 
                 }
                 /*}}}*/
+                void Shift(doubletype shift){/*{{{*/
+                        for(int i=0;i<this->m;i++)this->vector[i]+=shift;
+                }
+                /*}}}*/
                 void Copy(IssmAbsVec<doubletype>* toin){/*{{{*/
 
                         int i;

../trunk-jpl/src/c/toolkits/issm/IssmSeqVec.h

@@ -198 +198 @@
 
                 }
                 /*}}}*/
+                void Shift(doubletype shift){/*{{{*/
+
+                        for(int i=0;i<this->M;i++)this->vector[i]+=shift;
+
+                }
+                /*}}}*/
                 void Copy(IssmAbsVec<doubletype>* toin){/*{{{*/
 
                         int i;

../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.h

@@ -49 +49 @@
                 void        AXPY(PetscVec* X, IssmDouble a);
                 void        AYPX(PetscVec* X, IssmDouble a);
                 IssmDouble* ToMPISerial(void);
+                void        Shift(IssmDouble shift);
                 void        Copy(PetscVec* to);
                 IssmDouble  Norm(NormMode norm_type);
                 IssmDouble  Max(void);

../trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp

@@ -162 +162 @@
 
 }
 /*}}}*/
+void PetscVec::Shift(IssmDouble shift){/*{{{*/
+
+        if(this->vector) VecShift(this->vector,shift);
+
+}
+/*}}}*/
 void PetscVec::Copy(PetscVec* to){/*{{{*/
 
         if(this->vector) VecCopy(this->vector,to->vector);