+// Copyright (C) 2007-2008 CEA/DEN, EDF R&D
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
#ifndef __INTERPOLATION3D_TXX__
#define __INTERPOLATION3D_TXX__
case 8:
- intersector = new IntersectorHexa<MyMeshType>(srcMesh, targetMesh, targetIdx,getSplittingPolicy());
+ intersector = new IntersectorHexa<MyMeshType>(srcMesh, targetMesh, targetIdx,getSplittingPolicy());
break;
default:
+// Copyright (C) 2007-2008 CEA/DEN, EDF R&D
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
#ifndef __INTERPOLATIONPLANAR_TXX__
#define __INTERPOLATIONPLANAR_TXX__
_dimCaracteristic=std::min(DimCaracteristic_S, DimCaracteristic_P);
if (InterpolationOptions::getPrintLevel()>=1)
- {
- std::cout << " - Characteristic size of the source mesh : " << DimCaracteristic_S << std::endl;
- std::cout << " - Characteristic size of the target mesh: " << DimCaracteristic_P << std::endl;
- std::cout << "InterpolationPlanar::computation of the intersections" << std::endl;
- }
-
+ {
+ std::cout << " - Characteristic size of the source mesh : " << DimCaracteristic_S << std::endl;
+ std::cout << " - Characteristic size of the target mesh: " << DimCaracteristic_P << std::endl;
+ std::cout << "InterpolationPlanar::computation of the intersections" << std::endl;
+ }
+
PlanarIntersector<MyMeshType>* intersector;
-
+
switch (InterpolationOptions::getIntersectionType())
- {
- case Triangulation:
- intersector=new TriangulationIntersector<MyMeshType>(
- myMesh_P,
- myMesh_S,
- _dimCaracteristic,
- InterpolationOptions::getPrecision(),
- InterpolationOptions::getMedianPlane(),
- InterpolationOptions::getPrintLevel());
- break;
- case Convex:
- intersector=new ConvexIntersector<MyMeshType>(
- myMesh_P,
- myMesh_S,
- _dimCaracteristic,
- InterpolationOptions::getPrecision(),
- InterpolationOptions::getDoRotate(),
- InterpolationOptions::getMedianPlane(),
- InterpolationOptions::getPrintLevel());
- break;
- case Geometric2D:
- intersector=new Geometric2DIntersector<MyMeshType>(myMesh_P, myMesh_S, _dimCaracteristic, InterpolationOptions::getPrecision());
- break;
- // case MEDMEM::Generic:
- //intersector=new GenericIntersector<SPACEDIM>(myMesh_P,myMesh_S, _DimCaracteristic,_Precision,
- // 0, 0, _PrintLevel);
- //break;
- }
+ {
+ case Triangulation:
+ intersector=new TriangulationIntersector<MyMeshType>(
+ myMesh_P,
+ myMesh_S,
+ _dimCaracteristic,
+ InterpolationOptions::getPrecision(),
+ InterpolationOptions::getMedianPlane(),
+ InterpolationOptions::getPrintLevel());
+ break;
+ case Convex:
+ intersector=new ConvexIntersector<MyMeshType>(
+ myMesh_P,
+ myMesh_S,
+ _dimCaracteristic,
+ InterpolationOptions::getPrecision(),
+ InterpolationOptions::getDoRotate(),
+ InterpolationOptions::getMedianPlane(),
+ InterpolationOptions::getPrintLevel());
+ break;
+ case Geometric2D:
+ intersector=new Geometric2DIntersector<MyMeshType>(myMesh_P, myMesh_S, _dimCaracteristic, InterpolationOptions::getPrecision());
+ break;
+ // case MEDMEM::Generic:
+ //intersector=new GenericIntersector<SPACEDIM>(myMesh_P,myMesh_S, _DimCaracteristic,_Precision,
+ // 0, 0, _PrintLevel);
+ //break;
+ }
/****************************************************************/
/* Create a search tree based on the bounding boxes */
/* Instanciate the intersector and initialise the result vector */
/****************************************************************/
-
+
long start_filtering=clock();
std::vector<double> bbox;