inclusion du code de Pierre-Yves, gestion de différents types de mesh avec...

inclusion du code de Pierre-Yves, gestion de différents types de mesh avec ElementNodeData  avec element order 1 et 2

test_bem/main.cpp

@@ -8,16 +8,45 @@
 # include <stdio.h>
 #include "code_kevin.hpp"
 
-
 using namespace std;
 using namespace gmsh;
 using namespace Eigen;
 
 int main(int argc, char **argv)
 {
-    //cout << setprecision(2) << fixed;
+    if (argc < 2)
+    {
+        std::cout << "Usage: " << argv[0] << " <geo_file>\n";
+        return 0;
+    }
+
+    initialize();
+
+    open(argv[1]);
+    model::mesh::generate(2);
+
+    vector<int> elementTypes;
+    vector<std::vector<size_t>> elementTags;
+    vector<std::vector<size_t>> nodeTags;
+    model::mesh::getElements(elementTypes, elementTags, nodeTags, 1);
+
+    int nb1dElements = 0;
+    for(int i = 0; i < elementTypes.size(); i++)
+    {
+        nb1dElements += elementTags[i].size();
+    }
+
+    vector<size_t> finalElementTags(nb1dElements);
+    vector<double> electrostaticPressure(nb1dElements);
+
+    solverBEMFun(argc, argv, finalElementTags, electrostaticPressure);
+
+//    for(int i = 0; i < finalElementTags.size(); i++)
+//    {
+//        cout << "Element tag '" << finalElementTags[i] << "'\t-> Electrostatic pressure = " << electrostaticPressure[i] << endl;
+//    }
 
-    code_py(argc,argv);
+    finalize();
 
     return 0;
 }