#include "MEDLoaderDefines.hxx"
+#include "MEDFileFieldInternal.hxx"
+#include "MEDFileFieldGlobs.hxx"
#include "MEDFileFieldOverView.hxx"
#include "MEDFileUtilities.txx"
#include "MEDFileEntities.hxx"
class MEDFileMesh;
class MEDFileFieldVisitor;
- class MEDFileGTKeeper
- {
- public:
- virtual MEDFileGTKeeper *deepCopy() const = 0;
- virtual INTERP_KERNEL::NormalizedCellType getGeoType() const = 0;
- virtual std::string getRepr() const = 0;
- virtual bool isEqual(const MEDFileGTKeeper *other) const = 0;
- virtual ~MEDFileGTKeeper();
- };
-
- class MEDFileGTKeeperSta : public MEDFileGTKeeper
- {
- public:
- MEDFileGTKeeperSta(INTERP_KERNEL::NormalizedCellType gt):_geo_type(gt) { }
- MEDFileGTKeeper *deepCopy() const;
- INTERP_KERNEL::NormalizedCellType getGeoType() const;
- std::string getRepr() const;
- bool isEqual(const MEDFileGTKeeper *other) const;
- private:
- INTERP_KERNEL::NormalizedCellType _geo_type;
- };
-
-
- class MEDFileGTKeeperDyn : public MEDFileGTKeeper
- {
- public:
- MEDFileGTKeeperDyn(const MEDFileUMesh *mesh, const MEDFileUMesh *section, const MEDFileStructureElement *se);
- MEDFileGTKeeper *deepCopy() const;
- INTERP_KERNEL::NormalizedCellType getGeoType() const;
- std::string getRepr() const;
- bool isEqual(const MEDFileGTKeeper *other) const;
- const MEDFileUMesh *getMesh() const { return _mesh; }
- const MEDFileUMesh *getSection() const { return _section; }
- const MEDFileStructureElement *getSE() const { return _se; }
- private:
- MCConstAuto<MEDFileUMesh> _mesh;
- MCConstAuto<MEDFileUMesh> _section;
- MCConstAuto<MEDFileStructureElement> _se;
- };
-
- class MEDFileFieldLoc : public RefCountObject
- {
- public:
- MEDLOADER_EXPORT void simpleRepr(std::ostream& oss) const;
- MEDLOADER_EXPORT std::string getName() const { return _name; }
- MEDLOADER_EXPORT void setName(const std::string& name);
- static MEDFileFieldLoc *New(med_idt fid, const std::string& locName);
- static MEDFileFieldLoc *New(med_idt fid, int i, const MEDFileEntities *entities);
- 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);
- std::size_t getHeapMemorySizeWithoutChildren() const;
- std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- MEDFileFieldLoc *deepCopy() const;
- bool isOnStructureElement() const;
- const MEDFileGTKeeper *getUndergroundGTKeeper() const { return _gt; }
- MEDLOADER_EXPORT int getNbOfGaussPtPerCell() const { return _nb_gauss_pt; }
- MEDLOADER_EXPORT void writeLL(med_idt fid) const;
- MEDLOADER_EXPORT std::string repr() const;
- MEDLOADER_EXPORT bool isName(const std::string& name) const { return _name==name; }
- MEDLOADER_EXPORT int getDimension() const { return _dim; }
- MEDLOADER_EXPORT int getNumberOfGaussPoints() const { return _nb_gauss_pt; }
- MEDLOADER_EXPORT int getNumberOfPointsInCells() const { return _nb_node_per_cell; }
- MEDLOADER_EXPORT const std::vector<double>& getRefCoords() const { return _ref_coo; }
- MEDLOADER_EXPORT const std::vector<double>& getGaussCoords() const { return _gs_coo; }
- MEDLOADER_EXPORT const std::vector<double>& getGaussWeights() const { return _w; }
- MEDLOADER_EXPORT INTERP_KERNEL::NormalizedCellType getGeoType() const { return _gt->getGeoType(); }
- MEDLOADER_EXPORT bool isEqual(const MEDFileFieldLoc& other, double eps) const;
- private:
- MEDFileFieldLoc(const MEDFileFieldLoc& other);
- MEDFileFieldLoc(med_idt fid, const std::string& locName);
- MEDFileFieldLoc(med_idt fid, int id, const MEDFileEntities *entities);
- MEDFileFieldLoc(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
- private:
- int _dim;
- int _nb_gauss_pt;
- INTERP_KERNEL::AutoCppPtr<MEDFileGTKeeper> _gt;
- int _nb_node_per_cell;
- std::string _name;
- std::vector<double> _ref_coo;
- std::vector<double> _gs_coo;
- std::vector<double> _w;
- };
-
- /// @cond INTERNAL
- class MEDFileAnyTypeField1TSWithoutSDA;
- class MEDFileFieldPerMeshPerTypeCommon;
- class MEDFileFieldPerMeshPerType;
- class MEDFileField1TSWithoutSDA;
- class MEDFileFieldNameScope;
- class MEDFileFieldGlobsReal;
- class MEDFileFieldPerMesh;
-
- class MEDFileFieldPerMeshPerTypePerDisc : public RefCountObject, public MEDFileWritable
- {
- public:
- static MEDFileFieldPerMeshPerTypePerDisc *NewOnRead(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
- static MEDFileFieldPerMeshPerTypePerDisc *New(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int locId);
- static MEDFileFieldPerMeshPerTypePerDisc *New(const MEDFileFieldPerMeshPerTypePerDisc& other);
- std::size_t getHeapMemorySizeWithoutChildren() const;
- std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- MEDFileFieldPerMeshPerTypePerDisc *deepCopy(MEDFileFieldPerMeshPerTypeCommon *father) const;
- void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldTemplate *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldTemplate *field, const DataArray *arrr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob);
- void getCoarseData(TypeOfField& type, std::pair<int,int>& dad, std::string& pfl, std::string& loc) const;
- void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
- const MEDFileFieldPerMeshPerTypeCommon *getFather() const;
- void loadOnlyStructureOfDataRecursively(med_idt fid, int& start, const MEDFileFieldNameScope& nasc);
- void loadBigArray(med_idt fid, const MEDFileFieldNameScope& nasc);
- void setNewStart(int newValueOfStart);
- int getIteration() const;
- int getOrder() const;
- double getTime() const;
- std::string getMeshName() const;
- TypeOfField getType() const;
- void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
- void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
- void setType(TypeOfField newType);
- INTERP_KERNEL::NormalizedCellType getGeoType() const;
- int getNumberOfComponents() const;
- int getNumberOfTuples() const;
- int getStart() const { return _start; }
- int getEnd() const { return _end; }
- void setEnd(int endd) { _end=endd; }
- int getNumberOfVals() const { return _nval; }
- DataArray *getOrCreateAndGetArray();
- const DataArray *getOrCreateAndGetArray() const;
- const std::vector<std::string>& getInfo() const;
- std::string getProfile() const;
- void setProfile(const std::string& newPflName);
- std::string getLocalization() const;
- void setLocalization(const std::string& newLocName);
- int getLocId() const { return _loc_id; }
- void setLocId(int newId) const { _loc_id=newId; }
- void setFather(MEDFileFieldPerMeshPerTypeCommon *newFather) { _father=newFather; }
- void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- void getFieldAtLevel(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;
- void fillValues(int discId, int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
- int fillEltIdsFromCode(int offset, const std::vector<int>& codeOfMesh, const MEDFileFieldGlobsReal& glob, int *ptToFill) const;
- int fillTupleIds(int *ptToFill) const;
- static int ConvertType(TypeOfField type, int locId);
- static std::vector< std::vector< const MEDFileFieldPerMeshPerTypePerDisc *> > SplitPerDiscretization(const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entries);
- static bool RenumberChunks(int offset, const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
- const DataArrayInt *explicitIdsInMesh, const std::vector<int>& newCode,
- MEDFileFieldGlobsReal& glob, DataArrayDouble *arr, std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> >& result);
- static MEDFileFieldPerMeshPerTypePerDisc *NewObjectOnSameDiscThanPool(TypeOfField typeF, INTERP_KERNEL::NormalizedCellType geoType, DataArrayInt *idsOfMeshElt,
- bool isPfl, int nbi, int offset, std::list< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
- MEDFileFieldGlobsReal& glob, bool ¬InExisting);
- 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);
- MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type):_type(type),_father(fath),_start(-1),_end(-1),_nval(-1),_loc_id(-5),_profile_it(-1) { }
- private:
- MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
- MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const std::string& dummy);
- MEDFileFieldPerMeshPerTypePerDisc(const MEDFileFieldPerMeshPerTypePerDisc& other);
- MEDFileFieldPerMeshPerTypePerDisc();
- private:
- 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);
- private:
- TypeOfField _type;
- MEDFileFieldPerMeshPerTypeCommon *_father;
- int _start;
- int _end;
- //! _nval is different than end-start in case of ON_GAUSS_PT and ON_GAUSS_NE ! (_nval=(_end-_start)/nbi)
- int _nval;
- std::string _profile;
- std::string _localization;
- //! only on assignement -3 : ON_NODES, -2 : ON_CELLS, -1 : ON_GAUSS_NE, 0..* : ON_GAUSS_PT
- mutable int _loc_id;
- mutable int _profile_it;
- MCAuto<PartDefinition> _pd;
- public:
- mutable int _tmp_work1;
- };
-
- class MEDFileFieldPerMeshPerTypeCommon : public RefCountObject, public MEDFileWritable
- {
- public:
- std::size_t getHeapMemorySizeWithoutChildren() const;
- std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldTemplate *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
- void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- const MEDFileFieldPerMesh *getFather() const;
- void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
- void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
- void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
- bool isUniqueLevel(int& dim) const;
- void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
- void fillFieldSplitedByType(std::vector< std::pair<int,int> >& dads, std::vector<TypeOfField>& types, std::vector<std::string>& pfls, std::vector<std::string>& locs) const;
- int getIteration() const;
- int getOrder() const;
- double getTime() const;
- std::string getMeshName() const;
- void getSizes(int& globalSz, int& nbOfEntries) const;
- int getNumberOfComponents() const;
- bool presenceOfMultiDiscPerGeoType() const;
- void pushDiscretization(MEDFileFieldPerMeshPerTypePerDisc *disc);
- DataArray *getOrCreateAndGetArray();
- const DataArray *getOrCreateAndGetArray() const;
- const std::vector<std::string>& getInfo() const;
- std::vector<std::string> getPflsReallyUsed() const;
- std::vector<std::string> getLocsReallyUsed() const;
- std::vector<std::string> getPflsReallyUsedMulti() const;
- std::vector<std::string> getLocsReallyUsedMulti() const;
- void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId);
- const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId) const;
- int getNumberOfLoc() const { return _field_pm_pt_pd.size(); }
- void fillValues(int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
- void setLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
- bool keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
- bool keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
- static med_entity_type ConvertIntoMEDFileType(TypeOfField ikType, INTERP_KERNEL::NormalizedCellType ikGeoType, med_geometry_type& medfGeoType);
- MEDFileFieldPerMeshPerTypeCommon(MEDFileFieldPerMesh *father):_father(father) { }
- void setFather(MEDFileFieldPerMesh *father);
- void accept(MEDFileFieldVisitor& visitor) const;
- public:
- virtual ~MEDFileFieldPerMeshPerTypeCommon();
- virtual void getDimension(int& dim) const = 0;
- virtual INTERP_KERNEL::NormalizedCellType getGeoType() const = 0;
- virtual void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const = 0;
- virtual void simpleRepr(int bkOffset, std::ostream& oss, int id) const = 0;
- virtual std::string getGeoTypeRepr() const = 0;
- virtual MEDFileFieldPerMeshPerTypeCommon *deepCopy(MEDFileFieldPerMesh *father) const = 0;
- 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;
- protected:
- void deepCopyElements();
- std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldTemplate *field, int offset, int nbOfCells);
- std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldTemplate *field, int offset, int nbOfCells);
- std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldTemplate *field, const DataArrayInt *subCells);
- std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldTemplate *field, const DataArrayInt *subCells);
- private:
- MEDFileFieldPerMesh *_father;
- protected:
- std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> > _field_pm_pt_pd;
- };
-
- class MEDFileFieldPerMeshPerType : public MEDFileFieldPerMeshPerTypeCommon
- {
- public:
- static MEDFileFieldPerMeshPerType *New(MEDFileFieldPerMesh *fath, INTERP_KERNEL::NormalizedCellType geoType);
- static MEDFileFieldPerMeshPerType *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
- 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);
- public:// overload of abstract methods
- void getDimension(int& dim) const;
- INTERP_KERNEL::NormalizedCellType getGeoType() const;
- void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
- void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
- std::string getGeoTypeRepr() const;
- MEDFileFieldPerMeshPerType *deepCopy(MEDFileFieldPerMesh *father) const;
- 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;
- private:
- MEDFileFieldPerMeshPerType(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
- MEDFileFieldPerMeshPerType(MEDFileFieldPerMesh *father, INTERP_KERNEL::NormalizedCellType gt);
- private:
- INTERP_KERNEL::NormalizedCellType _geo_type;
- };
-
- class MEDFileFieldPerMeshPerTypeDyn : public MEDFileFieldPerMeshPerTypeCommon
- {
- public:
- static MEDFileFieldPerMeshPerTypeDyn *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileEntities *entities, int idGT, const MEDFileFieldNameScope& nasc);
- int getDynGT() const;
- std::string getModelName() const;
- public:
- void getDimension(int& dim) const;
- INTERP_KERNEL::NormalizedCellType getGeoType() const;
- void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
- void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
- std::string getGeoTypeRepr() const;
- MEDFileFieldPerMeshPerTypeDyn *deepCopy(MEDFileFieldPerMesh *father) const;
- 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;
- private:
- MEDFileFieldPerMeshPerTypeDyn(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileStructureElement *se, const MEDFileFieldNameScope& nasc);
- private:
- MCConstAuto<MEDFileStructureElement> _se;
- };
-
- class MEDFileMesh;
-
- class MEDFileFieldPerMesh : public RefCountObject, public MEDFileWritable
- {
- public:
- static MEDFileFieldPerMesh *New(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
- static MEDFileFieldPerMesh *NewOnRead(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
- std::size_t getHeapMemorySizeWithoutChildren() const;
- std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- MEDFileFieldPerMesh *deepCopy(MEDFileAnyTypeField1TSWithoutSDA *father) const;
- void simpleRepr(int bkOffset,std::ostream& oss, int id) const;
- void copyTinyInfoFrom(const MEDCouplingMesh *mesh);
- 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 MEDCouplingFieldTemplate *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignFieldNoProfileNoRenum(int& start, const std::vector<int>& code, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
- void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
- void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
- void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
- void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
- void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
- 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;
- void accept(MEDFileFieldVisitor& visitor) const;
- void getDimension(int& dim) const;
- bool isUniqueLevel(int& dim) const;
- double getTime() const;
- int getIteration() const;
- int getOrder() const;
- int getMeshIteration() const { return _mesh_iteration; }
- int getMeshOrder() const { return _mesh_order; }
- std::string getMeshName() const;
- void setMeshName(const std::string& meshName);
- int getNumberOfComponents() const;
- bool presenceOfMultiDiscPerGeoType() const;
- bool presenceOfStructureElements() const;
- bool onlyStructureElements() const;
- void killStructureElements();
- void keepOnlyStructureElements();
- void keepOnlyOnSE(const std::string& seName);
- void getMeshSENames(std::vector< std::pair<std::string,std::string> >& ps) const;
- DataArray *getOrCreateAndGetArray();
- const DataArray *getOrCreateAndGetArray() const;
- const std::vector<std::string>& getInfo() const;
- std::vector<std::string> getPflsReallyUsed() const;
- std::vector<std::string> getLocsReallyUsed() const;
- std::vector<std::string> getPflsReallyUsedMulti() const;
- std::vector<std::string> getLocsReallyUsedMulti() const;
- void convertMedBallIntoClassic();
- bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
- bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N, MEDFileFieldGlobsReal& glob);
- void keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
- void keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
- void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, const MEDFileFieldGlobsReal *glob, const MEDCouplingMesh *mesh, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
- DataArray *getFieldOnMeshAtLevelWithPfl(TypeOfField type, const MEDCouplingMesh *mesh, DataArrayInt *&pfl, const MEDFileFieldGlobsReal *glob, const MEDFileFieldNameScope& nasc) const;
- void getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
- MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId);
- const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId) const;
- 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);
- private:
- int addNewEntryIfNecessary(INTERP_KERNEL::NormalizedCellType type);
- MEDCouplingFieldDouble *finishField(TypeOfField type, const MEDFileFieldGlobsReal *glob,
- 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;
- MEDCouplingFieldDouble *finishField2(TypeOfField type, const MEDFileFieldGlobsReal *glob,
- const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
- const std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes,
- const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
- MEDCouplingFieldDouble *finishFieldNode2(const MEDFileFieldGlobsReal *glob,
- const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
- const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
- DataArray *finishField4(const std::vector< std::pair<int,int> >& dads, const DataArrayInt *pflIn, int nbOfElems, DataArrayInt *&pflOut) const;
- void assignNewLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
- static void SortArraysPerType(const MEDFileFieldGlobsReal *glob, TypeOfField type,
- 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,
- std::vector<int>& code, std::vector<DataArrayInt *>& notNullPfls);
- 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);
- MEDFileFieldPerMesh(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
- MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
- MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const std::string& meshName, int meshIt, int meshOrd):_father(fath),_mesh_iteration(meshIt),_mesh_order(meshOrd) { }
- private:
- int _mesh_iteration;
- int _mesh_order;
- MEDFileAnyTypeField1TSWithoutSDA *_father;
- std::vector< MCAuto< MEDFileFieldPerMeshPerTypeCommon > > _field_pm_pt;
- };
-
- class MEDFileFieldGlobsReal;
-
- class MEDFileFieldGlobs : public RefCountObject
- {
- public:
- static MEDFileFieldGlobs *New(med_idt fid);
- static MEDFileFieldGlobs *New();
- std::size_t getHeapMemorySizeWithoutChildren() const;
- std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- MEDFileFieldGlobs *deepCopy() const;
- MEDFileFieldGlobs *shallowCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
- MEDFileFieldGlobs *deepCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
- void simpleRepr(std::ostream& oss) const;
- void appendGlobs(const MEDFileFieldGlobs& other, double eps);
- void checkGlobsPflsPartCoherency(const std::vector<std::string>& pflsUsed) const;
- void checkGlobsLocsPartCoherency(const std::vector<std::string>& locsUsed) const;
- void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
- void loadProfileInFile(med_idt fid, int id);
- void loadGlobals(med_idt fid, const MEDFileFieldGlobsReal& real);
- void loadAllGlobals(med_idt fid, const MEDFileEntities *entities);
- void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
- std::vector<std::string> getPfls() const;
- std::vector<std::string> getLocs() const;
- bool existsPfl(const std::string& pflName) const;
- bool existsLoc(const std::string& locName) const;
- std::string createNewNameOfPfl() const;
- std::string createNewNameOfLoc() const;
- std::vector< std::vector<int> > whichAreEqualProfiles() const;
- std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
- void setFileName(const std::string& fileName) { _file_name=fileName; }
- void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- int getNbOfGaussPtPerCell(int locId) const;
- int getLocalizationId(const std::string& loc) const;
- std::string getFileName() const { return _file_name; }
- const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
- const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
- const DataArrayInt *getProfileFromId(int pflId) const;
- const DataArrayInt *getProfile(const std::string& pflName) const;
- MEDFileFieldLoc& getLocalizationFromId(int locId);
- MEDFileFieldLoc& getLocalization(const std::string& locName);
- DataArrayInt *getProfile(const std::string& pflName);
- DataArrayInt *getProfileFromId(int pflId);
- void killProfileIds(const std::vector<int>& pflIds);
- void killLocalizationIds(const std::vector<int>& locIds);
- void killStructureElementsInGlobs();
- //
- void appendProfile(DataArrayInt *pfl);
- 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);
- //
- static std::string CreateNewNameNotIn(const std::string& prefix, const std::vector<std::string>& namesToAvoid);
- protected:
- MEDFileFieldGlobs(med_idt fid);
- MEDFileFieldGlobs();
- ~MEDFileFieldGlobs();
- protected:
- std::vector< MCAuto<DataArrayInt> > _pfls;
- std::vector< MCAuto<MEDFileFieldLoc> > _locs;
- std::string _file_name;
- };
-
- /// @endcond INTERNAL
-
- class MEDFileFieldGlobsReal
- {
- public:
- MEDLOADER_EXPORT MEDFileFieldGlobsReal(med_idt fid);
- MEDLOADER_EXPORT MEDFileFieldGlobsReal();
- MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
- MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
- MEDLOADER_EXPORT void simpleReprGlobs(std::ostream& oss) const;
- MEDLOADER_EXPORT void resetContent();
- MEDLOADER_EXPORT void killStructureElementsInGlobs();
- MEDLOADER_EXPORT void shallowCpyGlobs(const MEDFileFieldGlobsReal& other);
- MEDLOADER_EXPORT void deepCpyGlobs(const MEDFileFieldGlobsReal& other);
- MEDLOADER_EXPORT void shallowCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
- MEDLOADER_EXPORT void deepCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
- MEDLOADER_EXPORT void appendGlobs(const MEDFileFieldGlobsReal& other, double eps);
- MEDLOADER_EXPORT void checkGlobsCoherency() const;
- MEDLOADER_EXPORT void checkGlobsPflsPartCoherency() const;
- MEDLOADER_EXPORT void checkGlobsLocsPartCoherency() const;
- MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsed() const = 0;
- MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsed() const = 0;
- MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsedMulti() const = 0;
- MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsedMulti() const = 0;
- MEDLOADER_EXPORT virtual void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
- MEDLOADER_EXPORT virtual void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
- MEDLOADER_EXPORT virtual ~MEDFileFieldGlobsReal();
- //
- MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
- MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id);
- MEDLOADER_EXPORT void loadGlobals(med_idt fid);
- MEDLOADER_EXPORT void loadAllGlobals(med_idt fid, const MEDFileEntities *entities=0);
- MEDLOADER_EXPORT void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
- MEDLOADER_EXPORT std::vector<std::string> getPfls() const;
- MEDLOADER_EXPORT std::vector<std::string> getLocs() const;
- MEDLOADER_EXPORT bool existsPfl(const std::string& pflName) const;
- MEDLOADER_EXPORT bool existsLoc(const std::string& locName) const;
- MEDLOADER_EXPORT std::string createNewNameOfPfl() const;
- MEDLOADER_EXPORT std::string createNewNameOfLoc() const;
- MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualProfiles() const;
- MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
- MEDLOADER_EXPORT void setFileName(const std::string& fileName);
- MEDLOADER_EXPORT void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDLOADER_EXPORT void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDLOADER_EXPORT void changePflsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDLOADER_EXPORT void changeLocsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
- MEDLOADER_EXPORT void changePflName(const std::string& oldName, const std::string& newName);
- MEDLOADER_EXPORT void changeLocName(const std::string& oldName, const std::string& newName);
- MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipPflsNames();
- MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipLocsNames(double eps);
- MEDLOADER_EXPORT int getNbOfGaussPtPerCell(int locId) const;
- MEDLOADER_EXPORT int getLocalizationId(const std::string& loc) const;
- MEDLOADER_EXPORT std::string getFileName() const;
- MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
- MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
- MEDLOADER_EXPORT MEDFileFieldLoc& getLocalizationFromId(int locId);
- MEDLOADER_EXPORT MEDFileFieldLoc& getLocalization(const std::string& locName);
- MEDLOADER_EXPORT const DataArrayInt *getProfile(const std::string& pflName) const;
- MEDLOADER_EXPORT const DataArrayInt *getProfileFromId(int pflId) const;
- MEDLOADER_EXPORT DataArrayInt *getProfile(const std::string& pflName);
- MEDLOADER_EXPORT DataArrayInt *getProfileFromId(int pflId);
- MEDLOADER_EXPORT void killProfileIds(const std::vector<int>& pflIds);
- MEDLOADER_EXPORT void killLocalizationIds(const std::vector<int>& locIds);
- //
- MEDLOADER_EXPORT void appendProfile(DataArrayInt *pfl);
- 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);
- protected:
- MEDFileFieldGlobs *contentNotNull();
- const MEDFileFieldGlobs *contentNotNull() const;
- protected:
- MCAuto< MEDFileFieldGlobs > _globals;
- };
-
- class MEDFileFieldNameScope
- {
- public:
- MEDLOADER_EXPORT MEDFileFieldNameScope();
- MEDLOADER_EXPORT MEDFileFieldNameScope(const std::string& fieldName, const std::string& meshName);
- MEDLOADER_EXPORT std::string getName() const;
- MEDLOADER_EXPORT void setName(const std::string& fieldName);
- MEDLOADER_EXPORT std::string getDtUnit() const;
- MEDLOADER_EXPORT void setDtUnit(const std::string& dtUnit);
- MEDLOADER_EXPORT void copyNameScope(const MEDFileFieldNameScope& other);
- MEDLOADER_EXPORT std::string getMeshName() const;
- MEDLOADER_EXPORT void setMeshName(const std::string& meshName);
- protected:
- std::string _name;
- std::string _dt_unit;
- std::string _mesh_name;
- };
-
class MEDFileMeshes;
/*!
--- /dev/null
+// Copyright (C) 2007-2016 CEA/DEN, EDF R&D
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+// Author : Anthony Geay (EDF R&D)
+
+#ifndef __MEDFILEFIELDGLOBS_HXX__
+#define __MEDFILEFIELDGLOBS_HXX__
+
+#include "MEDLoaderDefines.hxx"
+
+#include "MCAuto.hxx"
+
+#include "med.h"
+
+namespace MEDCoupling
+{
+ class MEDFileFieldGlobsReal;
+
+ class MEDFileFieldGlobs : public RefCountObject
+ {
+ public:
+ static MEDFileFieldGlobs *New(med_idt fid);
+ static MEDFileFieldGlobs *New();
+ std::size_t getHeapMemorySizeWithoutChildren() const;
+ std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ MEDFileFieldGlobs *deepCopy() const;
+ MEDFileFieldGlobs *shallowCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
+ MEDFileFieldGlobs *deepCpyPart(const std::vector<std::string>& pfls, const std::vector<std::string>& locs) const;
+ void simpleRepr(std::ostream& oss) const;
+ void appendGlobs(const MEDFileFieldGlobs& other, double eps);
+ void checkGlobsPflsPartCoherency(const std::vector<std::string>& pflsUsed) const;
+ void checkGlobsLocsPartCoherency(const std::vector<std::string>& locsUsed) const;
+ void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
+ void loadProfileInFile(med_idt fid, int id);
+ void loadGlobals(med_idt fid, const MEDFileFieldGlobsReal& real);
+ void loadAllGlobals(med_idt fid, const MEDFileEntities *entities);
+ void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
+ std::vector<std::string> getPfls() const;
+ std::vector<std::string> getLocs() const;
+ bool existsPfl(const std::string& pflName) const;
+ bool existsLoc(const std::string& locName) const;
+ std::string createNewNameOfPfl() const;
+ std::string createNewNameOfLoc() const;
+ std::vector< std::vector<int> > whichAreEqualProfiles() const;
+ std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
+ void setFileName(const std::string& fileName) { _file_name=fileName; }
+ void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ int getNbOfGaussPtPerCell(int locId) const;
+ int getLocalizationId(const std::string& loc) const;
+ std::string getFileName() const { return _file_name; }
+ const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
+ const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
+ const DataArrayInt *getProfileFromId(int pflId) const;
+ const DataArrayInt *getProfile(const std::string& pflName) const;
+ MEDFileFieldLoc& getLocalizationFromId(int locId);
+ MEDFileFieldLoc& getLocalization(const std::string& locName);
+ DataArrayInt *getProfile(const std::string& pflName);
+ DataArrayInt *getProfileFromId(int pflId);
+ void killProfileIds(const std::vector<int>& pflIds);
+ void killLocalizationIds(const std::vector<int>& locIds);
+ void killStructureElementsInGlobs();
+ //
+ void appendProfile(DataArrayInt *pfl);
+ 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);
+ //
+ static std::string CreateNewNameNotIn(const std::string& prefix, const std::vector<std::string>& namesToAvoid);
+ protected:
+ MEDFileFieldGlobs(med_idt fid);
+ MEDFileFieldGlobs();
+ ~MEDFileFieldGlobs();
+ protected:
+ std::vector< MCAuto<DataArrayInt> > _pfls;
+ std::vector< MCAuto<MEDFileFieldLoc> > _locs;
+ std::string _file_name;
+ };
+
+ /// @endcond INTERNAL
+
+ class MEDFileFieldGlobsReal
+ {
+ public:
+ MEDLOADER_EXPORT MEDFileFieldGlobsReal(med_idt fid);
+ MEDLOADER_EXPORT MEDFileFieldGlobsReal();
+ MEDLOADER_EXPORT std::size_t getHeapMemorySizeWithoutChildren() const;
+ MEDLOADER_EXPORT std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ MEDLOADER_EXPORT void simpleReprGlobs(std::ostream& oss) const;
+ MEDLOADER_EXPORT void resetContent();
+ MEDLOADER_EXPORT void killStructureElementsInGlobs();
+ MEDLOADER_EXPORT void shallowCpyGlobs(const MEDFileFieldGlobsReal& other);
+ MEDLOADER_EXPORT void deepCpyGlobs(const MEDFileFieldGlobsReal& other);
+ MEDLOADER_EXPORT void shallowCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
+ MEDLOADER_EXPORT void deepCpyOnlyUsedGlobs(const MEDFileFieldGlobsReal& other);
+ MEDLOADER_EXPORT void appendGlobs(const MEDFileFieldGlobsReal& other, double eps);
+ MEDLOADER_EXPORT void checkGlobsCoherency() const;
+ MEDLOADER_EXPORT void checkGlobsPflsPartCoherency() const;
+ MEDLOADER_EXPORT void checkGlobsLocsPartCoherency() const;
+ MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsed() const = 0;
+ MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsed() const = 0;
+ MEDLOADER_EXPORT virtual std::vector<std::string> getPflsReallyUsedMulti() const = 0;
+ MEDLOADER_EXPORT virtual std::vector<std::string> getLocsReallyUsedMulti() const = 0;
+ MEDLOADER_EXPORT virtual void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
+ MEDLOADER_EXPORT virtual void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif) = 0;
+ MEDLOADER_EXPORT virtual ~MEDFileFieldGlobsReal();
+ //
+ MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id, const std::string& pflName);
+ MEDLOADER_EXPORT void loadProfileInFile(med_idt fid, int id);
+ MEDLOADER_EXPORT void loadGlobals(med_idt fid);
+ MEDLOADER_EXPORT void loadAllGlobals(med_idt fid, const MEDFileEntities *entities=0);
+ MEDLOADER_EXPORT void writeGlobals(med_idt fid, const MEDFileWritable& opt) const;
+ MEDLOADER_EXPORT std::vector<std::string> getPfls() const;
+ MEDLOADER_EXPORT std::vector<std::string> getLocs() const;
+ MEDLOADER_EXPORT bool existsPfl(const std::string& pflName) const;
+ MEDLOADER_EXPORT bool existsLoc(const std::string& locName) const;
+ MEDLOADER_EXPORT std::string createNewNameOfPfl() const;
+ MEDLOADER_EXPORT std::string createNewNameOfLoc() const;
+ MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualProfiles() const;
+ MEDLOADER_EXPORT std::vector< std::vector<int> > whichAreEqualLocs(double eps) const;
+ MEDLOADER_EXPORT void setFileName(const std::string& fileName);
+ MEDLOADER_EXPORT void changePflsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDLOADER_EXPORT void changeLocsNamesInStruct(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDLOADER_EXPORT void changePflsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDLOADER_EXPORT void changeLocsNames(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDLOADER_EXPORT void changePflName(const std::string& oldName, const std::string& newName);
+ MEDLOADER_EXPORT void changeLocName(const std::string& oldName, const std::string& newName);
+ MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipPflsNames();
+ MEDLOADER_EXPORT std::vector< std::pair<std::vector<std::string>, std::string > > zipLocsNames(double eps);
+ MEDLOADER_EXPORT int getNbOfGaussPtPerCell(int locId) const;
+ MEDLOADER_EXPORT int getLocalizationId(const std::string& loc) const;
+ MEDLOADER_EXPORT std::string getFileName() const;
+ MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalizationFromId(int locId) const;
+ MEDLOADER_EXPORT const MEDFileFieldLoc& getLocalization(const std::string& locName) const;
+ MEDLOADER_EXPORT MEDFileFieldLoc& getLocalizationFromId(int locId);
+ MEDLOADER_EXPORT MEDFileFieldLoc& getLocalization(const std::string& locName);
+ MEDLOADER_EXPORT const DataArrayInt *getProfile(const std::string& pflName) const;
+ MEDLOADER_EXPORT const DataArrayInt *getProfileFromId(int pflId) const;
+ MEDLOADER_EXPORT DataArrayInt *getProfile(const std::string& pflName);
+ MEDLOADER_EXPORT DataArrayInt *getProfileFromId(int pflId);
+ MEDLOADER_EXPORT void killProfileIds(const std::vector<int>& pflIds);
+ MEDLOADER_EXPORT void killLocalizationIds(const std::vector<int>& locIds);
+ //
+ MEDLOADER_EXPORT void appendProfile(DataArrayInt *pfl);
+ 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);
+ protected:
+ MEDFileFieldGlobs *contentNotNull();
+ const MEDFileFieldGlobs *contentNotNull() const;
+ protected:
+ MCAuto< MEDFileFieldGlobs > _globals;
+ };
+
+ class MEDFileFieldNameScope
+ {
+ public:
+ MEDLOADER_EXPORT MEDFileFieldNameScope();
+ MEDLOADER_EXPORT MEDFileFieldNameScope(const std::string& fieldName, const std::string& meshName);
+ MEDLOADER_EXPORT std::string getName() const;
+ MEDLOADER_EXPORT void setName(const std::string& fieldName);
+ MEDLOADER_EXPORT std::string getDtUnit() const;
+ MEDLOADER_EXPORT void setDtUnit(const std::string& dtUnit);
+ MEDLOADER_EXPORT void copyNameScope(const MEDFileFieldNameScope& other);
+ MEDLOADER_EXPORT std::string getMeshName() const;
+ MEDLOADER_EXPORT void setMeshName(const std::string& meshName);
+ protected:
+ std::string _name;
+ std::string _dt_unit;
+ std::string _mesh_name;
+ };
+}
+
+#endif
--- /dev/null
+// Copyright (C) 2007-2016 CEA/DEN, EDF R&D
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
+// Author : Anthony Geay (EDF R&D)
+
+#ifndef __MEDFILEFIELDINTERNAL_HXX__
+#define __MEDFILEFIELDINTERNAL_HXX__
+
+#include "MEDLoaderDefines.hxx"
+#include "MEDFileUtilities.hxx"
+#include "NormalizedGeometricTypes"
+#include "InterpKernelAutoPtr.hxx"
+#include "MCAuto.hxx"
+
+#include "med.h"
+
+#include <string>
+#include <list>
+
+namespace MEDCoupling
+{
+ class MEDFileGTKeeper
+ {
+ public:
+ virtual MEDFileGTKeeper *deepCopy() const = 0;
+ virtual INTERP_KERNEL::NormalizedCellType getGeoType() const = 0;
+ virtual std::string getRepr() const = 0;
+ virtual bool isEqual(const MEDFileGTKeeper *other) const = 0;
+ virtual ~MEDFileGTKeeper();
+ };
+
+ class MEDFileGTKeeperSta : public MEDFileGTKeeper
+ {
+ public:
+ MEDFileGTKeeperSta(INTERP_KERNEL::NormalizedCellType gt):_geo_type(gt) { }
+ MEDFileGTKeeper *deepCopy() const;
+ INTERP_KERNEL::NormalizedCellType getGeoType() const;
+ std::string getRepr() const;
+ bool isEqual(const MEDFileGTKeeper *other) const;
+ private:
+ INTERP_KERNEL::NormalizedCellType _geo_type;
+ };
+
+
+ class MEDFileStructureElement;
+ class MEDFileUMesh;
+
+ class MEDFileGTKeeperDyn : public MEDFileGTKeeper
+ {
+ public:
+ MEDFileGTKeeperDyn(const MEDFileUMesh *mesh, const MEDFileUMesh *section, const MEDFileStructureElement *se);
+ MEDFileGTKeeper *deepCopy() const;
+ INTERP_KERNEL::NormalizedCellType getGeoType() const;
+ std::string getRepr() const;
+ bool isEqual(const MEDFileGTKeeper *other) const;
+ const MEDFileUMesh *getMesh() const { return _mesh; }
+ const MEDFileUMesh *getSection() const { return _section; }
+ const MEDFileStructureElement *getSE() const { return _se; }
+ private:
+ MCConstAuto<MEDFileUMesh> _mesh;
+ MCConstAuto<MEDFileUMesh> _section;
+ MCConstAuto<MEDFileStructureElement> _se;
+ };
+
+ class MEDFileEntities;
+
+ class MEDFileFieldLoc : public RefCountObject
+ {
+ public:
+ MEDLOADER_EXPORT void simpleRepr(std::ostream& oss) const;
+ MEDLOADER_EXPORT std::string getName() const { return _name; }
+ MEDLOADER_EXPORT void setName(const std::string& name);
+ static MEDFileFieldLoc *New(med_idt fid, const std::string& locName);
+ static MEDFileFieldLoc *New(med_idt fid, int i, const MEDFileEntities *entities);
+ 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);
+ std::size_t getHeapMemorySizeWithoutChildren() const;
+ std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ MEDFileFieldLoc *deepCopy() const;
+ bool isOnStructureElement() const;
+ const MEDFileGTKeeper *getUndergroundGTKeeper() const { return _gt; }
+ MEDLOADER_EXPORT int getNbOfGaussPtPerCell() const { return _nb_gauss_pt; }
+ MEDLOADER_EXPORT void writeLL(med_idt fid) const;
+ MEDLOADER_EXPORT std::string repr() const;
+ MEDLOADER_EXPORT bool isName(const std::string& name) const { return _name==name; }
+ MEDLOADER_EXPORT int getDimension() const { return _dim; }
+ MEDLOADER_EXPORT int getNumberOfGaussPoints() const { return _nb_gauss_pt; }
+ MEDLOADER_EXPORT int getNumberOfPointsInCells() const { return _nb_node_per_cell; }
+ MEDLOADER_EXPORT const std::vector<double>& getRefCoords() const { return _ref_coo; }
+ MEDLOADER_EXPORT const std::vector<double>& getGaussCoords() const { return _gs_coo; }
+ MEDLOADER_EXPORT const std::vector<double>& getGaussWeights() const { return _w; }
+ MEDLOADER_EXPORT INTERP_KERNEL::NormalizedCellType getGeoType() const { return _gt->getGeoType(); }
+ MEDLOADER_EXPORT bool isEqual(const MEDFileFieldLoc& other, double eps) const;
+ private:
+ MEDFileFieldLoc(const MEDFileFieldLoc& other);
+ MEDFileFieldLoc(med_idt fid, const std::string& locName);
+ MEDFileFieldLoc(med_idt fid, int id, const MEDFileEntities *entities);
+ MEDFileFieldLoc(const std::string& locName, INTERP_KERNEL::NormalizedCellType geoType, const std::vector<double>& refCoo, const std::vector<double>& gsCoo, const std::vector<double>& w);
+ private:
+ int _dim;
+ int _nb_gauss_pt;
+ INTERP_KERNEL::AutoCppPtr<MEDFileGTKeeper> _gt;
+ int _nb_node_per_cell;
+ std::string _name;
+ std::vector<double> _ref_coo;
+ std::vector<double> _gs_coo;
+ std::vector<double> _w;
+ };
+
+ /// @cond INTERNAL
+ class MEDFileFieldPerMeshPerTypeCommon;
+ class MEDFileFieldPerMeshPerType;
+ class MEDCouplingFieldTemplate;
+ class MEDFileFieldNameScope;
+ class MEDFileFieldGlobsReal;
+ class MEDCouplingMesh;
+
+ class MEDFileFieldPerMeshPerTypePerDisc : public RefCountObject, public MEDFileWritable
+ {
+ public:
+ static MEDFileFieldPerMeshPerTypePerDisc *NewOnRead(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
+ static MEDFileFieldPerMeshPerTypePerDisc *New(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int locId);
+ static MEDFileFieldPerMeshPerTypePerDisc *New(const MEDFileFieldPerMeshPerTypePerDisc& other);
+ std::size_t getHeapMemorySizeWithoutChildren() const;
+ std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ MEDFileFieldPerMeshPerTypePerDisc *deepCopy(MEDFileFieldPerMeshPerTypeCommon *father) const;
+ void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldTemplate *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldTemplate *field, const DataArray *arrr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arrr, MEDFileFieldGlobsReal& glob);
+ void getCoarseData(TypeOfField& type, std::pair<int,int>& dad, std::string& pfl, std::string& loc) const;
+ void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
+ const MEDFileFieldPerMeshPerTypeCommon *getFather() const;
+ void loadOnlyStructureOfDataRecursively(med_idt fid, int& start, const MEDFileFieldNameScope& nasc);
+ void loadBigArray(med_idt fid, const MEDFileFieldNameScope& nasc);
+ void setNewStart(int newValueOfStart);
+ int getIteration() const;
+ int getOrder() const;
+ double getTime() const;
+ std::string getMeshName() const;
+ TypeOfField getType() const;
+ void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
+ void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
+ void setType(TypeOfField newType);
+ INTERP_KERNEL::NormalizedCellType getGeoType() const;
+ int getNumberOfComponents() const;
+ int getNumberOfTuples() const;
+ int getStart() const { return _start; }
+ int getEnd() const { return _end; }
+ void setEnd(int endd) { _end=endd; }
+ int getNumberOfVals() const { return _nval; }
+ DataArray *getOrCreateAndGetArray();
+ const DataArray *getOrCreateAndGetArray() const;
+ const std::vector<std::string>& getInfo() const;
+ std::string getProfile() const;
+ void setProfile(const std::string& newPflName);
+ std::string getLocalization() const;
+ void setLocalization(const std::string& newLocName);
+ int getLocId() const { return _loc_id; }
+ void setLocId(int newId) const { _loc_id=newId; }
+ void setFather(MEDFileFieldPerMeshPerTypeCommon *newFather) { _father=newFather; }
+ void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ void getFieldAtLevel(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;
+ void fillValues(int discId, int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
+ int fillEltIdsFromCode(int offset, const std::vector<int>& codeOfMesh, const MEDFileFieldGlobsReal& glob, int *ptToFill) const;
+ int fillTupleIds(int *ptToFill) const;
+ static int ConvertType(TypeOfField type, int locId);
+ static std::vector< std::vector< const MEDFileFieldPerMeshPerTypePerDisc *> > SplitPerDiscretization(const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entries);
+ static bool RenumberChunks(int offset, const std::vector< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
+ const DataArrayInt *explicitIdsInMesh, const std::vector<int>& newCode,
+ MEDFileFieldGlobsReal& glob, DataArrayDouble *arr, std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> >& result);
+ static MEDFileFieldPerMeshPerTypePerDisc *NewObjectOnSameDiscThanPool(TypeOfField typeF, INTERP_KERNEL::NormalizedCellType geoType, DataArrayInt *idsOfMeshElt,
+ bool isPfl, int nbi, int offset, std::list< const MEDFileFieldPerMeshPerTypePerDisc *>& entriesOnSameDisc,
+ MEDFileFieldGlobsReal& glob, bool ¬InExisting);
+ 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);
+ MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type):_type(type),_father(fath),_start(-1),_end(-1),_nval(-1),_loc_id(-5),_profile_it(-1) { }
+ private:
+ MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const PartDefinition *pd);
+ MEDFileFieldPerMeshPerTypePerDisc(MEDFileFieldPerMeshPerTypeCommon *fath, TypeOfField type, int profileIt, const std::string& dummy);
+ MEDFileFieldPerMeshPerTypePerDisc(const MEDFileFieldPerMeshPerTypePerDisc& other);
+ MEDFileFieldPerMeshPerTypePerDisc();
+ private:
+ 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);
+ private:
+ TypeOfField _type;
+ MEDFileFieldPerMeshPerTypeCommon *_father;
+ int _start;
+ int _end;
+ //! _nval is different than end-start in case of ON_GAUSS_PT and ON_GAUSS_NE ! (_nval=(_end-_start)/nbi)
+ int _nval;
+ std::string _profile;
+ std::string _localization;
+ //! only on assignement -3 : ON_NODES, -2 : ON_CELLS, -1 : ON_GAUSS_NE, 0..* : ON_GAUSS_PT
+ mutable int _loc_id;
+ mutable int _profile_it;
+ MCAuto<PartDefinition> _pd;
+ public:
+ mutable int _tmp_work1;
+ };
+
+ class MEDFileFieldVisitor;
+ class MEDFileFieldPerMesh;
+
+ class MEDFileFieldPerMeshPerTypeCommon : public RefCountObject, public MEDFileWritable
+ {
+ public:
+ std::size_t getHeapMemorySizeWithoutChildren() const;
+ std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ void assignFieldNoProfile(int& start, int offset, int nbOfCells, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignFieldProfile(bool isPflAlone, int& start, const DataArrayInt *multiTypePfl, const DataArrayInt *idsInPfl, DataArrayInt *locIds, int nbOfEltsInWholeMesh, const MEDCouplingFieldTemplate *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
+ void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ const MEDFileFieldPerMesh *getFather() const;
+ void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
+ void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
+ void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
+ bool isUniqueLevel(int& dim) const;
+ void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
+ void fillFieldSplitedByType(std::vector< std::pair<int,int> >& dads, std::vector<TypeOfField>& types, std::vector<std::string>& pfls, std::vector<std::string>& locs) const;
+ int getIteration() const;
+ int getOrder() const;
+ double getTime() const;
+ std::string getMeshName() const;
+ void getSizes(int& globalSz, int& nbOfEntries) const;
+ int getNumberOfComponents() const;
+ bool presenceOfMultiDiscPerGeoType() const;
+ void pushDiscretization(MEDFileFieldPerMeshPerTypePerDisc *disc);
+ DataArray *getOrCreateAndGetArray();
+ const DataArray *getOrCreateAndGetArray() const;
+ const std::vector<std::string>& getInfo() const;
+ std::vector<std::string> getPflsReallyUsed() const;
+ std::vector<std::string> getLocsReallyUsed() const;
+ std::vector<std::string> getPflsReallyUsedMulti() const;
+ std::vector<std::string> getLocsReallyUsedMulti() const;
+ void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId);
+ const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenLocId(int locId) const;
+ int getNumberOfLoc() const { return _field_pm_pt_pd.size(); }
+ void fillValues(int& startEntryId, std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
+ void setLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
+ bool keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
+ bool keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
+ static med_entity_type ConvertIntoMEDFileType(TypeOfField ikType, INTERP_KERNEL::NormalizedCellType ikGeoType, med_geometry_type& medfGeoType);
+ MEDFileFieldPerMeshPerTypeCommon(MEDFileFieldPerMesh *father):_father(father) { }
+ void setFather(MEDFileFieldPerMesh *father);
+ void accept(MEDFileFieldVisitor& visitor) const;
+ public:
+ virtual ~MEDFileFieldPerMeshPerTypeCommon();
+ virtual void getDimension(int& dim) const = 0;
+ virtual INTERP_KERNEL::NormalizedCellType getGeoType() const = 0;
+ virtual void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const = 0;
+ virtual void simpleRepr(int bkOffset, std::ostream& oss, int id) const = 0;
+ virtual std::string getGeoTypeRepr() const = 0;
+ virtual MEDFileFieldPerMeshPerTypeCommon *deepCopy(MEDFileFieldPerMesh *father) const = 0;
+ 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;
+ protected:
+ void deepCopyElements();
+ std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldTemplate *field, int offset, int nbOfCells);
+ std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldTemplate *field, int offset, int nbOfCells);
+ std::vector<int> addNewEntryIfNecessary(const MEDCouplingFieldTemplate *field, const DataArrayInt *subCells);
+ std::vector<int> addNewEntryIfNecessaryGauss(const MEDCouplingFieldTemplate *field, const DataArrayInt *subCells);
+ private:
+ MEDFileFieldPerMesh *_father;
+ protected:
+ std::vector< MCAuto<MEDFileFieldPerMeshPerTypePerDisc> > _field_pm_pt_pd;
+ };
+
+ class MEDFileFieldPerMeshPerType : public MEDFileFieldPerMeshPerTypeCommon
+ {
+ public:
+ static MEDFileFieldPerMeshPerType *New(MEDFileFieldPerMesh *fath, INTERP_KERNEL::NormalizedCellType geoType);
+ static MEDFileFieldPerMeshPerType *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
+ 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);
+ public:// overload of abstract methods
+ void getDimension(int& dim) const;
+ INTERP_KERNEL::NormalizedCellType getGeoType() const;
+ void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
+ void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
+ std::string getGeoTypeRepr() const;
+ MEDFileFieldPerMeshPerType *deepCopy(MEDFileFieldPerMesh *father) const;
+ 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;
+ private:
+ MEDFileFieldPerMeshPerType(med_idt fid, MEDFileFieldPerMesh *fath, TypeOfField type, INTERP_KERNEL::NormalizedCellType geoType, const MEDFileFieldNameScope& nasc, const PartDefinition *pd);
+ MEDFileFieldPerMeshPerType(MEDFileFieldPerMesh *father, INTERP_KERNEL::NormalizedCellType gt);
+ private:
+ INTERP_KERNEL::NormalizedCellType _geo_type;
+ };
+
+ class MEDFileFieldPerMeshPerTypeDyn : public MEDFileFieldPerMeshPerTypeCommon
+ {
+ public:
+ static MEDFileFieldPerMeshPerTypeDyn *NewOnRead(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileEntities *entities, int idGT, const MEDFileFieldNameScope& nasc);
+ int getDynGT() const;
+ std::string getModelName() const;
+ public:
+ void getDimension(int& dim) const;
+ INTERP_KERNEL::NormalizedCellType getGeoType() const;
+ void entriesForMEDfile(TypeOfField mct, med_geometry_type& gt, med_entity_type& ent) const;
+ void simpleRepr(int bkOffset, std::ostream& oss, int id) const;
+ std::string getGeoTypeRepr() const;
+ MEDFileFieldPerMeshPerTypeDyn *deepCopy(MEDFileFieldPerMesh *father) const;
+ 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;
+ private:
+ MEDFileFieldPerMeshPerTypeDyn(med_idt fid, MEDFileFieldPerMesh *fath, const MEDFileStructureElement *se, const MEDFileFieldNameScope& nasc);
+ private:
+ MCConstAuto<MEDFileStructureElement> _se;
+ };
+
+ class MEDFileMesh;
+ class MEDFileAnyTypeField1TSWithoutSDA;
+ class MEDFileField1TSWithoutSDA;
+
+ class MEDFileFieldPerMesh : public RefCountObject, public MEDFileWritable
+ {
+ public:
+ static MEDFileFieldPerMesh *New(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
+ static MEDFileFieldPerMesh *NewOnRead(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
+ std::size_t getHeapMemorySizeWithoutChildren() const;
+ std::vector<const BigMemoryObject *> getDirectChildrenWithNull() const;
+ MEDFileFieldPerMesh *deepCopy(MEDFileAnyTypeField1TSWithoutSDA *father) const;
+ void simpleRepr(int bkOffset,std::ostream& oss, int id) const;
+ void copyTinyInfoFrom(const MEDCouplingMesh *mesh);
+ 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 MEDCouplingFieldTemplate *field, const DataArray *arr, const MEDCouplingMesh *mesh, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignFieldNoProfileNoRenum(int& start, const std::vector<int>& code, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void assignNodeFieldNoProfile(int& start, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob);
+ void assignNodeFieldProfile(int& start, const DataArrayInt *pfl, const MEDCouplingFieldTemplate *field, const DataArray *arr, MEDFileFieldGlobsReal& glob, const MEDFileFieldNameScope& nasc);
+ void loadOnlyStructureOfDataRecursively(med_idt fid, int &start, const MEDFileFieldNameScope& nasc);
+ void loadBigArraysRecursively(med_idt fid, const MEDFileFieldNameScope& nasc);
+ void writeLL(med_idt fid, const MEDFileFieldNameScope& nasc) const;
+ void fillTypesOfFieldAvailable(std::set<TypeOfField>& types) const;
+ 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;
+ void accept(MEDFileFieldVisitor& visitor) const;
+ void getDimension(int& dim) const;
+ bool isUniqueLevel(int& dim) const;
+ double getTime() const;
+ int getIteration() const;
+ int getOrder() const;
+ int getMeshIteration() const { return _mesh_iteration; }
+ int getMeshOrder() const { return _mesh_order; }
+ std::string getMeshName() const;
+ void setMeshName(const std::string& meshName);
+ int getNumberOfComponents() const;
+ bool presenceOfMultiDiscPerGeoType() const;
+ bool presenceOfStructureElements() const;
+ bool onlyStructureElements() const;
+ void killStructureElements();
+ void keepOnlyStructureElements();
+ void keepOnlyOnSE(const std::string& seName);
+ void getMeshSENames(std::vector< std::pair<std::string,std::string> >& ps) const;
+ DataArray *getOrCreateAndGetArray();
+ const DataArray *getOrCreateAndGetArray() const;
+ const std::vector<std::string>& getInfo() const;
+ std::vector<std::string> getPflsReallyUsed() const;
+ std::vector<std::string> getLocsReallyUsed() const;
+ std::vector<std::string> getPflsReallyUsedMulti() const;
+ std::vector<std::string> getLocsReallyUsedMulti() const;
+ void convertMedBallIntoClassic();
+ bool changeMeshNames(const std::vector< std::pair<std::string,std::string> >& modifTab);
+ bool renumberEntitiesLyingOnMesh(const std::string& meshName, const std::vector<int>& oldCode, const std::vector<int>& newCode, const DataArrayInt *renumO2N, MEDFileFieldGlobsReal& glob);
+ void keepOnlySpatialDiscretization(TypeOfField tof, int &globalNum, std::vector< std::pair<int,int> >& its);
+ void keepOnlyGaussDiscretization(std::size_t idOfDisc, int &globalNum, std::vector< std::pair<int,int> >& its);
+ void changePflsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ void changeLocsRefsNamesGen(const std::vector< std::pair<std::vector<std::string>, std::string > >& mapOfModif);
+ MEDCouplingFieldDouble *getFieldOnMeshAtLevel(TypeOfField type, const MEDFileFieldGlobsReal *glob, const MEDCouplingMesh *mesh, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
+ DataArray *getFieldOnMeshAtLevelWithPfl(TypeOfField type, const MEDCouplingMesh *mesh, DataArrayInt *&pfl, const MEDFileFieldGlobsReal *glob, const MEDFileFieldNameScope& nasc) const;
+ void getUndergroundDataArrayExt(std::vector< std::pair<std::pair<INTERP_KERNEL::NormalizedCellType,int>,std::pair<int,int> > >& entries) const;
+ MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId);
+ const MEDFileFieldPerMeshPerTypePerDisc *getLeafGivenTypeAndLocId(INTERP_KERNEL::NormalizedCellType typ, int locId) const;
+ 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);
+ private:
+ int addNewEntryIfNecessary(INTERP_KERNEL::NormalizedCellType type);
+ MEDCouplingFieldDouble *finishField(TypeOfField type, const MEDFileFieldGlobsReal *glob,
+ 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;
+ MEDCouplingFieldDouble *finishField2(TypeOfField type, const MEDFileFieldGlobsReal *glob,
+ const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
+ const std::vector<INTERP_KERNEL::NormalizedCellType>& geoTypes,
+ const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
+ MEDCouplingFieldDouble *finishFieldNode2(const MEDFileFieldGlobsReal *glob,
+ const std::vector< std::pair<int,int> >& dads, const std::vector<int>& locs,
+ const MEDCouplingMesh *mesh, const DataArrayInt *da, bool& isPfl, MCAuto<DataArray> &arrOut, const MEDFileFieldNameScope& nasc) const;
+ DataArray *finishField4(const std::vector< std::pair<int,int> >& dads, const DataArrayInt *pflIn, int nbOfElems, DataArrayInt *&pflOut) const;
+ void assignNewLeaves(const std::vector< MCAuto< MEDFileFieldPerMeshPerTypePerDisc > >& leaves);
+ static void SortArraysPerType(const MEDFileFieldGlobsReal *glob, TypeOfField type,
+ 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,
+ std::vector<int>& code, std::vector<DataArrayInt *>& notNullPfls);
+ 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);
+ MEDFileFieldPerMesh(med_idt fid, MEDFileAnyTypeField1TSWithoutSDA *fath, int meshCsit, int meshIteration, int meshOrder, const MEDFileFieldNameScope& nasc, const MEDFileMesh *mm, const MEDFileEntities *entities);
+ MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const MEDCouplingMesh *mesh);
+ MEDFileFieldPerMesh(MEDFileAnyTypeField1TSWithoutSDA *fath, const std::string& meshName, int meshIt, int meshOrd):_father(fath),_mesh_iteration(meshIt),_mesh_order(meshOrd) { }
+ private:
+ int _mesh_iteration;
+ int _mesh_order;
+ MEDFileAnyTypeField1TSWithoutSDA *_father;
+ std::vector< MCAuto< MEDFileFieldPerMeshPerTypeCommon > > _field_pm_pt;
+ };
+}
+
+#endif
