diff --git a/EOS/helmholtz/actual_eos.H b/EOS/helmholtz/actual_eos.H index 62e8b03a3a..2a7eb64cd3 100644 --- a/EOS/helmholtz/actual_eos.H +++ b/EOS/helmholtz/actual_eos.H @@ -686,8 +686,8 @@ void apply_coulomb_corrections (T& state) amrex::Real dsdd = z * dxnidd; [[maybe_unused]] amrex::Real dsda = z * dxnida; - amrex::Real lami = 1.0e0_rt / std::pow(s, onethird); - amrex::Real inv_lami = 1.0e0_rt / lami; + amrex::Real inv_lami = std::cbrt(s); + amrex::Real lami = 1.0e0_rt / inv_lami; z = -onethird * lami; amrex::Real lamidd = z * dsdd / s; [[maybe_unused]] amrex::Real lamida = z * dsda / s; diff --git a/EOS/ztwd/actual_eos.H b/EOS/ztwd/actual_eos.H index 533b92d71a..222e4fe2db 100644 --- a/EOS/ztwd/actual_eos.H +++ b/EOS/ztwd/actual_eos.H @@ -176,7 +176,7 @@ void actual_eos (I input, T& state) // Solve for the pressure and energy: - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); pres = pressure(x); eint = enth - pres / dens; @@ -188,7 +188,7 @@ void actual_eos (I input, T& state) // Solve for the pressure, energy and enthalpy: - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); pres = pressure(x); enth = enthalpy(x, B); eint = enth - pres / dens; @@ -203,7 +203,7 @@ void actual_eos (I input, T& state) pres_iter(pres, dens, B); - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); enth = enthalpy(x, B); eint = enth - pres / dens; @@ -215,7 +215,7 @@ void actual_eos (I input, T& state) // Solve for the enthalpy and energy: - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); enth = enthalpy(x, B); eint = enth - pres / dens; @@ -227,7 +227,7 @@ void actual_eos (I input, T& state) // Solve for the pressure and enthalpy: - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); pres = pressure(x); enth = enthalpy(x, B); @@ -241,7 +241,7 @@ void actual_eos (I input, T& state) pres_iter(pres, dens, B); - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); enth = enthalpy(x, B); eint = enth - pres / dens; @@ -318,7 +318,7 @@ void actual_eos (I input, T& state) // Density derivatives are computed using the chain rule, e.g. dpdr = dpdx * dxdr. - x = std::pow(dens / B, 1.0_rt / 3.0_rt); + x = std::cbrt(dens / B); dxdr = (1.0_rt / 3.0_rt) * x / dens; Real dpdr = dxdr * dpdx(x); diff --git a/conductivity/stellar/actual_conductivity.H b/conductivity/stellar/actual_conductivity.H index d1da23e9fa..2649d14f38 100644 --- a/conductivity/stellar/actual_conductivity.H +++ b/conductivity/stellar/actual_conductivity.H @@ -299,7 +299,7 @@ actual_conductivity (T& state) // from yakovlev & urpin soviet astro 1980 24 303 and // potekhin et al. 1997 aa 323 415 for degenerate regimes if (dlog10 > drel10) { - amrex::Real xmas = meff * std::pow(state.xne, third); + amrex::Real xmas = meff * std::cbrt(state.xne); amrex::Real ymas = std::sqrt(1.0_rt + xmas*xmas); amrex::Real wfac = weid * state.T/ymas * state.xne; amrex::Real cint = 1.0_rt; @@ -311,10 +311,10 @@ actual_conductivity (T& state) // electron-electron collision frequency and thermal conductivity amrex::Real tpe = xec * std::sqrt(state.xne/ymas); amrex::Real yg = rt3 * tpe/state.T; - amrex::Real xrel = 1.009_rt * std::pow(zbar/abar * state.rho * 1.0e-6_rt, third); + amrex::Real xrel = 1.009_rt * std::cbrt(zbar/abar * state.rho * 1.0e-6_rt); amrex::Real beta2 = xrel*xrel/(1.0_rt + xrel*xrel); - amrex::Real jy = (1.0_rt + 6.0_rt/(5.0_rt*xrel*xrel) + 2.0_rt/(5.0_rt*xrel*xrel*xrel*xrel)) - * ( yg*yg*yg / (3.0_rt * amrex::Math::powi<3>(1.0_rt + 0.07414_rt * yg)) + amrex::Real jy = (1.0_rt + 6.0_rt/(5.0_rt*xrel*xrel) + 2.0_rt/(5.0_rt*amrex::Math::powi<4>(xrel))) + * ( amrex::Math::powi<3>(yg) / (3.0_rt * amrex::Math::powi<3>(1.0_rt + 0.07414_rt * yg)) * std::log((2.81_rt - 0.810_rt*beta2 + yg)/yg) + std::pow(M_PI, 5/6.0_rt) * amrex::Math::powi<4>(yg/(13.91_rt + yg))); amrex::Real vee = 0.511_rt * state.T*state.T * xmas/(ymas*ymas) * std::sqrt(xmas/ymas) * jy; diff --git a/integration/BackwardEuler/be_integrator.H b/integration/BackwardEuler/be_integrator.H index f82cab4711..0c63a2f274 100644 --- a/integration/BackwardEuler/be_integrator.H +++ b/integration/BackwardEuler/be_integrator.H @@ -271,7 +271,7 @@ int be_integrator (BurnT& state, BeT& be) // can we potentially increase the timestep? // backward-Euler has a local truncation error of dt**2 - amrex::Real dt_new = dt_sub * std::pow(1.0_rt / rel_error, 0.5_rt); + amrex::Real dt_new = dt_sub * std::sqrt(1.0_rt / rel_error); dt_sub = std::clamp(dt_new, dt_sub / 2.0, 2.0 * dt_sub); } else { diff --git a/integration/RKC/rkc.H b/integration/RKC/rkc.H index 9598534f33..62973d03e1 100644 --- a/integration/RKC/rkc.H +++ b/integration/RKC/rkc.H @@ -126,8 +126,6 @@ int rkclow (BurnT& state, RkcT& rstate) const amrex::Real onep1 = 1.1_rt; const amrex::Real p4 = 0.4_rt; const amrex::Real p8 = 0.8_rt; - const amrex::Real one3rd = 1.0_rt/3.0_rt; - const amrex::Real two3rd = 2.0_rt/3.0_rt; // Initialize on the first call. @@ -265,7 +263,7 @@ int rkclow (BurnT& state, RkcT& rstate) if (err > 1.0_rt) { // Step is rejected. rstate.nrejct++; - absh = p8 * absh / std::pow(err, one3rd); + absh = p8 * absh / std::cbrt(err); if (absh < hmin) { return IERR_DT_UNDERFLOW; } else { @@ -296,13 +294,13 @@ int rkclow (BurnT& state, RkcT& rstate) } amrex::Real fac = 10.0_rt; if (rstate.naccpt == 1) { - amrex::Real temp2 = std::pow(err, one3rd); + amrex::Real temp2 = std::cbrt(err); if (p8 < fac * temp2) { fac = p8/temp2; } } else { - amrex::Real temp1 = p8 * absh * std::pow(errold, one3rd); - amrex::Real temp2 = std::abs(hold) * std::pow(err, two3rd); + amrex::Real temp1 = p8 * absh * std::cbrt(errold); + amrex::Real temp2 = std::abs(hold) * amrex::Math::powi<2>(std::cbrt(err)); if (temp1 < fac * temp2) { fac = temp1 / temp2; } 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 9da551205f..3d2073ec0e 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.08-206-g798ca3356d16)... -AMReX (23.08-206-g798ca3356d16) initialized +Initializing AMReX (23.05-658-g0165b6743355)... +AMReX (23.05-658-g0165b6743355) initialized starting the single zone burn... Maximum Time (s): 0.01 State Density (g/cm^3): 1000000 @@ -33,8 +33,8 @@ RHS at t = 0 Ni56 3.938787868e-40 ------------------------------------ successful? 1 - - Hnuc = 8.55839024e+18 - - added e = 8.55839024e+16 + - Hnuc = 8.558390237e+18 + - added e = 8.558390237e+16 - final T = 1516425860 ------------------------------------ e initial = 1.284393683e+17 @@ -43,31 +43,31 @@ e final = 2.140232707e+17 new mass fractions: He4 0.8760723967 C12 0.1064099566 -O16 0.0001403204362 +O16 0.0001403204361 Ne20 8.129701234e-05 -Mg24 0.0002972369573 +Mg24 0.0002972369574 Si28 0.01113151728 -S32 0.005297014805 -Ar36 0.0005641483217 -Ca40 6.109847214e-06 -Ti44 2.004306753e-09 -Cr48 3.2721052e-14 -Fe52 7.126143425e-20 -Ni56 1.224820426e-26 +S32 0.005297014784 +Ar36 0.0005641483118 +Ca40 6.109846937e-06 +Ti44 2.004306574e-09 +Cr48 3.272104685e-14 +Fe52 7.126141581e-20 +Ni56 1.22481994e-26 ------------------------------------ species creation rates: omegadot(He4): -12.39276033 omegadot(C12): 10.64099566 -omegadot(O16): 0.01403204362 +omegadot(O16): 0.01403204361 omegadot(Ne20): 0.008129701234 -omegadot(Mg24): 0.02972369573 +omegadot(Mg24): 0.02972369574 omegadot(Si28): 1.113151728 -omegadot(S32): 0.5297014805 -omegadot(Ar36): 0.05641483217 -omegadot(Ca40): 0.0006109847214 -omegadot(Ti44): 2.004306753e-07 -omegadot(Cr48): 3.2721052e-12 -omegadot(Fe52): 7.126143424e-18 -omegadot(Ni56): 1.224720426e-24 +omegadot(S32): 0.5297014784 +omegadot(Ar36): 0.05641483118 +omegadot(Ca40): 0.0006109846937 +omegadot(Ti44): 2.004306574e-07 +omegadot(Cr48): 3.272104685e-12 +omegadot(Fe52): 7.126141581e-18 +omegadot(Ni56): 1.22471994e-24 number of steps taken: 255 -AMReX (23.08-206-g798ca3356d16) finalized +AMReX (23.05-658-g0165b6743355) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/chamulak_VODE_unit_test.out b/unit_test/burn_cell/ci-benchmarks/chamulak_VODE_unit_test.out index 2f73a7431c..e8ab3c56a7 100644 --- a/unit_test/burn_cell/ci-benchmarks/chamulak_VODE_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/chamulak_VODE_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (24.07-16-gdcb9cc0383dc)... -AMReX (24.07-16-gdcb9cc0383dc) initialized +Initializing AMReX (23.05-658-g0165b6743355)... +AMReX (23.05-658-g0165b6743355) initialized starting the single zone burn... Maximum Time (s): 0.01585 State Density (g/cm^3): 1000000000 @@ -13,21 +13,21 @@ RHS at t = 0 ash 0.01230280576 ------------------------------------ successful? 1 - - Hnuc = 5.277400893e+17 - - added e = 8.364680415e+15 - - final T = 1433712612 + - Hnuc = 5.277036926e+17 + - added e = 8.364103527e+15 + - final T = 1433684587 ------------------------------------ e initial = 1.253426044e+18 -e final = 1.261790725e+18 +e final = 1.261790148e+18 ------------------------------------ new mass fractions: -C12 0.9657895158 +C12 0.9657918752 O16 1e-30 -ash 0.03421048417 +ash 0.03420812477 ------------------------------------ species creation rates: -omegadot(C12): -2.158390168 -omegadot(O16): 8.840999775e-44 -omegadot(ash): 2.158390168 +omegadot(C12): -2.15824131 +omegadot(O16): 7.735874803e-44 +omegadot(ash): 2.15824131 number of steps taken: 381 -AMReX (24.07-16-gdcb9cc0383dc) finalized +AMReX (23.05-658-g0165b6743355) finalized 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 4688ffd4fc..0501599a29 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 (23.11-4-ga7afcba3cffd)... -AMReX (23.11-4-ga7afcba3cffd) initialized +Initializing AMReX (23.05-658-g0165b6743355)... +AMReX (23.05-658-g0165b6743355) 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.111466233e+19 - - added e = 4.111466233e+17 - - final T = 6536235789 + - Hnuc = 4.520819292e+19 + - added e = 4.520819292e+17 + - final T = 6714483735 ------------------------------------ e initial = 5.995956082e+17 -e final = 1.010742231e+18 +e final = 1.051677537e+18 ------------------------------------ new mass fractions: -H1 0.007813581548 -He4 9.986703417e-31 -O16 7.812412175e-07 -O20 0.009524501615 -F20 0.00952319448 -Ne20 9.986703417e-31 -Mg24 8.489446236e-12 -Al27 9.986703417e-31 -Si28 0.2115906158 -P31 9.986703417e-31 -S32 0.7615473253 +H1 0.009195917106 +He4 9.999999989e-31 +O16 5.670739329e-07 +O20 0.009949514518 +F20 0.009948148298 +Ne20 9.076266252e-14 +Mg24 1.62767428e-09 +Al27 9.999999989e-31 +Si28 0.2969684506 +P31 9.999999989e-31 +S32 0.6739374007 ------------------------------------ species creation rates: -omegadot(H1): -0.2186418452 +omegadot(H1): -0.08040828943 omegadot(He4): -3 -omegadot(O16): -49.99992188 -omegadot(O20): -0.04754983851 -omegadot(F20): -0.04768055199 +omegadot(O16): -49.99994329 +omegadot(O20): -0.005048548176 +omegadot(F20): -0.005185170248 omegadot(Ne20): -30 -omegadot(Mg24): -9.999999999 +omegadot(Mg24): -9.999999837 omegadot(Al27): -1 -omegadot(Si28): 20.15906158 +omegadot(Si28): 28.69684506 omegadot(P31): -1 -omegadot(S32): 75.15473253 -number of steps taken: 9639 -AMReX (23.11-4-ga7afcba3cffd) finalized +omegadot(S32): 66.39374007 +number of steps taken: 45595 +AMReX (23.05-658-g0165b6743355) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/subch_simple_BE_unit_test.out b/unit_test/burn_cell/ci-benchmarks/subch_simple_BE_unit_test.out index f81d63f2e8..4215e4b70e 100644 --- a/unit_test/burn_cell/ci-benchmarks/subch_simple_BE_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/subch_simple_BE_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (24.10-20-gb9d549bcf4a6)... -AMReX (24.10-20-gb9d549bcf4a6) initialized +Initializing AMReX (23.05-658-g0165b6743355)... +AMReX (23.05-658-g0165b6743355) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... Maximum Time (s): 1e-05 @@ -52,59 +52,59 @@ RHS at t = 0 Ni56 2.574225196e-29 ------------------------------------ successful? 1 - - Hnuc = 7.242722161e+22 - - added e = 7.242722161e+17 - - final T = 3332747822 + - Hnuc = 7.242765333e+22 + - added e = 7.242765333e+17 + - final T = 3332749496 ------------------------------------ e initial = 1.396711859e+18 -e final = 2.120984075e+18 +e final = 2.120988393e+18 ------------------------------------ new mass fractions: H1 0.1076931362 -He4 0.08267440359 -C12 1.785165175e-05 -N13 3.39344959e-08 -N14 3.303771022e-07 -O16 0.1086084686 -F18 1.41796077e-07 -Ne20 0.007002255133 -Ne21 5.733598864e-06 -Na22 0.157136158 -Na23 4.819931671e-07 -Mg24 0.004455038438 -Al27 6.194114228e-06 -Si28 0.061469832 -P31 9.539520569e-06 -S32 0.1152093538 -Ar36 0.09902713352 -Ca40 0.250750628 -Ti44 0.005256054521 -Cr48 0.0006643440612 -Fe52 1.284023015e-05 -Ni56 4.684404848e-08 +He4 0.08267200873 +C12 1.785233903e-05 +N13 3.393565547e-08 +N14 3.304210894e-07 +O16 0.1086069817 +F18 1.418025805e-07 +Ne20 0.007002135899 +Ne21 5.733670948e-06 +Na22 0.1571361579 +Na23 4.819730053e-07 +Mg24 0.004454958807 +Al27 6.193841873e-06 +Si28 0.06146904694 +P31 9.539154561e-06 +S32 0.1152098618 +Ar36 0.0990289131 +Ca40 0.2507533613 +Ti44 0.005255952422 +Cr48 0.0006642930459 +Fe52 1.283819153e-05 +Ni56 4.68317219e-08 ------------------------------------ species creation rates: -omegadot(H1): 769.3136205 -omegadot(He4): -41732.55964 -omegadot(C12): -9998.214835 -omegadot(N13): -9999.996607 -omegadot(N14): -9999.966962 -omegadot(O16): 860.8468635 -omegadot(F18): 0.0141796077 -omegadot(Ne20): 700.2255133 -omegadot(Ne21): 0.5733598864 -omegadot(Na22): 15713.6158 -omegadot(Na23): 0.04819931671 -omegadot(Mg24): 445.5038438 -omegadot(Al27): 0.6194114228 -omegadot(Si28): 6146.9832 -omegadot(P31): 0.9539520569 -omegadot(S32): 11520.93538 -omegadot(Ar36): 9902.713352 -omegadot(Ca40): 25075.0628 -omegadot(Ti44): 525.6054521 -omegadot(Cr48): 66.43440612 -omegadot(Fe52): 1.284023015 -omegadot(Ni56): 0.004684404848 -number of steps taken: 9552 -AMReX (24.10-20-gb9d549bcf4a6) finalized +omegadot(H1): 769.3136185 +omegadot(He4): -41732.79913 +omegadot(C12): -9998.214766 +omegadot(N13): -9999.996606 +omegadot(N14): -9999.966958 +omegadot(O16): 860.6981667 +omegadot(F18): 0.01418025805 +omegadot(Ne20): 700.2135899 +omegadot(Ne21): 0.5733670948 +omegadot(Na22): 15713.61579 +omegadot(Na23): 0.04819730053 +omegadot(Mg24): 445.4958807 +omegadot(Al27): 0.6193841873 +omegadot(Si28): 6146.904694 +omegadot(P31): 0.9539154561 +omegadot(S32): 11520.98618 +omegadot(Ar36): 9902.89131 +omegadot(Ca40): 25075.33613 +omegadot(Ti44): 525.5952422 +omegadot(Cr48): 66.42930459 +omegadot(Fe52): 1.283819153 +omegadot(Ni56): 0.00468317219 +number of steps taken: 9520 +AMReX (23.05-658-g0165b6743355) finalized diff --git a/unit_test/burn_cell/ci-benchmarks/subch_simple_unit_test.out b/unit_test/burn_cell/ci-benchmarks/subch_simple_unit_test.out index 7417850737..91943fbdb9 100644 --- a/unit_test/burn_cell/ci-benchmarks/subch_simple_unit_test.out +++ b/unit_test/burn_cell/ci-benchmarks/subch_simple_unit_test.out @@ -1,5 +1,5 @@ -Initializing AMReX (24.10-20-gb9d549bcf4a6)... -AMReX (24.10-20-gb9d549bcf4a6) initialized +Initializing AMReX (23.05-658-g0165b6743355)... +AMReX (23.05-658-g0165b6743355) initialized starting the single zone burn... reading in network electron-capture / beta-decay tables... Maximum Time (s): 1e-05 @@ -52,59 +52,59 @@ RHS at t = 0 Ni56 2.574225196e-29 ------------------------------------ successful? 1 - - Hnuc = 7.241796955e+22 - - added e = 7.241796955e+17 - - final T = 3332712347 + - Hnuc = 7.2421617e+22 + - added e = 7.2421617e+17 + - final T = 3332726322 ------------------------------------ e initial = 1.396711859e+18 -e final = 2.120891555e+18 +e final = 2.120928029e+18 ------------------------------------ new mass fractions: -H1 0.1076931293 -He4 0.08271515334 -C12 1.784227067e-05 -N13 3.391967633e-08 -N14 3.298055329e-07 -O16 0.1086399603 -F18 1.416783115e-07 -Ne20 0.007004688063 -Ne21 5.732083096e-06 -Na22 0.1571361384 -Na23 4.823584771e-07 -Mg24 0.004456669669 -Al27 6.199112576e-06 -Si28 0.06147978562 -P31 9.545147785e-06 -S32 0.1151792002 -Ar36 0.09897516874 -Ca40 0.2507395336 -Ti44 0.005261227135 -Cr48 0.0006660893333 -Fe52 1.290260885e-05 -Ni56 4.720489636e-08 +H1 0.1076931297 +He4 0.08269937758 +C12 1.784557982e-05 +N13 3.392480906e-08 +N14 3.29915611e-07 +O16 0.1086304808 +F18 1.417242761e-07 +Ne20 0.007003921313 +Ne21 5.732677863e-06 +Na22 0.1571361405 +Na23 4.822287628e-07 +Mg24 0.004456132113 +Al27 6.197297936e-06 +Si28 0.06147467977 +P31 9.542765771e-06 +S32 0.1151815145 +Ar36 0.09898521539 +Ca40 0.2507594706 +Ti44 0.005260845954 +Cr48 0.0006658458928 +Fe52 1.289260143e-05 +Ni56 4.714390576e-08 ------------------------------------ species creation rates: -omegadot(H1): 769.3129335 -omegadot(He4): -41728.48467 -omegadot(C12): -9998.215773 +omegadot(H1): 769.3129719 +omegadot(He4): -41730.06224 +omegadot(C12): -9998.215442 omegadot(N13): -9999.996608 -omegadot(N14): -9999.967019 -omegadot(O16): 863.9960253 -omegadot(F18): 0.01416783115 -omegadot(Ne20): 700.4688063 -omegadot(Ne21): 0.5732083096 -omegadot(Na22): 15713.61384 -omegadot(Na23): 0.04823584771 -omegadot(Mg24): 445.6669669 -omegadot(Al27): 0.6199112576 -omegadot(Si28): 6147.978562 -omegadot(P31): 0.9545147785 -omegadot(S32): 11517.92002 -omegadot(Ar36): 9897.516874 -omegadot(Ca40): 25073.95336 -omegadot(Ti44): 526.1227135 -omegadot(Cr48): 66.60893333 -omegadot(Fe52): 1.290260885 -omegadot(Ni56): 0.004720489636 -number of steps taken: 564 -AMReX (24.10-20-gb9d549bcf4a6) finalized +omegadot(N14): -9999.967008 +omegadot(O16): 863.0480784 +omegadot(F18): 0.01417242761 +omegadot(Ne20): 700.3921313 +omegadot(Ne21): 0.5732677863 +omegadot(Na22): 15713.61405 +omegadot(Na23): 0.04822287628 +omegadot(Mg24): 445.6132113 +omegadot(Al27): 0.6197297936 +omegadot(Si28): 6147.467977 +omegadot(P31): 0.9542765771 +omegadot(S32): 11518.15145 +omegadot(Ar36): 9898.521539 +omegadot(Ca40): 25075.94706 +omegadot(Ti44): 526.0845954 +omegadot(Cr48): 66.58458928 +omegadot(Fe52): 1.289260143 +omegadot(Ni56): 0.004714390576 +number of steps taken: 568 +AMReX (23.05-658-g0165b6743355) finalized 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 69d51e0bea..2b4ca43587 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 @@ -6,8 +6,8 @@ 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.000510727 3.96209e-06 2.88476e-08 1.15228e-10 2.13986e-14 1.28346e-18 4.13423e-23 7.01099e-28 1.05157e-30 2.45097e-27 2.46121e-27 -3.16228e-07 1e+06 3.07032e+09 0.704166 6.23619e-07 0.191035 1.67132e-11 0.0110708 0.024113 0.0345257 0.0260972 0.00875392 0.000232005 5.60651e-06 7.35601e-08 4.44762e-11 8.74683e-15 9.19961e-19 2.00534e-23 3.69734e-27 9.94337e-27 7.52646e-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.58907e-19 1.73607e-22 4.89224e-26 3.55151e-20 + 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-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_NSE_state_over_time.txt b/unit_test/burn_cell_sdc/ci-benchmarks/aprox19_NSE_state_over_time.txt index f769696d24..acc229b06e 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.37008e+09 1e-30 2.88682e-10 0.131281 1.05261e-06 8.12352e-11 2.72153e-06 4.16582e-08 9.4534e-06 0.00582325 0.00649833 0.00424371 0.00488588 0.00129866 0.0372919 0.0580338 0.670922 0.0532456 1.44201e-06 0.0264607 -0.00115544 1e+07 5.37008e+09 1e-30 2.88682e-10 0.131281 1.05261e-06 8.12351e-11 2.72152e-06 4.16582e-08 9.4534e-06 0.00582325 0.00649833 0.00424371 0.00488588 0.00129865 0.0372919 0.0580338 0.670922 0.0532457 1.44201e-06 0.0264607 -0.00146854 1e+07 5.37008e+09 1e-30 2.88682e-10 0.131281 1.05261e-06 8.1235e-11 2.72152e-06 4.16581e-08 9.45339e-06 0.00582325 0.00649833 0.00424371 0.00488588 0.00129865 0.0372919 0.0580338 0.670922 0.0532457 1.44201e-06 0.0264607 -0.00186649 1e+07 5.37008e+09 1e-30 2.88682e-10 0.131281 1.0526e-06 8.12349e-11 2.72152e-06 4.16581e-08 9.45339e-06 0.00582325 0.00649834 0.00424371 0.00488588 0.00129865 0.0372919 0.0580338 0.670922 0.0532458 1.442e-06 0.0264607 -0.00237227 1e+07 5.37008e+09 1e-30 2.88681e-10 0.131281 1.0526e-06 8.12347e-11 2.72152e-06 4.1658e-08 9.45338e-06 0.00582325 0.00649834 0.00424372 0.00488589 0.00129865 0.0372919 0.0580338 0.670922 0.0532459 1.442e-06 0.0264607 -0.00301511 1e+07 5.37008e+09 1e-30 2.8868e-10 0.131281 1.0526e-06 8.12345e-11 2.72151e-06 4.1658e-08 9.45337e-06 0.00582325 0.00649834 0.00424372 0.00488589 0.00129865 0.0372919 0.0580338 0.670922 0.053246 1.44199e-06 0.0264607 -0.00383215 1e+07 5.37008e+09 1e-30 2.8868e-10 0.13128 1.0526e-06 8.12342e-11 2.72151e-06 4.16579e-08 9.45336e-06 0.00582325 0.00649834 0.00424372 0.0048859 0.00129865 0.0372919 0.0580338 0.670923 0.0532462 1.44199e-06 0.0264607 - 0.0048706 1e+07 5.37008e+09 1e-30 2.88679e-10 0.13128 1.05259e-06 8.12338e-11 2.7215e-06 4.16577e-08 9.45334e-06 0.00582325 0.00649834 0.00424372 0.0048859 0.00129865 0.0372919 0.0580338 0.670923 0.0532464 1.44198e-06 0.0264606 -0.00619044 1e+07 5.37008e+09 1e-30 2.88677e-10 0.13128 1.05259e-06 8.12333e-11 2.72149e-06 4.16576e-08 9.45332e-06 0.00582325 0.00649835 0.00424373 0.00488591 0.00129865 0.0372918 0.0580338 0.670923 0.0532467 1.44197e-06 0.0264606 -0.00786793 1e+07 5.37008e+09 1e-30 2.88676e-10 0.131279 1.05259e-06 8.12327e-11 2.72148e-06 4.16574e-08 9.45329e-06 0.00582325 0.00649835 0.00424373 0.00488593 0.00129865 0.0372918 0.0580338 0.670923 0.053247 1.44195e-06 0.0264606 - 0.01 1e+07 5.37008e+09 1e-30 2.88674e-10 0.131279 1.05258e-06 8.12319e-11 2.72146e-06 4.16571e-08 9.45326e-06 0.00582325 0.00649836 0.00424374 0.00488594 0.00129864 0.0372918 0.0580337 0.670923 0.0532475 1.44194e-06 0.0264606 +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