Salome HOME
4356a1c214c933e852501aca9d3c73d019c721b7
[tools/medcoupling.git] / src / INTERP_KERNEL / Planar2D1DIntersectorP0P0.hxx
1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D
2 //
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, or (at your option) any later version.
7 //
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.
12 //
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
16 //
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
18 //
19 // Author : Anthony Geay (CEA/DEN)
20
21 #ifndef __PLANAR2D1DINTERSECTORP0P0_HXX__
22 #define __PLANAR2D1DINTERSECTORP0P0_HXX__
23
24 #include "PlanarIntersector.hxx"
25
26 namespace INTERP_KERNEL
27 {
28   template<class MyMeshType, class MyMatrix, class ConcreteP0P0Intersector>
29   class Planar2D1DIntersectorP0P0 : public PlanarIntersector<MyMeshType,MyMatrix>
30   {
31   public:
32     static const int SPACEDIM=MyMeshType::MY_SPACEDIM;
33     static const int MESHDIM=MyMeshType::MY_MESHDIM;
34     typedef typename MyMeshType::MyConnType ConnType;
35     static const NumberingPolicy numPol=MyMeshType::My_numPol;
36   protected:
37     Planar2D1DIntersectorP0P0(const MyMeshType& meshT, const MyMeshType& meshS,
38                               double dimCaracteristic, double precision, double md3DSurf, double minDot3DSurf, double medianPlane, bool doRotate, int orientation, int printLevel);
39   public:
40     int getNumberOfRowsOfResMatrix() const;
41     int getNumberOfColsOfResMatrix() const;
42     const typename PlanarIntersector<MyMeshType,MyMatrix>::DuplicateFacesType* getIntersectFaces() const
43     {
44       return &_intersect_faces;
45     }
46     void intersectCells(ConnType icellT, const std::vector<ConnType>& icellsS, MyMatrix& res);
47     /*!
48      * Contrary to intersectCells method here icellS and icellT are \b not in \b C mode but in mode of MyMeshType.
49      */
50     double intersectGeometry1D(ConnType icellT, ConnType icellS, ConnType nbNodesT, ConnType nbNodesS,
51                                bool& isColinear)
52       { return asLeaf().intersectGeometry1D(icellT,icellS,nbNodesT,nbNodesS, isColinear); }
53   protected:
54     ConcreteP0P0Intersector& asLeaf() { return static_cast<ConcreteP0P0Intersector&>(*this); }
55   private:
56     typename PlanarIntersector<MyMeshType,MyMatrix>::DuplicateFacesType _intersect_faces;
57   };
58 }
59
60 #endif