[WIP] Use monotonic cubic Hermite interpolation for the NSE table

.github/workflows/test_nse_interp.yml

@@ -26,22 +26,22 @@ jobs:
           sudo apt-get update -y -qq
           sudo apt-get -qq -y install curl cmake jq clang g++>=9.3.0
 
-      - name: Compile, test_nse_interp (NSE, aprox19)
+      - name: Compile, nse_interp_check (NSE, aprox19)
         run: |
-          cd unit_test/test_nse_interp
+          cd unit_test/nse_interp_check
           make realclean
           make -j 4
 
-      - name: Run test_nse_interp (NSE, aprox19)
+      - name: Run nse_interp_check (NSE, aprox19)
         run: |
-          cd unit_test/test_nse_interp
+          cd unit_test/nse_interp_check
           ./main3d.gnu.ex amrex.fpe_trap_{invalid,zero,overflow}=1 > test.out
 
       - name: Print backtrace
-        if: ${{ failure() && hashFiles('unit_test/test_nse_interp/Backtrace.0') != '' }}
-        run: cat unit_test/test_nse_interp/Backtrace.0
+        if: ${{ failure() && hashFiles('unit_test/nse_interp_check/Backtrace.0') != '' }}
+        run: cat unit_test/nse_interp_check/Backtrace.0
 
       - name: Compare to stored output (NSE, aprox19)
         run: |
-          cd unit_test/test_nse_interp
+          cd unit_test/nse_interp_check
           diff -I "^Initializing AMReX" -I "^AMReX" -I "^reading in reaclib rates" test.out ci-benchmarks/aprox19.out

Docs/source/unit_tests.rst

@@ -174,17 +174,17 @@ One-zone tests
 Infrastructure tests
 ====================
 
-.. index:: test_linear_algebra, test_nse_interp, test_parameters, test_sdc_vode_rhs
+.. index:: test_linear_algebra, nse_interp_check, test_parameters, test_sdc_vode_rhs
 
 * ``test_linear_algebra`` :
 
   create a diagonally dominant matrix, multiply it by a test vector, $x$,
   to get $b = Ax$, and then call the linear algebra routines to see if we
   we recover $x$ from $b$.
 
-* ``test_nse_interp`` :
+* ``nse_interp_check`` :
 
-  run various tests of the NSE interpolation routines.
+  run various tests of the tabular NSE interpolation routines.
 
 * ``test_parameters`` :
 

nse_tabular/_parameters

@@ -29,3 +29,6 @@ nse_relax_factor  real   1.0
 
 # do we do tri-linear or tri-cubic interpolation on the table?
 nse_table_interp_linear   bool    0
+
+# do we always do tri-linear for species, even when doing cubic for other quantities?
+nse_force_linear_interp_species   bool   0