src/MEDLoader/MEDFileMesh.hxx

   1 // Copyright (C) 2007-2014  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 __MEDFILEMESH_HXX__
  22 #define __MEDFILEMESH_HXX__
  23
  24 #include "MEDLoaderDefines.hxx"
  25 #include "MEDFileMeshLL.hxx"
  26 #include "MEDFileUtilities.hxx"
  27 #include "MEDFileMeshReadSelector.hxx"
  28
  29 #include <map>
  30 #include <list>
  31
  32 namespace ParaMEDMEM
  33 {
  34   class MEDFileFieldGlobsReal;
  35   class MEDFileField1TSStructItem;
  36
  37   class MEDFileMesh : public RefCountObject, public MEDFileWritable
  38   {
  39   public:
  40     MEDLOADER_EXPORT static MEDFileMesh *New(const std::string& fileName, MEDFileMeshReadSelector *mrs=0);
  41     MEDLOADER_EXPORT static MEDFileMesh *New(const std::string& fileName, const std::string& mName, int dt=-1, int it=-1, MEDFileMeshReadSelector *mrs=0);
  42     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
  43     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
  44     MEDLOADER_EXPORT virtual MEDFileMesh *createNewEmpty() const = 0;
  45     MEDLOADER_EXPORT virtual MEDFileMesh *deepCpy() const = 0;
  46     MEDLOADER_EXPORT virtual MEDFileMesh *shallowCpy() const = 0;
  47     MEDLOADER_EXPORT virtual bool isEqual(const MEDFileMesh *other, double eps, std::string& what) const;
  48     MEDLOADER_EXPORT virtual void clearNonDiscrAttributes() const;
  49     MEDLOADER_EXPORT virtual void setName(const std::string& name);
  50     MEDLOADER_EXPORT bool changeNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
  51     MEDLOADER_EXPORT std::string getName() const { return _name; }
  52     MEDLOADER_EXPORT std::string getUnivName() const { return _univ_name; }
  53     MEDLOADER_EXPORT bool getUnivNameWrStatus() const { return _univ_wr_status; }
  54     MEDLOADER_EXPORT void setUnivNameWrStatus(bool newStatus) { _univ_wr_status=newStatus; }
  55     MEDLOADER_EXPORT void setDescription(const std::string& name) { _desc_name=name; }
  56     MEDLOADER_EXPORT std::string getDescription() const { return _desc_name; }
  57     MEDLOADER_EXPORT void setOrder(int order) { _order=order; }
  58     MEDLOADER_EXPORT int getOrder() const { return _order; }
  59     MEDLOADER_EXPORT void setIteration(int it) { _iteration=it; }
  60     MEDLOADER_EXPORT int getIteration() const { return _iteration; }
  61     MEDLOADER_EXPORT void setTimeValue(double time) { _time=time; }
  62     MEDLOADER_EXPORT void setTime(int dt, int it, double time) { _time=time; _iteration=dt; _order=it; }
  63     MEDLOADER_EXPORT double getTime(int& dt, int& it) { dt=_iteration; it=_order; return _time; }
  64     MEDLOADER_EXPORT double getTimeValue() const { return _time; }
  65     MEDLOADER_EXPORT void setTimeUnit(const std::string& unit) { _dt_unit=unit; }
  66     MEDLOADER_EXPORT std::string getTimeUnit() const { return _dt_unit; }
  67     MEDLOADER_EXPORT std::vector<INTERP_KERNEL::NormalizedCellType> getAllGeoTypes() const;
  68     MEDLOADER_EXPORT virtual int getNumberOfNodes() const = 0;
  69     MEDLOADER_EXPORT virtual bool hasImplicitPart() const = 0;
  70     MEDLOADER_EXPORT virtual int buildImplicitPartIfAny(INTERP_KERNEL::NormalizedCellType gt) const = 0;
  71     MEDLOADER_EXPORT virtual void releaseImplicitPartIfAny() const = 0;
  72     MEDLOADER_EXPORT virtual std::vector<INTERP_KERNEL::NormalizedCellType> getGeoTypesAtLevel(int meshDimRelToMax) const = 0;
  73     MEDLOADER_EXPORT virtual std::vector<int> getNonEmptyLevels() const = 0;
  74     MEDLOADER_EXPORT virtual std::vector<int> getNonEmptyLevelsExt() const = 0;
  75     MEDLOADER_EXPORT virtual std::vector<int> getFamArrNonEmptyLevelsExt() const = 0;
  76     MEDLOADER_EXPORT virtual std::vector<int> getNumArrNonEmptyLevelsExt() const = 0;
  77     MEDLOADER_EXPORT virtual std::vector<int> getNameArrNonEmptyLevelsExt() const = 0;
  78     MEDLOADER_EXPORT virtual void write(const std::string& fileName, int mode) const;
  79     MEDLOADER_EXPORT virtual void write(med_idt fid) const;
  80     MEDLOADER_EXPORT virtual int getSizeAtLevel(int meshDimRelToMaxExt) const = 0;
  81     MEDLOADER_EXPORT virtual MEDCouplingMesh *getGenMeshAtLevel(int meshDimRelToMax, bool renum=false) const = 0;
  82     MEDLOADER_EXPORT virtual std::vector<int> getDistributionOfTypes(int meshDimRelToMax) const;
  83     MEDLOADER_EXPORT virtual void whichAreNodesFetched(const MEDFileField1TSStructItem& st, const MEDFileFieldGlobsReal *globs, std::vector<bool>& nodesFetched) const = 0;
  84     //
  85     MEDLOADER_EXPORT bool areFamsEqual(const MEDFileMesh *other, std::string& what) const;
  86     MEDLOADER_EXPORT bool areGrpsEqual(const MEDFileMesh *other, std::string& what) const;
  87     MEDLOADER_EXPORT bool existsGroup(const std::string& groupName) const;
  88     MEDLOADER_EXPORT bool existsFamily(int famId) const;
  89     MEDLOADER_EXPORT bool existsFamily(const std::string& familyName) const;
  90     MEDLOADER_EXPORT void setFamilyId(const std::string& familyName, int id);
  91     MEDLOADER_EXPORT void setFamilyIdUnique(const std::string& familyName, int id);
  92     MEDLOADER_EXPORT virtual void addFamily(const std::string& familyName, int id);
  93     MEDLOADER_EXPORT virtual void createGroupOnAll(int meshDimRelToMaxExt, const std::string& groupName);
  94     MEDLOADER_EXPORT virtual bool keepFamIdsOnlyOnLevs(const std::vector<int>& famIds, const std::vector<int>& levs);
  95     MEDLOADER_EXPORT void addFamilyOnGrp(const std::string& grpName, const std::string& famName);
  96     MEDLOADER_EXPORT std::string findOrCreateAndGiveFamilyWithId(int id, bool& created);
  97     MEDLOADER_EXPORT void setFamilyInfo(const std::map<std::string,int>& info);
  98     MEDLOADER_EXPORT void setGroupInfo(const std::map<std::string, std::vector<std::string> >&info);
  99     MEDLOADER_EXPORT void copyFamGrpMapsFrom(const MEDFileMesh& other);
 100     MEDLOADER_EXPORT void clearGrpMap();
 101     MEDLOADER_EXPORT void clearFamMap();
 102     MEDLOADER_EXPORT void clearFamGrpMaps();
 103     MEDLOADER_EXPORT const std::map<std::string,int>& getFamilyInfo() const { return _families; }
 104     MEDLOADER_EXPORT const std::map<std::string, std::vector<std::string> >& getGroupInfo() const { return _groups; }
 105     MEDLOADER_EXPORT std::vector<std::string> getFamiliesOnGroup(const std::string& name) const;
 106     MEDLOADER_EXPORT std::vector<std::string> getFamiliesOnGroups(const std::vector<std::string>& grps) const;
 107     MEDLOADER_EXPORT std::vector<int> getFamiliesIdsOnGroup(const std::string& name) const;
 108     MEDLOADER_EXPORT void setFamiliesOnGroup(const std::string& name, const std::vector<std::string>& fams);
 109     MEDLOADER_EXPORT void setFamiliesIdsOnGroup(const std::string& name, const std::vector<int>& famIds);
 110     MEDLOADER_EXPORT std::vector<std::string> getGroupsOnFamily(const std::string& name) const;
 111     MEDLOADER_EXPORT void setGroupsOnFamily(const std::string& famName, const std::vector<std::string>& grps);
 112     MEDLOADER_EXPORT std::vector<std::string> getGroupsNames() const;
 113     MEDLOADER_EXPORT std::vector<std::string> getFamiliesNames() const;
 114     MEDLOADER_EXPORT void assignFamilyNameWithGroupName();
 115     MEDLOADER_EXPORT std::vector<std::string> removeEmptyGroups();
 116     MEDLOADER_EXPORT void removeGroup(const std::string& name);
 117     MEDLOADER_EXPORT void removeFamily(const std::string& name);
 118     MEDLOADER_EXPORT std::vector<std::string> removeOrphanGroups();
 119     MEDLOADER_EXPORT std::vector<std::string> removeOrphanFamilies();
 120     MEDLOADER_EXPORT void changeGroupName(const std::string& oldName, const std::string& newName);
 121     MEDLOADER_EXPORT void changeFamilyName(const std::string& oldName, const std::string& newName);
 122     MEDLOADER_EXPORT void changeFamilyId(int oldId, int newId);
 123     MEDLOADER_EXPORT void changeAllGroupsContainingFamily(const std::string& familyNameToChange, const std::vector<std::string>& newFamiliesNames);
 124     MEDLOADER_EXPORT int getFamilyId(const std::string& name) const;
 125     MEDLOADER_EXPORT int getMaxAbsFamilyId() const;
 126     MEDLOADER_EXPORT int getMaxFamilyId() const;
 127     MEDLOADER_EXPORT int getMinFamilyId() const;
 128     MEDLOADER_EXPORT int getTheMaxAbsFamilyId() const;
 129     MEDLOADER_EXPORT int getTheMaxFamilyId() const;
 130     MEDLOADER_EXPORT int getTheMinFamilyId() const;
 131     MEDLOADER_EXPORT virtual int getMaxAbsFamilyIdInArrays() const = 0;
 132     MEDLOADER_EXPORT virtual int getMaxFamilyIdInArrays() const = 0;
 133     MEDLOADER_EXPORT virtual int getMinFamilyIdInArrays() const = 0;
 134     MEDLOADER_EXPORT DataArrayInt *getAllFamiliesIdsReferenced() const;
 135     MEDLOADER_EXPORT DataArrayInt *computeAllFamilyIdsInUse() const;
 136     MEDLOADER_EXPORT std::vector<int> getFamiliesIds(const std::vector<std::string>& famNames) const;
 137     MEDLOADER_EXPORT std::string getFamilyNameGivenId(int id) const;
 138     MEDLOADER_EXPORT bool ensureDifferentFamIdsPerLevel();
 139     MEDLOADER_EXPORT void normalizeFamIdsTrio();
 140     MEDLOADER_EXPORT void normalizeFamIdsMEDFile();
 141     MEDLOADER_EXPORT virtual int getMeshDimension() const = 0;
 142     MEDLOADER_EXPORT virtual std::string simpleRepr() const;
 143     MEDLOADER_EXPORT virtual std::string advancedRepr() const = 0;
 144     //
 145     MEDLOADER_EXPORT virtual void setGroupsAtLevel(int meshDimRelToMaxExt, const std::vector<const DataArrayInt *>& grps, bool renum=false);
 146     MEDLOADER_EXPORT virtual void setFamilyFieldArr(int meshDimRelToMaxExt, DataArrayInt *famArr) = 0;
 147     MEDLOADER_EXPORT virtual void setRenumFieldArr(int meshDimRelToMaxExt, DataArrayInt *renumArr) = 0;
 148     MEDLOADER_EXPORT virtual void setNameFieldAtLevel(int meshDimRelToMaxExt, DataArrayAsciiChar *nameArr) = 0;
 149     MEDLOADER_EXPORT virtual const DataArrayInt *getFamilyFieldAtLevel(int meshDimRelToMaxExt) const = 0;
 150     MEDLOADER_EXPORT virtual const DataArrayInt *getNumberFieldAtLevel(int meshDimRelToMaxExt) const = 0;
 151     MEDLOADER_EXPORT virtual const DataArrayInt *getRevNumberFieldAtLevel(int meshDimRelToMaxExt) const = 0;
 152     MEDLOADER_EXPORT virtual const DataArrayAsciiChar *getNameFieldAtLevel(int meshDimRelToMaxExt) const = 0;
 153     MEDLOADER_EXPORT virtual DataArrayInt *getFamiliesArr(int meshDimRelToMaxExt, const std::vector<std::string>& fams, bool renum=false) const = 0;
 154     MEDLOADER_EXPORT virtual DataArrayInt *getGroupsArr(int meshDimRelToMaxExt, const std::vector<std::string>& grps, bool renum=false) const;
 155     MEDLOADER_EXPORT virtual DataArrayInt *getGroupArr(int meshDimRelToMaxExt, const std::string& grp, bool renum=false) const;
 156     MEDLOADER_EXPORT virtual DataArrayInt *getFamilyArr(int meshDimRelToMaxExt, const std::string& fam, bool renum=false) const;
 157     MEDLOADER_EXPORT virtual DataArrayInt *getNodeGroupArr(const std::string& grp, bool renum=false) const;
 158     MEDLOADER_EXPORT virtual DataArrayInt *getNodeGroupsArr(const std::vector<std::string>& grps, bool renum=false) const;
 159     MEDLOADER_EXPORT virtual DataArrayInt *getNodeFamilyArr(const std::string& fam, bool renum=false) const;
 160     MEDLOADER_EXPORT virtual DataArrayInt *getNodeFamiliesArr(const std::vector<std::string>& fams, bool renum=false) const;
 161     // tools
 162     MEDLOADER_EXPORT virtual bool unPolyze(std::vector<int>& oldCode, std::vector<int>& newCode, DataArrayInt *& o2nRenumCell) = 0;
 163   protected:
 164     MEDFileMesh();
 165     //! protected because no way in MED file API to specify this name
 166     void setUnivName(const std::string& name) { _univ_name=name; }
 167     void addFamilyOnAllGroupsHaving(const std::string& famName, const std::string& otherFamName);
 168     virtual void writeLL(med_idt fid) const = 0;
 169     void dealWithTinyInfo(const MEDCouplingMesh *m);
 170     virtual void synchronizeTinyInfoOnLeaves() const = 0;
 171     void getFamilyRepr(std::ostream& oss) const;
 172     virtual void appendFamilyEntries(const DataArrayInt *famIds, const std::vector< std::vector<int> >& fidsOfGrps, const std::vector<std::string>& grpNames);
 173     virtual void changeFamilyIdArr(int oldId, int newId) = 0;
 174     static void TranslateFamilyIds(int offset, DataArrayInt *famArr, std::vector< std::vector<int> >& famIdsPerGrp);
 175     static void ChangeAllGroupsContainingFamily(std::map<std::string, std::vector<std::string> >& groups, const std::string& familyNameToChange, const std::vector<std::string>& newFamiliesNames);
 176     static std::string FindOrCreateAndGiveFamilyWithId(std::map<std::string,int>& families, int id, bool& created);
 177     static std::string CreateNameNotIn(const std::string& nameTry, const std::vector<std::string>& namesToAvoid);
 178     static int PutInThirdComponentOfCodeOffset(std::vector<int>& code, int strt);
 179   protected:
 180     int _order;
 181     int _iteration;
 182     double _time;
 183     std::string _dt_unit;
 184     std::string _name;
 185     //! this attribute do not impact the state of instance -> mutable
 186     mutable std::string _univ_name;
 187     bool _univ_wr_status;
 188     std::string _desc_name;
 189   protected:
 190     std::map<std::string, std::vector<std::string> > _groups;
 191     std::map<std::string,int> _families;
 192   public:
 193     MEDLOADER_EXPORT static const char DFT_FAM_NAME[];
 194   };
 195
 196   class MEDFileUMesh : public MEDFileMesh
 197   {
 198     friend class MEDFileMesh;
 199   public:
 200     MEDLOADER_EXPORT static MEDFileUMesh *New(const std::string& fileName, const std::string& mName, int dt=-1, int it=-1, MEDFileMeshReadSelector *mrs=0);
 201     MEDLOADER_EXPORT static MEDFileUMesh *New(const std::string& fileName, MEDFileMeshReadSelector *mrs=0);
 202     MEDLOADER_EXPORT static MEDFileUMesh *New();
 203     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 204     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 205     MEDLOADER_EXPORT MEDFileMesh *createNewEmpty() const;
 206     MEDLOADER_EXPORT MEDFileMesh *deepCpy() const;
 207     MEDLOADER_EXPORT MEDFileMesh *shallowCpy() const;
 208     MEDLOADER_EXPORT bool isEqual(const MEDFileMesh *other, double eps, std::string& what) const;
 209     MEDLOADER_EXPORT void clearNonDiscrAttributes() const;
 210     MEDLOADER_EXPORT void setName(const std::string& name);
 211     MEDLOADER_EXPORT ~MEDFileUMesh();
 212     //
 213     MEDLOADER_EXPORT int getMaxAbsFamilyIdInArrays() const;
 214     MEDLOADER_EXPORT int getMaxFamilyIdInArrays() const;
 215     MEDLOADER_EXPORT int getMinFamilyIdInArrays() const;
 216     MEDLOADER_EXPORT int getMeshDimension() const;
 217     MEDLOADER_EXPORT int getSpaceDimension() const;
 218     MEDLOADER_EXPORT std::string simpleRepr() const;
 219     MEDLOADER_EXPORT std::string advancedRepr() const;
 220     MEDLOADER_EXPORT int getSizeAtLevel(int meshDimRelToMaxExt) const;
 221     MEDLOADER_EXPORT const DataArrayInt *getFamilyFieldAtLevel(int meshDimRelToMaxExt) const;
 222     MEDLOADER_EXPORT const DataArrayInt *getNumberFieldAtLevel(int meshDimRelToMaxExt) const;
 223     MEDLOADER_EXPORT const DataArrayInt *getRevNumberFieldAtLevel(int meshDimRelToMaxExt) const;
 224     MEDLOADER_EXPORT const DataArrayAsciiChar *getNameFieldAtLevel(int meshDimRelToMaxExt) const;
 225     MEDLOADER_EXPORT int getNumberOfNodes() const;
 226     MEDLOADER_EXPORT bool hasImplicitPart() const;
 227     MEDLOADER_EXPORT int buildImplicitPartIfAny(INTERP_KERNEL::NormalizedCellType gt) const;
 228     MEDLOADER_EXPORT void releaseImplicitPartIfAny() const;
 229     MEDLOADER_EXPORT std::vector<INTERP_KERNEL::NormalizedCellType> getGeoTypesAtLevel(int meshDimRelToMax) const;
 230     MEDLOADER_EXPORT void whichAreNodesFetched(const MEDFileField1TSStructItem& st, const MEDFileFieldGlobsReal *globs, std::vector<bool>& nodesFetched) const;
 231     MEDLOADER_EXPORT std::vector<int> getNonEmptyLevels() const;
 232     MEDLOADER_EXPORT std::vector<int> getNonEmptyLevelsExt() const;
 233     MEDLOADER_EXPORT std::vector<int> getFamArrNonEmptyLevelsExt() const;
 234     MEDLOADER_EXPORT std::vector<int> getNumArrNonEmptyLevelsExt() const;
 235     MEDLOADER_EXPORT std::vector<int> getNameArrNonEmptyLevelsExt() const;
 236     MEDLOADER_EXPORT std::vector<int> getGrpNonEmptyLevels(const std::string& grp) const;
 237     MEDLOADER_EXPORT std::vector<int> getGrpNonEmptyLevelsExt(const std::string& grp) const;
 238     MEDLOADER_EXPORT std::vector<int> getFamNonEmptyLevels(const std::string& fam) const;
 239     MEDLOADER_EXPORT std::vector<int> getFamNonEmptyLevelsExt(const std::string& fam) const;
 240     MEDLOADER_EXPORT std::vector<int> getGrpsNonEmptyLevels(const std::vector<std::string>& grps) const;
 241     MEDLOADER_EXPORT std::vector<int> getGrpsNonEmptyLevelsExt(const std::vector<std::string>& grps) const;
 242     MEDLOADER_EXPORT std::vector<int> getFamsNonEmptyLevels(const std::vector<std::string>& fams) const;
 243     MEDLOADER_EXPORT std::vector<int> getFamsNonEmptyLevelsExt(const std::vector<std::string>& fams) const;
 244     MEDLOADER_EXPORT std::vector<std::string> getGroupsOnSpecifiedLev(int meshDimRelToMaxExt) const;
 245     MEDLOADER_EXPORT DataArrayDouble *getCoords() const;
 246     MEDLOADER_EXPORT MEDCouplingUMesh *getGroup(int meshDimRelToMaxExt, const std::string& grp, bool renum=false) const;
 247     MEDLOADER_EXPORT MEDCouplingUMesh *getGroups(int meshDimRelToMaxExt, const std::vector<std::string>& grps, bool renum=false) const;
 248     MEDLOADER_EXPORT MEDCouplingUMesh *getFamily(int meshDimRelToMaxExt, const std::string& fam, bool renum=false) const;
 249     MEDLOADER_EXPORT MEDCouplingUMesh *getFamilies(int meshDimRelToMaxExt, const std::vector<std::string>& fams, bool renum=false) const;
 250     MEDLOADER_EXPORT DataArrayInt *getFamiliesArr(int meshDimRelToMaxExt, const std::vector<std::string>& fams, bool renum=false) const;
 251     MEDLOADER_EXPORT MEDCouplingUMesh *getMeshAtLevel(int meshDimRelToMaxExt, bool renum=false) const;
 252     MEDLOADER_EXPORT MEDCouplingMesh *getGenMeshAtLevel(int meshDimRelToMax, bool renum=false) const;
 253     MEDLOADER_EXPORT std::vector<int> getDistributionOfTypes(int meshDimRelToMax) const;
 254     MEDLOADER_EXPORT MEDCouplingUMesh *getLevel0Mesh(bool renum=false) const;
 255     MEDLOADER_EXPORT MEDCouplingUMesh *getLevelM1Mesh(bool renum=false) const;
 256     MEDLOADER_EXPORT MEDCouplingUMesh *getLevelM2Mesh(bool renum=false) const;
 257     MEDLOADER_EXPORT MEDCouplingUMesh *getLevelM3Mesh(bool renum=false) const;
 258     MEDLOADER_EXPORT void forceComputationOfParts() const;
 259     MEDLOADER_EXPORT std::vector<MEDCoupling1GTUMesh *> getDirectUndergroundSingleGeoTypeMeshes(int meshDimRelToMax) const;
 260     MEDLOADER_EXPORT MEDCoupling1GTUMesh *getDirectUndergroundSingleGeoTypeMesh(INTERP_KERNEL::NormalizedCellType gt) const;
 261     MEDLOADER_EXPORT DataArrayInt *extractFamilyFieldOnGeoType(INTERP_KERNEL::NormalizedCellType gt) const;
 262     MEDLOADER_EXPORT DataArrayInt *extractNumberFieldOnGeoType(INTERP_KERNEL::NormalizedCellType gt) const;
 263     //
 264     MEDLOADER_EXPORT void setFamilyNameAttachedOnId(int id, const std::string& newFamName);
 265     MEDLOADER_EXPORT void setCoords(DataArrayDouble *coords);
 266     MEDLOADER_EXPORT void eraseGroupsAtLevel(int meshDimRelToMaxExt);
 267     MEDLOADER_EXPORT void setFamilyFieldArr(int meshDimRelToMaxExt, DataArrayInt *famArr);
 268     MEDLOADER_EXPORT void setRenumFieldArr(int meshDimRelToMaxExt, DataArrayInt *renumArr);
 269     MEDLOADER_EXPORT void setNameFieldAtLevel(int meshDimRelToMaxExt, DataArrayAsciiChar *nameArr);
 270     MEDLOADER_EXPORT void addNodeGroup(const DataArrayInt *ids);
 271     MEDLOADER_EXPORT void addGroup(int meshDimRelToMaxExt, const DataArrayInt *ids);
 272     MEDLOADER_EXPORT void removeMeshAtLevel(int meshDimRelToMax);
 273     MEDLOADER_EXPORT void setMeshAtLevel(int meshDimRelToMax, MEDCoupling1GTUMesh *m);
 274     MEDLOADER_EXPORT void setMeshAtLevel(int meshDimRelToMax, MEDCouplingUMesh *m, bool newOrOld=false);
 275     MEDLOADER_EXPORT void setMeshes(const std::vector<const MEDCouplingUMesh *>& ms, bool renum=false);
 276     MEDLOADER_EXPORT void setGroupsFromScratch(int meshDimRelToMax, const std::vector<const MEDCouplingUMesh *>& ms, bool renum=false);
 277     MEDLOADER_EXPORT void setGroupsOnSetMesh(int meshDimRelToMax, const std::vector<const MEDCouplingUMesh *>& ms, bool renum=false);
 278     MEDLOADER_EXPORT void optimizeFamilies();
 279     // tools
 280     MEDLOADER_EXPORT void duplicateNodesOnM1Group(const std::string& grpNameM1, DataArrayInt *&nodesDuplicated, DataArrayInt *&cellsModified, DataArrayInt *&cellsNotModified);
 281     MEDLOADER_EXPORT bool unPolyze(std::vector<int>& oldCode, std::vector<int>& newCode, DataArrayInt *& o2nRenumCell);
 282     MEDLOADER_EXPORT DataArrayInt *zipCoords();
 283   private:
 284     void writeLL(med_idt fid) const;
 285     MEDFileUMesh();
 286     MEDFileUMesh(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 287     void loadUMeshFromFile(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 288     const MEDFileUMeshSplitL1 *getMeshAtLevSafe(int meshDimRelToMaxExt) const;
 289     MEDFileUMeshSplitL1 *getMeshAtLevSafe(int meshDimRelToMaxExt);
 290     void checkMeshDimCoherency(int meshDim, int meshDimRelToMax) const;
 291     DataArrayDouble *checkMultiMesh(const std::vector<const MEDCouplingUMesh *>& ms) const;
 292     void computeRevNum() const;
 293     void synchronizeTinyInfoOnLeaves() const;
 294     void changeFamilyIdArr(int oldId, int newId);
 295     std::list< MEDCouplingAutoRefCountObjectPtr<DataArrayInt> > getAllNonNullFamilyIds() const;
 296     void addGroupUnderground(bool isNodeGroup, const DataArrayInt *ids, DataArrayInt *famArr);
 297     MEDCouplingAutoRefCountObjectPtr<MEDFileUMeshSplitL1>& checkAndGiveEntryInSplitL1(int meshDimRelToMax, MEDCouplingPointSet *m);
 298   private:
 299     std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileUMeshSplitL1> > _ms;
 300     MEDCouplingAutoRefCountObjectPtr<DataArrayDouble> _coords;
 301     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _fam_coords;
 302     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _num_coords;
 303     MEDCouplingAutoRefCountObjectPtr<DataArrayAsciiChar> _name_coords;
 304     mutable MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _rev_num_coords;
 305   };
 306
 307   class MEDFileStructuredMesh : public MEDFileMesh
 308   {
 309     friend class MEDFileMesh;
 310   public:
 311     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 312     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 313     MEDLOADER_EXPORT int getMaxAbsFamilyIdInArrays() const;
 314     MEDLOADER_EXPORT int getMaxFamilyIdInArrays() const;
 315     MEDLOADER_EXPORT int getMinFamilyIdInArrays() const;
 316     MEDLOADER_EXPORT bool isEqual(const MEDFileMesh *other, double eps, std::string& what) const;
 317     MEDLOADER_EXPORT void clearNonDiscrAttributes() const;
 318     MEDLOADER_EXPORT DataArrayInt *getFamiliesArr(int meshDimRelToMaxExt, const std::vector<std::string>& fams, bool renum=false) const;
 319     MEDLOADER_EXPORT const DataArrayInt *getFamilyFieldAtLevel(int meshDimRelToMaxExt) const;
 320     MEDLOADER_EXPORT void setFamilyFieldArr(int meshDimRelToMaxExt, DataArrayInt *famArr);
 321     MEDLOADER_EXPORT void setRenumFieldArr(int meshDimRelToMaxExt, DataArrayInt *renumArr);
 322     MEDLOADER_EXPORT void setNameFieldAtLevel(int meshDimRelToMaxExt, DataArrayAsciiChar *nameArr);
 323     MEDLOADER_EXPORT const DataArrayInt *getNumberFieldAtLevel(int meshDimRelToMaxExt) const;
 324     MEDLOADER_EXPORT const DataArrayInt *getRevNumberFieldAtLevel(int meshDimRelToMaxExt) const;
 325     MEDLOADER_EXPORT const DataArrayAsciiChar *getNameFieldAtLevel(int meshDimRelToMaxExt) const;
 326     MEDLOADER_EXPORT std::vector<int> getNonEmptyLevels() const;
 327     MEDLOADER_EXPORT std::vector<int> getNonEmptyLevelsExt() const;
 328     MEDLOADER_EXPORT std::vector<int> getFamArrNonEmptyLevelsExt() const;
 329     MEDLOADER_EXPORT std::vector<int> getNumArrNonEmptyLevelsExt() const;
 330     MEDLOADER_EXPORT std::vector<int> getNameArrNonEmptyLevelsExt() const;
 331     MEDCouplingMesh *getGenMeshAtLevel(int meshDimRelToMax, bool renum=false) const;
 332     MEDLOADER_EXPORT int getSizeAtLevel(int meshDimRelToMaxExt) const;
 333     MEDLOADER_EXPORT int getNumberOfNodes() const;
 334     MEDLOADER_EXPORT bool hasImplicitPart() const;
 335     MEDLOADER_EXPORT int buildImplicitPartIfAny(INTERP_KERNEL::NormalizedCellType gt) const;
 336     MEDLOADER_EXPORT void releaseImplicitPartIfAny() const;
 337     MEDLOADER_EXPORT MEDCoupling1SGTUMesh *getImplicitFaceMesh() const;
 338     MEDLOADER_EXPORT std::vector<INTERP_KERNEL::NormalizedCellType> getGeoTypesAtLevel(int meshDimRelToMax) const;
 339     MEDLOADER_EXPORT void whichAreNodesFetched(const MEDFileField1TSStructItem& st, const MEDFileFieldGlobsReal *globs, std::vector<bool>& nodesFetched) const;
 340     MEDLOADER_EXPORT virtual const MEDCouplingStructuredMesh *getStructuredMesh() const = 0;
 341     // tools
 342     MEDLOADER_EXPORT bool unPolyze(std::vector<int>& oldCode, std::vector<int>& newCode, DataArrayInt *& o2nRenumCell);
 343   protected:
 344     void changeFamilyIdArr(int oldId, int newId);
 345     void deepCpyAttributes();
 346     void loadStrMeshFromFile(MEDFileStrMeshL2 *strm, med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 347     void writeStructuredLL(med_idt fid, const std::string& maa) const;
 348     void buildImplicitPart() const;
 349     void buildMinusOneImplicitPartIfNeeded() const;
 350     static med_geometry_type GetGeoTypeFromMeshDim(int meshDim);
 351   private:
 352     static void LoadStrMeshDAFromFile(med_idt fid, int meshDim, int dt, int it, const std::string& mName, MEDFileMeshReadSelector *mrs,
 353                                       MEDCouplingAutoRefCountObjectPtr<DataArrayInt>& famCells, MEDCouplingAutoRefCountObjectPtr<DataArrayInt>& numCells, MEDCouplingAutoRefCountObjectPtr<DataArrayAsciiChar>& namesCells);
 354   private:
 355     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _fam_nodes;
 356     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _num_nodes;
 357     MEDCouplingAutoRefCountObjectPtr<DataArrayAsciiChar> _names_nodes;
 358     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _fam_cells;
 359     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _num_cells;
 360     MEDCouplingAutoRefCountObjectPtr<DataArrayAsciiChar> _names_cells;
 361     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _fam_faces;
 362     MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _num_faces;
 363     MEDCouplingAutoRefCountObjectPtr<DataArrayAsciiChar> _names_faces;
 364     mutable MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _rev_num_nodes;
 365     mutable MEDCouplingAutoRefCountObjectPtr<DataArrayInt> _rev_num_cells;
 366     mutable MEDCouplingAutoRefCountObjectPtr<MEDCoupling1SGTUMesh> _faces_if_necessary;
 367   };
 368
 369   class MEDFileCMesh : public MEDFileStructuredMesh
 370   {
 371     friend class MEDFileMesh;
 372   public:
 373     MEDLOADER_EXPORT static MEDFileCMesh *New();
 374     MEDLOADER_EXPORT static MEDFileCMesh *New(const std::string& fileName, MEDFileMeshReadSelector *mrs=0);
 375     MEDLOADER_EXPORT static MEDFileCMesh *New(const std::string& fileName, const std::string& mName, int dt=-1, int it=-1, MEDFileMeshReadSelector *mrs=0);
 376     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 377     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 378     MEDLOADER_EXPORT MEDFileMesh *createNewEmpty() const;
 379     MEDLOADER_EXPORT MEDFileMesh *deepCpy() const;
 380     MEDLOADER_EXPORT MEDFileMesh *shallowCpy() const;
 381     MEDLOADER_EXPORT bool isEqual(const MEDFileMesh *other, double eps, std::string& what) const;
 382     MEDLOADER_EXPORT int getMeshDimension() const;
 383     MEDLOADER_EXPORT int getSpaceDimension() const;
 384     MEDLOADER_EXPORT std::string simpleRepr() const;
 385     MEDLOADER_EXPORT std::string advancedRepr() const;
 386     MEDLOADER_EXPORT void clearNonDiscrAttributes() const;
 387     MEDLOADER_EXPORT const MEDCouplingCMesh *getMesh() const;
 388     MEDLOADER_EXPORT void setMesh(MEDCouplingCMesh *m);
 389   private:
 390     const MEDCouplingStructuredMesh *getStructuredMesh() const;
 391     void writeLL(med_idt fid) const;
 392     MEDFileCMesh();
 393     void synchronizeTinyInfoOnLeaves() const;
 394     MEDFileCMesh(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 395     void loadCMeshFromFile(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 396   private:
 397     MEDCouplingAutoRefCountObjectPtr<MEDCouplingCMesh> _cmesh;
 398   };
 399
 400   class MEDFileCurveLinearMesh : public MEDFileStructuredMesh
 401   {
 402     friend class MEDFileMesh;
 403   public:
 404     MEDLOADER_EXPORT static MEDFileCurveLinearMesh *New();
 405     MEDLOADER_EXPORT static MEDFileCurveLinearMesh *New(const std::string& fileName, MEDFileMeshReadSelector *mrs=0);
 406     MEDLOADER_EXPORT static MEDFileCurveLinearMesh *New(const std::string& fileName, const std::string& mName, int dt=-1, int it=-1, MEDFileMeshReadSelector *mrs=0);
 407     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 408     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 409     MEDLOADER_EXPORT MEDFileMesh *createNewEmpty() const;
 410     MEDLOADER_EXPORT MEDFileMesh *deepCpy() const;
 411     MEDLOADER_EXPORT MEDFileMesh *shallowCpy() const;
 412     MEDLOADER_EXPORT bool isEqual(const MEDFileMesh *other, double eps, std::string& what) const;
 413     MEDLOADER_EXPORT int getMeshDimension() const;
 414     MEDLOADER_EXPORT std::string simpleRepr() const;
 415     MEDLOADER_EXPORT std::string advancedRepr() const;
 416     MEDLOADER_EXPORT void clearNonDiscrAttributes() const;
 417     MEDLOADER_EXPORT const MEDCouplingCurveLinearMesh *getMesh() const;
 418     MEDLOADER_EXPORT void setMesh(MEDCouplingCurveLinearMesh *m);
 419   private:
 420     MEDFileCurveLinearMesh();
 421     MEDFileCurveLinearMesh(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);
 422     const MEDCouplingStructuredMesh *getStructuredMesh() const;
 423     void synchronizeTinyInfoOnLeaves() const;
 424     void writeLL(med_idt fid) const;
 425     void loadCLMeshFromFile(med_idt fid, const std::string& mName, int dt, int it, MEDFileMeshReadSelector *mrs);//to imp
 426   private:
 427     MEDCouplingAutoRefCountObjectPtr<MEDCouplingCurveLinearMesh> _clmesh;
 428   };
 429
 430   class MEDFileMeshMultiTS : public RefCountObject, public MEDFileWritable
 431   {
 432   public:
 433     MEDLOADER_EXPORT static MEDFileMeshMultiTS *New();
 434     MEDLOADER_EXPORT static MEDFileMeshMultiTS *New(const std::string& fileName);
 435     MEDLOADER_EXPORT static MEDFileMeshMultiTS *New(const std::string& fileName, const std::string& mName);
 436     MEDLOADER_EXPORT MEDFileMeshMultiTS *deepCpy() const;
 437     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 438     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 439     MEDLOADER_EXPORT std::string getName() const;
 440     MEDLOADER_EXPORT void setName(const std::string& newMeshName);
 441     MEDLOADER_EXPORT bool changeNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 442     MEDLOADER_EXPORT MEDFileMesh *getOneTimeStep() const;
 443     MEDLOADER_EXPORT void write(med_idt fid) const;
 444     MEDLOADER_EXPORT void write(const std::string& fileName, int mode) const;
 445     MEDLOADER_EXPORT void setOneTimeStep(MEDFileMesh *mesh1TimeStep);
 446   private:
 447     void loadFromFile(const std::string& fileName, const std::string& mName);
 448     MEDFileMeshMultiTS();
 449     MEDFileMeshMultiTS(const std::string& fileName);
 450     MEDFileMeshMultiTS(const std::string& fileName, const std::string& mName);
 451   private:
 452     std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileMesh> > _mesh_one_ts;
 453   };
 454
 455   class MEDFileMeshesIterator;
 456
 457   class MEDFileMeshes : public RefCountObject, public MEDFileWritable
 458   {
 459   public:
 460     MEDLOADER_EXPORT static MEDFileMeshes *New();
 461     MEDLOADER_EXPORT static MEDFileMeshes *New(const std::string& fileName);
 462     MEDLOADER_EXPORT MEDFileMeshes *deepCpy() const;
 463     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 464     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 465     MEDLOADER_EXPORT std::string simpleRepr() const;
 466     MEDLOADER_EXPORT void simpleReprWithoutHeader(std::ostream& oss) const;
 467     MEDLOADER_EXPORT void write(const std::string& fileName, int mode) const;
 468     MEDLOADER_EXPORT void write(med_idt fid) const;
 469     MEDLOADER_EXPORT int getNumberOfMeshes() const;
 470     MEDLOADER_EXPORT MEDFileMeshesIterator *iterator();
 471     MEDLOADER_EXPORT MEDFileMesh *getMeshAtPos(int i) const;
 472     MEDLOADER_EXPORT MEDFileMesh *getMeshWithName(const std::string& mname) const;
 473     MEDLOADER_EXPORT std::vector<std::string> getMeshesNames() const;
 474     MEDLOADER_EXPORT bool changeNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 475     //
 476     MEDLOADER_EXPORT void resize(int newSize);
 477     MEDLOADER_EXPORT void pushMesh(MEDFileMesh *mesh);
 478     MEDLOADER_EXPORT void setMeshAtPos(int i, MEDFileMesh *mesh);
 479     MEDLOADER_EXPORT void destroyMeshAtPos(int i);
 480   private:
 481     void checkCoherency() const;
 482     void loadFromFile(const std::string& fileName);
 483     MEDFileMeshes();
 484     MEDFileMeshes(const std::string& fileName);
 485   private:
 486     std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileMeshMultiTS> > _meshes;
 487   };
 488
 489   class MEDFileMeshesIterator
 490   {
 491   public:
 492     MEDLOADER_EXPORT MEDFileMeshesIterator(MEDFileMeshes *ms);
 493     MEDLOADER_EXPORT ~MEDFileMeshesIterator();
 494     MEDLOADER_EXPORT MEDFileMesh *nextt();
 495   private:
 496     MEDCouplingAutoRefCountObjectPtr<MEDFileMeshes> _ms;
 497     int _iter_id;
 498     int _nb_iter;
 499   };
 500 }
 501
 502 #endif