From 893b4579692243e7a2ba8b33ea5242508c562130 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Mon, 2 Oct 2023 10:43:24 -0400 Subject: [PATCH] fix the coefficients for 2nd order Radau SDC --- Source/sdc/sdc_util.cpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/Source/sdc/sdc_util.cpp b/Source/sdc/sdc_util.cpp index 8c91596abc..68b06805d2 100644 --- a/Source/sdc/sdc_util.cpp +++ b/Source/sdc/sdc_util.cpp @@ -222,8 +222,8 @@ Castro::ca_sdc_compute_C2_radau(const Box& bx, C(i,j,k,n) = -R_1_old(i,j,k,n) + (A_m(i,j,k,n) - A_0_old(i,j,k,n)) + (dt/dt_m) * (1.0_rt/12.0_rt) * - 5.0_rt*(A_1_old(i,j,k,n) + R_1_old(i,j,k,n)) - - (A_2_old(i,j,k,n) + R_2_old(i,j,k,n)); + (5.0_rt*(A_1_old(i,j,k,n) + R_1_old(i,j,k,n)) - + (A_2_old(i,j,k,n) + R_2_old(i,j,k,n))); }); } else if (m_start == 1) @@ -233,8 +233,8 @@ Castro::ca_sdc_compute_C2_radau(const Box& bx, C(i,j,k,n) = -R_2_old(i,j,k,n) + (A_m(i,j,k,n) - A_1_old(i,j,k,n)) + (dt/dt_m) * (1.0_rt/3.0_rt) * - (A_1_old(i,j,k,n) + R_1_old(i,j,k,n)) + - (A_2_old(i,j,k,n) + R_2_old(i,j,k,n)); + ((A_1_old(i,j,k,n) + R_1_old(i,j,k,n)) + + (A_2_old(i,j,k,n) + R_2_old(i,j,k,n))); }); } } // end ca_sdc_compute_C2_radau @@ -308,7 +308,7 @@ Castro::ca_sdc_compute_C4_radau(const Box& bx, // compute the 'C' term for the 4th-order solve with reactions // note: this 'C' is cell-averages - // Gauss-Lobatto (Simpsons) + // Radau if (m_start == 0) {