From 7352f04249df810e479da67b5b72ddd236b59907 Mon Sep 17 00:00:00 2001 From: William Throwe Date: Wed, 11 Sep 2024 18:57:13 -0400 Subject: [PATCH] Change non-monotonic AM predictor coefficients Either version works, but this version is simpler and is what will be in the paper. --- src/Time/TimeSteppers/AdamsMoultonPc.cpp | 11 +---------- 1 file changed, 1 insertion(+), 10 deletions(-) diff --git a/src/Time/TimeSteppers/AdamsMoultonPc.cpp b/src/Time/TimeSteppers/AdamsMoultonPc.cpp index 2d64b4a4f1db..904aec9aafe9 100644 --- a/src/Time/TimeSteppers/AdamsMoultonPc.cpp +++ b/src/Time/TimeSteppers/AdamsMoultonPc.cpp @@ -396,26 +396,17 @@ void AdamsMoultonPc::add_boundary_delta_impl( } else { adams_lts::AdamsScheme scheme{adams_lts::SchemeType::Implicit, current_order}; - auto remote_scheme = scheme; if (local_times.number_of_substeps(local_times.size() - 1) == 1) { // Predictor scheme = {adams_lts::SchemeType::Explicit, current_order - 1}; ASSERT(remote_times.back() <= local_times.back(), "Unexpected remote values available."); - // If the sides are not aligned, we use the predictor data - // available from the neighbor. If they are, that data has not - // been received. - remote_scheme = {remote_times.back() == local_times.back() - ? adams_lts::SchemeType::Explicit - : adams_lts::SchemeType::Implicit, - current_order - 1}; } const auto lts_coefficients = adams_lts::lts_coefficients( local_times, remote_times, local_times.back().step_time(), - local_times.back().step_time() + time_step, scheme, remote_scheme, - scheme); + local_times.back().step_time() + time_step, scheme, scheme, scheme); adams_lts::apply_coefficients(result, lts_coefficients, coupling); } }