factor 2 elec pressure

srcs/BEM/solverBEM.cpp

@@ -235,7 +235,7 @@ int solverBEM(int argc, char **argv, std::map<int, double> &electrostaticPressur
     double epsilon = 8.8541878128e-12;
     for(size_t i = 0; i < elementVector.size(); i++)
     {
-        electrostaticPressure[elementVector[i].tag] = epsilon * pow(elementVector[i].bcValue[1], 2.0);
+        electrostaticPressure[elementVector[i].tag] = epsilon * pow(elementVector[i].bcValue[1], 2.0) / 2;
     }
 
     end_total = std::chrono::system_clock::now();