add method to retreive updated position from IBImplicitStrategy

IBAMR · boyceg · Oct 3, 2019 · Oct 3, 2019 · Oct 3, 2019 · Dec 30, 2019
commit 5e587ec86cedf9c4295eb3f14bd65670b2d82f68
diff --git a/include/ibamr/IBImplicitStrategy.h b/include/ibamr/IBImplicitStrategy.h
@@ -83,6 +83,11 @@ class IBImplicitStrategy : public IBStrategy
      */
     virtual void setUpdatedPosition(Vec& X_new_vec) = 0;
 
+    /*!
+     * Get the value of the updated position vector.
+     */
+    virtual void getUpdatedPosition(Vec& X_new_vec) = 0;
+
     /*!
      * Compute the nonlinear residual for backward Euler time stepping.
      */

diff --git a/include/ibamr/IBMethod.h b/include/ibamr/IBMethod.h
@@ -203,6 +203,11 @@ class IBMethod : public IBImplicitStrategy
      */
     void setUpdatedPosition(Vec& X_new_vec) override;
 
+    /*!
+     * Get the value of the updated position vector.
+     */
+    void getUpdatedPosition(Vec& X_new_vec) override;
+
     /*!
      * Compute the nonlinear residual for backward Euler time stepping.
      */

diff --git a/src/IB/IBMethod.cpp b/src/IB/IBMethod.cpp
@@ -524,6 +524,16 @@ IBMethod::setUpdatedPosition(Vec& X_new_vec)
     return;
 } // setUpdatedPosition
 
+void
+IBMethod::getUpdatedPosition(Vec& X_new_vec)
+{
+    PetscErrorCode ierr;
+    const int level_num = d_hierarchy->getFinestLevelNumber();
+    ierr = VecCopy(d_X_new_data[level_num]->getVec(), X_new_vec);
+    IBTK_CHKERRQ(ierr);
+    return;
+} // getUpdatedPosition
+
 void
 IBMethod::computeResidualBackwardEuler(Vec& R_vec)
 {