From: eap <eap@opencascade.com>
Date: Wed, 14 Jan 2009 10:24:42 +0000 (+0000)
Subject: MEDMEM Industrialization 2008
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=6bbd7063452f77318cffe388dee3a469044f6f9a;p=tools%2Fmedcoupling.git

MEDMEM Industrialization 2008
  Form interpolation matrix in intersectors

     void interpolateMeshes(const MyMeshType& mesh1, const MyMeshType& mesh2, MatrixType& result,
+                           const std::string& method)
---

diff --git a/src/INTERP_KERNEL/Interpolation.hxx b/src/INTERP_KERNEL/Interpolation.hxx
index fac260389..ad7b7e44b 100644
--- a/src/INTERP_KERNEL/Interpolation.hxx
+++ b/src/INTERP_KERNEL/Interpolation.hxx
@@ -10,6 +10,8 @@
 #include "INTERPKERNEL_defines.hxx"
 #include "InterpolationOptions.hxx"
 
+#include <string>
+
 namespace INTERP_KERNEL
 {
   template<class TrueMainInterpolator>
@@ -28,8 +30,9 @@ namespace INTERP_KERNEL
 		Interpolation(const InterpolationOptions& io) :InterpolationOptions(io){}
     //interpolation of two triangular meshes.
     template<class MatrixType, class MyMeshType>
-    void interpolateMeshes(const MyMeshType& mesh1, const MyMeshType& mesh2, MatrixType& result)
-    { return asLeaf().interpolateMeshes(mesh1,mesh2,result); }
+    void interpolateMeshes(const MyMeshType& mesh1, const MyMeshType& mesh2, MatrixType& result,
+                           const std::string& method)
+    { return asLeaf().interpolateMeshes(mesh1,mesh2,result,method); }
   protected:
     TrueMainInterpolator& asLeaf() { return static_cast<TrueMainInterpolator&>(*this); }
   };
diff --git a/src/INTERP_KERNEL/Interpolation3D.hxx b/src/INTERP_KERNEL/Interpolation3D.hxx
index 89d820d44..05a0a12d0 100644
--- a/src/INTERP_KERNEL/Interpolation3D.hxx
+++ b/src/INTERP_KERNEL/Interpolation3D.hxx
@@ -13,7 +13,8 @@ namespace INTERP_KERNEL
     Interpolation3D();
     Interpolation3D(const InterpolationOptions& io);
     template<class MatrixType, class MyMeshType>
-    void interpolateMeshes(const MyMeshType& srcMesh, const MyMeshType& targetMesh, MatrixType& result);
+    void interpolateMeshes(const MyMeshType& srcMesh, const MyMeshType& targetMesh, MatrixType& result,
+                           const std::string& method);
   private:
     SplittingPolicy _splitting_policy;
   };
diff --git a/src/INTERP_KERNEL/InterpolationPlanar.hxx b/src/INTERP_KERNEL/InterpolationPlanar.hxx
index 111ade9e6..094c7c3d6 100755
--- a/src/INTERP_KERNEL/InterpolationPlanar.hxx
+++ b/src/INTERP_KERNEL/InterpolationPlanar.hxx
@@ -25,7 +25,8 @@ namespace INTERP_KERNEL
 
     // Main function to interpolate triangular and quadratic meshes
     template<class MatrixType, class MyMeshType>
-    void interpolateMeshes(const MyMeshType& mesh1, const MyMeshType& mesh2, MatrixType& result);
+    void interpolateMeshes(const MyMeshType& mesh1, const MyMeshType& mesh2, MatrixType& result,
+                           const std::string& method);
 
   public:
     bool doRotate() const { return asLeafInterpPlanar().doRotate(); }
diff --git a/src/INTERP_KERNEL/InterpolationPlanar.txx b/src/INTERP_KERNEL/InterpolationPlanar.txx
index f9061110f..b41546430 100644
--- a/src/INTERP_KERNEL/InterpolationPlanar.txx
+++ b/src/INTERP_KERNEL/InterpolationPlanar.txx
@@ -249,15 +249,9 @@ namespace INTERP_KERNEL
                 getBarycentricCoordinates<MyMeshType,3>( baryCentre, OTT<ConnType,numPol>::indFC(i_P),
                                                          myMesh_P, baryCoords );
                 if ( nodeCoeff ) { // 1 coeff per target node
-                  cout << "P elem " << OTT<ConnType,numPol>::indFC(i_P)
-                       << ", S elem " << OTT<ConnType,numPol>::indFC(i_S)
-                       << ", surf: " << surf
-                       << endl;
                   for ( int n = 0; n < 3; ++n ) {
                     ConnType node = getGlobalNumberOfNode(n, OTT<ConnType,numPol>::indFC(i_P), myMesh_P);
                     result[node].insert(make_pair(OTT<ConnType,numPol>::indFC(i_S),surf*baryCoords[n]));
-                    cout << "node " << node << " baryCoord[" << n << "]=" << baryCoords[n]
-                         << ", surf*bc = " << surf*baryCoords[n] << endl;
                   }
                 }
                 else { // several coeffs per target node