src/MULTIPR/MULTIPR_Obj.hxx

   1 // Project MULTIPR, IOLS WP1.2.1 - EDF/CS
   2 // Partitioning/decimation module for the SALOME v3.2 platform
   3
   4 /**
   5  * \file    MULTIPR_GaussLoc.hxx
   6  *
   7  * \brief   Class MULTIPR_Obj.
   8  *          This class is used as an interface to implement MULTIPR services in the salome MODULE
   9  *          as described in the spec/design doc.
  10  *
  11  * \author  Olivier LE ROUX - CS, Virtual Reality Dpt
  12  *
  13  * \date    01/2007
  14  */
  15
  16 #ifndef MULTIPR_OBJ_HXX
  17 #define MULTIPR_OBJ_HXX
  18
  19 //*****************************************************************************
  20 // Includes section
  21 //*****************************************************************************
  22
  23 #include <iostream>
  24 #include <string>
  25 #include <vector>
  26
  27
  28 namespace multipr
  29 {
  30
  31 class MeshDis;
  32
  33
  34 //*****************************************************************************
  35 // Class Obj
  36 //*****************************************************************************
  37
  38 enum ObjState
  39 {
  40         MULTIPR_OBJ_STATE_ERROR,
  41         MULTIPR_OBJ_STATE_RESET,
  42         MULTIPR_OBJ_STATE_SEQ_INIT,
  43         MULTIPR_OBJ_STATE_SEQ,
  44         MULTIPR_OBJ_STATE_DIS,
  45         MULTIPR_OBJ_STATE_DIS_MEM
  46 }; // enum ObjState
  47
  48
  49 class Obj
  50 {
  51 public:
  52
  53         /**
  54          * Builds an empty Gauss reference (default constructor).
  55          */
  56         Obj();
  57
  58         /**
  59          * Destructor. Removes everything.
  60          */
  61         ~Obj();
  62
  63         /**
  64          * Resets this object in its state by default (empty). Cleans memory.
  65          */
  66         void reset();
  67
  68         /**
  69          * Associate a MED file (sequential or distributed) with this object.
  70          * This method also:
  71          * - reset everything before starting
  72          * - determine if the given file is a sequential (SEQ) or a distributed (DIS) MED file
  73          * - read the ASCII master file if it is a distributed MED file
  74          * - set state
  75          */
  76         void create(const char* pMEDfilename);
  77
  78         //---------------------------------------------------------------------
  79         // Basic accessors/mutators
  80         //--------------------------------------------------------------------
  81
  82         /**
  83          * Returns true iff this obj represents a valid sequential MED file.
  84          * \return true iff this obj represents a valid sequential MED file.
  85          */
  86         bool isValidSequentialMEDFile() const { return (mState == MULTIPR_OBJ_STATE_SEQ) || (mState == MULTIPR_OBJ_STATE_SEQ_INIT); }
  87
  88         /**
  89          * Returns true iff this obj represents a valid distributed MED file.
  90          * \return true iff this obj represents a valid distributed  MED file.
  91          */
  92         bool isValidDistributedMEDFile() const { return (mState == MULTIPR_OBJ_STATE_DIS) || (mState == MULTIPR_OBJ_STATE_DIS_MEM); }
  93
  94         /**
  95          * Defines the mesh to be processed.
  96          * \param  pMeshName name of the mesh to be partitionned.
  97          */
  98         void setMesh(const char* pMeshName);
  99
 100         /**
 101          * Returns the list of meshes contained in the sequential MED file.
 102          * Assumes this object encapsulates a sequential MED file.
 103          * \return the list of meshes contained in the sequential MED file.
 104          */
 105         std::vector<std::string> getMeshes() const;
 106
 107         /**
 108          * Returns the list of fields contained in the sequential MED file.
 109          * Assumes this object encapsulates a sequential MED file.
 110          * \return the list of fields contained in the sequential MED file.
 111          */
 112         std::vector<std::string> getFields() const;
 113
 114         /**
 115          * Returns the number of timestamps for a given field.
 116          * Assumes this object encapsulates a sequential MED file.
 117          * \param  pFieldName name of any field.
 118          * \return the number of timestamps for a given field; 0 if field not found.
 119          */
 120         int getTimeStamps(const char* pFieldName) const;
 121
 122         /**
 123          * Returns the name of all partitions.
 124          * Assumes this object encapsulates a distributed MED file.
 125          * \return the name of all partitions.
 126          */
 127         std::vector<std::string> getParts() const;
 128
 129         /**
 130          * Returns all information about a part.
 131          * \param  pPartName name of the part.
 132          * \return all information about a part.
 133          */
 134         std::string getPartInfo(const char* pPartName) const;
 135
 136         /**
 137          * Removes all the part beginning by pPrefixPartName from the distributed MED file.
 138          * Assume this object encapsulates a distributed MED file.
 139          * \param  pPrefixPartName name of the part.
 140          */
 141         void removeParts(const char* pPrefixPartName);
 142
 143
 144         //---------------------------------------------------------------------
 145         // Algorithms
 146         //--------------------------------------------------------------------
 147
 148         /**
 149          * Creates a distributed MED file (v2.3) by extracting all the groups from the current mesh of the current MED sequential MED file.
 150          *         Assumes:
 151          *         - the file is in MED format and can be read using MED file v2.3.
 152          *         - the file is sequential (not a distributed MED).
 153          *         - the file only contains TETRA10 elements (dimension of space and mesh is 3).
 154          *         - the file have no profil.
 155          * \return the name of each part.
 156          * \throw  RuntimeException if any error occurs.
 157          */
 158         std::vector<std::string> partitionneDomaine();
 159
 160         /**
 161          * Creates a distributed MED file (V2.3) by splitting a group of a MED file previously created by partitionneDomaine.
 162          *         Assumes:
 163          *         - the file is a distributed MED file, previously created by partitionneDomaine()
 164          *           (=> each part only contain 1 mesh, TETRA10 elements only)
 165          *         - nbPart > 1
 166          * \param  pPartName     name of the part to be splitted.
 167          * \param  pNbParts      number of parts; must be > 1.
 168          * \param  pPartitionner use value 0=MULTIPR_METIS for Metis or 1=MULTIPR_SCOTCH for Scotch.
 169          * \return the name of each part.
 170          * \throw  RuntimeException if any error occurs.
 171          */
 172         std::vector<std::string> partitionneGrain(
 173                 const char* pPartName,
 174                 int         pNbParts,
 175                 int         pPartitionner=0);
 176
 177         /**
 178          * Creates 3 resolutions of the given part of a distributed MED file (V2.3).
 179          *         Assumes:
 180          *         - the file is a distributed MED file, previously created by partitionneDomaine() or partitionneGrain()
 181          *           (=> each part only contain 1 mesh, TETRA10 elements only)
 182          * \param  pPartName    name of the part to be decimated.
 183          * \param  pFieldName   name of the field used for decimation.
 184          * \param  pFieldIt     iteration (time step) of the field.
 185          * \param  pFilterName  name of the filter to be used.
 186          * \param  pTmed        threshold used for medium resolution.
 187          * \param  pTlow        threshold used for low resolution; tmed must be less than tlow
 188          * \param  pTadius      radius used to determine the neighbourhood.
 189          * \param  pBoxing      number of cells along each axis; must be >= 1; e.g. if 100 then acceleration grid will have 100*100*100 = 10**6 cells.
 190          * \return the name of each part.
 191          * \throw  RuntimeException if any error occurs.
 192          */
 193         std::vector<std::string> decimePartition(
 194                 const char* pPartName,
 195                 const char* pFieldName,
 196                 int         pFieldIt,
 197                 const char* pFilterName,
 198                 double      pTmed,
 199                 double      pTlow,
 200                 double      pRadius,
 201                 int         pBoxing);
 202
 203         //---------------------------------------------------------------------
 204         // I/O
 205         //---------------------------------------------------------------------
 206
 207         /**
 208          * Saves the associated MED file if necessary.
 209          * \throw  IOException if any i/o error occurs.
 210          */
 211         void save();
 212
 213         /**
 214          * Dumps any Obj to the given output stream.
 215          * \param  pOs any output stream.
 216          * \param  pO  any Obj.
 217          * \return the output stream pOs.
 218          */
 219         friend std::ostream& operator<<(std::ostream& pOs, Obj& pO);
 220
 221 private:
 222
 223         /**
 224          * Returns the list of parts contained in the current distributed mesh.
 225          * \return the list of parts contained in the current distributed mesh.
 226          * \throw  IllegalStateException if the distributed mesh does not still exist.
 227          */
 228         std::vector<std::string> getListParts() const;
 229
 230 private:
 231
 232         std::string  mMEDfilename;     /**< Name of the MED file: sequential or distributed. */
 233         std::string  mMeshName;        /**< Mesh to be partitionned. */
 234         ObjState     mState;           /**< State of this object. */
 235         MeshDis*     mMeshDis;         /**< Distributed mesh. */
 236
 237
 238 private:
 239
 240         // do not allow copy constructor
 241         Obj(const Obj&);
 242
 243         // do not allow copy
 244         Obj& operator=(const Obj&);
 245
 246         // do not allow operator ==
 247         bool operator==(const Obj&);
 248
 249 }; // class Obj
 250
 251
 252 } // namespace MULTIPR
 253
 254
 255 #endif // MULTIPR_OBJ_HXX
 256
 257 // EOF