src/MEDLoader/MEDFileFieldOverView.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 __MEDFILEFIELDOVERVIEW_HXX__
  22 #define __MEDFILEFIELDOVERVIEW_HXX__
  23
  24 #include "MEDCouplingAutoRefCountObjectPtr.hxx"
  25 #include "MEDCouplingRefCountObject.hxx"
  26 #include "MEDCoupling1GTUMesh.hxx"
  27
  28 #include "NormalizedUnstructuredMesh.hxx"
  29 #include "InterpKernelException.hxx"
  30
  31 #include <vector>
  32
  33 namespace ParaMEDMEM
  34 {
  35   class DataArrayInt;
  36   class MEDCouplingMesh;
  37   class MEDFileMesh;
  38   class MEDFileUMesh;
  39   class MEDFileCMesh;
  40   class MEDFileCurveLinearMesh;
  41   class MEDFileFieldGlobs;
  42   class MEDFileFieldGlobsReal;
  43   class MEDFileAnyTypeField1TS;
  44   class MEDFileAnyTypeFieldMultiTS;
  45
  46   class MEDFileMeshStruct : public RefCountObject
  47   {
  48   public:
  49     static MEDFileMeshStruct *New(const MEDFileMesh *mesh);
  50     std::size_t getHeapMemorySize() const;
  51     const MEDFileMesh *getTheMesh() const { return _mesh; }
  52     int getNumberOfNodes() const { return _nb_nodes; }
  53     int getNumberOfElemsOfGeoType(INTERP_KERNEL::NormalizedCellType t) const throw(INTERP_KERNEL::Exception);
  54     int getLevelOfGeoType(INTERP_KERNEL::NormalizedCellType t) const throw(INTERP_KERNEL::Exception);
  55     int getNumberOfLevs() const;
  56     int getNumberOfGeoTypesInLev(int relativeLev) const throw(INTERP_KERNEL::Exception);
  57   private:
  58     MEDFileMeshStruct(const MEDFileMesh *mesh);
  59   private:
  60     const MEDFileMesh *_mesh;
  61     std::string _name;
  62     int _nb_nodes;
  63     std::vector< std::vector<int> > _geo_types_distrib;
  64   };
  65
  66   class MEDFileField1TSStructItem;
  67
  68   class MEDMeshMultiLev : public RefCountObject
  69   {
  70   public:
  71     std::size_t getHeapMemorySize() const;
  72   public:
  73     static MEDMeshMultiLev *New(const MEDFileMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities) throw(INTERP_KERNEL::Exception);
  74     static MEDMeshMultiLev *New(const MEDFileMesh *m, const std::vector<int>& levs) throw(INTERP_KERNEL::Exception);
  75     static MEDMeshMultiLev *NewOnlyOnNode(const MEDFileMesh *m, const DataArrayInt *pflOnNode) throw(INTERP_KERNEL::Exception);
  76     void setNodeReduction(const DataArrayInt *nr);
  77     bool isFastlyTheSameStruct(const MEDFileField1TSStructItem& fst, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
  78     DataArray *buildDataArray(const MEDFileField1TSStructItem& fst, const MEDFileFieldGlobsReal *globs, const DataArray *vals) const throw(INTERP_KERNEL::Exception);
  79     virtual void selectPartOfNodes(const DataArrayInt *pflNodes) throw(INTERP_KERNEL::Exception) = 0;
  80     virtual MEDMeshMultiLev *prepare() const throw(INTERP_KERNEL::Exception) = 0;
  81   protected:
  82     std::string getPflNameOfId(int id) const;
  83     DataArray *constructDataArray(const MEDFileField1TSStructItem& fst, const MEDFileFieldGlobsReal *globs, const DataArray *vals) const throw(INTERP_KERNEL::Exception);
  84   protected:
  85     MEDMeshMultiLev();
  86     MEDMeshMultiLev(const MEDMeshMultiLev& other);
  87     MEDMeshMultiLev(const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities);
  88   protected:
  89     std::vector< MEDCouplingAutoRefCountObjectPtr<DataArrayInt> > _pfls;
  90     std::vector< INTERP_KERNEL::NormalizedCellType > _geo_types;
  91     std::vector<int> _nb_entities;
  92     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _node_reduction;
  93   };
  94
  95   class MEDStructuredMeshMultiLev;
  96
  97   class MEDUMeshMultiLev : public MEDMeshMultiLev
  98   {
  99   public:
 100     static MEDUMeshMultiLev *New(const MEDFileUMesh *m, const std::vector<int>& levs) throw(INTERP_KERNEL::Exception);
 101     static MEDUMeshMultiLev *New(const MEDFileUMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities) throw(INTERP_KERNEL::Exception);
 102     void selectPartOfNodes(const DataArrayInt *pflNodes) throw(INTERP_KERNEL::Exception);
 103     MEDMeshMultiLev *prepare() const throw(INTERP_KERNEL::Exception);
 104     MEDUMeshMultiLev(const MEDStructuredMeshMultiLev& other, const MEDCouplingAutoRefCountObjectPtr<MEDCoupling1GTUMesh>& part);
 105   private:
 106     MEDUMeshMultiLev(const MEDUMeshMultiLev& other);
 107     MEDUMeshMultiLev(const MEDFileUMesh *m, const std::vector<int>& levs);
 108     MEDUMeshMultiLev(const MEDFileUMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities);
 109   private:
 110     std::vector< MEDCouplingAutoRefCountObjectPtr<MEDCoupling1GTUMesh> > _parts;
 111   };
 112
 113   class MEDStructuredMeshMultiLev : public MEDMeshMultiLev
 114   {
 115   public:
 116     void selectPartOfNodes(const DataArrayInt *pflNodes) throw(INTERP_KERNEL::Exception);
 117     virtual std::vector<int> getNodeGridStructure() const throw(INTERP_KERNEL::Exception) = 0;
 118   protected:
 119     MEDStructuredMeshMultiLev();
 120     MEDStructuredMeshMultiLev(const MEDStructuredMeshMultiLev& other);
 121     MEDStructuredMeshMultiLev(const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities);
 122   };
 123
 124   class MEDCMeshMultiLev : public MEDStructuredMeshMultiLev
 125   {
 126   public:
 127     static MEDCMeshMultiLev *New(const MEDFileCMesh *m, const std::vector<int>& levs) throw(INTERP_KERNEL::Exception);
 128     static MEDCMeshMultiLev *New(const MEDFileCMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities) throw(INTERP_KERNEL::Exception);
 129     std::vector<int> getNodeGridStructure() const throw(INTERP_KERNEL::Exception);
 130     MEDMeshMultiLev *prepare() const throw(INTERP_KERNEL::Exception);
 131   private:
 132     MEDCMeshMultiLev(const MEDCMeshMultiLev& other);
 133     MEDCMeshMultiLev(const MEDFileCMesh *m, const std::vector<int>& levs);
 134     MEDCMeshMultiLev(const MEDFileCMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities);
 135   private:
 136     std::vector< MEDCouplingAutoRefCountObjectPtr<DataArrayDouble> > _coords;
 137   };
 138
 139   class MEDCurveLinearMeshMultiLev : public MEDStructuredMeshMultiLev
 140   {
 141   public:
 142     static MEDCurveLinearMeshMultiLev *New(const MEDFileCurveLinearMesh *m, const std::vector<int>& levs) throw(INTERP_KERNEL::Exception);
 143     static MEDCurveLinearMeshMultiLev *New(const MEDFileCurveLinearMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls , const std::vector<int>& nbEntities) throw(INTERP_KERNEL::Exception);
 144     std::vector<int> getNodeGridStructure() const throw(INTERP_KERNEL::Exception);
 145     MEDMeshMultiLev *prepare() const throw(INTERP_KERNEL::Exception);
 146   private:
 147     MEDCurveLinearMeshMultiLev(const MEDCurveLinearMeshMultiLev& other);
 148     MEDCurveLinearMeshMultiLev(const MEDFileCurveLinearMesh *m, const std::vector<int>& levs);
 149     MEDCurveLinearMeshMultiLev(const MEDFileCurveLinearMesh *m, const std::vector<INTERP_KERNEL::NormalizedCellType>& gts, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& nbEntities);
 150   private:
 151     MEDCouplingAutoRefCountObjectPtr<DataArrayDouble> _coords;
 152     std::vector<int> _structure;
 153   };
 154
 155   class MEDFileField1TSStructItem2 : public RefCountObject
 156   {
 157   public:
 158     MEDFileField1TSStructItem2();
 159     MEDFileField1TSStructItem2(INTERP_KERNEL::NormalizedCellType a, const std::pair<int,int>& b, const std::string& pfl, const std::string& loc);
 160     void checkWithMeshStructForCells(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 161     void checkWithMeshStructForGaussNE(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 162     void checkWithMeshStructForGaussPT(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 163     //
 164     std::size_t getHeapMemorySize() const;
 165     //
 166     const DataArrayInt *getPfl(const MEDFileFieldGlobsReal *globs) const;
 167     INTERP_KERNEL::NormalizedCellType getGeo() const { return _geo_type; }
 168     int getNbEntity() const { return _nb_of_entity; }
 169     const std::pair<int,int>& getStartStop() const { return _start_end; }
 170     std::string getPflName() const;
 171     int getNbOfIntegrationPts(const MEDFileFieldGlobsReal *globs) const;
 172     //! warning this method also set _nb_of_entity attribute !
 173     void checkInRange(int nbOfEntity, int nip, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 174     bool isFastlyEqual(int& startExp, INTERP_KERNEL::NormalizedCellType gt, const char *pflName) const;
 175     bool operator==(const MEDFileField1TSStructItem2& other) const throw(INTERP_KERNEL::Exception);
 176     bool isCellSupportEqual(const MEDFileField1TSStructItem2& other, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 177     bool isNodeSupportEqual(const MEDFileField1TSStructItem2& other, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 178     static MEDFileField1TSStructItem2 BuildAggregationOf(const std::vector<const MEDFileField1TSStructItem2 *>& objs, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 179   public:
 180     static const char NEWLY_CREATED_PFL_NAME[];
 181   private:
 182     INTERP_KERNEL::NormalizedCellType _geo_type;
 183     std::pair<int,int> _start_end;
 184     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _pfl;
 185     std::string _loc;
 186     int _nb_of_entity;
 187   };
 188
 189   class MEDFileField1TSStructItem : public RefCountObject
 190   {
 191   public:
 192     MEDFileField1TSStructItem(TypeOfField a, const std::vector< MEDFileField1TSStructItem2 >& b);
 193     void checkWithMeshStruct(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 194     bool operator==(const MEDFileField1TSStructItem& other) const throw(INTERP_KERNEL::Exception);
 195     std::size_t getHeapMemorySize() const;
 196     bool isEntityCell() const;
 197     bool isComputed() const { return _computed; }
 198     TypeOfField getType() const { return _type; }
 199     std::size_t getNumberOfItems() const { return _items.size(); }
 200     const MEDFileField1TSStructItem2& operator[](std::size_t i) const throw(INTERP_KERNEL::Exception);
 201     //
 202     bool isCellSupportEqual(const MEDFileField1TSStructItem& other, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 203     bool isNodeSupportEqual(const MEDFileField1TSStructItem& other, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 204     MEDFileField1TSStructItem simplifyMeOnCellEntity(const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 205     bool isCompatibleWithNodesDiscr(const MEDFileField1TSStructItem& other, const MEDFileMeshStruct *meshSt, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 206     bool isFullyOnOneLev(const MEDFileMeshStruct *meshSt, int& theFirstLevFull) const throw(INTERP_KERNEL::Exception);
 207     MEDMeshMultiLev *buildFromScratchDataSetSupportOnCells(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 208   private:
 209     bool _computed;
 210     TypeOfField _type;
 211     std::vector< MEDFileField1TSStructItem2 > _items;
 212   };
 213
 214   class MEDFileField1TSStruct : public RefCountObject
 215   {
 216   public:
 217     static MEDFileField1TSStruct *New(const MEDFileAnyTypeField1TS *ref, MEDFileMeshStruct *mst) throw(INTERP_KERNEL::Exception);
 218     void checkWithMeshStruct(MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) throw(INTERP_KERNEL::Exception);
 219     std::size_t getHeapMemorySize() const;
 220     bool isEqualConsideringThePast(const MEDFileAnyTypeField1TS *other, const MEDFileMeshStruct *mst) const throw(INTERP_KERNEL::Exception);
 221     bool isSupportSameAs(const MEDFileAnyTypeField1TS *other, const MEDFileMeshStruct *meshSt) throw(INTERP_KERNEL::Exception);
 222     bool isCompatibleWithNodesDiscr(const MEDFileAnyTypeField1TS *other, const MEDFileMeshStruct *meshSt) throw(INTERP_KERNEL::Exception);
 223     MEDMeshMultiLev *buildFromScratchDataSetSupport(const MEDFileMeshStruct *mst, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 224     bool isDataSetSupportFastlyEqualTo(const MEDFileField1TSStruct& other, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 225     static MEDFileField1TSStructItem BuildItemFrom(const MEDFileAnyTypeField1TS *ref, const MEDFileMeshStruct *meshSt);
 226   private:
 227     MEDFileField1TSStruct(const MEDFileAnyTypeField1TS *ref, MEDFileMeshStruct *mst);
 228     bool presenceOfCellDiscr(int& pos) const throw(INTERP_KERNEL::Exception);
 229     bool presenceOfPartialNodeDiscr(int& pos) const throw(INTERP_KERNEL::Exception);
 230   private:
 231     std::vector<MEDFileField1TSStructItem> _already_checked;
 232   };
 233
 234   class MEDFileFastCellSupportComparator : public RefCountObject
 235   {
 236   public:
 237     static MEDFileFastCellSupportComparator *New(const MEDFileMeshStruct *m, const MEDFileAnyTypeFieldMultiTS *ref) throw(INTERP_KERNEL::Exception);
 238     MEDMeshMultiLev *buildFromScratchDataSetSupport(int timeStepId, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 239     bool isDataSetSupportEqualToThePreviousOne(int timeStepId, const MEDFileFieldGlobsReal *globs) const throw(INTERP_KERNEL::Exception);
 240     bool isEqual(const MEDFileAnyTypeFieldMultiTS *other) throw(INTERP_KERNEL::Exception);
 241     bool isCompatibleWithNodesDiscr(const MEDFileAnyTypeFieldMultiTS *other) throw(INTERP_KERNEL::Exception);
 242     std::size_t getHeapMemorySize() const;
 243   private:
 244     MEDFileFastCellSupportComparator(const MEDFileMeshStruct *m, const MEDFileAnyTypeFieldMultiTS *ref);
 245   private:
 246     MEDCouplingAutoRefCountObjectPtr<MEDFileMeshStruct> _mesh_comp;
 247     std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileField1TSStruct> > _f1ts_cmps;
 248   };
 249 }
 250
 251 #endif