From 68ec508460bcc25a61c10566917d7045744f7a42 Mon Sep 17 00:00:00 2001 From: David Pecoraro Date: Thu, 19 Dec 2024 17:32:39 +0900 Subject: [PATCH 1/2] time_integrator: add in sadd for parabolic vector --- source/time_integrator.template.h | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/source/time_integrator.template.h b/source/time_integrator.template.h index 20ab34c7..6524e485 100644 --- a/source/time_integrator.template.h +++ b/source/time_integrator.template.h @@ -11,10 +11,6 @@ namespace ryujin { using namespace dealii; - - /** - * TODO: clear out precomputed vector and also scale add V. - */ template void sadd(StateVector &dst, const Number s, const Number b, const StateVector &src) @@ -22,6 +18,10 @@ namespace ryujin 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); } From 18a366b436163b6a9fc3f5d8b9a3a10207b33d6a Mon Sep 17 00:00:00 2001 From: David Pecoraro Date: Thu, 19 Dec 2024 17:35:15 +0900 Subject: [PATCH 2/2] hyperbolic_module: chain through the parabolic state vector --- source/hyperbolic_module.template.h | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/source/hyperbolic_module.template.h b/source/hyperbolic_module.template.h index dfb329e3..cce9eb7e 100644 --- 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? */