Merge pull request #205 from ChrisPec27/chris/pass_through_parabolic_…

…data Chris/pass through parabolic data
conservation-laws · Dec 31, 2024 · 1b24c75 · 1b24c75
2 parents 0348cbb + 18a366b
commit 1b24c75
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 4 deletions.
diff --git a/source/hyperbolic_module.template.h b/source/hyperbolic_module.template.h
@@ -1190,6 +1190,13 @@ namespace ryujin
 
     CALLGRIND_STOP_INSTRUMENTATION;
 
+    /*
+     * Pass through the parabolic state vector
+     */
+    const auto &old_V = std::get<2>(old_state_vector);
+    auto &new_V = std::get<2>(new_state_vector);
+    new_V = old_V;
+
     /*
      * Do we have to restart?
      */

diff --git a/source/time_integrator.template.h b/source/time_integrator.template.h
@@ -11,17 +11,17 @@ namespace ryujin
 {
   using namespace dealii;
 
-
-  /**
-   * TODO: clear out precomputed vector and also scale add V.
-   */
   template <typename StateVector, typename Number>
   void
   sadd(StateVector &dst, const Number s, const Number b, const StateVector &src)
   {
     auto &dst_U = std::get<0>(dst);
     auto &src_U = std::get<0>(src);
     dst_U.sadd(s, b, src_U);
+
+    auto &dst_V = std::get<2>(dst);
+    auto &src_V = std::get<2>(src);
+    dst_V.sadd(s, b, src_V);
   }