From a50eada7f8026baf1a8fb57ee429f336e14b8fc8 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Wed, 14 Jan 2026 19:14:08 -0500 Subject: [PATCH 01/15] disable integrator.do_species_clip by default it seems to cause more problems then it fixes --- Docs/source/ode_integrators.rst | 6 ++++-- integration/_parameters | 2 +- 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/Docs/source/ode_integrators.rst b/Docs/source/ode_integrators.rst index 4644a347a5..c180bfaa42 100644 --- a/Docs/source/ode_integrators.rst +++ b/Docs/source/ode_integrators.rst @@ -275,9 +275,11 @@ constraint on the intermediate states during the integration. The default is ``1.e-30``. * ``integrator.do_species_clip`` : this enforces that the mass fractions - all in $[\mathtt{SMALL\_X\_SAFE}, 1.0]$. + all in $[\mathtt{SMALL\_X\_SAFE}, 1.0]$ before calling the network righthand + side function. - This is enabled by default. + This is off by default. Turning this on can sometimes make the integrator + work a lot harder. diff --git a/integration/_parameters b/integration/_parameters index a11b5f3ea8..bf697b44f7 100644 --- a/integration/_parameters +++ b/integration/_parameters @@ -91,7 +91,7 @@ retry_atol_enuc real -1 # in the clean_state process, do we clip the species such that they # are in [0, 1]? -do_species_clip bool 1 +do_species_clip bool 0 # flag for turning on the use of number densities for all species use_number_densities bool 0 From 290d278dcba63ef621d92601a7ab822d130128bd Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 15:48:38 -0500 Subject: [PATCH 02/15] VODE: bail out of the nonlinear solve if we predict negative mass fractions --- integration/VODE/vode_dvnlsd.H | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) diff --git a/integration/VODE/vode_dvnlsd.H b/integration/VODE/vode_dvnlsd.H index 46ede5d86c..60fbc401d4 100644 --- a/integration/VODE/vode_dvnlsd.H +++ b/integration/VODE/vode_dvnlsd.H @@ -112,6 +112,32 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) (vstate.RL1 * vstate.yh(i,2) + vstate.acor(i)); } + // sometime VODE goes way off tangent. If these mass fractions + // are really bad, then let's just bail now + + { + amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; +#ifdef SDC + amrex::Real thresh = species_failure_tolerance * rho_current; +#else + amrex::Real thresh = species_failure_tolerance; +#endif + bool fail{}; + for (int i = 1; i <= NumSpec; ++i) { + if (vstate.y(i) < -thresh) { + fail = true; + break; + } + } + + if (fail) { + NFLAG = -1; + vstate.ICF = 2; + vstate.IPUP = 1; + return ACNRM; + } + } + #ifdef NEW_NETWORK_IMPLEMENTATION RHS::dgesl(vstate.jac, vstate.y); #else From 779fa534822ba9c48ead0a8f842fbadbeaf0ac1e Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 17:16:49 -0500 Subject: [PATCH 03/15] add param --- integration/VODE/vode_dvnlsd.H | 8 ++++---- integration/_parameters | 4 ++++ 2 files changed, 8 insertions(+), 4 deletions(-) diff --git a/integration/VODE/vode_dvnlsd.H b/integration/VODE/vode_dvnlsd.H index 60fbc401d4..ad956f8ff7 100644 --- a/integration/VODE/vode_dvnlsd.H +++ b/integration/VODE/vode_dvnlsd.H @@ -115,12 +115,12 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) // sometime VODE goes way off tangent. If these mass fractions // are really bad, then let's just bail now - { - amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; + if (integrator_rp::do_corrector_validation) { #ifdef SDC - amrex::Real thresh = species_failure_tolerance * rho_current; + const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; + const amrex::Real thresh = species_failure_tolerance * rho_current; #else - amrex::Real thresh = species_failure_tolerance; + const amrex::Real thresh = species_failure_tolerance; #endif bool fail{}; for (int i = 1; i <= NumSpec; ++i) { diff --git a/integration/_parameters b/integration/_parameters index a11b5f3ea8..679e554598 100644 --- a/integration/_parameters +++ b/integration/_parameters @@ -93,6 +93,10 @@ retry_atol_enuc real -1 # are in [0, 1]? do_species_clip bool 1 +# in the corrector loop, do we check if the predicted state is +# valid (X > 0) before calling the RHS? +do_corrector_validation bool 0 + # flag for turning on the use of number densities for all species use_number_densities bool 0 From 21785072caaf61c5af5ebb413339389a9a406430 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 18:55:25 -0500 Subject: [PATCH 04/15] move the check --- integration/VODE/vode_dvnlsd.H | 46 +++++++++++++++++----------------- 1 file changed, 23 insertions(+), 23 deletions(-) diff --git a/integration/VODE/vode_dvnlsd.H b/integration/VODE/vode_dvnlsd.H index ad956f8ff7..11baf5563d 100644 --- a/integration/VODE/vode_dvnlsd.H +++ b/integration/VODE/vode_dvnlsd.H @@ -115,29 +115,6 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) // sometime VODE goes way off tangent. If these mass fractions // are really bad, then let's just bail now - if (integrator_rp::do_corrector_validation) { -#ifdef SDC - const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; - const amrex::Real thresh = species_failure_tolerance * rho_current; -#else - const amrex::Real thresh = species_failure_tolerance; -#endif - bool fail{}; - for (int i = 1; i <= NumSpec; ++i) { - if (vstate.y(i) < -thresh) { - fail = true; - break; - } - } - - if (fail) { - NFLAG = -1; - vstate.ICF = 2; - vstate.IPUP = 1; - return ACNRM; - } - } - #ifdef NEW_NETWORK_IMPLEMENTATION RHS::dgesl(vstate.jac, vstate.y); #else @@ -171,6 +148,29 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) vstate.y(i) = vstate.yh(i,1) + vstate.acor(i); } + if (integrator_rp::do_corrector_validation) { +#ifdef SDC + const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; + const amrex::Real thresh = species_failure_tolerance * rho_current; +#else + const amrex::Real thresh = species_failure_tolerance; +#endif + bool fail{}; + for (int i = 1; i <= NumSpec; ++i) { + if (vstate.y(i) < -thresh) { + fail = true; + break; + } + } + + if (fail) { + NFLAG = -1; + vstate.ICF = 2; + vstate.IPUP = 1; + return ACNRM; + } + } + // Test for convergence. If M > 0, an estimate of the convergence // rate constant is stored in CRATE, and this is used in the test. From b662c7f6887a9f11374d11731802f76b262fa24a Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 19:07:33 -0500 Subject: [PATCH 05/15] move the comment too --- integration/VODE/vode_dvnlsd.H | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/integration/VODE/vode_dvnlsd.H b/integration/VODE/vode_dvnlsd.H index 11baf5563d..0af50f9ec5 100644 --- a/integration/VODE/vode_dvnlsd.H +++ b/integration/VODE/vode_dvnlsd.H @@ -112,9 +112,6 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) (vstate.RL1 * vstate.yh(i,2) + vstate.acor(i)); } - // sometime VODE goes way off tangent. If these mass fractions - // are really bad, then let's just bail now - #ifdef NEW_NETWORK_IMPLEMENTATION RHS::dgesl(vstate.jac, vstate.y); #else @@ -148,6 +145,9 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) vstate.y(i) = vstate.yh(i,1) + vstate.acor(i); } + // sometime VODE goes way off tangent. If these mass fractions + // are really bad, then let's just bail now + if (integrator_rp::do_corrector_validation) { #ifdef SDC const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; From 7bb7d8c793ced815dbc70a7eb44fda48575209e4 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 19:16:31 -0500 Subject: [PATCH 06/15] update docs --- Docs/source/ode_integrators.rst | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/Docs/source/ode_integrators.rst b/Docs/source/ode_integrators.rst index 4644a347a5..c2a59fe16d 100644 --- a/Docs/source/ode_integrators.rst +++ b/Docs/source/ode_integrators.rst @@ -248,7 +248,7 @@ fractions. Looser than this can produce large errors. Controlling Species $\sum_k X_k = 1$ ==================================== -.. index:: integrator.renormalize_abundances, integrator.SMALL_X_SAFE, integrator.do_species_clip +.. index:: integrator.renormalize_abundances, integrator.SMALL_X_SAFE, integrator.do_species_clip, integrator.do_corrector_validation The ODE integrators don't know about the constraint that @@ -279,6 +279,12 @@ constraint on the intermediate states during the integration. This is enabled by default. +* ``integrator.do_corrector_validation`` : in the nonlinear solve + corrector loop, when we get a corrected integration state, do we + check to make sure the mass fractions are valid before calling the + righthand function? This is needed in some cases if + ``integrator.do_species_clip`` is disabled. Note: this is not + implemented for every integrator. Retry Mechanism From 97de4a31f5a80af339e25fcec80a136a4c80548a Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 15 Jan 2026 19:18:42 -0500 Subject: [PATCH 07/15] enable the new validation --- integration/_parameters | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/integration/_parameters b/integration/_parameters index 679e554598..dcfe82bdb7 100644 --- a/integration/_parameters +++ b/integration/_parameters @@ -95,7 +95,7 @@ do_species_clip bool 1 # in the corrector loop, do we check if the predicted state is # valid (X > 0) before calling the RHS? -do_corrector_validation bool 0 +do_corrector_validation bool 1 # flag for turning on the use of number densities for all species use_number_densities bool 0 From df57315140e0f12d3f39059c03ad79fdd8d685ea Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:19:16 -0500 Subject: [PATCH 08/15] don't do this for number densities also add some comments --- integration/VODE/vode_dvnlsd.H | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/integration/VODE/vode_dvnlsd.H b/integration/VODE/vode_dvnlsd.H index 0af50f9ec5..f080958a90 100644 --- a/integration/VODE/vode_dvnlsd.H +++ b/integration/VODE/vode_dvnlsd.H @@ -148,7 +148,7 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) // sometime VODE goes way off tangent. If these mass fractions // are really bad, then let's just bail now - if (integrator_rp::do_corrector_validation) { + if (integrator_rp::do_corrector_validation && !integrator_rp::use_number_densities) { #ifdef SDC const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO]; const amrex::Real thresh = species_failure_tolerance * rho_current; @@ -163,6 +163,10 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate) } } + // this resets the flags in the same fashion as is done above + // for a singular matrix. ICF = 2 means "unrecoverable error", + // IPUP = 1 forces a Jacobian reevaluation. The return value, + // ACNRM is ignored in this case (since NFLAG != 0) if (fail) { NFLAG = -1; vstate.ICF = 2; From 13b02e1323912b096881582d9ba3c342f27d3506 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:49:17 -0500 Subject: [PATCH 09/15] some udpates --- .../ci-benchmarks/ecsn_unit_test.out | 52 +++++------ .../ci-benchmarks/he-burn-19am_unit_test.out | 90 +++++++++---------- .../triple_alpha_plus_cago_FE_unit_test.out | 5 +- 3 files changed, 74 insertions(+), 73 deletions(-) diff --git a/unit_test/burn_cell/ci-benchmarks/ecsn_unit_test.out b/unit_test/burn_cell/ci-benchmarks/ecsn_unit_test.out index 75357037c4..cc5c01c849 100644 --- a/unit_test/burn_cell/ci-benchmarks/ecsn_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/ecsn_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (25.12-25-g10b67a795031)... -AMReX (25.12-25-g10b67a795031) initialized +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... Maximum Time (s): 0.01 @@ -30,37 +30,37 @@ RHS at t = 0 S32 1323519.174 ------------------------------------ successful? 1 - - Hnuc = 4.572458836e+19 - - added e = 4.572458836e+17 - - final T = 6740499652 + - Hnuc = 4.397909277e+19 + - added e = 4.397909277e+17 + - final T = 6656291303 ------------------------------------ e initial = 5.995956082e+17 -e final = 1.056841492e+18 +e final = 1.039386536e+18 ------------------------------------ new mass fractions: -H1 0.008875172349 -He4 9.999999998e-31 -O16 6.258606287e-07 -O20 0.009862665409 -F20 0.009861311457 -Ne20 8.532947059e-14 -Mg24 1.247632983e-06 -Al27 1.79296072e-18 -Si28 0.3055707459 -P31 9.999999998e-31 -S32 0.6658282314 +H1 0.009307048983 +He4 3.598820253e-13 +O16 6.945042998e-07 +O20 0.009997687144 +F20 0.009996314357 +Ne20 1.027891736e-13 +Mg24 1.24621466e-06 +Al27 1.645253981e-18 +Si28 0.3262434305 +P31 8.175490055e-14 +S32 0.6444535783 ------------------------------------ species creation rates: -omegadot(H1): -0.1124827651 +omegadot(H1): -0.06929510166 omegadot(He4): -3 -omegadot(O16): -49.99993741 -omegadot(O20): -0.01373345906 -omegadot(F20): -0.01386885428 +omegadot(O16): -49.99993055 +omegadot(O20): -0.0002312855859 +omegadot(F20): -0.0003685642773 omegadot(Ne20): -30 -omegadot(Mg24): -9.999875237 +omegadot(Mg24): -9.999875379 omegadot(Al27): -1 -omegadot(Si28): 29.55707459 +omegadot(Si28): 31.62434305 omegadot(P31): -1 -omegadot(S32): 65.58282314 -number of steps taken: 16498 -AMReX (25.12-25-g10b67a795031) finalized +omegadot(S32): 63.44535783 +number of steps taken: 598 +AMReX (25.11-47-gc4fcda1fb2ea) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/he-burn-19am_unit_test.out b/unit_test/burn_cell/ci-benchmarks/he-burn-19am_unit_test.out index 60af806226..ef2be38b72 100644 --- a/unit_test/burn_cell/ci-benchmarks/he-burn-19am_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/he-burn-19am_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (23.05-1038-gbd922c6216e0)... -AMReX (23.05-1038-gbd922c6216e0) initialized +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... Maximum Time (s): 1e-05 @@ -46,53 +46,53 @@ RHS at t = 0 Ni56 2.574224319e-29 ------------------------------------ successful? 1 - - Hnuc = 1.324552639e+23 - - added e = 1.324552639e+18 - - final T = 3539825699 + - Hnuc = 1.324544896e+23 + - added e = 1.324544896e+18 + - final T = 3539823486 ------------------------------------ e initial = 1.396711859e+18 -e final = 2.721264498e+18 +e final = 2.721256755e+18 ------------------------------------ new mass fractions: -H1 0.06140805092 -He4 0.001406317534 -C12 1.188197119e-07 -N13 7.909402052e-11 -N14 0.001127053623 -O16 6.931646535e-06 -Ne20 0.001267408673 -Na23 5.838358901e-09 -Mg24 0.001846400897 -Al27 1.108967039e-07 -Si28 0.3024943162 -P31 2.222585402e-06 -S32 0.3485154437 -Ar36 0.1639550483 -Ca40 0.1170346558 -Ti44 0.0008747967244 -Cr48 6.048749462e-05 -Fe52 6.288550582e-07 -Ni56 1.355779724e-09 +H1 0.06140746232 +He4 0.001405438558 +C12 1.188941922e-07 +N13 7.914334437e-11 +N14 0.001127796195 +O16 6.936214993e-06 +Ne20 0.001268241136 +Na23 5.834920863e-09 +Mg24 0.001847642975 +Al27 1.10886506e-07 +Si28 0.3024939412 +P31 2.221172998e-06 +S32 0.3484994557 +Ar36 0.1639525261 +Ca40 0.1170518309 +Ti44 0.0008751143826 +Cr48 6.0526545e-05 +Fe52 6.29481464e-07 +Ni56 1.357689393e-09 ------------------------------------ species creation rates: -omegadot(H1): -3859.194908 -omegadot(He4): -49859.36825 -omegadot(C12): -9999.988118 +omegadot(H1): -3859.253768 +omegadot(He4): -49859.45614 +omegadot(C12): -9999.988111 omegadot(N13): -9999.999992 -omegadot(N14): -9887.294638 -omegadot(O16): -9999.306835 -omegadot(Ne20): 126.7408673 -omegadot(Na23): 0.0005838358901 -omegadot(Mg24): 184.6400897 -omegadot(Al27): 0.01108967039 -omegadot(Si28): 30249.43162 -omegadot(P31): 0.2222585402 -omegadot(S32): 34851.54437 -omegadot(Ar36): 16395.50483 -omegadot(Ca40): 11703.46558 -omegadot(Ti44): 87.47967244 -omegadot(Cr48): 6.048749462 -omegadot(Fe52): 0.06288550582 -omegadot(Ni56): 0.0001355779724 -number of steps taken: 733 -AMReX (23.05-1038-gbd922c6216e0) finalized +omegadot(N14): -9887.220381 +omegadot(O16): -9999.306379 +omegadot(Ne20): 126.8241136 +omegadot(Na23): 0.0005834920863 +omegadot(Mg24): 184.7642975 +omegadot(Al27): 0.0110886506 +omegadot(Si28): 30249.39412 +omegadot(P31): 0.2221172998 +omegadot(S32): 34849.94557 +omegadot(Ar36): 16395.25261 +omegadot(Ca40): 11705.18309 +omegadot(Ti44): 87.51143826 +omegadot(Cr48): 6.0526545 +omegadot(Fe52): 0.0629481464 +omegadot(Ni56): 0.0001357689393 +number of steps taken: 679 +AMReX (25.11-47-gc4fcda1fb2ea) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/triple_alpha_plus_cago_FE_unit_test.out b/unit_test/burn_cell/ci-benchmarks/triple_alpha_plus_cago_FE_unit_test.out index d985a36dc4..74392a57b7 100644 --- a/unit_test/burn_cell/ci-benchmarks/triple_alpha_plus_cago_FE_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/triple_alpha_plus_cago_FE_unit_test.out @@ -1,4 +1,5 @@ -AMReX (23.07-7-g88f03408f18a) initialized +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized starting the single zone burn... Maximum Time (s): 0.001 State Density (g/cm^3): 1000000 @@ -33,4 +34,4 @@ omegadot(C12): 0.01497617246 omegadot(O16): 0.5494210118 omegadot(Fe56): 0 number of steps taken: 85669 -AMReX (23.07-7-g88f03408f18a) finalized +AMReX (25.11-47-gc4fcda1fb2ea) finalized From 0cf314bd5148e4dd7b597f7fb42d364060b2ab3a Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:51:19 -0500 Subject: [PATCH 10/15] more updates --- .../ci-benchmarks/aprox19_state_over_time.txt | 64 +++++++++++++++---- 1 file changed, 51 insertions(+), 13 deletions(-) diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt index 059c202f5b..5d9047dd8a 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt @@ -1,13 +1,51 @@ - # Time Density Temperature H1 He3 He4 C12 N14 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Fe54 Ni56 n p - 0 1e+06 3e+09 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 - 1e-10 1e+06 3.02858e+09 0.706587 0.00187808 0.215408 0.018277 0.0328286 0.0249997 2.17031e-05 3.97686e-09 1.44287e-11 1.21807e-15 1.84286e-20 1.24859e-25 1.03449e-30 1e-30 1e-30 1e-30 1e-30 2.32606e-30 2.32606e-30 -3.16228e-10 1e+06 3.03478e+09 0.705507 0.000621921 0.216227 0.0092854 0.043275 0.024999 8.50466e-05 4.81558e-08 4.59635e-11 5.40157e-15 2.46986e-19 5.24295e-24 5.63193e-30 1e-30 1e-30 1e-30 1e-30 5.98338e-30 5.98338e-30 - 1e-09 1e+06 3.03971e+09 0.704282 0.000198228 0.216432 0.00109041 0.0526569 0.0249968 0.000342898 6.19859e-07 5.31673e-10 5.54672e-14 5.58299e-18 3.12982e-22 7.99114e-28 1.00063e-30 1e-30 1e-30 1e-30 1.95324e-29 1.95324e-29 -3.16228e-09 1e+06 3.04074e+09 0.704146 6.28602e-05 0.216257 1.73328e-06 0.0533156 0.0249897 0.00121971 7.35522e-06 1.86066e-08 4.22894e-12 9.22256e-16 1.12401e-19 6.48983e-25 2.22129e-30 1.00001e-30 1e-30 1e-30 6.51843e-29 6.51843e-29 - 1e-08 1e+06 3.04176e+09 0.70416 2.04828e-05 0.215457 1e-30 0.0514256 0.0249673 0.00389241 7.69853e-05 6.38368e-07 4.73546e-10 3.40222e-13 1.3847e-16 2.72235e-21 1.7909e-26 1.06638e-30 1.00001e-30 1e-30 2.12345e-28 2.12345e-28 -3.16228e-08 1e+06 3.04471e+09 0.704164 6.34237e-06 0.212969 2.41781e-11 0.0459028 0.0248981 0.0113145 0.000725596 1.98421e-05 4.89394e-08 1.18275e-10 1.63701e-13 1.10302e-17 2.50276e-22 3.20374e-27 1.00897e-30 1e-30 6.99475e-28 6.99475e-28 - 1e-07 1e+06 3.05276e+09 0.704166 1.97272e-06 0.206101 2.10567e-11 0.0322273 0.0246918 0.0266706 0.00563202 0.00050569 3.92045e-06 2.8844e-08 1.18231e-10 2.29592e-14 1.4715e-18 5.18729e-23 3.90644e-28 1.07265e-30 2.4494e-27 2.46383e-27 -3.16228e-07 1e+06 3.07034e+09 0.704166 6.16355e-07 0.191015 1.67081e-11 0.0110228 0.0241129 0.0345787 0.0261422 0.00872611 0.000229534 5.49795e-06 7.15038e-08 4.2921e-11 8.40418e-15 8.83784e-19 1.93663e-23 3.5714e-27 9.93669e-27 7.27207e-25 - 1e-06 1e+06 3.09568e+09 0.704166 1.90977e-07 0.169097 1.16129e-11 0.000479845 0.0224668 0.00930122 0.0309807 0.0570159 0.00594684 0.000520989 2.41528e-05 5.01772e-08 3.40108e-11 1.22121e-14 9.05195e-19 1.63882e-22 4.89062e-26 3.35258e-20 -3.16228e-06 1e+06 3.11505e+09 0.704166 5.90985e-08 0.150118 8.16827e-12 5.39796e-08 0.0180596 0.000891114 0.00273936 0.0735111 0.0352666 0.0127917 0.00243636 1.92951e-05 4.85823e-08 6.32194e-11 1.6726e-14 2.98487e-18 3.33563e-25 6.19482e-16 - 1e-05 1e+06 3.14015e+09 0.704166 1.8682e-08 0.120635 4.27245e-12 4.83365e-09 0.00987352 0.000461095 0.00086113 0.0268146 0.036848 0.0509385 0.0478247 0.00156184 1.52568e-05 7.37586e-08 7.06632e-11 1.23934e-14 2.01895e-21 2.61715e-12 +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized +starting the single zone burn... +Maximum Time (s): 1e-05 +State Density (g/cm^3): 1000000 +State Temperature (K): 3000000000 +State rhoe (erg/cm^3): -1 +Mass Fraction (H1): 0.7 +Mass Fraction (He3): 0.025 +Mass Fraction (He4): 0.2 +Mass Fraction (C12): 0.025 +Mass Fraction (N14): 0.025 +Mass Fraction (O16): 0.025 +Mass Fraction (Ne20): 0 +Mass Fraction (Mg24): 0 +Mass Fraction (Si28): 0 +Mass Fraction (S32): 0 +Mass Fraction (Ar36): 0 +Mass Fraction (Ca40): 0 +Mass Fraction (Ti44): 0 +Mass Fraction (Cr48): 0 +Mass Fraction (Fe52): 0 +Mass Fraction (Fe54): 0 +Mass Fraction (Ni56): 0 +Mass Fraction (n): 0 +Mass Fraction (p): 0 +rho = 1000000 +T = 3000000000 +xn = 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 +(i, j, k) = 0 0 0 +y[SRHO] = 1000000 +y[SEINT] = 1.65414035784e+24 +y[SFS:] = 700000 25000 200000 25000 25000 25000 0 0 0 0 0 0 0 0 0 0 0 0 0 +ydot_a[SRHO] = 0 +ydot_a[SEINT] = 0 +ydot_a[SFS:] = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + +burning for dt = 1e-10 +burning for dt = 2.16227766017e-10 +burning for dt = 6.83772233983e-10 +burning for dt = 2.16227766017e-09 +burning for dt = 6.83772233983e-09 +burning for dt = 2.16227766017e-08 +burning for dt = 6.83772233983e-08 +burning for dt = 2.16227766017e-07 +burning for dt = 6.83772233983e-07 +burning for dt = 2.16227766017e-06 +burning for dt = 6.83772233983e-06 +successful? 1 +number of steps taken: 257 +AMReX (25.11-47-gc4fcda1fb2ea) finalized From fef048e71eee092981e184b60c79b02ace091c40 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:53:11 -0500 Subject: [PATCH 11/15] update bench --- .../aprox19_NSE_state_over_time.txt | 22 +++++++++---------- 1 file changed, 11 insertions(+), 11 deletions(-) diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_NSE_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_NSE_state_over_time.txt index acc229b06e..d8bc2940d4 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_NSE_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_NSE_state_over_time.txt @@ -1,13 +1,13 @@ # Time Density Temperature H1 He3 He4 C12 N14 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Fe54 Ni56 n p 0 1e+07 5e+09 0.1 0.025 0.8 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000909091 1e+07 5.36923e+09 1e-30 2.87875e-10 0.131064 1.04986e-06 8.09251e-11 2.71575e-06 4.15584e-08 9.43974e-06 0.00582282 0.00650057 0.00424621 0.00489184 0.00129753 0.0372752 0.058024 0.670993 0.0534183 1.43513e-06 0.0264521 -0.00115544 1e+07 5.36923e+09 1e-30 2.87875e-10 0.131064 1.04986e-06 8.0925e-11 2.71575e-06 4.15584e-08 9.43973e-06 0.00582282 0.00650057 0.00424621 0.00489184 0.00129753 0.0372752 0.058024 0.670993 0.0534183 1.43512e-06 0.0264521 -0.00146854 1e+07 5.36923e+09 1e-30 2.87874e-10 0.131064 1.04986e-06 8.09249e-11 2.71574e-06 4.15583e-08 9.43973e-06 0.00582282 0.00650057 0.00424621 0.00489185 0.00129753 0.0372752 0.058024 0.670993 0.0534184 1.43512e-06 0.0264521 -0.00186649 1e+07 5.36923e+09 1e-30 2.87874e-10 0.131064 1.04986e-06 8.09248e-11 2.71574e-06 4.15583e-08 9.43972e-06 0.00582282 0.00650057 0.00424621 0.00489185 0.00129753 0.0372752 0.058024 0.670993 0.0534184 1.43512e-06 0.0264521 -0.00237227 1e+07 5.36923e+09 1e-30 2.87873e-10 0.131064 1.04986e-06 8.09246e-11 2.71574e-06 4.15582e-08 9.43972e-06 0.00582282 0.00650057 0.00424621 0.00489185 0.00129753 0.0372752 0.058024 0.670993 0.0534186 1.43511e-06 0.0264521 -0.00301511 1e+07 5.36923e+09 1e-30 2.87873e-10 0.131063 1.04986e-06 8.09244e-11 2.71573e-06 4.15582e-08 9.4397e-06 0.00582282 0.00650058 0.00424621 0.00489186 0.00129753 0.0372751 0.058024 0.670993 0.0534187 1.43511e-06 0.0264521 -0.00383215 1e+07 5.36923e+09 1e-30 2.87872e-10 0.131063 1.04985e-06 8.09241e-11 2.71573e-06 4.15581e-08 9.43969e-06 0.00582282 0.00650058 0.00424622 0.00489186 0.00129753 0.0372751 0.058024 0.670993 0.0534189 1.4351e-06 0.0264521 - 0.0048706 1e+07 5.36923e+09 1e-30 2.87871e-10 0.131063 1.04985e-06 8.09237e-11 2.71572e-06 4.15579e-08 9.43967e-06 0.00582282 0.00650058 0.00424622 0.00489187 0.00129753 0.0372751 0.058024 0.670993 0.0534191 1.43509e-06 0.0264521 -0.00619044 1e+07 5.36923e+09 1e-30 2.8787e-10 0.131063 1.04985e-06 8.09232e-11 2.71571e-06 4.15578e-08 9.43965e-06 0.00582282 0.00650058 0.00424622 0.00489188 0.00129752 0.0372751 0.058024 0.670993 0.0534193 1.43508e-06 0.0264521 -0.00786793 1e+07 5.36923e+09 1e-30 2.87868e-10 0.131062 1.04984e-06 8.09226e-11 2.7157e-06 4.15576e-08 9.43963e-06 0.00582282 0.00650059 0.00424623 0.00489189 0.00129752 0.037275 0.0580239 0.670993 0.0534197 1.43507e-06 0.026452 - 0.01 1e+07 5.36923e+09 1e-30 2.87866e-10 0.131062 1.04983e-06 8.09218e-11 2.71569e-06 4.15573e-08 9.43959e-06 0.00582282 0.00650059 0.00424624 0.00489191 0.00129752 0.037275 0.0580239 0.670994 0.0534201 1.43505e-06 0.026452 +0.000909091 1e+07 5.16542e+09 0.10647 9.90099e-31 0.788419 1.34585e-05 0.0164549 0.0274406 0.0173537 0.0125208 0.0259179 0.00434822 0.00094121 0.000118497 1.87907e-06 1.24233e-07 5.57196e-09 9.31398e-11 7.88442e-15 3.34056e-14 3.48303e-12 +0.00115544 1e+07 5.38284e+09 1e-30 3.03736e-10 0.135293 1.10373e-06 8.70507e-11 2.82863e-06 4.35126e-08 9.70415e-06 0.00582948 0.00645566 0.00419711 0.00477662 0.00131915 0.0375917 0.0582282 0.669506 0.0501732 1.57018e-06 0.0266143 +0.00146854 1e+07 5.38284e+09 1e-30 3.03735e-10 0.135293 1.10373e-06 8.70506e-11 2.82862e-06 4.35126e-08 9.70415e-06 0.00582948 0.00645566 0.00419711 0.00477663 0.00131915 0.0375917 0.0582282 0.669506 0.0501733 1.57018e-06 0.0266143 +0.00186649 1e+07 5.38284e+09 1e-30 3.03735e-10 0.135293 1.10373e-06 8.70504e-11 2.82862e-06 4.35125e-08 9.70414e-06 0.00582948 0.00645566 0.00419711 0.00477663 0.00131915 0.0375917 0.0582281 0.669506 0.0501733 1.57018e-06 0.0266143 +0.00237227 1e+07 5.38284e+09 1e-30 3.03735e-10 0.135293 1.10373e-06 8.70502e-11 2.82862e-06 4.35124e-08 9.70413e-06 0.00582948 0.00645566 0.00419711 0.00477663 0.00131915 0.0375917 0.0582281 0.669507 0.0501734 1.57017e-06 0.0266143 +0.00301511 1e+07 5.38284e+09 1e-30 3.03734e-10 0.135292 1.10373e-06 8.705e-11 2.82861e-06 4.35124e-08 9.70412e-06 0.00582948 0.00645566 0.00419712 0.00477664 0.00131915 0.0375917 0.0582281 0.669507 0.0501736 1.57017e-06 0.0266143 +0.00383215 1e+07 5.38284e+09 1e-30 3.03733e-10 0.135292 1.10373e-06 8.70497e-11 2.82861e-06 4.35123e-08 9.70411e-06 0.00582948 0.00645566 0.00419712 0.00477664 0.00131915 0.0375917 0.0582281 0.669507 0.0501737 1.57016e-06 0.0266143 + 0.0048706 1e+07 5.38284e+09 1e-30 3.03732e-10 0.135292 1.10372e-06 8.70493e-11 2.8286e-06 4.35121e-08 9.70409e-06 0.00582948 0.00645567 0.00419712 0.00477665 0.00131915 0.0375916 0.0582281 0.669507 0.0501739 1.57015e-06 0.0266143 +0.00619044 1e+07 5.38284e+09 1e-30 3.03731e-10 0.135292 1.10372e-06 8.70487e-11 2.82859e-06 4.3512e-08 9.70407e-06 0.00582947 0.00645567 0.00419712 0.00477666 0.00131915 0.0375916 0.0582281 0.669507 0.0501742 1.57014e-06 0.0266143 +0.00786793 1e+07 5.38284e+09 1e-30 3.03729e-10 0.135291 1.10371e-06 8.70481e-11 2.82858e-06 4.35118e-08 9.70404e-06 0.00582947 0.00645567 0.00419713 0.00477667 0.00131915 0.0375916 0.0582281 0.669507 0.0501745 1.57013e-06 0.0266142 + 0.01 1e+07 5.38284e+09 1e-30 3.03727e-10 0.135291 1.1037e-06 8.70473e-11 2.82856e-06 4.35115e-08 9.704e-06 0.00582947 0.00645568 0.00419714 0.00477669 0.00131914 0.0375915 0.058228 0.669507 0.0501749 1.57011e-06 0.0266142 From f61b3dc7337050945f2fd2112257ac38a1801409 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:56:07 -0500 Subject: [PATCH 12/15] remove whitespace --- .../burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt index 5d9047dd8a..7c367a5037 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt @@ -26,14 +26,14 @@ Mass Fraction (n): 0 Mass Fraction (p): 0 rho = 1000000 T = 3000000000 -xn = 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 +xn = 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 (i, j, k) = 0 0 0 y[SRHO] = 1000000 y[SEINT] = 1.65414035784e+24 -y[SFS:] = 700000 25000 200000 25000 25000 25000 0 0 0 0 0 0 0 0 0 0 0 0 0 +y[SFS:] = 700000 25000 200000 25000 25000 25000 0 0 0 0 0 0 0 0 0 0 0 0 0 ydot_a[SRHO] = 0 ydot_a[SEINT] = 0 -ydot_a[SFS:] = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +ydot_a[SFS:] = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 burning for dt = 1e-10 burning for dt = 2.16227766017e-10 From 99118a67004d43ed177bf3d0e2f15b5574388257 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 10:58:41 -0500 Subject: [PATCH 13/15] fix benchmark --- .../ci-benchmarks/aprox19_state_over_time.txt | 64 ++++--------------- 1 file changed, 13 insertions(+), 51 deletions(-) diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt index 7c367a5037..e778a293b7 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_state_over_time.txt @@ -1,51 +1,13 @@ -Initializing AMReX (25.11-47-gc4fcda1fb2ea)... -AMReX (25.11-47-gc4fcda1fb2ea) initialized -starting the single zone burn... -Maximum Time (s): 1e-05 -State Density (g/cm^3): 1000000 -State Temperature (K): 3000000000 -State rhoe (erg/cm^3): -1 -Mass Fraction (H1): 0.7 -Mass Fraction (He3): 0.025 -Mass Fraction (He4): 0.2 -Mass Fraction (C12): 0.025 -Mass Fraction (N14): 0.025 -Mass Fraction (O16): 0.025 -Mass Fraction (Ne20): 0 -Mass Fraction (Mg24): 0 -Mass Fraction (Si28): 0 -Mass Fraction (S32): 0 -Mass Fraction (Ar36): 0 -Mass Fraction (Ca40): 0 -Mass Fraction (Ti44): 0 -Mass Fraction (Cr48): 0 -Mass Fraction (Fe52): 0 -Mass Fraction (Fe54): 0 -Mass Fraction (Ni56): 0 -Mass Fraction (n): 0 -Mass Fraction (p): 0 -rho = 1000000 -T = 3000000000 -xn = 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 -(i, j, k) = 0 0 0 -y[SRHO] = 1000000 -y[SEINT] = 1.65414035784e+24 -y[SFS:] = 700000 25000 200000 25000 25000 25000 0 0 0 0 0 0 0 0 0 0 0 0 0 -ydot_a[SRHO] = 0 -ydot_a[SEINT] = 0 -ydot_a[SFS:] = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 - -burning for dt = 1e-10 -burning for dt = 2.16227766017e-10 -burning for dt = 6.83772233983e-10 -burning for dt = 2.16227766017e-09 -burning for dt = 6.83772233983e-09 -burning for dt = 2.16227766017e-08 -burning for dt = 6.83772233983e-08 -burning for dt = 2.16227766017e-07 -burning for dt = 6.83772233983e-07 -burning for dt = 2.16227766017e-06 -burning for dt = 6.83772233983e-06 -successful? 1 -number of steps taken: 257 -AMReX (25.11-47-gc4fcda1fb2ea) finalized + # Time Density Temperature H1 He3 He4 C12 N14 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Fe54 Ni56 n p + 0 1e+06 3e+09 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 + 1e-10 1e+06 3.02858e+09 0.706587 0.00187808 0.215408 0.018277 0.0328286 0.0249997 2.17031e-05 3.97686e-09 1.44287e-11 1.21807e-15 1.84286e-20 1.24858e-25 1e-30 1e-30 1e-30 1e-30 1e-30 1.32606e-30 1.32606e-30 +3.16228e-10 1e+06 3.03478e+09 0.705507 0.000621921 0.216227 0.0092854 0.043275 0.024999 8.50466e-05 4.81558e-08 4.59635e-11 5.40157e-15 2.46986e-19 5.24295e-24 5.59745e-30 1e-30 1e-30 1e-30 1e-30 4.98338e-30 4.98338e-30 + 1e-09 1e+06 3.03971e+09 0.704282 0.000198228 0.216432 0.00109041 0.0526569 0.0249968 0.000342898 6.19859e-07 5.31673e-10 5.54672e-14 5.58299e-18 3.12982e-22 7.9908e-28 1.00063e-30 1e-30 1e-30 1e-30 1.85324e-29 1.85324e-29 +3.16228e-09 1e+06 3.04074e+09 0.704146 6.28602e-05 0.216257 1.73328e-06 0.0533156 0.0249897 0.00121971 7.35522e-06 1.86066e-08 4.22894e-12 9.22256e-16 1.12401e-19 6.48983e-25 2.22129e-30 1.00001e-30 1e-30 1e-30 6.41843e-29 6.41843e-29 + 1e-08 1e+06 3.04176e+09 0.70416 1.99373e-05 0.215457 1.27512e-09 0.0514254 0.0249673 0.00389243 7.6966e-05 6.42582e-07 4.83588e-10 3.54601e-13 1.4788e-16 2.98481e-21 2.01631e-26 1.07665e-30 1.00001e-30 1e-30 2.11345e-28 2.11345e-28 +3.16228e-08 1e+06 3.04471e+09 0.704165 6.27825e-06 0.212969 2.33419e-11 0.0459026 0.0248981 0.0113143 0.000725634 1.98576e-05 4.90582e-08 1.18872e-10 1.65113e-13 1.11748e-17 2.54889e-22 3.28266e-27 1.00926e-30 1e-30 6.98481e-28 6.98481e-28 + 1e-07 1e+06 3.05275e+09 0.704166 1.96625e-06 0.206101 2.10567e-11 0.0322271 0.0246918 0.0266704 0.005632 0.000505708 3.92091e-06 2.88513e-08 1.18284e-10 2.29751e-14 1.47295e-18 5.19417e-23 3.91304e-28 1.07278e-30 2.44839e-27 2.46285e-27 +3.16228e-07 1e+06 3.07034e+09 0.704166 6.15727e-07 0.191015 1.67081e-11 0.0110228 0.0241129 0.0345785 0.0261421 0.00872608 0.000229535 5.498e-06 7.15053e-08 4.29225e-11 8.40464e-15 8.83852e-19 1.93683e-23 3.57177e-27 9.93564e-27 7.27279e-25 + 1e-06 1e+06 3.09568e+09 0.704167 1.90917e-07 0.169097 1.16129e-11 0.000479842 0.0224668 0.00930117 0.0309806 0.0570157 0.00594682 0.000520987 2.41527e-05 5.01772e-08 3.40108e-11 1.22121e-14 9.05196e-19 1.63883e-22 4.89048e-26 3.35258e-20 +3.16228e-06 1e+06 3.11505e+09 0.704167 5.90928e-08 0.150119 8.16832e-12 5.39794e-08 0.0180596 0.000891114 0.00273935 0.0735108 0.0352664 0.0127917 0.00243635 1.92951e-05 4.85822e-08 6.32192e-11 1.6726e-14 2.98487e-18 3.33559e-25 6.1948e-16 + 1e-05 1e+06 3.14015e+09 0.704167 1.86815e-08 0.120635 4.27248e-12 4.83368e-09 0.0098735 0.000461094 0.00086113 0.0268144 0.0368478 0.0509383 0.0478246 0.00156184 1.52568e-05 7.37585e-08 7.06631e-11 1.23934e-14 2.01893e-21 2.61715e-12 From 18cae4df40270963f96c06d44389e2bcc77c0549 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 11:45:52 -0500 Subject: [PATCH 14/15] more updates --- .../ci-benchmarks/aprox13_RKC_unit_test.out | 6 +- .../he-burn-19am_BE_unit_test.out | 90 +++++++++---------- 2 files changed, 48 insertions(+), 48 deletions(-) diff --git a/unit_test/burn_cell/ci-benchmarks/aprox13_RKC_unit_test.out b/unit_test/burn_cell/ci-benchmarks/aprox13_RKC_unit_test.out index 3d2073ec0e..968c557982 100644 --- a/unit_test/burn_cell/ci-benchmarks/aprox13_RKC_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/aprox13_RKC_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (23.05-658-g0165b6743355)... -AMReX (23.05-658-g0165b6743355) initialized +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized starting the single zone burn... Maximum Time (s): 0.01 State Density (g/cm^3): 1000000 @@ -70,4 +70,4 @@ omegadot(Cr48): 3.272104685e-12 omegadot(Fe52): 7.126141581e-18 omegadot(Ni56): 1.22471994e-24 number of steps taken: 255 -AMReX (23.05-658-g0165b6743355) finalized +AMReX (25.11-47-gc4fcda1fb2ea) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/he-burn-19am_BE_unit_test.out b/unit_test/burn_cell/ci-benchmarks/he-burn-19am_BE_unit_test.out index 77836e824e..cecc7cec15 100644 --- a/unit_test/burn_cell/ci-benchmarks/he-burn-19am_BE_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/he-burn-19am_BE_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (23.05-1038-gbd922c6216e0)... -AMReX (23.05-1038-gbd922c6216e0) initialized +Initializing AMReX (25.11-47-gc4fcda1fb2ea)... +AMReX (25.11-47-gc4fcda1fb2ea) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... Maximum Time (s): 1e-05 @@ -46,53 +46,53 @@ RHS at t = 0 Ni56 2.574224319e-29 ------------------------------------ successful? 1 - - Hnuc = 1.324591682e+23 - - added e = 1.324591682e+18 - - final T = 3539838349 + - Hnuc = 1.324591986e+23 + - added e = 1.324591986e+18 + - final T = 3539838445 ------------------------------------ e initial = 1.396711859e+18 -e final = 2.721303541e+18 +e final = 2.721303845e+18 ------------------------------------ new mass fractions: -H1 0.0614041578 -He4 0.001390967515 -C12 1.194150996e-07 -N13 7.948320928e-11 -N14 0.001132572545 -O16 6.939462514e-06 -Ne20 0.001268664041 -Na23 5.77734127e-09 -Mg24 0.001851928121 -Al27 1.100081874e-07 -Si28 0.3024482452 -P31 2.19817828e-06 -S32 0.3484511353 -Ar36 0.1639686458 -Ca40 0.1171367255 -Ti44 0.0008763032297 -Cr48 6.064947631e-05 -Fe52 6.31118656e-07 -Ni56 1.361868236e-09 +H1 0.06140416809 +He4 0.001390777362 +C12 1.194102635e-07 +N13 7.947998734e-11 +N14 0.001132527207 +O16 6.938893863e-06 +Ne20 0.001268559515 +Na23 5.77607584e-09 +Mg24 0.001851806271 +Al27 1.099859099e-07 +Si28 0.3024474385 +P31 2.197872323e-06 +S32 0.3484525859 +Ar36 0.1639692923 +Ca40 0.117135911 +Ti44 0.0008762828478 +Cr48 6.064662425e-05 +Fe52 6.310696052e-07 +Ni56 1.36171123e-09 ------------------------------------ species creation rates: -omegadot(H1): -3859.58422 -omegadot(He4): -49860.90325 -omegadot(C12): -9999.988058 +omegadot(H1): -3859.583191 +omegadot(He4): -49860.92226 +omegadot(C12): -9999.988059 omegadot(N13): -9999.999992 -omegadot(N14): -9886.742745 -omegadot(O16): -9999.306054 -omegadot(Ne20): 126.8664041 -omegadot(Na23): 0.000577734127 -omegadot(Mg24): 185.1928121 -omegadot(Al27): 0.01100081874 -omegadot(Si28): 30244.82452 -omegadot(P31): 0.219817828 -omegadot(S32): 34845.11353 -omegadot(Ar36): 16396.86458 -omegadot(Ca40): 11713.67255 -omegadot(Ti44): 87.63032297 -omegadot(Cr48): 6.064947631 -omegadot(Fe52): 0.0631118656 -omegadot(Ni56): 0.0001361868236 -number of steps taken: 9720 -AMReX (23.05-1038-gbd922c6216e0) finalized +omegadot(N14): -9886.747279 +omegadot(O16): -9999.306111 +omegadot(Ne20): 126.8559515 +omegadot(Na23): 0.000577607584 +omegadot(Mg24): 185.1806271 +omegadot(Al27): 0.01099859099 +omegadot(Si28): 30244.74385 +omegadot(P31): 0.2197872323 +omegadot(S32): 34845.25859 +omegadot(Ar36): 16396.92923 +omegadot(Ca40): 11713.5911 +omegadot(Ti44): 87.62828478 +omegadot(Cr48): 6.064662425 +omegadot(Fe52): 0.06310696052 +omegadot(Ni56): 0.000136171123 +number of steps taken: 9678 +AMReX (25.11-47-gc4fcda1fb2ea) finalized From ff6af3243f8657d9f03896cade50b126009c1fbf Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 16 Jan 2026 11:48:10 -0500 Subject: [PATCH 15/15] more updates --- .../aprox19_BE_state_over_time.txt | 20 +++++++++---------- .../aprox19_RKC_state_over_time.txt | 20 +++++++++---------- 2 files changed, 20 insertions(+), 20 deletions(-) diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_BE_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_BE_state_over_time.txt index 2b4ca43587..d450cbb94d 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_BE_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_BE_state_over_time.txt @@ -1,13 +1,13 @@ # Time Density Temperature H1 He3 He4 C12 N14 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Fe54 Ni56 n p 0 1e+06 3e+09 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 - 1e-10 1e+06 3.02856e+09 0.706582 0.00189497 0.215397 0.0182827 0.0328219 0.0249997 2.17185e-05 4.00707e-09 1.44235e-11 1.21931e-15 1.85797e-20 1.27593e-25 1.03595e-30 1e-30 1e-30 1e-30 1e-30 2.32624e-30 2.32624e-30 -3.16228e-10 1e+06 3.03476e+09 0.705509 0.000628581 0.216222 0.00931244 0.0432434 0.024999 8.512e-05 4.85269e-08 4.61301e-11 5.41599e-15 2.49264e-19 5.35054e-24 5.79192e-30 1e-30 1e-30 1e-30 1e-30 5.98461e-30 5.98461e-30 - 1e-09 1e+06 3.03969e+09 0.704286 0.000200506 0.21643 0.00111892 0.0526236 0.0249968 0.000342966 6.23243e-07 5.38214e-10 5.62735e-14 5.68043e-18 3.2007e-22 8.23321e-28 1.00066e-30 1e-30 4.24516e-28 1e-30 1.95328e-29 1.95328e-29 -3.16228e-09 1e+06 3.04074e+09 0.704146 6.42104e-05 0.216256 2.27397e-06 0.0533151 0.0249897 0.00121946 7.46561e-06 1.94541e-08 4.61129e-12 1.05971e-15 1.37192e-19 8.46684e-25 2.71189e-30 1.00001e-30 4.24516e-28 1e-30 6.51863e-29 6.51863e-29 - 1e-08 1e+06 3.04176e+09 0.70416 2.09463e-05 0.215457 2.67096e-11 0.0514278 0.0249673 0.00388704 7.92663e-05 6.95194e-07 5.53831e-10 4.29729e-13 1.88314e-16 3.94438e-21 2.71834e-26 1.10323e-30 4.24516e-28 1e-30 2.12447e-28 2.12447e-28 -3.16228e-08 1e+06 3.04471e+09 0.704165 6.50597e-06 0.212971 2.32879e-11 0.0459107 0.0248981 0.0112961 0.000733137 2.01775e-05 4.919e-08 1.14434e-10 1.47598e-13 8.93693e-18 1.75295e-22 1.86294e-27 4.24521e-28 1e-30 6.99839e-28 6.9984e-28 - 1e-07 1e+06 3.05275e+09 0.704166 2.0132e-06 0.206109 2.10593e-11 0.032257 0.0246919 0.0266086 0.00565053 0.000510725 3.96206e-06 2.88471e-08 1.15225e-10 2.13979e-14 1.2834e-18 4.13398e-23 7.01078e-28 1.05156e-30 2.45097e-27 2.46121e-27 -3.16228e-07 1e+06 3.07032e+09 0.704166 6.23618e-07 0.191035 1.67132e-11 0.0110708 0.024113 0.0345257 0.0260972 0.00875392 0.000232005 5.6065e-06 7.35599e-08 4.44761e-11 8.74679e-15 9.19956e-19 2.00533e-23 3.69731e-27 9.94337e-27 7.52641e-25 - 1e-06 1e+06 3.09565e+09 0.704167 1.92199e-07 0.169122 1.16179e-11 0.000490962 0.022467 0.00934994 0.0309577 0.0569389 0.00595735 0.000524875 2.45272e-05 5.14426e-08 3.52397e-11 1.27949e-14 9.58906e-19 1.73607e-22 4.89224e-26 3.55151e-20 -3.16228e-06 1e+06 3.11504e+09 0.704167 5.937e-08 0.150126 8.16958e-12 7.84999e-08 0.0180619 0.000893683 0.00277257 0.0734901 0.0352149 0.0128016 0.00245239 1.95824e-05 4.98184e-08 6.56227e-11 1.7601e-14 3.14095e-18 3.3984e-25 6.51888e-16 + 1e-10 1e+06 3.02856e+09 0.706582 0.00189497 0.215397 0.0182827 0.0328219 0.0249997 2.17185e-05 4.00707e-09 1.44235e-11 1.21931e-15 1.85797e-20 1.27592e-25 1e-30 1e-30 1e-30 1e-30 1e-30 1.32624e-30 1.32624e-30 +3.16228e-10 1e+06 3.03476e+09 0.705509 0.000628581 0.216222 0.00931244 0.0432434 0.024999 8.512e-05 4.85269e-08 4.61301e-11 5.41599e-15 2.49264e-19 5.35054e-24 5.75597e-30 1e-30 1e-30 1e-30 1e-30 4.98461e-30 4.98461e-30 + 1e-09 1e+06 3.03969e+09 0.704286 0.000200506 0.21643 0.00111892 0.0526236 0.0249968 0.000342966 6.23243e-07 5.38214e-10 5.62735e-14 5.68043e-18 3.2007e-22 8.23285e-28 1.00066e-30 1e-30 4.24516e-28 1e-30 1.85328e-29 1.85328e-29 +3.16228e-09 1e+06 3.04074e+09 0.704146 6.42104e-05 0.216256 2.27397e-06 0.0533151 0.0249897 0.00121946 7.46561e-06 1.94541e-08 4.61129e-12 1.05971e-15 1.37192e-19 8.46684e-25 2.71189e-30 1.00001e-30 4.24516e-28 1e-30 6.41863e-29 6.41863e-29 + 1e-08 1e+06 3.04176e+09 0.70416 2.09463e-05 0.215457 2.67096e-11 0.0514278 0.0249673 0.00388704 7.92663e-05 6.95194e-07 5.53831e-10 4.29729e-13 1.88314e-16 3.94438e-21 2.71834e-26 1.10323e-30 4.24516e-28 1e-30 2.11447e-28 2.11447e-28 +3.16228e-08 1e+06 3.04471e+09 0.704165 6.5068e-06 0.212971 2.32879e-11 0.0459108 0.0248981 0.011296 0.000733175 2.01812e-05 4.92086e-08 1.14507e-10 1.47741e-13 8.94913e-18 1.75615e-22 1.86732e-27 4.2452e-28 1e-30 6.98842e-28 6.98842e-28 + 1e-07 1e+06 3.05275e+09 0.704166 2.01328e-06 0.206109 2.10593e-11 0.0322571 0.0246919 0.0266085 0.00565055 0.000510734 3.96223e-06 2.88493e-08 1.15238e-10 2.14014e-14 1.28369e-18 4.13518e-23 7.01178e-28 1.05158e-30 2.44997e-27 2.46022e-27 +3.16228e-07 1e+06 3.07032e+09 0.704166 6.23626e-07 0.191035 1.67132e-11 0.0110708 0.024113 0.0345257 0.0260972 0.00875393 0.000232006 5.60654e-06 7.35607e-08 4.44768e-11 8.74697e-15 9.19981e-19 2.0054e-23 3.69744e-27 9.94237e-27 7.52666e-25 + 1e-06 1e+06 3.09565e+09 0.704167 1.922e-07 0.169122 1.16179e-11 0.000490963 0.022467 0.00934994 0.0309577 0.0569389 0.00595735 0.000524875 2.45272e-05 5.14427e-08 3.52398e-11 1.2795e-14 9.58909e-19 1.73608e-22 4.89214e-26 3.55152e-20 +3.16228e-06 1e+06 3.11504e+09 0.704167 5.93701e-08 0.150126 8.16958e-12 7.85e-08 0.0180619 0.000893683 0.00277257 0.0734901 0.0352149 0.0128016 0.00245239 1.95824e-05 4.98184e-08 6.56227e-11 1.7601e-14 3.14095e-18 3.39839e-25 6.51888e-16 1e-05 1e+06 3.14013e+09 0.704167 1.8749e-08 0.120661 4.27519e-12 4.8354e-09 0.00988207 0.000461506 0.000861871 0.0268741 0.0368383 0.0508397 0.0478301 0.00156971 1.54457e-05 7.53675e-08 7.2984e-11 1.28002e-14 2.10934e-21 2.70311e-12 diff --git a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_RKC_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_RKC_state_over_time.txt index 8fe4376100..3c58f9f115 100644 --- a/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_RKC_state_over_time.txt +++ b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_RKC_state_over_time.txt @@ -1,13 +1,13 @@ # Time Density Temperature H1 He3 He4 C12 N14 O16 Ne20 Mg24 Si28 S32 Ar36 Ca40 Ti44 Cr48 Fe52 Fe54 Ni56 n p 0 1e+06 3e+09 0.7 0.025 0.2 0.025 0.025 0.025 0 0 0 0 0 0 0 0 0 0 0 0 0 - 1e-10 1e+06 3.02857e+09 0.706585 0.00188489 0.215404 0.0182772 0.0328284 0.0249997 2.17031e-05 3.97592e-09 1.44272e-11 1.21802e-15 1.84217e-20 1.24506e-25 1.03418e-30 1e-30 1e-30 1e-30 1e-30 2.32582e-30 2.32582e-30 -3.16228e-10 1e+06 3.03477e+09 0.705506 0.000626042 0.216224 0.00928939 0.0432703 0.024999 8.50563e-05 4.81208e-08 4.58849e-11 5.39715e-15 2.46511e-19 5.20229e-24 5.5281e-30 1e-30 1e-30 1e-30 1e-30 5.983e-30 5.983e-30 - 1e-09 1e+06 3.0397e+09 0.704284 0.000199673 0.216431 0.00110358 0.0526416 0.0249968 0.000342932 6.19714e-07 5.29897e-10 5.50639e-14 5.52076e-18 3.07859e-22 7.78601e-28 1.00061e-30 1e-30 1e-30 1e-30 1.95319e-29 1.95319e-29 -3.16228e-09 1e+06 3.04074e+09 0.704146 6.3294e-05 0.216257 2.38193e-06 0.0533148 0.0249897 0.00121971 7.35457e-06 1.85457e-08 4.17808e-12 8.96377e-16 1.06551e-19 5.94637e-25 2.0718e-30 1.00001e-30 1e-30 1e-30 6.51819e-29 6.51819e-29 - 1e-08 1e+06 3.04176e+09 0.70416 2.05221e-05 0.215457 2.39436e-07 0.0514251 0.0249673 0.00389247 7.69714e-05 6.27047e-07 4.41924e-10 2.8774e-13 1.00103e-16 1.57253e-21 7.67348e-27 1.01949e-30 1.00001e-30 1e-30 2.1234e-28 2.1234e-28 -3.16228e-08 1e+06 3.04471e+09 0.704164 6.68022e-06 0.212968 1.03434e-08 0.0459013 0.0248981 0.0113164 0.000725848 1.91918e-05 4.34088e-08 9.00858e-11 9.93813e-14 4.9269e-18 7.56624e-23 6.02336e-28 1.00099e-30 1e-30 6.99356e-28 6.99356e-28 - 1e-07 1e+06 3.05275e+09 0.704166 2.03484e-06 0.206102 3.52636e-11 0.0322322 0.0246918 0.0266555 0.00564306 0.000503599 3.81369e-06 2.67434e-08 1.01326e-10 1.7554e-14 9.64649e-19 2.79266e-23 1.65534e-28 1.03067e-30 2.44921e-27 2.4553e-27 -3.16228e-07 1e+06 3.07033e+09 0.704166 6.2411e-07 0.191017 1.67296e-11 0.0110301 0.0241129 0.034563 0.026145 0.00873 0.000229225 5.45954e-06 7.01956e-08 4.13428e-11 7.86914e-15 7.95554e-19 1.65483e-23 3.05149e-27 9.93618e-27 6.22838e-25 - 1e-06 1e+06 3.09567e+09 0.704167 1.91871e-07 0.1691 1.16134e-11 0.000483324 0.0224667 0.00930451 0.0309675 0.0570185 0.00594801 0.000520809 2.41023e-05 4.98998e-08 3.36261e-11 1.19663e-14 8.7559e-19 1.58524e-22 4.89047e-26 3.24293e-20 -3.16228e-06 1e+06 3.11504e+09 0.704167 5.92136e-08 0.15012 8.16857e-12 7.3465e-08 0.0180596 0.00089227 0.00274651 0.0734999 0.0352657 0.0127937 0.00243611 1.92743e-05 4.84326e-08 6.28099e-11 1.65304e-14 2.95002e-18 3.31773e-25 6.12237e-16 + 1e-10 1e+06 3.02857e+09 0.706585 0.00188489 0.215404 0.0182772 0.0328284 0.0249997 2.17031e-05 3.97592e-09 1.44272e-11 1.21802e-15 1.84217e-20 1.24505e-25 1e-30 1e-30 1e-30 1e-30 1e-30 1.32586e-30 1.32586e-30 +3.16228e-10 1e+06 3.03477e+09 0.705506 0.000626042 0.216224 0.00928939 0.0432703 0.024999 8.50563e-05 4.81208e-08 4.58849e-11 5.39715e-15 2.46511e-19 5.20229e-24 5.49392e-30 1e-30 1e-30 1e-30 1e-30 4.98305e-30 4.98305e-30 + 1e-09 1e+06 3.0397e+09 0.704284 0.000199673 0.216431 0.00110358 0.0526416 0.0249968 0.000342932 6.19714e-07 5.29897e-10 5.50639e-14 5.52076e-18 3.07859e-22 7.78566e-28 1.00061e-30 1e-30 1e-30 1e-30 1.8532e-29 1.8532e-29 +3.16228e-09 1e+06 3.04074e+09 0.704146 6.3294e-05 0.216257 2.38193e-06 0.0533148 0.0249897 0.00121971 7.35457e-06 1.85457e-08 4.17808e-12 8.96377e-16 1.06551e-19 5.94637e-25 2.0718e-30 1.00001e-30 1e-30 1e-30 6.4182e-29 6.4182e-29 + 1e-08 1e+06 3.04176e+09 0.70416 2.05221e-05 0.215457 2.39436e-07 0.0514251 0.0249673 0.00389247 7.69714e-05 6.27047e-07 4.41924e-10 2.8774e-13 1.00103e-16 1.57253e-21 7.67348e-27 1.01949e-30 1.00001e-30 1e-30 2.1134e-28 2.1134e-28 +3.16228e-08 1e+06 3.04471e+09 0.704164 6.68022e-06 0.212968 1.03434e-08 0.0459013 0.0248981 0.0113164 0.000725848 1.91918e-05 4.34088e-08 9.00858e-11 9.93813e-14 4.9269e-18 7.56624e-23 6.02336e-28 1.00099e-30 1e-30 6.98356e-28 6.98356e-28 + 1e-07 1e+06 3.05275e+09 0.704166 2.03484e-06 0.206102 3.52636e-11 0.0322322 0.0246918 0.0266555 0.00564306 0.000503599 3.81369e-06 2.67434e-08 1.01326e-10 1.7554e-14 9.64649e-19 2.79266e-23 1.65534e-28 1.03066e-30 2.44821e-27 2.4543e-27 +3.16228e-07 1e+06 3.07033e+09 0.704166 6.2411e-07 0.191017 1.67296e-11 0.0110301 0.0241129 0.034563 0.026145 0.00873 0.000229225 5.45954e-06 7.01956e-08 4.13428e-11 7.86914e-15 7.95554e-19 1.65483e-23 3.05149e-27 9.93518e-27 6.22837e-25 + 1e-06 1e+06 3.09567e+09 0.704167 1.91871e-07 0.1691 1.16134e-11 0.000483324 0.0224667 0.00930451 0.0309675 0.0570185 0.00594801 0.000520809 2.41023e-05 4.98998e-08 3.36261e-11 1.19663e-14 8.7559e-19 1.58524e-22 4.89037e-26 3.24293e-20 +3.16228e-06 1e+06 3.11504e+09 0.704167 5.92136e-08 0.15012 8.16857e-12 7.3465e-08 0.0180596 0.00089227 0.00274651 0.0734999 0.0352657 0.0127937 0.00243611 1.92743e-05 4.84326e-08 6.28099e-11 1.65304e-14 2.95002e-18 3.31772e-25 6.12237e-16 1e-05 1e+06 3.14015e+09 0.704167 1.8699e-08 0.120635 4.2725e-12 4.83322e-09 0.00987394 0.000461115 0.000861121 0.0268198 0.0368385 0.0509335 0.0478329 0.00156179 1.52423e-05 7.35358e-08 7.01785e-11 1.23088e-14 1.99059e-21 2.5992e-12