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         // Algorithms
 138         //--------------------------------------------------------------------
 139
 140         /**
 141          * Creates a distributed MED file (v2.3) by extracting all the groups from the current mesh of the current MED sequential MED file.
 142          *         Assumes:
 143          *         - the file is in MED format and can be read using MED file v2.3.
 144          *         - the file is sequential (not a distributed MED).
 145          *         - the file only contains TETRA10 elements (dimension of space and mesh is 3).
 146          *         - the file have no profil.
 147          * \return the name of each part.
 148          * \throw  RuntimeException if any error occurs.
 149          */
 150         std::vector<std::string> partitionneDomaine();
 151
 152         /**
 153          * Creates a distributed MED file (V2.3) by splitting a group of a MED file previously created by partitionneDomaine.
 154          *         Assumes:
 155          *         - the file is a distributed MED file, previously created by partitionneDomaine()
 156          *           (=> each part only contain 1 mesh, TETRA10 elements only)
 157          *         - nbPart > 1
 158          * \param  pPartName     name of the part to be splitted.
 159          * \param  pNbParts      number of parts; must be > 1.
 160          * \param  pPartitionner use value 0=MULTIPR_METIS for Metis or 1=MULTIPR_SCOTCH for Scotch.
 161          * \return the name of each part.
 162          * \throw  RuntimeException if any error occurs.
 163          */
 164         std::vector<std::string> partitionneGrain(
 165                 const char* pPartName,
 166                 int         pNbParts,
 167                 int         pPartitionner=0);
 168
 169         /**
 170          * Creates 3 resolutions of the given part of a distributed MED file (V2.3).
 171          *         Assumes:
 172          *         - the file is a distributed MED file, previously created by partitionneDomaine() or partitionneGrain()
 173          *           (=> each part only contain 1 mesh, TETRA10 elements only)
 174          * \param  pPartName    name of the part to be decimated.
 175          * \param  pFieldName   name of the field used for decimation.
 176          * \param  pFieldIt     iteration (time step) of the field.
 177          * \param  pFilterName  name of the filter to be used.
 178          * \param  pTmed        threshold used for medium resolution.
 179          * \param  pTlow        threshold used for low resolution; tmed must be less than tlow
 180          * \param  pTadius      radius used to determine the neighbourhood.
 181          * \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.
 182          * \return the name of each part.
 183          * \throw  RuntimeException if any error occurs.
 184          */
 185         std::vector<std::string> decimePartition(
 186                 const char* pPartName,
 187                 const char* pFieldName,
 188                 int         pFieldIt,
 189                 const char* pFilterName,
 190                 double      pTmed,
 191                 double      pTlow,
 192                 double      pRadius,
 193                 int         pBoxing);
 194
 195         //---------------------------------------------------------------------
 196         // I/O
 197         //---------------------------------------------------------------------
 198
 199         /**
 200          * Saves the associated MED file if necessary.
 201          * \throw  IOException if any i/o error occurs.
 202          */
 203         void save();
 204
 205         /**
 206          * Dumps any Obj to the given output stream.
 207          * \param  pOs any output stream.
 208          * \param  pO  any Obj.
 209          * \return the output stream pOs.
 210          */
 211         friend std::ostream& operator<<(std::ostream& pOs, Obj& pO);
 212
 213 private:
 214
 215         /**
 216          * Returns the list of parts contained in the current distributed mesh.
 217          * \return the list of parts contained in the current distributed mesh.
 218          * \throw  IllegalStateException if the distributed mesh does not still exist.
 219          */
 220         std::vector<std::string> getListParts() const;
 221
 222 private:
 223
 224         std::string  mMEDfilename;     /**< Name of the MED file: sequential or distributed. */
 225         std::string  mMeshName;        /**< Mesh to be partitionned. */
 226         ObjState     mState;           /**< State of this object. */
 227         MeshDis*     mMeshDis;         /**< Distributed mesh. */
 228
 229
 230 private:
 231
 232         // do not allow copy constructor
 233         Obj(const Obj&);
 234
 235         // do not allow copy
 236         Obj& operator=(const Obj&);
 237
 238         // do not allow operator ==
 239         bool operator==(const Obj&);
 240
 241 }; // class Obj
 242
 243
 244 } // namespace MULTIPR
 245
 246
 247 #endif // MULTIPR_OBJ_HXX
 248
 249 // EOF