adds unit test

CMakeLists.txt

@@ -416,6 +416,10 @@ if (BUILD_TESTING)
         NAME fluid_solver_step_check
         COMMAND ${PYTHON_EXECUTABLE} ${PROJECT_SOURCE_DIR}/examples/solver_comparison/sanity_check.py
         WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
+    add_test(
+        NAME fluid_solver_acceleration_check
+        COMMAND ${PYTHON_EXECUTABLE} ${PROJECT_SOURCE_DIR}/examples/Lagrangian_acceleration/acc_check.py
+        WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
     ### basic particle functionality
     add_test(
         NAME test_particle_set_init

examples/Lagrangian_acceleration/acc_check.py

@@ -84,14 +84,15 @@ class acc_test(TurTLE.PP):
         return pars
 
 def main():
-    if not os.path.exists('test.h5'):
+    if not os.path.exists('acceleration_check_data.h5'):
         cc = TurTLE.DNS()
-        cc.launch(['NSVE', '--np', '2'])
+        cc.launch(['NSVE', '--np', '2', '--simname', 'acceleration_check_data'])
     bla = acc_test()
     bla.launch(
             ['acceleration_check',
              '--np', '1',
-             '--ntpp', '1']
+             '--ntpp', '1',
+             '--simname', 'acceleration_check_data']
             + sys.argv[1:])
     return None
 

examples/solver_comparison/fluid_solver_check_analytic.cpp

@@ -164,9 +164,11 @@ int check_solvers(
             energy_NSE,
             energy_NSVE,
             analytic_energy);
-    const bool difference_small =
-        (energy_NSE - energy_NSVE) / std::sqrt(energy_NSE*energy_NSVE) < 1e-5;
-    if (difference_small)
+    const bool difference_small_NSE =
+        (energy_NSE - analytic_energy) / std::sqrt(energy_NSE*analytic_energy) < 1e-5;
+    const bool difference_small_NSVE =
+        (energy_NSVE - analytic_energy) / std::sqrt(energy_NSVE*analytic_energy) < 1e-5;
+    if (difference_small_NSE && difference_small_NSVE)
         return EXIT_SUCCESS;
     else
         return EXIT_FAILURE;

examples/solver_comparison/sanity_check.py

@@ -97,14 +97,14 @@ def main():
             ['fluid_solver_check',
              '--np', '1',
              '--ntpp', '1',
-             '--simname', 'fluid_solver_onestep_sanity_check']
+             '--simname', 'fluid_solver_onestep_analytic_check']
             + sys.argv[1:])
     bla = sanity_check(check_type = 'general')
     bla.launch(
             ['fluid_solver_check',
              '--np', '1',
              '--ntpp', '1',
-             '--simname', 'fluid_solver_onestep_sanity_check']
+             '--simname', 'fluid_solver_onestep_general_check']
             + sys.argv[1:])
     return None