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