1 // Copyright (C) 2007-2012 CEA/DEN, EDF R&D
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 #ifndef __PlanarIntersectorP1P0Bary_HXX__
21 #define __PlanarIntersectorP1P0Bary_HXX__
23 #include "PlanarIntersector.hxx"
25 namespace INTERP_KERNEL
27 template<class MyMeshType, class MyMatrix, class ConcreteP1P0Intersector>
28 class PlanarIntersectorP1P0Bary : public PlanarIntersector<MyMeshType,MyMatrix>
31 static const int SPACEDIM=MyMeshType::MY_SPACEDIM;
32 static const int MESHDIM=MyMeshType::MY_MESHDIM;
33 typedef typename MyMeshType::MyConnType ConnType;
34 static const NumberingPolicy numPol=MyMeshType::My_numPol;
36 PlanarIntersectorP1P0Bary(const MyMeshType& meshT, const MyMeshType& meshS, double dimCaracteristic, double precision, double md3DSurf, double medianPlane, bool doRotate, int orientation, int printLevel);
38 void intersectCells(ConnType icellT, const std::vector<ConnType>& icellsS, MyMatrix& res);
39 int getNumberOfRowsOfResMatrix() const;
40 int getNumberOfColsOfResMatrix() const;
42 * Contrary to intersectCells method here icellS and icellT are \b not in \b C mode but in mode of MyMeshType.
44 double intersectGeoBary(const std::vector<double>& targetCell,
45 bool targetCellQuadratic,
46 const double * sourceTria,
47 std::vector<double>& res)
48 { return asLeaf().intersectGeoBary(targetCell,targetCellQuadratic,sourceTria,res); }
50 ConcreteP1P0Intersector& asLeaf() { return static_cast<ConcreteP1P0Intersector&>(*this); }