Changeset 11704

Timestamp:

03/14/12 09:09:05 (13 years ago)

Author:

Mathieu Morlighem

Message:

Added constructor from PETSc vector

Location:

issm/trunk-jpl/src/c/objects/Numerics

Files:

• Vector.cpp (modified) (3 diffs)
• Vector.h (modified) (1 diff)

issm/trunk-jpl/src/c/objects/Numerics/Vector.cpp

@@ -78 +78 @@
 }
 /*}}}*/
+#ifdef _HAVE_PETSC_
+/*FUNCTION Vector::Vector(Vec petsc_vec){{{1*/
+Vector::Vector(Vec petsc_vec){
+
+        /*Get Vector size*/
+        VecGetSize(petsc_vec,&this->M);
+
+        /*copy vector*/
+        VecDuplicate(petsc_vec,&this->vector);
+        VecCopy(petsc_vec,this->vector);
+
+}
+/*}}}*/
+#endif
                 /*FUNCTION Vector::~Vector(){{{1*/
 Vector::~Vector(){
@@ -91 +105 @@
 }
 /*}}}*/
-
 
 /*Vector specific routines: */
@@ -240 +253 @@
                 _error_("not implemented yet!");
         #endif
-
 }
 /*}}}*/

issm/trunk-jpl/src/c/objects/Numerics/Vector.h

@@ -46 +46 @@
                 Vector(int M,bool fromlocalsize=false);
                 Vector(double* serial_vec,int pM);
+                #ifdef _HAVE_PETSC_
+                Vector(Vec petsc_vec);
+                #endif
                 ~Vector();
                 /*}}}*/