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 __INTEGRALUNIFORMINTERSECTOR_HXX__
21 #define __INTEGRALUNIFORMINTERSECTOR_HXX__
23 #include "TargetIntersector.hxx"
27 namespace INTERP_KERNEL
29 template<class MyMeshType, class MyMatrix>
30 class IntegralUniformIntersector : public TargetIntersector<MyMeshType,MyMatrix>
33 typedef typename MyMeshType::MyConnType ConnType;
35 IntegralUniformIntersector(const MyMeshType& mesh, bool isAbs);
36 double performNormalization(double val) const { if(_is_abs) return fabs(val); else return val; }
37 void setFromTo(bool val) { _from_to=val; }
38 void putValueIn(ConnType i, double val, MyMatrix& res) const;
40 const MyMeshType& _mesh;
41 //! if false means fromIntegralUniform if true means toIntegralUniform
46 template<class MyMeshType, class MyMatrix>
47 class IntegralUniformIntersectorP0 : public IntegralUniformIntersector<MyMeshType,MyMatrix>
50 typedef typename MyMeshType::MyConnType ConnType;
52 IntegralUniformIntersectorP0(const MyMeshType& mesh, bool isAbs);
53 int getNumberOfRowsOfResMatrix() const;
54 int getNumberOfColsOfResMatrix() const;
55 void intersectCells(ConnType targetCell, const std::vector<ConnType>& srcCells, MyMatrix& res);
58 template<class MyMeshType, class MyMatrix>
59 class IntegralUniformIntersectorP1 : public IntegralUniformIntersector<MyMeshType,MyMatrix>
62 typedef typename MyMeshType::MyConnType ConnType;
64 IntegralUniformIntersectorP1(const MyMeshType& mesh, bool isAbs);
65 int getNumberOfRowsOfResMatrix() const;
66 int getNumberOfColsOfResMatrix() const;
67 void intersectCells(ConnType targetCell, const std::vector<ConnType>& srcCells, MyMatrix& res);