Added single precision test. Examples now work in both double and single precision.

ikorotkin · ikorotkin · commit 3bb9d48b635a · 2019-05-17T18:20:18.000+01:00
diff --git a/examples/diffusion_2d/diffusion_2d.cpp b/examples/diffusion_2d/diffusion_2d.cpp
@@ -80,9 +80,7 @@ int main()
     // parameters defined in solver_options.h
     dae::SolverOptions opt;
 
-    opt.dt_init         = 5.0e-5;  // Change initial time step
-    opt.fact_every_iter = false;   // Gain some speed. The matrices will be
-                                   // factorized only once each time step.
+    opt.dt_init = 5.0e-5;  // Change initial time step
 
     // We can override Jacobian class from dae-cpp library and provide
     // analytical Jacobian. But we will use numerically estimated one.
@@ -211,7 +209,11 @@ int solution_check(dae::state_type &x, MKL_INT N, double t, double D)
               << "% deviation from the analytical value)\n";
     std::cout << "Maximum relative error: " << err_max << "%\n";
 
-    if(err_max < 1.0 && err_conc < 1.0e-6)
+#ifdef DAE_SINGLE
+    if(err_max < 1.0 && err_conc < 1.0e-5)
+#else
+    if(err_max < 1.0 && err_conc < 1.0e-10)
+#endif
         return 0;
     else
         return 1;
diff --git a/examples/perovskite/perovskite.cpp b/examples/perovskite/perovskite.cpp
@@ -99,9 +99,7 @@ int main()
     // parameters defined in solver_options.h
     dae::SolverOptions opt;
 
-    opt.time_stepping   = 1;      // Choose Stability-based time stepper
-    opt.fact_every_iter = false;  // Gain some speed. The matrices will be
-                                  // factorized only once each time step.
+    opt.time_stepping = 1;  // Choose Stability-based time stepper
 
     // Create an instance of the solver with particular RHS, Mass matrix,
     // Jacobian and solver options
diff --git a/test_single.sh b/test_single.sh
@@ -0,0 +1,6 @@
+rm -r build/
+mkdir build
+cd build
+cmake -DDAE_SINGLE=ON ..
+make -j 4
+ctest
