diff --git a/example/Grid/3_gray_scott_3d/main.cpp b/example/Grid/3_gray_scott_3d/main.cpp
index a72029fc694b2ae24f3b560e45c1a78d6e305b40..82ea1fd77b2cde05b6954fc8e63d608791104a02 100644
--- a/example/Grid/3_gray_scott_3d/main.cpp
+++ b/example/Grid/3_gray_scott_3d/main.cpp
@@ -230,7 +230,7 @@ int main(int argc, char* argv[])
 		// visualization
 		if (i % 500 == 0)
 		{
-			Old.save("output_" + std::to_string(count));
+//			Old.save("output_" + std::to_string(count));
 			count++;
 		}
 	}
diff --git a/example/Numerics/Vortex_in_cell/main_vic_petsc_opt.cpp b/example/Numerics/Vortex_in_cell/main_vic_petsc_opt.cpp
index 5c4d29a81579e50baa687071f2d78f281ba18cbd..0eb8c9d43ab591f2ce1eee4304ed1090f279edc5 100644
--- a/example/Numerics/Vortex_in_cell/main_vic_petsc_opt.cpp
+++ b/example/Numerics/Vortex_in_cell/main_vic_petsc_opt.cpp
@@ -505,9 +505,9 @@ template<typename grid> void calc_rhs(grid & g_vort, grid & g_vel, grid & g_dwp)
 
 	// calculate several pre-factors for the stencil finite
 	// difference
-	float fac1 = 1.0f*nu/(g_vort.spacing(0)*g_vort.spacing(0));
-	float fac2 = 1.0f*nu/(g_vort.spacing(1)*g_vort.spacing(1));
-	float fac3 = 1.0f*nu/(g_vort.spacing(2)*g_vort.spacing(2));
+	float fac1 = 2.0f*nu/(g_vort.spacing(0)*g_vort.spacing(0));
+	float fac2 = 2.0f*nu/(g_vort.spacing(1)*g_vort.spacing(1));
+	float fac3 = 2.0f*nu/(g_vort.spacing(2)*g_vort.spacing(2));
 
 	float fac4 = 0.5f/(g_vort.spacing(0));
 	float fac5 = 0.5f/(g_vort.spacing(1));
diff --git a/example/SparseGrid/4_gray_scott_3d_sparse_surface_cs/main.cpp b/example/SparseGrid/4_gray_scott_3d_sparse_surface_cs/main.cpp
index 33eea4a727ef25971bafc929996a3dbddaca924e..ab932e019c638d5203b0804982e50a8b7084ad87 100644
--- a/example/SparseGrid/4_gray_scott_3d_sparse_surface_cs/main.cpp
+++ b/example/SparseGrid/4_gray_scott_3d_sparse_surface_cs/main.cpp
@@ -188,13 +188,13 @@ int main(int argc, char* argv[])
 	// Diffusion constant for specie V
 	double dv = 0.5*1e-5;
 
-//#ifdef TEST_RUN
+#ifdef TEST_RUN
         // Number of timesteps
-//        size_t timeSteps = 200;
-//#else
+        size_t timeSteps = 200;
+#else
 	// Number of timesteps
         size_t timeSteps = 150000;
-//#endif
+#endif
 
 	// K and F (Physical constant in the equation)
         double K = 0.053;
diff --git a/example/SparseGrid/5_gray_scott_3d_surface_cs_opt/main.cpp b/example/SparseGrid/5_gray_scott_3d_surface_cs_opt/main.cpp
index 8dba1fcdf1bd4697ab4b3bfb9e5ba825c3299ea4..53e1c15f7b83437b2a37cca40479a2f14052e729 100644
--- a/example/SparseGrid/5_gray_scott_3d_surface_cs_opt/main.cpp
+++ b/example/SparseGrid/5_gray_scott_3d_surface_cs_opt/main.cpp
@@ -202,13 +202,13 @@ int main(int argc, char* argv[])
 	// Diffusion constant for specie V
 	double dv = 0.5*1e-5;
 
-//#ifdef TEST_RUN
+#ifdef TEST_RUN
         // Number of timesteps
-//        size_t timeSteps = 200;
-//#else
+        size_t timeSteps = 200;
+#else
 	// Number of timesteps
         size_t timeSteps = 100000;
-//#endif
+#endif
 
 	// K and F (Physical constant in the equation)
         double K = 0.053;