bug fixed

srcs/FEM/mainFEM.cpp

@@ -27,8 +27,6 @@ int main(int argc, char **argv)
 
     solverFEM(argc, argv, electrostaticPressure,0);
 
-    // first commit
-
     gmsh::fltk::run();
     
     gmsh::finalize();

srcs/FEM/my_geo.geo

@@ -19,7 +19,7 @@ Transfinite Surface {1};
 
 Recombine Surface {1}; // quads instead of triangles
 
-Mesh.ElementOrder = 3;
+Mesh.ElementOrder = 1;
 
 Physical Curve("left_edge", 5) = {4,5};
 Physical Surface("domain", 6) = {1};
@@ -35,11 +35,11 @@ SetNumber("Materials/domain/Poisson", 0.3);
 SetNumber("Materials/domain/rho",7800); //volumic mass of acier
 SetNumber("Boundary Conditions/top_edge/tx", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
 SetNumber("Boundary Conditions/top_edge/ty", 0.); //set to some non-zero value to induce vertical deflection
-//SetNumber("Boundary Conditions/right_edge/tx", 21e3); // for simple tension conditions
+SetNumber("Boundary Conditions/right_edge/tx", 21e3); // for simple tension conditions
 SetNumber("Boundary Conditions/right_edge/ty", 0.);
 SetNumber("Volumic Forces/domain/bx",0.);
 SetNumber("Volumic Forces/domain/by",0.); //set to -9.81 for gravity
-//SetNumber("Boundary Conditions/fixed_node/uy",2.);
-SetNumber("Boundary Conditions/fixed_node/uy",0.);
+SetNumber("Boundary Conditions/fixed_node/uy",2.);
+//SetNumber("Boundary Conditions/fixed_node/uy",0.);
 
 Physical Curve("BEM_FEM_boundary", 10) = {1,2,3};
\ No newline at end of file

srcs/FEM/solverFEM.cpp

@@ -157,7 +157,8 @@ int solverFEM(int argc, char **argv, std::map<int, double> &electrostaticPressur
     //std::cout << "time for neumann: " << start3 - start2 << "\n";
 
     /* ELECTROSTATIC PRESSURE */
-    applyElecPressure(integration_rule, node_map, groups, electrostaticPressure, full_f_vector);
+    if(electrostaticPressure.size() != 0)
+        applyElecPressure(integration_rule, node_map, groups, electrostaticPressure, full_f_vector);
     //double start4 = omp_get_wtime();
     //std::cout << "time for elec pressure: " << start4 - start3 << "\n";