src/MEDCoupling/MEDCouplingGaussLocalization.hxx

   1 // Copyright (C) 2007-2013  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.
   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 __PARAMEDMEM_MEDCOUPLINGGAUSSLOCALIZATION_HXX__
  22 #define __PARAMEDMEM_MEDCOUPLINGGAUSSLOCALIZATION_HXX__
  23
  24 #include "MEDCoupling.hxx"
  25 #include "NormalizedUnstructuredMesh.hxx"
  26 #include "InterpKernelException.hxx"
  27
  28 #include <vector>
  29
  30 namespace ParaMEDMEM
  31 {
  32   class MEDCouplingMesh;
  33
  34   class MEDCOUPLING_EXPORT MEDCouplingGaussLocalization
  35   {
  36   public:
  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);
  52     //
  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);
  65     //
  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);
  68   private:
  69     int checkCoherencyOfRequest(int gaussPtIdInCell, int comp) const throw(INTERP_KERNEL::Exception);
  70   private:
  71     INTERP_KERNEL::NormalizedCellType _type;
  72     std::vector<double> _ref_coord;
  73     std::vector<double> _gauss_coord;
  74     std::vector<double> _weight;
  75   };
  76 }
  77
  78 #endif