src/MEDLoader/MEDFileField.hxx

   1 // Copyright (C) 2007-2016  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 __MEDFILEFIELD_HXX__
  22 #define __MEDFILEFIELD_HXX__
  23
  24 #include "MEDLoaderDefines.hxx"
  25
  26 #include "MEDFileFieldOverView.hxx"
  27 #include "MEDFileUtilities.txx"
  28 #include "MEDFileEntities.hxx"
  29
  30 #include "MCAuto.hxx"
  31 #include "MEDLoaderTraits.hxx"
  32 #include "MEDCouplingTraits.hxx"
  33 #include "MEDCouplingRefCountObject.hxx"
  34 #include "MEDCouplingFieldInt.hxx"
  35 #include "MEDCouplingMemArray.hxx"
  36
  37 #include "NormalizedUnstructuredMesh.hxx"
  38 #include "InterpKernelException.hxx"
  39
  40 #include <vector>
  41 #include <string>
  42 #include <list>
  43 #include <set>
  44
  45 #include "med.h"
  46
  47 namespace MEDCoupling
  48 {
  49   class MEDFileFieldGlobs;
  50   class MEDCouplingMesh;
  51   class MEDCouplingFieldDouble;
  52   class MEDFileMesh;
  53
  54   class MEDFileFieldLoc : public RefCountObject
  55   {
  56   public:
  57     MEDLOADER_EXPORT void simpleRepr(std::ostream& oss) const;
  58     MEDLOADER_EXPORT std::string getName() const { return _name; }
  59     MEDLOADER_EXPORT void setName(const std::string& name);
  60     static MEDFileFieldLoc *New(med_idt fid, const std::string& locName);
  61     static MEDFileFieldLoc *New(med_idt fid, int id);
  62     static MEDFileFieldLoc *New(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
  63     std::size_t getHeapMemorySizeWithoutChildren() const;
  64     std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
  65     MEDFileFieldLoc *deepCopy() const;
  66     MEDLOADER_EXPORT int getNbOfGaussPtPerCell() const { return _nb_gauss_pt; }
  67     MEDLOADER_EXPORT void writeLL(med_idt fid) const;
  68     MEDLOADER_EXPORT std::string repr() const;
  69     MEDLOADER_EXPORT bool isName(const std::string& name) const { return _name==name; }
  70     MEDLOADER_EXPORT int getDimension() const { return _dim; }
  71     MEDLOADER_EXPORT int getNumberOfGaussPoints() const { return _nb_gauss_pt; }
  72     MEDLOADER_EXPORT int getNumberOfPointsInCells() const { return _nb_node_per_cell; }
  73     MEDLOADER_EXPORT const std::vector<double>& getRefCoords() const { return _ref_coo; }
  74     MEDLOADER_EXPORT const std::vector<double>& getGaussCoords() const { return _gs_coo; }
  75     MEDLOADER_EXPORT const std::vector<double>& getGaussWeights() const { return _w; }
  76     MEDLOADER_EXPORT INTERP_KERNEL::NormalizedCellType getGeoType() const { return _geo_type; }
  77     MEDLOADER_EXPORT bool isEqual(const MEDFileFieldLoc& other, double eps) const;
  78   private:
  79     MEDFileFieldLoc(med_idt fid, const std::string& locName);
  80     MEDFileFieldLoc(med_idt fid, int id);
  81     MEDFileFieldLoc(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
  82   private:
  83     int _dim;
  84     int _nb_gauss_pt;
  85     int _nb_node_per_cell;
  86     std::string _name;
  87     INTERP_KERNEL::NormalizedCellType _geo_type;
  88     std::vector<double> _ref_coo;
  89     std::vector<double> _gs_coo;
  90     std::vector<double> _w;
  91   };
  92
  93   /// @cond INTERNAL
  94   class MEDFileAnyTypeField1TSWithoutSDA;
  95   class MEDFileFieldPerMeshPerTypeCommon;
  96   class MEDFileFieldPerMeshPerType;
  97   class MEDFileField1TSWithoutSDA;
  98   class MEDFileFieldNameScope;
  99   class MEDFileFieldGlobsReal;
 100   class MEDFileFieldPerMesh;
 101   class PartDefinition;
 102
 103   class MEDFileFieldPerMeshPerTypePerDisc : public RefCountObject, public MEDFileWritable
 104   {
 105   public:
 106     static MEDFileFieldPerMeshPerTypePerDisc *NewOnRead(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
 107     static MEDFileFieldPerMeshPerTypePerDisc *New(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int locId);
 108     static MEDFileFieldPerMeshPerTypePerDisc *New(const MEDFileFieldPerMeshPerTypePerDisc& other);
 109     std::size_t getHeapMemorySizeWithoutChildren() const;
 110     std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 111     MEDFileFieldPerMeshPerTypePerDisc *deepCopy(MEDFileFieldPerMeshPerTypeCommon *father) const;
 112     void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldDouble *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 113     void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldDouble *field, const DataArray *arrr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 114     void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldDouble *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob);
 115     void getCoarseData(TypeOfField& type, std::pair<int,int>& dad, std::string& pfl, std::string& loc) const;
 116     void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
 117     const MEDFileFieldPerMeshPerTypeCommon *getFather() const;
 118     void loadOnlyStructureOfDataRecursively(med_idt fid, int& start, const MEDFileFieldNameScope& nasc);
 119     void loadBigArray(med_idt fid, const MEDFileFieldNameScope& nasc);
 120     void setNewStart(int newValueOfStart);
 121     int getIteration() const;
 122     int getOrder() const;
 123     double getTime() const;
 124     std::string getMeshName() const;
 125     TypeOfField getType() const;
 126     void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
 127     void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
 128     void setType(TypeOfField newType);
 129     INTERP_KERNEL::NormalizedCellType getGeoType() const;
 130     int getNumberOfComponents() const;
 131     int getNumberOfTuples() const;
 132     int getStart() const { return _start; }
 133     int getEnd() const { return _end; }
 134     int getNumberOfVals() const { return _nval; }
 135     DataArray *getOrCreateAndGetArray();
 136     const DataArray *getOrCreateAndGetArray() const;
 137     const std::vector<std::string>& getInfo() const;
 138     std::string getProfile() const;
 139     void setProfile(const std::string& newPflName);
 140     std::string getLocalization() const;
 141     void setLocalization(const std::string& newLocName);
 142     int getLocId() const { return _loc_id; }
 143     void setLocId(int newId) const { _loc_id=newId; }
 144     void setFather(MEDFileFieldPerMeshPerTypeCommon *newFather) { _father=newFather; }
 145     void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 146     void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 147     void getFieldAtLevel(TypeOfField type, const MEDFileFieldGlobsReal *glob, std::vector< std::pair<int,int> >& dads, std::vector<const DataArrayInt *>& pfls, std::vector<int>& locs,
 148                          std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes) const;
 149     void fillValues(int discId, int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 150     int fillEltIdsFromCode(int offset, const std::vector<int>& codeOfMesh, const MEDFileFieldGlobsReal& glob, int *ptToFill) const;
 151     int fillTupleIds(int *ptToFill) const;
 152     static int ConvertType(TypeOfField type, int locId);
 153     static std::vector< std::vector< const MEDFileFieldPerMeshPerTypePerDisc *> > SplitPerDiscretization(const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entries);
 154     static bool RenumberChunks(int offset, const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
 155                                const DataArrayInt *explicitIdsInMesh, const std::vector<int>& newCode,
 156                                MEDFileFieldGlobsReal& glob, DataArrayDouble *arr, std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> >& result);
 157     static MEDFileFieldPerMeshPerTypePerDisc *NewObjectOnSameDiscThanPool(TypeOfField typeF, INTERP_KERNEL::NormalizedCellType geoType, DataArrayInt *idsOfMeshElt,
 158                                                                           bool isPfl, int nbi, int offset, std::list< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
 159                                                                           MEDFileFieldGlobsReal& glob, bool &notInExisting);
 160     static MCAuto<MEDFileFieldPerMeshPerTypePerDisc> Aggregate(int &start, const std::vector<std::pair<int,const MEDFileFieldPerMeshPerTypePerDisc *> >& pms, const std::vector< std::vector< std::pair<int,int> > >& dts, TypeOfField tof, MEDFileFieldPerMeshPerType *father, std::vector<std::pair< int, std::pair<int,int> > >& extractInfo);
 161     MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type):_type(type),_father(fath),_start(-1),_end(-1),_nval(-1),_loc_id(-5),_profile_it(-1) { }
 162   private:
 163     MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
 164     MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const std::string& dummy);
 165     MEDFileFieldPerMeshPerTypePerDisc(const MEDFileFieldPerMeshPerTypePerDisc& other);
 166     MEDFileFieldPerMeshPerTypePerDisc();
 167   private:
 168     void goReadZeValuesInFile(med_idt fid, const std::string& fieldName, int nbOfCompo, int iteration, int order, med_entity_type menti, med_geometry_type mgeoti, unsigned char *startFeedingPtr);
 169   private:
 170     TypeOfField _type;
 171     MEDFileFieldPerMeshPerTypeCommon *_father;
 172     int _start;
 173     int _end;
 174     //! _nval is different than end-start in case of ON_GAUSS_PT and ON_GAUSS_NE ! (_nval=(_end-_start)/nbi)
 175     int _nval;
 176     std::string _profile;
 177     std::string _localization;
 178     //! only on assignement -3 : ON_NODES, -2 : ON_CELLS, -1 : ON_GAUSS_NE, 0..* : ON_GAUSS_PT
 179     mutable int _loc_id;
 180     mutable int _profile_it;
 181     MCAuto<PartDefinition> _pd;
 182   public:
 183     mutable int _tmp_work1;
 184   };
 185
 186   class MEDFileFieldPerMeshPerTypeCommon : public RefCountObject, public MEDFileWritable
 187   {
 188   public:
 189     std::size_t getHeapMemorySizeWithoutChildren() const;
 190     std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 191     void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 192     void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldDouble *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 193     void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
 194     void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 195     const MEDFileFieldPerMesh *getFather() const;
 196     void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
 197     void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
 198     void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
 199     bool isUniqueLevel(int& dim) const;
 200     void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
 201     void fillFieldSplitedByType(std::vector< std::pair<int,int> >& dads, std::vector<TypeOfField>& types, std::vector<std::string>& pfls, std::vector<std::string>& locs) const;
 202     int getIteration() const;
 203     int getOrder() const;
 204     double getTime() const;
 205     std::string getMeshName() const;
 206     void getSizes(int& globalSz, int& nbOfEntries) const;
 207     int getNumberOfComponents() const;
 208     bool presenceOfMultiDiscPerGeoType() const;
 209     DataArray *getOrCreateAndGetArray();
 210     const DataArray *getOrCreateAndGetArray() const;
 211     const std::vector<std::string>& getInfo() const;
 212     std::vector<std::string> getPflsReallyUsed() const;
 213     std::vector<std::string> getLocsReallyUsed() const;
 214     std::vector<std::string> getPflsReallyUsedMulti() const;
 215     std::vector<std::string> getLocsReallyUsedMulti() const;
 216     void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 217     void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 218     MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId);
 219     const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId) const;
 220     int getNumberOfLoc() const { return _field_pm_pt_pd.size(); }
 221     void fillValues(int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 222     void setLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
 223     bool keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
 224     bool keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
 225     static med_entity_type ConvertIntoMEDFileType(TypeOfField ikType, INTERP_KERNEL::NormalizedCellType ikGeoType, med_geometry_type& medfGeoType);
 226     MEDFileFieldPerMeshPerTypeCommon(MEDFileFieldPerMesh *father):_father(father) { }
 227   public:
 228     virtual ~MEDFileFieldPerMeshPerTypeCommon();
 229     virtual void getDimension(int& dim) const = 0;
 230     virtual INTERP_KERNEL::NormalizedCellType getGeoType() const = 0;
 231     virtual void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const = 0;
 232     virtual void simpleRepr(int bkOffset, std::ostream& oss, int id) const = 0;
 233     virtual std::string getGeoTypeRepr() const = 0;
 234     virtual MEDFileFieldPerMeshPerTypeCommon *deepCopy(MEDFileFieldPerMesh *father) const = 0;
 235     virtual void getFieldAtLevel(int meshDim, TypeOfField type, const MEDFileFieldGlobsReal *glob, std::vector< std::pair<int,int> >& dads, std::vector<const DataArrayInt *>& pfls, std::vector<int>& locs, std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes) const = 0;
 236   protected:
 237     void setFather(MEDFileFieldPerMesh *father);
 238     void deepCopyElements();
 239     std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldDouble *field, int offset, int nbOfCells);
 240     std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldDouble *field, int offset, int nbOfCells);
 241     std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldDouble *field, const DataArrayInt *subCells);
 242     std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldDouble *field, const DataArrayInt *subCells);
 243   private:
 244     MEDFileFieldPerMesh *_father;
 245   protected:
 246     std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> > _field_pm_pt_pd;
 247   };
 248
 249   class MEDFileFieldPerMeshPerType : public MEDFileFieldPerMeshPerTypeCommon
 250   {
 251   public:
 252     static MEDFileFieldPerMeshPerType *New(MEDFileFieldPerMesh *fath, INTERP_KERNEL::NormalizedCellType geoType);
 253     static MEDFileFieldPerMeshPerType *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
 254     static MCAuto<MEDFileFieldPerMeshPerType> Aggregate(int &start, const std::vector< std::pair<int,const MEDFileFieldPerMeshPerType *> >& pms, const std::vector< std::vector< std::pair<int,int> > >& dts, INTERP_KERNEL::NormalizedCellType gt, MEDFileFieldPerMesh *father, std::vector<std::pair< int, std::pair<int,int> > >& extractInfo);
 255   public:// overload of abstract methods
 256     void getDimension(int& dim) const;
 257     INTERP_KERNEL::NormalizedCellType getGeoType() const;
 258     void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
 259     void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
 260     std::string getGeoTypeRepr() const;
 261     MEDFileFieldPerMeshPerType *deepCopy(MEDFileFieldPerMesh *father) const;
 262     void getFieldAtLevel(int meshDim, TypeOfField type, const MEDFileFieldGlobsReal *glob, std::vector< std::pair<int,int> >& dads, std::vector<const DataArrayInt *>& pfls, std::vector<int>& locs, std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes) const;
 263   private:
 264     MEDFileFieldPerMeshPerType(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
 265     MEDFileFieldPerMeshPerType(MEDFileFieldPerMesh *father, INTERP_KERNEL::NormalizedCellType gt);
 266   private:
 267     INTERP_KERNEL::NormalizedCellType _geo_type;
 268   };
 269
 270   class MEDFileFieldPerMeshPerTypeDyn : public MEDFileFieldPerMeshPerTypeCommon
 271   {
 272   public:
 273     static MEDFileFieldPerMeshPerTypeDyn *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileEntities *entities, int idGT, const MEDFileFieldNameScope& nasc);
 274     int getDynGT() const;
 275   public:
 276     void getDimension(int& dim) const;
 277     INTERP_KERNEL::NormalizedCellType getGeoType() const;
 278     void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
 279     void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
 280     std::string getGeoTypeRepr() const;
 281     MEDFileFieldPerMeshPerTypeDyn *deepCopy(MEDFileFieldPerMesh *father) const;
 282     void getFieldAtLevel(int meshDim, TypeOfField type, const MEDFileFieldGlobsReal *glob, std::vector< std::pair<int,int> >& dads, std::vector<const DataArrayInt *>& pfls, std::vector<int>& locs, std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes) const;
 283   private:
 284     MEDFileFieldPerMeshPerTypeDyn(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileStructureElement *se, const MEDFileFieldNameScope& nasc);
 285   private:
 286     MCConstAuto<MEDFileStructureElement> _se;
 287   };
 288
 289   class MEDFileMesh;
 290
 291   class MEDFileFieldPerMesh : public RefCountObject, public MEDFileWritable
 292   {
 293   public:
 294     static MEDFileFieldPerMesh *New(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
 295     static MEDFileFieldPerMesh *NewOnRead(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
 296     std::size_t getHeapMemorySizeWithoutChildren() const;
 297     std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 298     MEDFileFieldPerMesh *deepCopy(MEDFileAnyTypeField1TSWithoutSDA *father) const;
 299     void simpleRepr(int bkOffset,std::ostream& oss, int id) const;
 300     void copyTinyInfoFrom(const MEDCouplingMesh *mesh);
 301     void assignFieldProfile(int& start, const DataArrayInt *multiTypePfl, const std::vector<int>& code, const std::vector<int>& code2, const std::vector<DataArrayInt *>& idsInPflPerType, const std::vector<DataArrayInt *>& idsPerType, const MEDCouplingFieldDouble *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 302     void assignFieldNoProfileNoRenum(int& start, const std::vector<int>& code, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 303     void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
 304     void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 305     void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
 306     void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
 307     void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
 308     void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
 309     std::vector< std::vector< std::pair<int,int> > > getFieldSplitedByType(std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 310     void getDimension(int& dim) const;
 311     bool isUniqueLevel(int& dim) const;
 312     double getTime() const;
 313     int getIteration() const;
 314     int getOrder() const;
 315     int getMeshIteration() const { return _mesh_iteration; }
 316     int getMeshOrder() const { return _mesh_order; }
 317     std::string getMeshName() const { return _mesh_name; }
 318     int getNumberOfComponents() const;
 319     bool presenceOfMultiDiscPerGeoType() const;
 320     bool presenceOfStructureElements() const;
 321     bool onlyStructureElements() const;
 322     void killStructureElements();
 323     void keepOnlyStructureElements();
 324     DataArray *getOrCreateAndGetArray();
 325     const DataArray *getOrCreateAndGetArray() const;
 326     const std::vector<std::string>& getInfo() const;
 327     std::vector<std::string> getPflsReallyUsed() const;
 328     std::vector<std::string> getLocsReallyUsed() const;
 329     std::vector<std::string> getPflsReallyUsedMulti() const;
 330     std::vector<std::string> getLocsReallyUsedMulti() const;
 331     bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 332     bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N, MEDFileFieldGlobsReal& glob);
 333     void keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
 334     void keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
 335     void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 336     void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 337     MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, const MEDFileFieldGlobsReal *glob, const MEDCouplingMesh *mesh, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 338     DataArray *getFieldOnMeshAtLevelWithPfl(TypeOfField type, const MEDCouplingMesh *mesh, DataArrayInt *&pfl, const MEDFileFieldGlobsReal *glob, const MEDFileFieldNameScope& nasc) const;
 339     void getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 340     MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId);
 341     const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId) const;
 342     static MCAuto<MEDFileFieldPerMesh> Aggregate(int &start, const std::vector<const MEDFileFieldPerMesh *>& pms, const std::vector< std::vector< std::pair<int,int> > >& dts, MEDFileAnyTypeField1TSWithoutSDA *father, std::vector<std::pair< int, std::pair<int,int> > >& extractInfo);
 343   private:
 344     int addNewEntryIfNecessary(INTERP_KERNEL::NormalizedCellType type);
 345     MEDCouplingFieldDouble *finishField(TypeOfField type, const MEDFileFieldGlobsReal *glob,
 346                                         const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs, const MEDCouplingMesh *mesh, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 347     MEDCouplingFieldDouble *finishField2(TypeOfField type, const MEDFileFieldGlobsReal *glob,
 348                                          const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
 349                                          const std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes,
 350                                          const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 351     MEDCouplingFieldDouble *finishFieldNode2(const MEDFileFieldGlobsReal *glob,
 352                                              const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
 353                                              const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 354     DataArray *finishField4(const std::vector< std::pair<int,int> >& dads, const DataArrayInt *pflIn, int nbOfElems, DataArrayInt *&pflOut) const;
 355     void assignNewLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
 356     static void SortArraysPerType(const MEDFileFieldGlobsReal *glob, TypeOfField type,
 357                                   const std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes, const std::vector< std::pair<int,int> >& dads, const std::vector<const DataArrayInt *>& pfls, const std::vector<int>& locs,
 358                                   std::vector<int>& code, std::vector<DataArrayInt *>& notNullPfls);
 359     static int ComputeNbOfElems(const MEDFileFieldGlobsReal *glob, TypeOfField type, const std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes, const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs);
 360     MEDFileFieldPerMesh(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
 361     MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
 362     MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const std::string& meshName, int meshIt, int meshOrd):_father(fath),_mesh_name(meshName),_mesh_iteration(meshIt),_mesh_order(meshOrd) { }
 363   private:
 364     std::string _mesh_name;
 365     int _mesh_iteration;
 366     int _mesh_order;
 367     MEDFileAnyTypeField1TSWithoutSDA *_father;
 368     std::vector< MCAuto< MEDFileFieldPerMeshPerTypeCommon > > _field_pm_pt;
 369   };
 370
 371   class MEDFileFieldGlobsReal;
 372
 373   class MEDFileFieldGlobs : public RefCountObject
 374   {
 375   public:
 376     static MEDFileFieldGlobs *New(med_idt fid);
 377     static MEDFileFieldGlobs *New();
 378     std::size_t getHeapMemorySizeWithoutChildren() const;
 379     std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 380     MEDFileFieldGlobs *deepCopy() const;
 381     MEDFileFieldGlobs *shallowCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
 382     MEDFileFieldGlobs *deepCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
 383     void simpleRepr(std::ostream& oss) const;
 384     void appendGlobs(const MEDFileFieldGlobs& other, double eps);
 385     void checkGlobsPflsPartCoherency(const std::vector<std::string>& pflsUsed) const;
 386     void checkGlobsLocsPartCoherency(const std::vector<std::string>& locsUsed) const;
 387     void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
 388     void loadProfileInFile(med_idt fid, int id);
 389     void loadGlobals(med_idt fid, const MEDFileFieldGlobsReal& real);
 390     void loadAllGlobals(med_idt fid);
 391     void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
 392     std::vector<std::string> getPfls() const;
 393     std::vector<std::string> getLocs() const;
 394     bool existsPfl(const std::string& pflName) const;
 395     bool existsLoc(const std::string& locName) const;
 396     std::string createNewNameOfPfl() const;
 397     std::string createNewNameOfLoc() const;
 398     std::vector< std::vector<int> > whichAreEqualProfiles() const;
 399     std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
 400     void setFileName(const std::string& fileName) { _file_name=fileName; }
 401     void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 402     void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 403     int getNbOfGaussPtPerCell(int locId) const;
 404     int getLocalizationId(const std::string& loc) const;
 405     std::string getFileName() const { return _file_name; }
 406     const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
 407     const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
 408     const DataArrayInt *getProfileFromId(int pflId) const;
 409     const DataArrayInt *getProfile(const std::string& pflName) const;
 410     MEDFileFieldLoc& getLocalizationFromId(int locId);
 411     MEDFileFieldLoc& getLocalization(const std::string& locName);
 412     DataArrayInt *getProfile(const std::string& pflName);
 413     DataArrayInt *getProfileFromId(int pflId);
 414     void killProfileIds(const std::vector<int>& pflIds);
 415     void killLocalizationIds(const std::vector<int>& locIds);
 416     //
 417     void appendProfile(DataArrayInt *pfl);
 418     void appendLoc(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
 419     //
 420     static std::string CreateNewNameNotIn(const std::string& prefix, const std::vector<std::string>& namesToAvoid);
 421   protected:
 422     MEDFileFieldGlobs(med_idt fid);
 423     MEDFileFieldGlobs();
 424     ~MEDFileFieldGlobs();
 425   protected:
 426     std::vector< MCAuto<DataArrayInt> > _pfls;
 427     std::vector< MCAuto<MEDFileFieldLoc> > _locs;
 428     std::string _file_name;
 429   };
 430
 431   /// @endcond INTERNAL
 432
 433   class MEDFileFieldGlobsReal
 434   {
 435   public:
 436     MEDLOADER_EXPORT MEDFileFieldGlobsReal(med_idt fid);
 437     MEDLOADER_EXPORT MEDFileFieldGlobsReal();
 438     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 439     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 440     MEDLOADER_EXPORT void simpleReprGlobs(std::ostream& oss) const;
 441     MEDLOADER_EXPORT void resetContent();
 442     MEDLOADER_EXPORT void shallowCpyGlobs(const MEDFileFieldGlobsReal& other);
 443     MEDLOADER_EXPORT void deepCpyGlobs(const MEDFileFieldGlobsReal& other);
 444     MEDLOADER_EXPORT void shallowCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
 445     MEDLOADER_EXPORT void deepCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
 446     MEDLOADER_EXPORT void appendGlobs(const MEDFileFieldGlobsReal& other, double eps);
 447     MEDLOADER_EXPORT void checkGlobsCoherency() const;
 448     MEDLOADER_EXPORT void checkGlobsPflsPartCoherency() const;
 449     MEDLOADER_EXPORT void checkGlobsLocsPartCoherency() const;
 450     MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsed() const = 0;
 451     MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsed() const = 0;
 452     MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsedMulti() const = 0;
 453     MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsedMulti() const = 0;
 454     MEDLOADER_EXPORT virtual void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
 455     MEDLOADER_EXPORT virtual void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
 456     MEDLOADER_EXPORT virtual ~MEDFileFieldGlobsReal();
 457     //
 458     MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
 459     MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id);
 460     MEDLOADER_EXPORT void loadGlobals(med_idt fid);
 461     MEDLOADER_EXPORT void loadAllGlobals(med_idt fid);
 462     MEDLOADER_EXPORT void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
 463     MEDLOADER_EXPORT std::vector<std::string> getPfls() const;
 464     MEDLOADER_EXPORT std::vector<std::string> getLocs() const;
 465     MEDLOADER_EXPORT bool existsPfl(const std::string& pflName) const;
 466     MEDLOADER_EXPORT bool existsLoc(const std::string& locName) const;
 467     MEDLOADER_EXPORT std::string createNewNameOfPfl() const;
 468     MEDLOADER_EXPORT std::string createNewNameOfLoc() const;
 469     MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualProfiles() const;
 470     MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
 471     MEDLOADER_EXPORT void setFileName(const std::string& fileName);
 472     MEDLOADER_EXPORT void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 473     MEDLOADER_EXPORT void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 474     MEDLOADER_EXPORT void changePflsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 475     MEDLOADER_EXPORT void changeLocsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 476     MEDLOADER_EXPORT void changePflName(const std::string& oldName, const std::string& newName);
 477     MEDLOADER_EXPORT void changeLocName(const std::string& oldName, const std::string& newName);
 478     MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipPflsNames();
 479     MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipLocsNames(double eps);
 480     MEDLOADER_EXPORT int getNbOfGaussPtPerCell(int locId) const;
 481     MEDLOADER_EXPORT int getLocalizationId(const std::string& loc) const;
 482     MEDLOADER_EXPORT std::string getFileName() const;
 483     MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
 484     MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
 485     MEDLOADER_EXPORT MEDFileFieldLoc& getLocalizationFromId(int locId);
 486     MEDLOADER_EXPORT MEDFileFieldLoc& getLocalization(const std::string& locName);
 487     MEDLOADER_EXPORT const DataArrayInt *getProfile(const std::string& pflName) const;
 488     MEDLOADER_EXPORT const DataArrayInt *getProfileFromId(int pflId) const;
 489     MEDLOADER_EXPORT DataArrayInt *getProfile(const std::string& pflName);
 490     MEDLOADER_EXPORT DataArrayInt *getProfileFromId(int pflId);
 491     MEDLOADER_EXPORT void killProfileIds(const std::vector<int>& pflIds);
 492     MEDLOADER_EXPORT void killLocalizationIds(const std::vector<int>& locIds);
 493     //
 494     MEDLOADER_EXPORT void appendProfile(DataArrayInt *pfl);
 495     MEDLOADER_EXPORT void appendLoc(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
 496   protected:
 497     MEDFileFieldGlobs *contentNotNull();
 498     const MEDFileFieldGlobs *contentNotNull() const;
 499   protected:
 500     MCAuto< MEDFileFieldGlobs > _globals;
 501   };
 502
 503   class MEDFileFieldNameScope
 504   {
 505   public:
 506     MEDLOADER_EXPORT MEDFileFieldNameScope();
 507     MEDLOADER_EXPORT MEDFileFieldNameScope(const std::string& fieldName);
 508     MEDLOADER_EXPORT std::string getName() const;
 509     MEDLOADER_EXPORT void setName(const std::string& fieldName);
 510     MEDLOADER_EXPORT std::string getDtUnit() const;
 511     MEDLOADER_EXPORT void setDtUnit(const std::string& dtUnit);
 512     MEDLOADER_EXPORT void copyNameScope(const MEDFileFieldNameScope& other);
 513   protected:
 514     std::string _name;
 515     std::string _dt_unit;
 516   };
 517
 518   class MEDFileMeshes;
 519
 520   /*!
 521    * SDA is for Shared Data Arrays such as profiles.
 522    */
 523   class MEDFileAnyTypeField1TSWithoutSDA : public RefCountObject, public MEDFileFieldNameScope
 524   {
 525   public:
 526     MEDLOADER_EXPORT MEDFileAnyTypeField1TSWithoutSDA();
 527     MEDLOADER_EXPORT MEDFileAnyTypeField1TSWithoutSDA(const std::string& fieldName, int csit, int iteration, int order);
 528     MEDLOADER_EXPORT int getIteration() const { return _iteration; }
 529     MEDLOADER_EXPORT int getOrder() const { return _order; }
 530     MEDLOADER_EXPORT double getTime(int& iteration, int& order) const { iteration=_iteration; order=_order; return _dt; }
 531     MEDLOADER_EXPORT void setTime(int iteration, int order, double val) { _dt=val; _iteration=iteration; _order=order; }
 532     MEDLOADER_EXPORT int getDimension() const;
 533     MEDLOADER_EXPORT std::string getMeshName() const;
 534     MEDLOADER_EXPORT void setMeshName(const std::string& newMeshName);
 535     MEDLOADER_EXPORT bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 536     MEDLOADER_EXPORT int getMeshIteration() const;
 537     MEDLOADER_EXPORT int getMeshOrder() const;
 538     MEDLOADER_EXPORT bool isDealingTS(int iteration, int order) const;
 539     MEDLOADER_EXPORT std::pair<int,int> getDtIt() const;
 540     MEDLOADER_EXPORT void fillIteration(std::pair<int,int>& p) const;
 541     MEDLOADER_EXPORT void fillTypesOfFieldAvailable(std::vector<TypeOfField>& types) const;
 542     MEDLOADER_EXPORT std::vector<TypeOfField> getTypesOfFieldAvailable() const;
 543     //
 544     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsed2() const;
 545     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsed2() const;
 546     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsedMulti2() const;
 547     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsedMulti2() const;
 548     MEDLOADER_EXPORT void changePflsRefsNamesGen2(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 549     MEDLOADER_EXPORT void changeLocsRefsNamesGen2(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 550     //
 551     MEDLOADER_EXPORT int getNonEmptyLevels(const std::string& mname, std::vector<int>& levs) const;
 552     MEDLOADER_EXPORT std::vector< std::vector<std::pair<int,int> > > getFieldSplitedByType(const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 553     //
 554     MEDLOADER_EXPORT MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenMeshAndTypeAndLocId(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId);
 555     MEDLOADER_EXPORT const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenMeshAndTypeAndLocId(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId) const;
 556     MEDLOADER_EXPORT void deepCpyLeavesFrom(const MEDFileAnyTypeField1TSWithoutSDA& other);
 557   public:
 558     MEDLOADER_EXPORT int getNumberOfComponents() const;
 559     MEDLOADER_EXPORT const std::vector<std::string>& getInfo() const;
 560     MEDLOADER_EXPORT std::vector<std::string>& getInfo();
 561     MEDLOADER_EXPORT bool presenceOfMultiDiscPerGeoType() const;
 562     MEDLOADER_EXPORT bool presenceOfStructureElements() const;
 563     MEDLOADER_EXPORT bool onlyStructureElements() const;
 564     MEDLOADER_EXPORT void killStructureElements();
 565     MEDLOADER_EXPORT void keepOnlyStructureElements();
 566     MEDLOADER_EXPORT void setInfo(const std::vector<std::string>& infos);
 567     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 568     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 569     MEDLOADER_EXPORT int copyTinyInfoFrom(const MEDCouplingFieldDouble *field, const DataArray *arr);
 570     MEDLOADER_EXPORT void setFieldNoProfileSBT(const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 571     MEDLOADER_EXPORT void setFieldProfile(const MEDCouplingFieldDouble *field, const DataArray *arrOfVals, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
 572     MEDLOADER_EXPORT virtual void simpleRepr(int bkOffset, std::ostream& oss, int f1tsId) const;
 573     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TSWithoutSDA *deepCopy() const = 0;
 574     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TSWithoutSDA *shallowCpy() const = 0;
 575     MEDLOADER_EXPORT virtual std::vector< MCAuto<MEDFileAnyTypeField1TSWithoutSDA> > splitComponents() const;
 576     MEDLOADER_EXPORT virtual const char *getTypeStr() const = 0;
 577     MEDLOADER_EXPORT virtual DataArray *getUndergroundDataArray() const = 0;
 578     MEDLOADER_EXPORT virtual DataArray *getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const = 0;
 579     MEDLOADER_EXPORT virtual void setArray(DataArray *arr) = 0;
 580     MEDLOADER_EXPORT virtual DataArray *createNewEmptyDataArrayInstance() const = 0;
 581     MEDLOADER_EXPORT virtual DataArray *getOrCreateAndGetArray() = 0;
 582     MEDLOADER_EXPORT virtual const DataArray *getOrCreateAndGetArray() const = 0;
 583   public:
 584     MEDLOADER_EXPORT MEDCouplingFieldDouble *fieldOnMesh(const MEDFileFieldGlobsReal *glob, const MEDFileMesh *mesh, MCAuto<DataArray>& arrOut, const MEDFileFieldNameScope& nasc) const;
 585     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtLevel(TypeOfField type, int meshDimRelToMax, const std::string& mName, int renumPol, const MEDFileFieldGlobsReal *glob, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 586     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, int meshDimRelToMax, int renumPol, const MEDFileFieldGlobsReal *glob, const MEDFileMesh *mesh, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 587     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtTopLevel(TypeOfField type, const std::string& mName, int renumPol, const MEDFileFieldGlobsReal *glob, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 588     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, int renumPol, const MEDFileFieldGlobsReal *glob, const MEDCouplingMesh *mesh, const DataArrayInt *cellRenum, const DataArrayInt *nodeRenum, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
 589     DataArray *getFieldWithProfile(TypeOfField type, int meshDimRelToMax, const MEDFileMesh *mesh, DataArrayInt *&pfl, const MEDFileFieldGlobsReal *glob, const MEDFileFieldNameScope& nasc) const;
 590   public:
 591     MEDLOADER_EXPORT bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N, MEDFileFieldGlobsReal& glob);
 592     MEDLOADER_EXPORT std::vector< MCAuto<MEDFileAnyTypeField1TSWithoutSDA> > splitDiscretizations() const;
 593     MEDLOADER_EXPORT std::vector< MCAuto<MEDFileAnyTypeField1TSWithoutSDA> > splitMultiDiscrPerGeoTypes() const;
 594     MEDLOADER_EXPORT int keepOnlySpatialDiscretization(TypeOfField tof, std::vector< std::pair<int,int> >& its);
 595     MEDLOADER_EXPORT int keepOnlyGaussDiscretization(std::size_t idOfDisc, std::vector< std::pair<int,int> >& its);
 596   public:
 597     MEDLOADER_EXPORT void allocNotFromFile(int newNbOfTuples);
 598     MEDLOADER_EXPORT bool allocIfNecessaryTheArrayToReceiveDataFromFile();
 599     MEDLOADER_EXPORT void loadOnlyStructureOfDataRecursively(med_idt fid, const MEDFileFieldNameScope& nasc, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 600     MEDLOADER_EXPORT void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
 601     MEDLOADER_EXPORT void loadBigArraysRecursivelyIfNecessary(med_idt fid, const MEDFileFieldNameScope& nasc);
 602     MEDLOADER_EXPORT void loadStructureAndBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 603     MEDLOADER_EXPORT void unloadArrays();
 604     MEDLOADER_EXPORT void writeLL(med_idt fid, const MEDFileWritable& opts, const MEDFileFieldNameScope& nasc) const;
 605   protected:
 606     int getMeshIdFromMeshName(const std::string& mName) const;
 607     int addNewEntryIfNecessary(const MEDCouplingMesh *mesh);
 608     void updateData(int newLgth, const std::vector< std::pair<int,int> >& oldStartStops);
 609   protected:
 610     std::vector< MCAuto< MEDFileFieldPerMesh > > _field_per_mesh;
 611     int _iteration;
 612     int _order;
 613     double _dt;
 614   public:
 615     //! only useable on reading
 616     mutable int _csit;
 617     // -3 means allocated and build from scratch
 618     // -2 means allocated and read from a file
 619     // -1 means not allocated and build from scratch
 620     // >=0 means not allocated and read from a file
 621     mutable int _nb_of_tuples_to_be_allocated;
 622   };
 623
 624   class MEDFileIntField1TSWithoutSDA;
 625
 626   template<class T>
 627   class MEDFileField1TSTemplateWithoutSDA : public MEDFileAnyTypeField1TSWithoutSDA
 628   {
 629   protected:
 630     MEDFileField1TSTemplateWithoutSDA(const std::string& fieldName, int csit, int iteration, int order):MEDFileAnyTypeField1TSWithoutSDA(fieldName,csit,iteration,order) { }
 631     MEDFileField1TSTemplateWithoutSDA():MEDFileAnyTypeField1TSWithoutSDA() { }
 632   public:
 633     MEDLOADER_EXPORT void setArray(DataArray *arr);
 634     MEDLOADER_EXPORT DataArray *createNewEmptyDataArrayInstance() const;
 635     MEDLOADER_EXPORT typename Traits<T>::ArrayType *getOrCreateAndGetArrayTemplate();
 636     MEDLOADER_EXPORT typename Traits<T>::ArrayType const *getOrCreateAndGetArrayTemplate() const;
 637     MEDLOADER_EXPORT typename Traits<T>::ArrayType *getUndergroundDataArrayTemplate() const;
 638     MEDLOADER_EXPORT DataArray *getOrCreateAndGetArray();
 639     MEDLOADER_EXPORT const DataArray *getOrCreateAndGetArray() const;
 640     MEDLOADER_EXPORT DataArray *getUndergroundDataArray() const;
 641     MEDLOADER_EXPORT void aggregate(const typename std::vector< typename MLFieldTraits<T>::F1TSWSDAType const * >& f1tss, const std::vector< std::vector< std::pair<int,int> > >& dts);
 642   protected:
 643     MCAuto< typename Traits<T>::ArrayType > _arr;
 644   };
 645
 646   /*!
 647    * SDA is for Shared Data Arrays such as profiles.
 648    */
 649   class MEDFileField1TSWithoutSDA : public MEDFileField1TSTemplateWithoutSDA<double>
 650   {
 651   public:
 652     MEDLOADER_EXPORT const char *getTypeStr() const;
 653     MEDLOADER_EXPORT DataArray *getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 654     MEDLOADER_EXPORT DataArrayDouble *getUndergroundDataArrayDoubleExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 655     MEDLOADER_EXPORT std::vector< std::vector<DataArrayDouble *> > getFieldSplitedByType2(const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 656     MEDLOADER_EXPORT static void CheckMeshDimRel(int meshDimRelToMax);
 657     MEDLOADER_EXPORT static std::vector<int> CheckSBTMesh(const MEDCouplingMesh *mesh);
 658     MEDLOADER_EXPORT static MEDFileField1TSWithoutSDA *New(const std::string& fieldName, int csit, int iteration, int order, const std::vector<std::string>& infos);
 659   public:
 660     MEDLOADER_EXPORT MEDFileField1TSWithoutSDA();
 661     MEDLOADER_EXPORT MEDFileField1TSWithoutSDA(const std::string& fieldName, int csit, int iteration, int order, const std::vector<std::string>& infos);
 662     MEDLOADER_EXPORT MEDFileField1TSWithoutSDA *shallowCpy() const;
 663     MEDLOADER_EXPORT MEDFileField1TSWithoutSDA *deepCopy() const;
 664     MEDLOADER_EXPORT MEDFileIntField1TSWithoutSDA *convertToInt() const;
 665   public:
 666     static const char TYPE_STR[];
 667   };
 668
 669   /*!
 670    * SDA is for Shared Data Arrays such as profiles.
 671    */
 672   class MEDFileIntField1TSWithoutSDA : public MEDFileField1TSTemplateWithoutSDA<int>
 673   {
 674   public:
 675     MEDLOADER_EXPORT MEDFileIntField1TSWithoutSDA();
 676     MEDLOADER_EXPORT static MEDFileIntField1TSWithoutSDA *New(const std::string& fieldName, int csit, int iteration, int order, const std::vector<std::string>& infos);
 677     MEDLOADER_EXPORT MEDFileIntField1TSWithoutSDA *deepCopy() const;
 678     MEDLOADER_EXPORT MEDFileIntField1TSWithoutSDA *shallowCpy() const;
 679     MEDLOADER_EXPORT const char *getTypeStr() const;
 680     MEDLOADER_EXPORT DataArray *getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 681     MEDLOADER_EXPORT DataArrayInt *getUndergroundDataArrayIntExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 682     MEDLOADER_EXPORT MEDFileField1TSWithoutSDA *convertToDouble() const;
 683   protected:
 684     MEDFileIntField1TSWithoutSDA(const std::string& fieldName, int csit, int iteration, int order, const std::vector<std::string>& infos);
 685   public:
 686     MEDLOADER_EXPORT static const char TYPE_STR[];
 687   };
 688
 689   /*!
 690    * User class.
 691    */
 692   class MEDFileAnyTypeField1TS : public RefCountObject, public MEDFileWritableStandAlone, public MEDFileFieldGlobsReal
 693   {
 694   protected:
 695     MEDLOADER_EXPORT MEDFileAnyTypeField1TS();
 696     MEDLOADER_EXPORT MEDFileAnyTypeField1TS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
 697     MEDLOADER_EXPORT MEDFileAnyTypeField1TS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms);
 698     MEDLOADER_EXPORT MEDFileAnyTypeField1TS(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll, const MEDFileMeshes *ms);
 699     MEDLOADER_EXPORT MEDFileAnyTypeField1TS(const MEDFileAnyTypeField1TSWithoutSDA& other, bool shallowCopyOfContent);
 700     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *BuildNewInstanceFromContent(MEDFileAnyTypeField1TSWithoutSDA *c);
 701     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *BuildNewInstanceFromContent(MEDFileAnyTypeField1TSWithoutSDA *c, med_idt fid);
 702     MEDLOADER_EXPORT static MEDFileAnyTypeField1TSWithoutSDA *BuildContentFrom(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
 703     MEDLOADER_EXPORT static MEDFileAnyTypeField1TSWithoutSDA *BuildContentFrom(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms);
 704     MEDLOADER_EXPORT static MEDFileAnyTypeField1TSWithoutSDA *BuildContentFrom(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll, const MEDFileMeshes *ms);
 705     MEDLOADER_EXPORT void writeLL(med_idt fid) const;
 706     // direct forwarding to MEDFileAnyTypeField1TSWithoutSDA instance _content
 707   public:
 708     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(const std::string& fileName, bool loadAll=true);
 709     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(med_idt fid, bool loadAll=true);
 710     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
 711     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
 712     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(const std::string& fileName, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 713     MEDLOADER_EXPORT static MEDFileAnyTypeField1TS *New(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 714     MEDLOADER_EXPORT int getDimension() const;
 715     MEDLOADER_EXPORT int getIteration() const;
 716     MEDLOADER_EXPORT int getOrder() const;
 717     MEDLOADER_EXPORT double getTime(int& iteration, int& order) const;
 718     MEDLOADER_EXPORT void setTime(int iteration, int order, double val);
 719     MEDLOADER_EXPORT std::string getName() const;
 720     MEDLOADER_EXPORT void setName(const std::string& name);
 721     MEDLOADER_EXPORT std::string simpleRepr() const;
 722     MEDLOADER_EXPORT void simpleRepr(int bkOffset, std::ostream& oss, int f1tsId) const;
 723     MEDLOADER_EXPORT std::string getDtUnit() const;
 724     MEDLOADER_EXPORT void setDtUnit(const std::string& dtUnit);
 725     MEDLOADER_EXPORT std::string getMeshName() const;
 726     MEDLOADER_EXPORT void setMeshName(const std::string& newMeshName);
 727     MEDLOADER_EXPORT bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 728     MEDLOADER_EXPORT int getMeshIteration() const;
 729     MEDLOADER_EXPORT int getMeshOrder() const;
 730     MEDLOADER_EXPORT int getNumberOfComponents() const;
 731     MEDLOADER_EXPORT bool isDealingTS(int iteration, int order) const;
 732     MEDLOADER_EXPORT std::pair<int,int> getDtIt() const;
 733     MEDLOADER_EXPORT void fillIteration(std::pair<int,int>& p) const;
 734     MEDLOADER_EXPORT void fillTypesOfFieldAvailable(std::vector<TypeOfField>& types) const;
 735     MEDLOADER_EXPORT void setInfo(const std::vector<std::string>& infos);
 736     MEDLOADER_EXPORT const std::vector<std::string>& getInfo() const;
 737     MEDLOADER_EXPORT std::vector<std::string>& getInfo();
 738     MEDLOADER_EXPORT bool presenceOfMultiDiscPerGeoType() const;
 739     MEDLOADER_EXPORT std::vector<TypeOfField> getTypesOfFieldAvailable() const;
 740     MEDLOADER_EXPORT std::vector< std::vector<std::pair<int,int> > > getFieldSplitedByType(const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF,
 741         std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 742     MEDLOADER_EXPORT MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenMeshAndTypeAndLocId(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId);
 743     MEDLOADER_EXPORT const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenMeshAndTypeAndLocId(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId) const;
 744     MEDLOADER_EXPORT int getNonEmptyLevels(const std::string& mname, std::vector<int>& levs) const;
 745   public:
 746     MEDLOADER_EXPORT void loadArrays();
 747     MEDLOADER_EXPORT void loadArraysIfNecessary();
 748     MEDLOADER_EXPORT void unloadArrays();
 749     MEDLOADER_EXPORT void unloadArraysWithoutDataLoss();
 750     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeField1TS > > splitComponents() const;
 751     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeField1TS > > splitDiscretizations() const;
 752     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeField1TS > > splitMultiDiscrPerGeoTypes() const;
 753     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *deepCopy() const;
 754     MEDLOADER_EXPORT int copyTinyInfoFrom(const MEDCouplingFieldDouble *field, const DataArray *arr);
 755     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TS *shallowCpy() const = 0;
 756   public:
 757     //! underground method see MEDFileField1TSWithoutSDA::setProfileNameOnLeaf
 758     MEDLOADER_EXPORT void setProfileNameOnLeaf(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId, const std::string& newPflName, bool forceRenameOnGlob=false);
 759     //! underground method see MEDFileField1TSWithoutSDA::setLocNameOnLeaf
 760     MEDLOADER_EXPORT void setLocNameOnLeaf(const std::string& mName, INTERP_KERNEL::NormalizedCellType typ, int locId, const std::string& newLocName, bool forceRenameOnGlob=false);
 761     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 762     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 763     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsed() const;
 764     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsed() const;
 765     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsedMulti() const;
 766     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsedMulti() const;
 767     MEDLOADER_EXPORT void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 768     MEDLOADER_EXPORT void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 769   public:
 770     MEDLOADER_EXPORT static int LocateField2(med_idt fid, int fieldIdCFormat, bool checkFieldId, std::string& fieldName, med_field_type& typcha, std::vector<std::string>& infos, std::string& dtunitOut);
 771     MEDLOADER_EXPORT static int LocateField(med_idt fid, const std::string& fieldName, int& posCFormat, med_field_type& typcha, std::vector<std::string>& infos, std::string& dtunitOut);
 772   public:
 773     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TS *extractPart(const std::map<int, MCAuto<DataArrayInt> >& extractDef, MEDFileMesh *mm) const = 0;
 774   public:
 775     MEDLOADER_EXPORT virtual med_field_type getMEDFileFieldType() const = 0;
 776     MEDLOADER_EXPORT MEDFileAnyTypeField1TSWithoutSDA *contentNotNullBase();
 777     MEDLOADER_EXPORT const MEDFileAnyTypeField1TSWithoutSDA *contentNotNullBase() const;
 778   protected:
 779     MCAuto<MEDFileAnyTypeField1TSWithoutSDA> _content;
 780   };
 781
 782   class MEDFileIntField1TS;
 783
 784   /*!
 785    * User class.
 786    */
 787   class MEDFileField1TS : public MEDFileAnyTypeField1TS
 788   {
 789   public:
 790     MEDLOADER_EXPORT static MEDFileField1TS *New(const std::string& fileName, bool loadAll=true);
 791     MEDLOADER_EXPORT static MEDFileField1TS *New(med_idt fid, bool loadAll=true);
 792     MEDLOADER_EXPORT static MEDFileField1TS *New(DataArrayByte *db) { return BuildFromMemoryChunk<MEDFileField1TS>(db); }
 793     MEDLOADER_EXPORT static MEDFileField1TS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
 794     MEDLOADER_EXPORT static MEDFileField1TS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
 795     MEDLOADER_EXPORT static MEDFileField1TS *New(const std::string& fileName, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 796     MEDLOADER_EXPORT static MEDFileField1TS *New(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 797     MEDLOADER_EXPORT static MEDFileField1TS *New(const MEDFileField1TSWithoutSDA& other, bool shallowCopyOfContent);
 798     MEDLOADER_EXPORT static MEDFileField1TS *New();
 799     MEDLOADER_EXPORT MEDFileIntField1TS *convertToInt(bool isDeepCpyGlobs=true) const;
 800     //
 801     MEDLOADER_EXPORT MEDCouplingFieldDouble *field(const MEDFileMesh *mesh) const;
 802     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtLevel(TypeOfField type, int meshDimRelToMax, int renumPol=0) const;
 803     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtTopLevel(TypeOfField type, int renumPol=0) const;
 804     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, int meshDimRelToMax, const MEDFileMesh *mesh, int renumPol=0) const;
 805     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, const MEDCouplingMesh *mesh, int renumPol=0) const;
 806     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtLevelOld(TypeOfField type, const std::string& mname, int meshDimRelToMax, int renumPol=0) const;
 807     MEDLOADER_EXPORT DataArrayDouble *getFieldWithProfile(TypeOfField type, int meshDimRelToMax, const MEDFileMesh *mesh, DataArrayInt *&pfl) const;
 808     //
 809     MEDLOADER_EXPORT void setFieldNoProfileSBT(const MEDCouplingFieldDouble *field);
 810     MEDLOADER_EXPORT void setFieldProfile(const MEDCouplingFieldDouble *field, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile);
 811     // direct forwarding to MEDFileField1TSWithoutSDA instance _content
 812   public:
 813     MEDLOADER_EXPORT MEDFileField1TS *shallowCpy() const;
 814     MEDLOADER_EXPORT DataArrayDouble *getUndergroundDataArray() const;
 815     MEDLOADER_EXPORT DataArrayDouble *getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 816
 817     MEDLOADER_EXPORT std::vector< std::vector<DataArrayDouble *> > getFieldSplitedByType2(const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF,
 818         std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 819   public:
 820     MEDLOADER_EXPORT static void SetDataArrayDoubleInField(MEDCouplingFieldDouble *f, MCAuto<DataArray>& arr);
 821     MEDLOADER_EXPORT static DataArrayDouble *ReturnSafelyDataArrayDouble(MCAuto<DataArray>& arr);
 822   public:
 823     MEDLOADER_EXPORT MEDFileField1TS *extractPart(const std::map<int, MCAuto<DataArrayInt> >& extractDef, MEDFileMesh *mm) const;
 824   private:
 825     med_field_type getMEDFileFieldType() const { return MED_FLOAT64; }
 826     const MEDFileField1TSWithoutSDA *contentNotNull() const;
 827     MEDFileField1TSWithoutSDA *contentNotNull();
 828   private:
 829     ~MEDFileField1TS() { }
 830     MEDFileField1TS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
 831     MEDFileField1TS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms);
 832     MEDFileField1TS(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll, const MEDFileMeshes *ms);
 833     MEDFileField1TS(const MEDFileField1TSWithoutSDA& other, bool shallowCopyOfContent);
 834     MEDFileField1TS();
 835   };
 836
 837   class MEDFileIntField1TS : public MEDFileAnyTypeField1TS
 838   {
 839   public:
 840     MEDLOADER_EXPORT static MEDFileIntField1TS *New();
 841     MEDLOADER_EXPORT static MEDFileIntField1TS *New(const std::string& fileName, bool loadAll=true);
 842     MEDLOADER_EXPORT static MEDFileIntField1TS *New(med_idt fid, bool loadAll=true);
 843     MEDLOADER_EXPORT static MEDFileIntField1TS *New(DataArrayByte *db) { return BuildFromMemoryChunk<MEDFileIntField1TS>(db); }
 844     MEDLOADER_EXPORT static MEDFileIntField1TS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
 845     MEDLOADER_EXPORT static MEDFileIntField1TS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
 846     MEDLOADER_EXPORT static MEDFileIntField1TS *New(const std::string& fileName, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 847     MEDLOADER_EXPORT static MEDFileIntField1TS *New(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll=true);
 848     MEDLOADER_EXPORT static MEDFileIntField1TS *New(const MEDFileIntField1TSWithoutSDA& other, bool shallowCopyOfContent);
 849     MEDLOADER_EXPORT MEDFileField1TS *convertToDouble(bool isDeepCpyGlobs=true) const;
 850     MEDLOADER_EXPORT MEDFileIntField1TS *shallowCpy() const;
 851     //
 852     MEDLOADER_EXPORT MEDCouplingFieldInt *field(const MEDFileMesh *mesh) const;
 853     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtLevel(TypeOfField type, int meshDimRelToMax, int renumPol=0) const;
 854     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtTopLevel(TypeOfField type, int renumPol=0) const;
 855     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldOnMeshAtLevel(TypeOfField type, int meshDimRelToMax, const MEDFileMesh *mesh, int renumPol=0) const;
 856     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldOnMeshAtLevel(TypeOfField type, const MEDCouplingMesh *mesh, int renumPol=0) const;
 857     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtLevelOld(TypeOfField type, const std::string& mname, int meshDimRelToMax, int renumPol=0) const;
 858     MEDLOADER_EXPORT DataArrayInt *getFieldWithProfile(TypeOfField type, int meshDimRelToMax, const MEDFileMesh *mesh, DataArrayInt *&pfl) const;
 859     //
 860     MEDLOADER_EXPORT void setFieldNoProfileSBT(const MEDCouplingFieldInt *field);
 861     MEDLOADER_EXPORT void setFieldProfile(const MEDCouplingFieldInt *field, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile);
 862     MEDLOADER_EXPORT DataArrayInt *getUndergroundDataArray() const;
 863   public:
 864     MEDLOADER_EXPORT static DataArrayInt *ReturnSafelyDataArrayInt(MCAuto<DataArray>& arr);
 865     MEDLOADER_EXPORT static MCAuto<MEDCouplingFieldInt> SetDataArrayDoubleInIntField(MEDCouplingFieldDouble *f, MCAuto<DataArray>& arr);
 866     MEDLOADER_EXPORT static MCAuto<MEDCouplingFieldDouble> ConvertFieldIntToFieldDouble(const MEDCouplingFieldInt *f);
 867   public:
 868     MEDLOADER_EXPORT MEDFileIntField1TS *extractPart(const std::map<int, MCAuto<DataArrayInt> >& extractDef, MEDFileMesh *mm) const;
 869   private:
 870     med_field_type getMEDFileFieldType() const { return MED_INT32; }
 871     const MEDFileIntField1TSWithoutSDA *contentNotNull() const;
 872     MEDFileIntField1TSWithoutSDA *contentNotNull();
 873   private:
 874     ~MEDFileIntField1TS() { }
 875     MEDFileIntField1TS();
 876     MEDFileIntField1TS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
 877     MEDFileIntField1TS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms);
 878     MEDFileIntField1TS(med_idt fid, const std::string& fieldName, int iteration, int order, bool loadAll, const MEDFileMeshes *ms);
 879     MEDFileIntField1TS(const MEDFileIntField1TSWithoutSDA& other, bool shallowCopyOfContent);
 880   };
 881
 882   class MEDFileAnyTypeFieldMultiTSWithoutSDA : public RefCountObject, public MEDFileFieldNameScope
 883   {
 884   protected:
 885     MEDFileAnyTypeFieldMultiTSWithoutSDA();
 886     MEDFileAnyTypeFieldMultiTSWithoutSDA(const std::string& fieldName);
 887     MEDFileAnyTypeFieldMultiTSWithoutSDA(med_idt fid, int fieldId, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 888     MEDFileAnyTypeFieldMultiTSWithoutSDA(med_idt fid, const std::string& fieldName, med_field_type fieldTyp, const std::vector<std::string>& infos, int nbOfStep, const std::string& dtunit, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 889   public:
 890     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
 891     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
 892     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTSWithoutSDA *deepCopy() const;
 893     MEDLOADER_EXPORT virtual std::vector< MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> > splitComponents() const;
 894     MEDLOADER_EXPORT virtual std::vector< MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> > splitDiscretizations() const;
 895     MEDLOADER_EXPORT virtual std::vector< MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> > splitMultiDiscrPerGeoTypes() const;
 896     MEDLOADER_EXPORT virtual const char *getTypeStr() const = 0;
 897     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTSWithoutSDA *shallowCpy() const = 0;
 898     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTSWithoutSDA *createNew() const = 0;
 899     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TSWithoutSDA *createNew1TSWithoutSDAEmptyInstance() const = 0;
 900     MEDLOADER_EXPORT virtual void checkCoherencyOfType(const MEDFileAnyTypeField1TSWithoutSDA *f1ts) const = 0;
 901     MEDLOADER_EXPORT const std::vector<std::string>& getInfo() const;
 902     MEDLOADER_EXPORT bool presenceOfMultiDiscPerGeoType() const;
 903     MEDLOADER_EXPORT void setInfo(const std::vector<std::string>& info);
 904     MEDLOADER_EXPORT int getTimeStepPos(int iteration, int order) const;
 905     MEDLOADER_EXPORT const MEDFileAnyTypeField1TSWithoutSDA& getTimeStepEntry(int iteration, int order) const;
 906     MEDLOADER_EXPORT MEDFileAnyTypeField1TSWithoutSDA& getTimeStepEntry(int iteration, int order);
 907     MEDLOADER_EXPORT std::string getMeshName() const;
 908     MEDLOADER_EXPORT void setMeshName(const std::string& newMeshName);
 909     MEDLOADER_EXPORT bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
 910     MEDLOADER_EXPORT int getNumberOfTS() const;
 911     MEDLOADER_EXPORT void eraseEmptyTS();
 912     MEDLOADER_EXPORT void eraseTimeStepIds(const int *startIds, const int *endIds);
 913     MEDLOADER_EXPORT void eraseTimeStepIds2(int bg, int end, int step);
 914     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *buildFromTimeStepIds(const int *startIds, const int *endIds) const;
 915     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *buildFromTimeStepIds2(int bg, int end, int step) const;
 916     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *partOfThisLyingOnSpecifiedTimeSteps(const std::vector< std::pair<int,int> >& timeSteps) const;
 917     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *partOfThisNotLyingOnSpecifiedTimeSteps(const std::vector< std::pair<int,int> >& timeSteps) const;
 918     MEDLOADER_EXPORT bool presenceOfStructureElements() const;
 919     MEDLOADER_EXPORT bool onlyStructureElements() const;
 920     MEDLOADER_EXPORT void killStructureElements();
 921     MEDLOADER_EXPORT void keepOnlyStructureElements();
 922     MEDLOADER_EXPORT int getPosOfTimeStep(int iteration, int order) const;
 923     MEDLOADER_EXPORT int getPosGivenTime(double time, double eps=1e-8) const;
 924     MEDLOADER_EXPORT std::vector< std::pair<int,int> > getIterations() const;
 925     MEDLOADER_EXPORT std::vector< std::pair<int,int> > getTimeSteps(std::vector<double>& ret1) const;
 926     MEDLOADER_EXPORT void pushBackTimeStep(MCAuto<MEDFileAnyTypeField1TSWithoutSDA>& tse);
 927     MEDLOADER_EXPORT void synchronizeNameScope();
 928     MEDLOADER_EXPORT void simpleRepr(int bkOffset, std::ostream& oss, int fmtsId) const;
 929     MEDLOADER_EXPORT int getNonEmptyLevels(int iteration, int order, const std::string& mname, std::vector<int>& levs) const;
 930     MEDLOADER_EXPORT void appendFieldNoProfileSBT(const MEDCouplingFieldDouble *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
 931     MEDLOADER_EXPORT void appendFieldProfile(const MEDCouplingFieldDouble *field, const DataArray *arr, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile, MEDFileFieldGlobsReal& glob);
 932     MEDLOADER_EXPORT std::vector< std::vector< std::pair<int,int> > > getFieldSplitedByType(int iteration, int order, const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 933     MEDLOADER_EXPORT std::vector< std::vector<TypeOfField> > getTypesOfFieldAvailable() const;
 934     MEDLOADER_EXPORT DataArray *getUndergroundDataArray(int iteration, int order) const;
 935     MEDLOADER_EXPORT DataArray *getUndergroundDataArrayExt(int iteration, int order, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
 936     MEDLOADER_EXPORT bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N, MEDFileFieldGlobsReal& glob);
 937     MEDLOADER_EXPORT void loadStructureOrStructureAndBigArraysRecursively(med_idt fid, int nbPdt, med_field_type fieldTyp, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 938     MEDLOADER_EXPORT void writeLL(med_idt fid, const MEDFileWritable& opts) const;
 939     MEDLOADER_EXPORT void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
 940     MEDLOADER_EXPORT void loadBigArraysRecursivelyIfNecessary(med_idt fid, const MEDFileFieldNameScope& nasc);
 941     MEDLOADER_EXPORT void unloadArrays();
 942   public:
 943     MEDLOADER_EXPORT const MEDFileAnyTypeField1TSWithoutSDA *getTimeStepAtPos2(int pos) const;
 944     MEDLOADER_EXPORT MEDFileAnyTypeField1TSWithoutSDA *getTimeStepAtPos2(int pos);
 945     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsed2() const;
 946     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsed2() const;
 947     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsedMulti2() const;
 948     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsedMulti2() const;
 949     MEDLOADER_EXPORT void changePflsRefsNamesGen2(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 950     MEDLOADER_EXPORT void changeLocsRefsNamesGen2(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
 951     MEDLOADER_EXPORT void setIteration(int i, MCAuto<MEDFileAnyTypeField1TSWithoutSDA> ts);
 952   protected:
 953     virtual med_field_type getMEDFileFieldType() const = 0;
 954     void copyTinyInfoFrom(const MEDCouplingFieldDouble *field, const DataArray *arr);
 955     void checkCoherencyOfTinyInfo(const MEDCouplingFieldDouble *field, const DataArray *arr) const;
 956     void checkThatComponentsMatch(const std::vector<std::string>& compos) const;
 957     void checkThatNbOfCompoOfTSMatchThis() const;
 958   protected:
 959     std::vector<std::string> _infos;
 960     std::vector< MCAuto<MEDFileAnyTypeField1TSWithoutSDA> > _time_steps;
 961   };
 962
 963   class MEDFileIntFieldMultiTSWithoutSDA;
 964
 965   class MEDFileFieldMultiTSWithoutSDA : public MEDFileAnyTypeFieldMultiTSWithoutSDA
 966   {
 967   public:
 968     MEDLOADER_EXPORT static MEDFileFieldMultiTSWithoutSDA *New(med_idt fid, const std::string& fieldName, med_field_type fieldTyp, const std::vector<std::string>& infos, int nbOfStep, const std::string& dtunit, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 969     MEDLOADER_EXPORT MEDFileFieldMultiTSWithoutSDA(med_idt fid, int fieldId, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 970     MEDLOADER_EXPORT const char *getTypeStr() const;
 971     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *shallowCpy() const;
 972     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *createNew() const;
 973     MEDLOADER_EXPORT std::vector< std::vector<DataArrayDouble *> > getFieldSplitedByType2(int iteration, int order, const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
 974     MEDLOADER_EXPORT MEDFileIntFieldMultiTSWithoutSDA *convertToInt() const;
 975   protected:
 976     MEDFileFieldMultiTSWithoutSDA(const std::string& fieldName);
 977     MEDFileFieldMultiTSWithoutSDA(med_idt fid, const std::string& fieldName, med_field_type fieldTyp, const std::vector<std::string>& infos, int nbOfStep, const std::string& dtunit, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 978     med_field_type getMEDFileFieldType() const { return MED_FLOAT64; }
 979     MEDFileAnyTypeField1TSWithoutSDA *createNew1TSWithoutSDAEmptyInstance() const;
 980     void checkCoherencyOfType(const MEDFileAnyTypeField1TSWithoutSDA *f1ts) const;
 981   public:
 982     MEDLOADER_EXPORT MEDFileFieldMultiTSWithoutSDA();
 983   };
 984
 985   class MEDFileIntFieldMultiTSWithoutSDA : public MEDFileAnyTypeFieldMultiTSWithoutSDA
 986   {
 987   public:
 988     MEDLOADER_EXPORT static MEDFileIntFieldMultiTSWithoutSDA *New(med_idt fid, const std::string& fieldName, med_field_type fieldTyp, const std::vector<std::string>& infos, int nbOfStep, const std::string& dtunit, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 989     MEDLOADER_EXPORT MEDFileIntFieldMultiTSWithoutSDA(med_idt fid, int fieldId, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 990     MEDLOADER_EXPORT const char *getTypeStr() const;
 991     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *shallowCpy() const;
 992     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSWithoutSDA *createNew() const;
 993     MEDLOADER_EXPORT MEDFileFieldMultiTSWithoutSDA *convertToDouble() const;
 994   protected:
 995     MEDFileIntFieldMultiTSWithoutSDA(const std::string& fieldName);
 996     MEDFileIntFieldMultiTSWithoutSDA(med_idt fid, const std::string& fieldName, med_field_type fieldTyp, const std::vector<std::string>& infos, int nbOfStep, const std::string& dtunit, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
 997     med_field_type getMEDFileFieldType() const { return MED_INT32; }
 998     MEDFileAnyTypeField1TSWithoutSDA *createNew1TSWithoutSDAEmptyInstance() const;
 999     void checkCoherencyOfType(const MEDFileAnyTypeField1TSWithoutSDA *f1ts) const;
1000   public:
1001     MEDLOADER_EXPORT MEDFileIntFieldMultiTSWithoutSDA();
1002   };
1003
1004   class MEDFileAnyTypeFieldMultiTSIterator;
1005   class MEDFileFastCellSupportComparator;
1006   /*!
1007    * User class.
1008    */
1009   class MEDFileAnyTypeFieldMultiTS : public RefCountObject, public MEDFileWritableStandAlone, public MEDFileFieldGlobsReal
1010   {
1011   protected:
1012     MEDFileAnyTypeFieldMultiTS();
1013     MEDFileAnyTypeFieldMultiTS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
1014     MEDFileAnyTypeFieldMultiTS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities=0);
1015     MEDFileAnyTypeFieldMultiTS(const MEDFileAnyTypeFieldMultiTSWithoutSDA& other, bool shallowCopyOfContent);
1016     static MEDFileAnyTypeFieldMultiTS *BuildNewInstanceFromContent(MEDFileAnyTypeFieldMultiTSWithoutSDA *c);
1017     static MEDFileAnyTypeFieldMultiTS *BuildNewInstanceFromContent(MEDFileAnyTypeFieldMultiTSWithoutSDA *c, med_idt fid);
1018     static MEDFileAnyTypeFieldMultiTSWithoutSDA *BuildContentFrom(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
1019     static MEDFileAnyTypeFieldMultiTSWithoutSDA *BuildContentFrom(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
1020   public:
1021     MEDLOADER_EXPORT static MEDFileAnyTypeFieldMultiTS *New(const std::string& fileName, bool loadAll=true);
1022     MEDLOADER_EXPORT static MEDFileAnyTypeFieldMultiTS *New(med_idt fid, bool loadAll=true);
1023     MEDLOADER_EXPORT static MEDFileAnyTypeFieldMultiTS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
1024     MEDLOADER_EXPORT static MEDFileAnyTypeFieldMultiTS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
1025     MEDLOADER_EXPORT void loadArrays();
1026     MEDLOADER_EXPORT void loadArraysIfNecessary();
1027     MEDLOADER_EXPORT void unloadArrays();
1028     MEDLOADER_EXPORT void unloadArraysWithoutDataLoss();
1029     MEDLOADER_EXPORT void writeLL(med_idt fid) const;
1030     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
1031     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
1032     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTS *deepCopy() const;
1033     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeFieldMultiTS > > splitComponents() const;
1034     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeFieldMultiTS > > splitDiscretizations() const;
1035     MEDLOADER_EXPORT std::vector< MCAuto< MEDFileAnyTypeFieldMultiTS > > splitMultiDiscrPerGeoTypes() const;
1036     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTS *shallowCpy() const = 0;
1037     MEDLOADER_EXPORT virtual void checkCoherencyOfType(const MEDFileAnyTypeField1TS *f1ts) const = 0;
1038     //
1039     MEDLOADER_EXPORT virtual MEDFileAnyTypeField1TS *getTimeStepAtPos(int pos) const = 0;
1040     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *getTimeStep(int iteration, int order) const;
1041     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *getTimeStepGivenTime(double time, double eps=1e-8) const;
1042     MEDLOADER_EXPORT static std::vector< std::vector<MEDFileAnyTypeFieldMultiTS *> > SplitIntoCommonTimeSeries(const std::vector<MEDFileAnyTypeFieldMultiTS *>& vectFMTS);
1043     MEDLOADER_EXPORT static std::vector< std::vector<MEDFileAnyTypeFieldMultiTS *> > SplitPerCommonSupport(const std::vector<MEDFileAnyTypeFieldMultiTS *>& vectFMTS, const MEDFileMesh *mesh, std::vector< MCAuto<MEDFileFastCellSupportComparator> >& fsc);
1044     MEDLOADER_EXPORT static int CheckSupportAcrossTime(MEDFileAnyTypeFieldMultiTS *f0, MEDFileAnyTypeFieldMultiTS *f1, const MEDFileMesh *mesh, TypeOfField& tof0, TypeOfField& tof1);
1045   public:// direct forwarding to MEDFileField1TSWithoutSDA instance _content
1046     MEDLOADER_EXPORT std::string getName() const;
1047     MEDLOADER_EXPORT void setName(const std::string& name);
1048     MEDLOADER_EXPORT std::string getDtUnit() const;
1049     MEDLOADER_EXPORT void setDtUnit(const std::string& dtUnit);
1050     MEDLOADER_EXPORT std::string getMeshName() const;
1051     MEDLOADER_EXPORT void setMeshName(const std::string& newMeshName);
1052     MEDLOADER_EXPORT std::string simpleRepr() const;
1053     MEDLOADER_EXPORT void simpleRepr(int bkOffset, std::ostream& oss, int fmtsId) const;
1054     MEDLOADER_EXPORT int getNumberOfTS() const;
1055     MEDLOADER_EXPORT void eraseEmptyTS();
1056     MEDLOADER_EXPORT void eraseTimeStepIds(const int *startIds, const int *endIds);
1057     MEDLOADER_EXPORT void eraseTimeStepIds2(int bg, int end, int step);
1058     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *buildSubPart(const int *startIds, const int *endIds) const;
1059     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *buildSubPartSlice(int bg, int end, int step) const;
1060     MEDLOADER_EXPORT std::vector< std::pair<int,int> > getTimeSteps(std::vector<double>& ret1) const;
1061     MEDLOADER_EXPORT std::vector< std::pair<int,int> > getIterations() const;
1062     MEDLOADER_EXPORT void pushBackTimeSteps(const std::vector<MEDFileAnyTypeField1TS *>& f1ts);
1063     MEDLOADER_EXPORT void pushBackTimeSteps(MEDFileAnyTypeFieldMultiTS *fmts);
1064     MEDLOADER_EXPORT void pushBackTimeStep(MEDFileAnyTypeField1TS *f1ts);
1065     MEDLOADER_EXPORT void synchronizeNameScope();
1066     MEDLOADER_EXPORT int getPosOfTimeStep(int iteration, int order) const;
1067     MEDLOADER_EXPORT int getPosGivenTime(double time, double eps=1e-8) const;
1068     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSIterator *iterator();
1069     MEDLOADER_EXPORT bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
1070     MEDLOADER_EXPORT const std::vector<std::string>& getInfo() const;
1071     MEDLOADER_EXPORT bool presenceOfMultiDiscPerGeoType() const;
1072     MEDLOADER_EXPORT void killStructureElements();
1073     MEDLOADER_EXPORT void keepOnlyStructureElements();
1074     MEDLOADER_EXPORT void setInfo(const std::vector<std::string>& info);
1075     MEDLOADER_EXPORT int getNumberOfComponents() const;
1076     MEDLOADER_EXPORT int getNonEmptyLevels(int iteration, int order, const std::string& mname, std::vector<int>& levs) const;
1077     MEDLOADER_EXPORT std::vector< std::vector<TypeOfField> > getTypesOfFieldAvailable() const;
1078     MEDLOADER_EXPORT std::vector< std::vector< std::pair<int,int> > > getFieldSplitedByType(int iteration, int order, const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
1079     MEDLOADER_EXPORT MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> getContent();
1080   public:
1081     MEDLOADER_EXPORT virtual MEDFileAnyTypeFieldMultiTS *buildNewEmpty() const = 0;
1082     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *extractPart(const std::map<int, MCAuto<DataArrayInt> >& extractDef, MEDFileMesh *mm) const;
1083     MEDLOADER_EXPORT static MCAuto<MEDFileAnyTypeFieldMultiTS> Aggregate(const std::vector<const MEDFileAnyTypeFieldMultiTS *>& fmtss, const std::vector< std::vector< std::pair<int,int> > >& dts);
1084   public:
1085     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsed() const;
1086     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsed() const;
1087     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsedMulti() const;
1088     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsedMulti() const;
1089     MEDLOADER_EXPORT void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
1090     MEDLOADER_EXPORT void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
1091   protected:
1092     MEDFileAnyTypeFieldMultiTSWithoutSDA *contentNotNullBase();
1093     const MEDFileAnyTypeFieldMultiTSWithoutSDA *contentNotNullBase() const;
1094   private:
1095     static std::vector< std::vector<MEDFileAnyTypeFieldMultiTS *> > SplitPerCommonSupportNotNodesAlg(const std::vector<MEDFileAnyTypeFieldMultiTS *>& vectFMTS, const MEDFileMesh *mesh, std::vector< MCAuto<MEDFileFastCellSupportComparator> >& cmps);
1096   protected:
1097     MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> _content;
1098   };
1099
1100   class MEDFileIntFieldMultiTS;
1101
1102   /*!
1103    * User class.
1104    */
1105   class MEDFileFieldMultiTS : public MEDFileAnyTypeFieldMultiTS
1106   {
1107   public:
1108     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New();
1109     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(const std::string& fileName, bool loadAll=true);
1110     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(med_idt fid, bool loadAll=true);
1111     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(DataArrayByte *db) { return BuildFromMemoryChunk<MEDFileFieldMultiTS>(db); }
1112     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
1113     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
1114     MEDLOADER_EXPORT static MEDFileFieldMultiTS *New(const MEDFileFieldMultiTSWithoutSDA& other, bool shallowCopyOfContent);
1115     MEDLOADER_EXPORT static MEDFileFieldMultiTS *LoadSpecificEntities(const std::string& fileName, const std::string& fieldName, const std::vector< std::pair<TypeOfField,INTERP_KERNEL::NormalizedCellType> >& entities, bool loadAll=true);
1116     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *shallowCpy() const;
1117     MEDLOADER_EXPORT void checkCoherencyOfType(const MEDFileAnyTypeField1TS *f1ts) const;
1118     MEDLOADER_EXPORT MEDFileIntFieldMultiTS *convertToInt(bool isDeepCpyGlobs=true) const;
1119     //
1120     MEDLOADER_EXPORT MEDFileField1TS *getTimeStepAtPos(int pos) const;
1121     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *getTimeStep(int iteration, int order) const;
1122     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *getTimeStepGivenTime(double time, double eps=1e-8) const;
1123     //
1124     MEDLOADER_EXPORT MEDCouplingFieldDouble *field(int iteration, int order, const MEDFileMesh *mesh) const;
1125     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtLevel(TypeOfField type, int iteration, int order, int meshDimRelToMax, int renumPol=0) const;
1126     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtTopLevel(TypeOfField type, int iteration, int order, int renumPol=0) const;
1127     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, int iteration, int order, int meshDimRelToMax, const MEDFileMesh *mesh, int renumPol=0) const;
1128     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, int iteration, int order, const MEDCouplingMesh *mesh, int renumPol=0) const;
1129     MEDLOADER_EXPORT MEDCouplingFieldDouble *getFieldAtLevelOld(TypeOfField type, const std::string& mname, int iteration, int order, int meshDimRelToMax, int renumPol=0) const;
1130     MEDLOADER_EXPORT DataArrayDouble *getFieldWithProfile(TypeOfField type, int iteration, int order, int meshDimRelToMax, const MEDFileMesh *mesh, DataArrayInt *&pfl) const;
1131     //
1132     MEDLOADER_EXPORT void appendFieldNoProfileSBT(const MEDCouplingFieldDouble *field);
1133     MEDLOADER_EXPORT void appendFieldProfile(const MEDCouplingFieldDouble *field, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile);
1134     MEDLOADER_EXPORT std::vector< std::vector<DataArrayDouble *> > getFieldSplitedByType2(int iteration, int order, const std::string& mname, std::vector<INTERP_KERNEL::NormalizedCellType>& types, std::vector< std::vector<TypeOfField> >& typesF, std::vector< std::vector<std::string> >& pfls, std::vector< std::vector<std::string> >& locs) const;
1135     MEDLOADER_EXPORT DataArrayDouble *getUndergroundDataArray(int iteration, int order) const;
1136     MEDLOADER_EXPORT DataArrayDouble *getUndergroundDataArrayExt(int iteration, int order, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
1137   public:
1138     MEDLOADER_EXPORT MEDFileFieldMultiTS *buildNewEmpty() const;
1139   private:
1140     const MEDFileFieldMultiTSWithoutSDA *contentNotNull() const;
1141     MEDFileFieldMultiTSWithoutSDA *contentNotNull();
1142   private:
1143     ~MEDFileFieldMultiTS() { }
1144     MEDFileFieldMultiTS();
1145     MEDFileFieldMultiTS(const MEDFileFieldMultiTSWithoutSDA& other, bool shallowCopyOfContent);
1146     MEDFileFieldMultiTS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
1147     MEDFileFieldMultiTS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities=0);
1148   };
1149
1150   /*!
1151    * User class.
1152    */
1153   class MEDFileIntFieldMultiTS : public MEDFileAnyTypeFieldMultiTS
1154   {
1155   public:
1156     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New();
1157     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(const std::string& fileName, bool loadAll=true);
1158     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(med_idt fid, bool loadAll=true);
1159     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(DataArrayByte *db) { return BuildFromMemoryChunk<MEDFileIntFieldMultiTS>(db); }
1160     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(const std::string& fileName, const std::string& fieldName, bool loadAll=true);
1161     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(med_idt fid, const std::string& fieldName, bool loadAll=true);
1162     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *New(const MEDFileIntFieldMultiTSWithoutSDA& other, bool shallowCopyOfContent);
1163     MEDLOADER_EXPORT static MEDFileIntFieldMultiTS *LoadSpecificEntities(const std::string& fileName, const std::string& fieldName, const std::vector< std::pair<TypeOfField,INTERP_KERNEL::NormalizedCellType> >& entities, bool loadAll=true);
1164     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *shallowCpy() const;
1165     MEDLOADER_EXPORT void checkCoherencyOfType(const MEDFileAnyTypeField1TS *f1ts) const;
1166     MEDLOADER_EXPORT MEDFileIntField1TS *getTimeStepAtPos(int pos) const;
1167     MEDLOADER_EXPORT MEDFileFieldMultiTS *convertToDouble(bool isDeepCpyGlobs=true) const;
1168     //
1169     MEDLOADER_EXPORT MEDCouplingFieldInt *field(int iteration, int order, const MEDFileMesh *mesh) const;
1170     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtLevel(TypeOfField type, int iteration, int order, int meshDimRelToMax, int renumPol=0) const;
1171     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtTopLevel(TypeOfField type, int iteration, int order, int renumPol=0) const;
1172     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldOnMeshAtLevel(TypeOfField type, int iteration, int order, int meshDimRelToMax, const MEDFileMesh *mesh, int renumPol=0) const;
1173     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldOnMeshAtLevel(TypeOfField type, int iteration, int order, const MEDCouplingMesh *mesh, int renumPol=0) const;
1174     MEDLOADER_EXPORT MEDCouplingFieldInt *getFieldAtLevelOld(TypeOfField type, int iteration, int order, const std::string& mname, int meshDimRelToMax, int renumPol=0) const;
1175     MEDLOADER_EXPORT DataArrayInt *getFieldWithProfile(TypeOfField type, int iteration, int order, int meshDimRelToMax, const MEDFileMesh *mesh, DataArrayInt *&pfl) const;
1176     //
1177     MEDLOADER_EXPORT void appendFieldNoProfileSBT(const MEDCouplingFieldInt *field);
1178     MEDLOADER_EXPORT void appendFieldProfile(const MEDCouplingFieldInt *field, const MEDFileMesh *mesh, int meshDimRelToMax, const DataArrayInt *profile);
1179     //
1180     MEDLOADER_EXPORT DataArrayInt *getUndergroundDataArray(int iteration, int order) const;
1181   public:
1182     MEDLOADER_EXPORT MEDFileIntFieldMultiTS *buildNewEmpty() const;
1183   private:
1184     const MEDFileIntFieldMultiTSWithoutSDA *contentNotNull() const;
1185     MEDFileIntFieldMultiTSWithoutSDA *contentNotNull();
1186   private:
1187     ~MEDFileIntFieldMultiTS() { }
1188     MEDFileIntFieldMultiTS();
1189     MEDFileIntFieldMultiTS(const MEDFileIntFieldMultiTSWithoutSDA& other, bool shallowCopyOfContent);
1190     MEDFileIntFieldMultiTS(med_idt fid, bool loadAll, const MEDFileMeshes *ms);
1191     MEDFileIntFieldMultiTS(med_idt fid, const std::string& fieldName, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities=0);
1192   };
1193
1194   class MEDFileAnyTypeFieldMultiTSIterator
1195   {
1196   public:
1197     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTSIterator(MEDFileAnyTypeFieldMultiTS *fmts);
1198     MEDLOADER_EXPORT ~MEDFileAnyTypeFieldMultiTSIterator();
1199     MEDLOADER_EXPORT MEDFileAnyTypeField1TS *nextt();
1200   private:
1201     MCAuto<MEDFileAnyTypeFieldMultiTS> _fmts;
1202     int _iter_id;
1203     int _nb_iter;
1204   };
1205
1206   class MEDFileFieldsIterator;
1207   class MEDFileStructureElements;
1208
1209   /*!
1210    * Use class.
1211    */
1212   class MEDFileFields : public RefCountObject, public MEDFileFieldGlobsReal, public MEDFileWritableStandAlone
1213   {
1214   public:
1215     MEDLOADER_EXPORT static MEDFileFields *New();
1216     MEDLOADER_EXPORT static MEDFileFields *New(const std::string& fileName, bool loadAll=true);
1217     MEDLOADER_EXPORT static MEDFileFields *New(med_idt fid, bool loadAll=true);
1218     MEDLOADER_EXPORT static MEDFileFields *NewWithDynGT(const std::string& fileName, const MEDFileStructureElements *se, bool loadAll=true);
1219     MEDLOADER_EXPORT static MEDFileFields *NewWithDynGT(med_idt fid, const MEDFileStructureElements *se, bool loadAll=true);
1220     MEDLOADER_EXPORT static MEDFileFields *New(DataArrayByte *db) { return BuildFromMemoryChunk<MEDFileFields>(db); }
1221     MEDLOADER_EXPORT static MEDFileFields *LoadPartOf(const std::string& fileName, bool loadAll=true, const MEDFileMeshes *ms=0);
1222     MEDLOADER_EXPORT static MEDFileFields *LoadSpecificEntities(const std::string& fileName, const std::vector< std::pair<TypeOfField,INTERP_KERNEL::NormalizedCellType> >& entities, bool loadAll=true);
1223     MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
1224     MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
1225     MEDLOADER_EXPORT MEDFileFields *deepCopy() const;
1226     MEDLOADER_EXPORT MEDFileFields *shallowCpy() const;
1227     MEDLOADER_EXPORT void writeLL(med_idt fid) const;
1228     MEDLOADER_EXPORT void loadArrays();
1229     MEDLOADER_EXPORT void loadArraysIfNecessary();
1230     MEDLOADER_EXPORT void unloadArrays();
1231     MEDLOADER_EXPORT void unloadArraysWithoutDataLoss();
1232     MEDLOADER_EXPORT int getNumberOfFields() const;
1233     MEDLOADER_EXPORT std::vector< std::pair<int,int> > getCommonIterations(bool& areThereSomeForgottenTS) const;
1234     MEDLOADER_EXPORT std::vector<std::string> getFieldsNames() const;
1235     MEDLOADER_EXPORT std::vector<std::string> getMeshesNames() const;
1236     MEDLOADER_EXPORT std::string simpleRepr() const;
1237     MEDLOADER_EXPORT void simpleRepr(int bkOffset, std::ostream& oss) const;
1238     //
1239     MEDLOADER_EXPORT void resize(int newSize);
1240     MEDLOADER_EXPORT void pushField(MEDFileAnyTypeFieldMultiTS *field);
1241     MEDLOADER_EXPORT void pushFields(const std::vector<MEDFileAnyTypeFieldMultiTS *>& fields);
1242     MEDLOADER_EXPORT void setFieldAtPos(int i, MEDFileAnyTypeFieldMultiTS *field);
1243     MEDLOADER_EXPORT int getPosFromFieldName(const std::string& fieldName) const;
1244     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *getFieldAtPos(int i) const;
1245     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *getFieldWithName(const std::string& fieldName) const;
1246     MEDLOADER_EXPORT MEDFileFields *buildSubPart(const int *startIds, const int *endIds) const;
1247     MEDLOADER_EXPORT bool removeFieldsWithoutAnyTimeStep();
1248     MEDLOADER_EXPORT MEDFileFields *partOfThisLyingOnSpecifiedMeshName(const std::string& meshName) const;
1249     MEDLOADER_EXPORT MEDFileFields *partOfThisLyingOnSpecifiedTimeSteps(const std::vector< std::pair<int,int> >& timeSteps) const;
1250     MEDLOADER_EXPORT MEDFileFields *partOfThisNotLyingOnSpecifiedTimeSteps(const std::vector< std::pair<int,int> >& timeSteps) const;
1251     MEDLOADER_EXPORT bool presenceOfStructureElements() const;
1252     MEDLOADER_EXPORT void killStructureElements();
1253     MEDLOADER_EXPORT void keepOnlyStructureElements();
1254     MEDLOADER_EXPORT MCAuto<MEDFileFields> partOfThisOnStructureElements() const;
1255     MEDLOADER_EXPORT MEDFileFieldsIterator *iterator();
1256     MEDLOADER_EXPORT void destroyFieldAtPos(int i);
1257     MEDLOADER_EXPORT void destroyFieldsAtPos(const int *startIds, const int *endIds);
1258     MEDLOADER_EXPORT void destroyFieldsAtPos2(int bg, int end, int step);
1259     MEDLOADER_EXPORT bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
1260     MEDLOADER_EXPORT bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N);
1261   public:
1262     MEDLOADER_EXPORT MEDFileFields *extractPart(const std::map<int, MCAuto<DataArrayInt> >& extractDef, MEDFileMesh *mm) const;
1263   public:
1264     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsed() const;
1265     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsed() const;
1266     MEDLOADER_EXPORT std::vector<std::string> getPflsReallyUsedMulti() const;
1267     MEDLOADER_EXPORT std::vector<std::string> getLocsReallyUsedMulti() const;
1268     MEDLOADER_EXPORT void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
1269     MEDLOADER_EXPORT void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
1270   private:
1271     ~MEDFileFields() { }
1272     MEDFileFields();
1273     MEDFileFields(med_idt fid, bool loadAll, const MEDFileMeshes *ms, const MEDFileEntities *entities);
1274   private:
1275     std::vector< MCAuto<MEDFileAnyTypeFieldMultiTSWithoutSDA> > _fields;
1276   };
1277
1278   class MEDFileFieldsIterator
1279   {
1280   public:
1281     MEDLOADER_EXPORT MEDFileFieldsIterator(MEDFileFields *fs);
1282     MEDLOADER_EXPORT ~MEDFileFieldsIterator();
1283     MEDLOADER_EXPORT MEDFileAnyTypeFieldMultiTS *nextt();
1284   private:
1285     MCAuto<MEDFileFields> _fs;
1286     int _iter_id;
1287     int _nb_iter;
1288   };
1289 }
1290
1291 #endif