create an NSE table with finer spacing in rho and Ye (#1427)

networks/aprox19/nse_table_size.H

@@ -12,20 +12,20 @@ namespace nse_table_size {
     const std::string table_name{"nse_aprox19.tbl"};
 
     constexpr int ntemp = 101;
-    constexpr int nden = 31;
-    constexpr int nye = 15;
+    constexpr int nden = 61;
+    constexpr int nye = 29;
 
     constexpr Real logT_min = 9.4;
     constexpr Real logT_max = 10.4;
     constexpr Real dlogT = 0.01;
 
     constexpr Real logrho_min = 7.0;
     constexpr Real logrho_max = 10.0;
-    constexpr Real dlogrho = 0.1;
+    constexpr Real dlogrho = 0.05;
 
     constexpr Real ye_min = 0.43;
     constexpr Real ye_max = 0.5;
-    constexpr Real dye = 0.005;
+    constexpr Real dye = 0.0025;
 
 }
 #endif

nse_tabular/make_nse_table.py

@@ -128,8 +128,8 @@ def generate_table():
     nse_net = make_nse_network()
 
     Ts = np.logspace(9.4, 10.4, 101)
-    rhos = np.logspace(7, 10, 31)
-    yes = np.linspace(0.43, 0.5, 15)
+    rhos = np.logspace(7, 10, 61)
+    yes = np.linspace(0.43, 0.5, 29)
 
     #Ts = np.logspace(9.4, 10.4, 26)
     #rhos = np.logspace(7, 10, 16)

unit_test/burn_cell/ci-benchmarks/aprox19_nse_unit_test.out

@@ -1,5 +1,5 @@
-Initializing AMReX (23.11-5-gd36463103dae)...
-AMReX (23.11-5-gd36463103dae) initialized
+Initializing AMReX (23.12-21-gef38229189e3)...
+AMReX (23.12-21-gef38229189e3) initialized
 starting the single zone burn...
 reading the NSE table (C++) ...
 Maximum Time (s): 0.01
@@ -48,53 +48,53 @@ burn entered NSE during integration (after 280 steps), zone = (0, 0, 0)
 recovering burn failure in NSE, zone = (0, 0, 0)
 ------------------------------------
 successful? 1
- - Hnuc = 2.576395456e+21
- - added e = 2.576395456e+19
- - final T = 1.063285386e+10
+ - Hnuc = 2.578034812e+21
+ - added e = 2.578034812e+19
+ - final T = 1.063446163e+10
 ------------------------------------
 e initial = 1.516256085e+18
-e final =   2.728021064e+19
+e final =   2.72966042e+19
 ------------------------------------
 new mass fractions: 
-H1 3.698300896e-05
-He3 7.832720283e-11
-He4 0.06033012727
-C12 9.999991743e-31
-N14 9.875775959e-14
-O16 2.423296446e-15
-Ne20 9.999991743e-31
-Mg24 3.915447813e-20
-Si28 2.173376783e-22
-S32 9.999991743e-31
-Ar36 9.999991743e-31
-Ca40 9.999991743e-31
-Ti44 4.768494533e-30
-Cr48 9.999991743e-31
-Fe52 2.410423559e-30
-Fe54 9.999991743e-31
-Ni56 9.999991743e-31
-n 0.4681726624
-p 0.4714602272
+H1 5.812961393e-07
+He3 9.821888481e-12
+He4 0.06008117992
+C12 9.999992174e-31
+N14 1.534220194e-15
+O16 2.310011368e-15
+Ne20 1.352206382e-19
+Mg24 5.363834622e-22
+Si28 4.61640554e-23
+S32 9.999992174e-31
+Ar36 9.999992174e-31
+Ca40 4.564517527e-29
+Ti44 9.999992174e-31
+Cr48 5.163875678e-30
+Fe52 9.999992174e-31
+Fe54 3.381529446e-28
+Ni56 9.999992174e-31
+n 0.4683168262
+p 0.4716014125
 ------------------------------------
 species creation rates: 
-omegadot(H1): -9.996301699
-omegadot(He3): -2.499999992
-omegadot(He4): -73.96698727
+omegadot(H1): -9.99994187
+omegadot(He3): -2.499999999
+omegadot(He4): -73.99188201
 omegadot(C12): -2.5
 omegadot(N14): -2.5
 omegadot(O16): -2.5
-omegadot(Ne20): -8.256539967e-35
-omegadot(Mg24): 3.915447813e-18
-omegadot(Si28): 2.173376773e-20
-omegadot(S32): -8.256539967e-35
-omegadot(Ar36): -8.256539967e-35
-omegadot(Ca40): -8.256539967e-35
-omegadot(Ti44): 3.768494533e-28
-omegadot(Cr48): -8.256539967e-35
-omegadot(Fe52): 1.410423559e-28
-omegadot(Fe54): -8.256539967e-35
-omegadot(Ni56): -8.256539967e-35
-omegadot(n): 46.81726624
-omegadot(p): 47.14602272
-number of steps taken: 38629
-AMReX (23.11-5-gd36463103dae) finalized
+omegadot(Ne20): 1.352206382e-17
+omegadot(Mg24): 5.363834612e-20
+omegadot(Si28): 4.61640544e-21
+omegadot(S32): -7.826200767e-35
+omegadot(Ar36): -7.826200767e-35
+omegadot(Ca40): 4.464517527e-27
+omegadot(Ti44): -7.826200767e-35
+omegadot(Cr48): 4.163875678e-28
+omegadot(Fe52): -7.826200767e-35
+omegadot(Fe54): 3.371529446e-26
+omegadot(Ni56): -7.826200767e-35
+omegadot(n): 46.83168262
+omegadot(p): 47.16014125
+number of steps taken: 27367
+AMReX (23.12-21-gef38229189e3) finalized

unit_test/test_nse/ci-benchmarks/aprox21_ci.out

@@ -1,6 +1,7 @@
-AMReX (22.10-25-g735c3513153f) initialized
+Initializing AMReX (23.12-21-gef38229189e3)...
+AMReX (23.12-21-gef38229189e3) initialized
 starting the single zone burn...
-0x7ffd73ad6d20
+0x7ffcb5676e40
 State Density (g/cm^3): 277355338.4
 State Temperature (K): 5197769252
 Mass Fraction (H1): 0.4
@@ -50,4 +51,4 @@ Fe56 : 1.686482626e-05
 Ni56 : 0.8982137413
 n : 9.442515318e-10
 p : 0.001666851761
-AMReX (22.10-25-g735c3513153f) finalized
+AMReX (23.12-21-gef38229189e3) finalized

unit_test/test_nse_interp/ci-benchmarks/aprox19.out

@@ -1,78 +1,78 @@
-Initializing AMReX (23.11-5-gd36463103dae)...
-AMReX (23.11-5-gd36463103dae) initialized
+Initializing AMReX (23.12-21-gef38229189e3)...
+AMReX (23.12-21-gef38229189e3) initialized
 starting the single zone burn...
 reading the NSE table (C++) ...
 rho, T, Ye = 1230000000 5180000000 0.472
-density value brackets: 9 < 9.089905111 < 9.1
+density value brackets: 9.05 < 9.089905111 < 9.1
 temperature value brackets: 9.71 < 9.71432976 < 9.72
-Ye value brackets: 0.475 < 0.472 < 0.47
+Ye value brackets: 0.4725 < 0.472 < 0.47
 
 4 rho values: 
-8.9 55.16599379
-9 55.28915823
-9.1 55.39796653
-9.2 55.49591413
+9 55.5223665
+9.05 55.57080116
+9.1 55.61707675
+9.15 55.66151541
 
-cubic interpolated value: 55.38753351
+cubic interpolated value: 55.6078922
 
 4 T values: 
-9.7 55.16599379
-9.71 54.97411678
-9.72 54.71072288
-9.73 54.35084601
+9.7 55.5223665
+9.71 55.39489494
+9.72 55.22387568
+9.73 54.99343581
 
-cubic interpolated value: 54.87031639
+cubic interpolated value: 55.32712406
 
 4 Ye values: 
-0.48 55.16599379
-0.475 55.41756427
-0.47 55.68822983
-0.465 55.93040523
+0.475 55.5223665
+0.4725 55.65227675
+0.47 55.78142906
+0.4675 55.90244371
 
-cubic interpolated value: 55.58071765
+cubic interpolated value: 55.678404
 
 tricubic interpolated values: 
-abar    = 55.60621897
-bea     = 8.755258739
-dyedt   = -0.004462626048
-dabardt = 2.472895091e-16
-dbeadt = 0.0001077443794
-e_nu    = 1.081686275e+16
+abar    = 55.60652462
+bea     = 8.755119523
+dyedt   = -0.004527880871
+dabardt = 1.310223024e-16
+dbeadt = 0.000111945508
+e_nu    = 1.101641642e+16
 X(H1) = 0
-X(He3) = 2.646187332e-14
-X(He4) = 0.0003164651992
-X(C12) = 7.949931897e-10
-X(N14) = 3.087926554e-14
-X(O16) = 1.98556354e-09
-X(Ne20) = 2.420691526e-11
-X(Mg24) = 1.01690623e-08
-X(Si28) = 7.957002233e-06
-X(S32) = 1.148035178e-05
-X(Ar36) = 1.188979307e-05
-X(Ca40) = 1.11725697e-05
-X(Ti44) = 0.0001435646298
-X(Cr48) = 0.01141230313
-X(Fe52) = 0.07731922235
-X(Fe54) = 0.9104458693
-X(Ni56) = 0.0003104223608
-X(n) = 2.004738472e-08
-X(p) = 9.620335872e-06
+X(He3) = 2.723457525e-14
+X(He4) = 0.0003164137294
+X(C12) = 7.975878903e-10
+X(N14) = 3.102372259e-14
+X(O16) = 1.998161464e-09
+X(Ne20) = 2.444604777e-11
+X(Mg24) = 1.029036116e-08
+X(Si28) = 8.173012779e-06
+X(S32) = 1.194076593e-05
+X(Ar36) = 1.243414613e-05
+X(Ca40) = 1.259531492e-05
+X(Ti44) = 0.0001438001413
+X(Cr48) = 0.01146726617
+X(Fe52) = 0.07746625779
+X(Fe54) = 0.9098424639
+X(Ni56) = 0.0007087774945
+X(n) = 2.078852742e-08
+X(p) = 9.843630258e-06
 
 testing temperature derivatives of cubic
 first finite-difference derivatives
-dAbar/dT = -1.070778042e-09
-dbea/dT = -6.886297725e-12
+dAbar/dT = -1.072562355e-09
+dbea/dT = -6.867059575e-12
 now using derivative of the interpolant
-dAbar/dT = -1.070778122e-09
-dbea/dT = -6.886272826e-12
+dAbar/dT = -1.072562604e-09
+dbea/dT = -6.867048589e-12
 
-EOS e consistency check (old method): 1.395273834e+18 1.388449367e+18
-updated T: 6394612134
-change in abar: 55.60621897 50.27196857
-EOS e consistency check (new method): 1.388449367e+18 1.388449367e+18
+EOS e consistency check (old method): 1.395278886e+18 1.38844906e+18
+updated T: 6394534499
+change in abar: 55.60652462 50.26831386
+EOS e consistency check (new method): 1.38844906e+18 1.38844906e+18
 
-EOS p consistency check (old method): 6.622132093e+26 6.57785215e+26
-updated T: 6466848772
-change in abar: 55.60621897 49.50670173
-EOS p consistency check (new method): 6.57785215e+26 6.57785215e+26
-AMReX (23.11-5-gd36463103dae) finalized
+EOS p consistency check (old method): 6.622159603e+26 6.577850616e+26
+updated T: 6466757500
+change in abar: 55.60652462 49.50320619
+EOS p consistency check (new method): 6.577850616e+26 6.577850616e+26
+AMReX (23.12-21-gef38229189e3) finalized