add table lookup to utility unit test (#601)

Author: baperry2

.github/workflows/ci.yml

@@ -190,7 +190,7 @@ jobs:
           ccache -z
           make -j ${{env.NPROCS}} TINY_PROFILE=TRUE USE_CCACHE=TRUE USE_AMREX_LIB=TRUE ${{matrix.amrex_build_args}}; \
           if [ "${{matrix.comp}}" == 'gnu' ] || [ "${{matrix.comp}}" == 'llvm' ]; then \
-            ./Pele2d.${{matrix.comp}}.TPROF.ex;
+            ./Pele2d.${{matrix.comp}}.TPROF.ex input;
           fi;
           make realclean;
           if [ $? -ne 0 ]; then exit 1; fi; \

Source/Utility/BlackBoxFunction/Table.H

@@ -4,6 +4,7 @@
 #define MAXD_TABLE 5
 #include <iostream>
 #include <fstream>
+#include <AMReX_ParmParse.H>
 #include "BlackBoxFunction.H"
 #include "Utilities.H"
 

Testing/Exec/UtilityUnitTest/input

@@ -0,0 +1 @@
+test.table.filename = ../../../Support/CMLM/peletable_2dim_10grid

Testing/Exec/UtilityUnitTest/main.cpp

@@ -1,6 +1,8 @@
 #include <AMReX.H>
 #include <AMReX_Print.H>
 #include "Utilities.H"
+#include "Table.H"
+#include "BlackBoxFunction.H"
 
 namespace m2c = pele::physics::utilities::mks2cgs;
 namespace c2m = pele::physics::utilities::cgs2mks;
@@ -10,6 +12,23 @@ main(int argc, char* argv[])
 {
   amrex::Initialize(argc, argv);
 
+  // Test table interpolation
+  pele::physics::TabulatedFunctionParams tabfunc_par;
+  tabfunc_par.host_only_parm().parm_parse_prefix = "test";
+  tabfunc_par.initialize();
+  pele::physics::TabulatedFunction tabfunc(
+    static_cast<const pele::physics::TabulatedFunctionData*>(
+      &tabfunc_par.host_parm()));
+  amrex::Real outval;
+  const amrex::GpuArray<amrex::Real, 2> invals{1.0, 1.0};
+  const std::string varname = "T";
+  const int idxT =
+    pele::physics::get_var_index(varname.c_str(), &tabfunc_par.host_parm());
+  tabfunc.get_value(idxT, invals.data(), outval);
+  AMREX_ALWAYS_ASSERT(amrex::almostEqual(outval, 3.60));
+  tabfunc_par.deallocate();
+  amrex::Print() << "TabulatedFunction Tests passed" << std::endl;
+
   // Test Unit Conversions
   AMREX_ALWAYS_ASSERT(m2c::Length(1.0) == 100.0);
   AMREX_ALWAYS_ASSERT(c2m::Length(1.0, 2) == 1.0e-4);