accélération du code BEM

Votre code BEM est anormalement lent. J'ai donc regardé où le code passait son temps. La majorité du CPU était consommée dans l'appel de computeGeometricParam qui utilise des arguments transférés par copie (dont une énorme map).

J'ai donc ajouté des références et le temps CPU total du test couplé clampedMicroBeam.geo est divisé par 2.

srcs/BEM/functionsBEM.cpp

@@ -659,7 +659,7 @@ void getBC(std::map<size_t, int> &allNodeIndex, const std::vector<double> &allCo
     }
 }
 
-void computeGeometricParam(elementStruct &element, std::map<std::size_t, int> allNodeIndex, const std::vector<double> &allCoord, std::pair<double, double> nodalDisp1, std::pair<double, double> nodalDisp2)
+void computeGeometricParam(elementStruct &element, std::map<std::size_t, int> &allNodeIndex, const std::vector<double> &allCoord, std::pair<double, double> &nodalDisp1, std::pair<double, double> &nodalDisp2)
 {
     int i1Tag = element.nodeTags[0];
     int i2Tag = element.nodeTags[1];