Fixing numerics on GCC6

src/FiniteDifference/eq_unit_test.hpp

@@ -242,15 +242,20 @@ template<typename solver_type,typename lid_nn> void lid_driven_cavity_2d()
 
 #else
 
-	#if __GNUC__ == 5
+	#if __GNUC__ == 6
 
-		std::string file1 = std::string("test/") + s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + "_test_GCC5.vtk";
-		std::string file2 = s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + ".vtk";
+	std::string file1 = std::string("test/") + s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + "_test_GCC6.vtk";
+	std::string file2 = s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + ".vtk";
+
+	#elif __GNUC__ == 5
+
+	std::string file1 = std::string("test/") + s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + "_test_GCC5.vtk";
+	std::string file2 = s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + ".vtk";
 
 	#else
 
-		std::string file1 = std::string("test/") + s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + "_test_GCC4.vtk";
-		std::string file2 = s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + ".vtk";
+	std::string file1 = std::string("test/") + s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + "_test_GCC4.vtk";
+	std::string file2 = s + "lid_driven_cavity_p" + std::to_string(v_cl.getProcessingUnits()) + "_grid_" + std::to_string(v_cl.getProcessUnitID()) + ".vtk";
 
 	#endif
 

src/Solvers/petsc_solver.hpp

@@ -1336,6 +1336,19 @@ public:
 		return true;
 	}
 
+	/*! \brief this function give you the possibility to set PETSC options
+	 *
+	 * this function call PetscOptionsSetValue
+	 *
+	 * \param name the name of the option
+	 * \param value the value of the option
+	 *
+	 */
+	void setPetscOption(const char * name, const char * value)
+	{
+		PetscOptionsSetValue(name,value);
+	}
+
 	/*! \brief Try to solve the system using all the solvers and generate a report
 	 *
 	 * In this mode the system will try different Solvers, Preconditioner and