element node data

test_bem/code_kevin.hpp

@@ -34,6 +34,6 @@ void compute_g(double &jMidX, double &jDeltaX, double &jMidY, double &jDeltaY, d
 
 void compute_internal_phi(vector<double> &allCoord,vector<size_t> &NodeTags2d, vector<elementStruct> &elementVector, vector<vector<double>> &data_boundary_nodes, vector<vector<double>> &data2d, map<int,double> &nodeMapDomain);
 
-void en_chantier();
+void element_node_data(vector<vector<vector<double>>> &data, vector<vector<size_t>> &elementTags2d, vector<int> &elementTypes2d, vector<vector<size_t>> &nodeTags2d, vector<elementStruct> &elementVector, map<int, double> &nodeMapDomain, bool &show_element_node_data);
 
 

test_bem/main.cpp

@@ -17,8 +17,6 @@ int main(int argc, char **argv)
 {
     //cout << setprecision(2) << fixed;
 
-    //code_kevin(argc,argv);
-
     code_py(argc,argv);
 
     return 0;

test_bem/my_geo.geo

 Lx = 1;
 Ly = 1;
-nx = 4;
-ny = 4;
+nx = 10;
+ny = 10;
 
 Point(1) = {0, 0, 0, 1.0};
 Point(2) = {Lx, 0, 0, 1.0};
@@ -17,7 +17,7 @@ Transfinite Curve {3, 1} = nx+1 Using Progression 1;
 Transfinite Curve {4, 2} = ny+1 Using Progression 1;
 Transfinite Surface {1};
 
-Recombine Surface {1}; // quads instead of triangles
+//Recombine Surface {1}; // quads instead of triangles
 
 //à gauche c'est le tag qu'on donne, il doit être différent pour chaque curve ou surface et à droite c'est l'élément définit correspondant 4 courbes et 1 surface
 Physical Curve("bottom_edge", 1) = {1};