1 // Copyright (C) 2007-2013 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
19 // Author : Anthony Geay (CEA/DEN)
21 #ifndef __PARAMEDMEM_MEDCOUPLINGGAUSSLOCALIZATION_HXX__
22 #define __PARAMEDMEM_MEDCOUPLINGGAUSSLOCALIZATION_HXX__
24 #include "MEDCoupling.hxx"
25 #include "NormalizedUnstructuredMesh.hxx"
26 #include "InterpKernelException.hxx"
32 class MEDCouplingMesh;
34 class MEDCOUPLING_EXPORT MEDCouplingGaussLocalization
37 MEDCouplingGaussLocalization(INTERP_KERNEL::NormalizedCellType type, const std::vector<double>& refCoo,
38 const std::vector<double>& gsCoo, const std::vector<double>& w) throw(INTERP_KERNEL::Exception);
39 MEDCouplingGaussLocalization(INTERP_KERNEL::NormalizedCellType typ) throw(INTERP_KERNEL::Exception);
40 INTERP_KERNEL::NormalizedCellType getType() const { return _type; }
41 void setType(INTERP_KERNEL::NormalizedCellType typ) throw(INTERP_KERNEL::Exception);
42 int getNumberOfGaussPt() const { return (int)_weight.size(); }
43 int getDimension() const;
44 int getNumberOfPtsInRefCell() const;
45 std::string getStringRepr() const;
46 std::size_t getHeapMemorySize() const;
47 void checkCoherency() const throw(INTERP_KERNEL::Exception);
48 bool isEqual(const MEDCouplingGaussLocalization& other, double eps) const;
49 void pushTinySerializationIntInfo(std::vector<int>& tinyInfo) const;
50 void pushTinySerializationDblInfo(std::vector<double>& tinyInfo) const;
51 const double *fillWithValues(const double *vals);
53 const std::vector<double>& getRefCoords() const { return _ref_coord; }
54 double getRefCoord(int ptIdInCell, int comp) const throw(INTERP_KERNEL::Exception);
55 const std::vector<double>& getGaussCoords() const { return _gauss_coord; }
56 double getGaussCoord(int gaussPtIdInCell, int comp) const throw(INTERP_KERNEL::Exception);
57 const std::vector<double>& getWeights() const { return _weight; }
58 double getWeight(int gaussPtIdInCell, double newVal) const throw(INTERP_KERNEL::Exception);
59 void setRefCoord(int ptIdInCell, int comp, double newVal) throw(INTERP_KERNEL::Exception);
60 void setGaussCoord(int gaussPtIdInCell, int comp, double newVal) throw(INTERP_KERNEL::Exception);
61 void setWeight(int gaussPtIdInCell, double newVal) throw(INTERP_KERNEL::Exception);
62 void setRefCoords(const std::vector<double>& refCoo) throw(INTERP_KERNEL::Exception);
63 void setGaussCoords(const std::vector<double>& gsCoo) throw(INTERP_KERNEL::Exception);
64 void setWeights(const std::vector<double>& w) throw(INTERP_KERNEL::Exception);
66 static MEDCouplingGaussLocalization BuildNewInstanceFromTinyInfo(int dim, const std::vector<int>& tinyData);
67 static bool AreAlmostEqual(const std::vector<double>& v1, const std::vector<double>& v2, double eps);
69 int checkCoherencyOfRequest(int gaussPtIdInCell, int comp) const throw(INTERP_KERNEL::Exception);
71 INTERP_KERNEL::NormalizedCellType _type;
72 std::vector<double> _ref_coord;
73 std::vector<double> _gauss_coord;
74 std::vector<double> _weight;