X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FINTERP_KERNEL%2FInterpolation2D1D.txx;h=7cc6ca44393db0bc5587e8751b409275c31ad112;hb=b219559763498c4bd10c730cd3d2c62b1eed45db;hp=1e644d7813c31a94a878a92aa89943d0411f4f47;hpb=48d2e4c6fa2057fcc6e9d8474eb73516886c861b;p=tools%2Fmedcoupling.git diff --git a/src/INTERP_KERNEL/Interpolation2D1D.txx b/src/INTERP_KERNEL/Interpolation2D1D.txx old mode 100644 new mode 100755 index 1e644d781..7cc6ca443 --- a/src/INTERP_KERNEL/Interpolation2D1D.txx +++ b/src/INTERP_KERNEL/Interpolation2D1D.txx @@ -46,19 +46,19 @@ namespace INTERP_KERNEL * */ template - int Interpolation2D1D::interpolateMeshes(const MyMeshType& myMeshS, const MyMeshType& myMeshT, MatrixType& result, const std::string& method) + typename MyMeshType::MyConnType Interpolation2D1D::interpolateMeshes(const MyMeshType& myMeshS, const MyMeshType& myMeshT, MatrixType& result, const std::string& method) { static const int SPACEDIM=MyMeshType::MY_SPACEDIM; typedef typename MyMeshType::MyConnType ConnType; static const NumberingPolicy numPol=MyMeshType::My_numPol; long global_start =clock(); - int counter=0; + std::size_t counter=0; /***********************************************************/ /* Check both meshes are made of triangles and quadrangles */ /***********************************************************/ - long nbMailleS=myMeshS.getNumberOfElements(); + ConnType nbMailleS=myMeshS.getNumberOfElements(); /**************************************************/ /* Search the characteristic size of the meshes */ @@ -114,12 +114,12 @@ namespace INTERP_KERNEL /* Loop on the target cells - core of the algorithm */ /****************************************************/ long start_intersection=clock(); - long nbelem_type=myMeshT.getNumberOfElements(); + ConnType nbelem_type=myMeshT.getNumberOfElements(); const ConnType *connIndxT=myMeshT.getConnectivityIndexPtr(); - for(int iT=0; iT intersecting_elems; + ConnType nb_nodesT=connIndxT[iT+1]-connIndxT[iT]; + std::vector intersecting_elems; double bb[2*SPACEDIM]; intersector->getElemBB(bb,myMeshT,OTT::indFC(iT),nb_nodesT); my_tree.getIntersectingElems(bb, intersecting_elems); @@ -127,7 +127,7 @@ namespace INTERP_KERNEL counter+=intersecting_elems.size(); intersecting_elems.clear(); } - int ret=intersector->getNumberOfColsOfResMatrix(); + ConnType ret=intersector->getNumberOfColsOfResMatrix(); const DuplicateFacesType& intersectFaces = *intersector->getIntersectFaces(); DuplicateFacesType::const_iterator iter;