Add CI coverage for ignition_chamulak (#1284)

.github/workflows/burn_cell.yml

@@ -53,11 +53,27 @@ jobs:
           cd unit_test/burn_cell
           ./main3d.gnu.ex inputs_ecsn > test.out
 
-      - name: Compare to stored output (ECSN)
+      - name: Compare to stored output (VODE, ECSN)
         run: |
           cd unit_test/burn_cell
           diff -I "^AMReX" -I "^reading in reaclib rates" test.out ci-benchmarks/ecsn_unit_test.out
 
+      - name: Compile, burn_cell (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/burn_cell
+          make realclean
+          make NETWORK_DIR=ignition_chamulak -j 4
+
+      - name: Run burn_cell (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/burn_cell
+          ./main3d.gnu.ex inputs_ignition_chamulak > test.out
+
+      - name: Compare to stored output (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/burn_cell
+          diff -I "^AMReX" -I "^reading in reaclib rates" test.out ci-benchmarks/chamulak_VODE_unit_test.out
+
       - name: Compile, burn_cell (ForwardEuler, triple_alpha_plus_cago)
         run: |
           cd unit_test/burn_cell
@@ -121,5 +137,3 @@ jobs:
         run: |
           cd unit_test/burn_cell
           diff -I "^AMReX" -I "^reading in reaclib rates" test.out ci-benchmarks/aprox13_RKC_unit_test.out
-
-

.github/workflows/test_rhs.yml

@@ -48,6 +48,23 @@ jobs:
           cd unit_test/test_rhs
           diff test.out ci-benchmarks/ignition_simple.out
 
+      - name: Compile, test_rhs (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/test_rhs
+          make realclean
+          make NETWORK_DIR=ignition_chamulak -j 4
+
+      - name: Run test_rhs (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/test_rhs
+          ./main3d.gnu.ex inputs_chamulak
+          ../../external/amrex/Tools/Plotfile/fextrema.gnu.ex react_chamulak_test_rhs.VODE > test.out
+
+      - name: Compare to stored output (VODE, ignition_chamulak)
+        run: |
+          cd unit_test/test_rhs
+          diff test.out ci-benchmarks/chamulak.out
+
       - name: Compile, test_rhs (VODE, iso7)
         run: |
           cd unit_test/test_rhs

unit_test/burn_cell/ci-benchmarks/chamulak_VODE_unit_test.out

@@ -0,0 +1,32 @@
+AMReX (23.07-24-g01b750d909b0) initialized
+starting the single zone burn...
+Maximum Time (s): 0.01585
+State Density (g/cm^3): 1000000000
+State Temperature (K): 1000000000
+Mass Fraction (C12): 1
+Mass Fraction (O16): 0
+Mass Fraction (ash): 0
+RHS at t = 0
+   C12 -0.01845420864
+   O16 0
+   ash 0.01230280576
+------------------------------------
+successful? 1
+ - Hnuc = 5.162791463e+17
+ - added e = 8.183024469e+15
+ - final T = 1424869952
+------------------------------------
+e initial = 1.253426044e+18
+e final =   1.261609069e+18
+------------------------------------
+new mass fractions: 
+C12 0.966535706
+O16 1.000003322e-30
+ash 0.033464294
+------------------------------------
+species creation rates: 
+omegadot(C12): -2.111311924
+omegadot(O16): 2.096110914e-34
+omegadot(ash): 2.111311924
+number of steps taken: 384
+AMReX (23.07-24-g01b750d909b0) finalized

unit_test/burn_cell/inputs_ignition_chamulak

@@ -0,0 +1,27 @@
+unit_test.run_prefix = "react_ignition_"
+
+unit_test.small_temp = 1e5
+unit_test.small_dens = 1e5
+
+integrator.burner_verbose = 0
+
+# Set which jacobian to use
+# 1 = analytic jacobian
+# 2 = numerical jacobian
+integrator.jacobian   = 1
+
+integrator.renormalize_abundances = 0
+
+integrator.rtol_spec = 1.0e-6
+integrator.rtol_enuc = 1.0e-6
+integrator.atol_spec = 1.0e-6
+integrator.atol_enuc = 1.0e-6
+
+unit_test.tmax    = 1.585e-2
+
+unit_test.density = 1.e9
+unit_test.temperature = 1.e9
+
+unit_test.X1  = 1.0
+unit_test.X2  = 0.0
+unit_test.X3  = 0.0

unit_test/test_rhs/ci-benchmarks/chamulak.out

@@ -0,0 +1,29 @@
+ plotfile = react_chamulak_test_rhs.VODE
+ time = 0
+              variables          minimum value          maximum value
+ density                                 10000              100000000
+ temperature                          50000000             5000000000
+ Ydot_carbon-12                  -2263658781.5      -3.2218541118e-71
+ Ydot_oxygen-16                              0                      0
+ Ydot_ash                     2.1479027412e-71           1509105854.4
+ Xold_carbon-12                            0.2                    0.9
+ Xold_oxygen-16                           0.05                    0.4
+ Xold_ash                                 0.05                    0.4
+ Edot                         9.3682885505e-53        6.584051211e+27
+ J_carbon-12_carbon-12           -5030352847.9      -3.2218541118e-70
+ J_oxygen-16_carbon-12                       0                      0
+ J_ash_carbon-12              2.1479027412e-70           3353568565.3
+ J_E_carbon-12                9.3682885505e-52       1.4631224913e+28
+ J_carbon-12_oxygen-16                       0                      0
+ J_oxygen-16_oxygen-16                       0                      0
+ J_ash_oxygen-16                             0                      0
+ J_E_oxygen-16                              -0                     -0
+ J_carbon-12_ash                             0                      0
+ J_oxygen-16_ash                             0                      0
+ J_ash_ash                                   0                      0
+ J_E_ash                                    -0                     -0
+ J_carbon-12_E               -4.4154660816e-08      -7.3754668622e-85
+ J_oxygen-16_E                               0                      0
+ J_ash_E                                     0                      0
+ J_E_E                        2.1445881583e-66       1.2842772523e+11
+