ReadIDEAS.hpp

Go to the documentation of this file.

    1 #ifndef READIDEAS_HPP
    2 #define READIDEAS_HPP
    3  
    4 #ifndef IS_BUILDING_MB
    5 #error "ReadIDEAS.hpp isn't supposed to be included into an application"
    6 #endif
    7  
    8 #include <iostream>
    9 #include <fstream>
   10 #include <vector>
   11  
   12 #include "moab/ReaderIface.hpp"
   13 #include "moab/Interface.hpp"
   14 #include "moab/RangeMap.hpp"
   15  
   16 #define MAT_PROP_TABLE_TAG  "mat_prop_table"
   17 #define PHYS_PROP_TABLE_TAG "phys_prop_table"
   18  
   19 namespace moab
   20 {
   21  
   22 class ReadUtilIface;
   23  
   24 class ReadIDEAS : public ReaderIface
   25 {
   26  
   27   public:
   28     static ReaderIface* factory( Interface* );
   29  
   30     ErrorCode load_file( const char* file_name,
   31                          const EntityHandle* file_set,
   32                          const FileOptions& opts,
   33                          const SubsetList* subset_list = 0,
   34                          const Tag* file_id_tag        = 0 );
   35  
   36     ErrorCode read_tag_values( const char* file_name,
   37                                const char* tag_name,
   38                                const FileOptions& opts,
   39                                std::vector< int >& tag_values_out,
   40                                const SubsetList* subset_list = 0 );
   41  
   42     //! Constructor
   43     ReadIDEAS( Interface* impl = NULL );
   44  
   45     //! Destructor
   46     virtual ~ReadIDEAS() {}
   47  
   48   protected:
   49     ErrorCode skip_header();
   50     ErrorCode create_vertices( EntityHandle& first_vertex, const Tag* file_id_tag );
   51     ErrorCode create_elements( EntityHandle first_vertex, const Tag* file_id_tag );
   52  
   53   private:
   54     std::ifstream file;
   55     RangeMap< int, EntityHandle > nodeIdMap;
   56  
   57     // Read mesh interface
   58     ReadUtilIface* readMeshIface;
   59  
   60     // MOAB Interface
   61     Interface* MBI;
   62  
   63     /* Universal dataset numbers
   64        An integer describes a chunk of information. These chunks include headers,
   65        units, nodes, elements, patches, etc... described in the OpenFOAM IDEAS
   66        reader.
   67  
   68        1) http://amira.zib.de/usersguide31/hxideas/HxFileFormat_IDEAS.html
   69        55,2414 data at nodes
   70  
   71        2) http://www.sdrl.uc.edu/universal-file-formats-for-modal-analysis-testing-1
   72        /file-format-storehouse/unv_0015.htm/
   73        15 nodes with single precision coordinates
   74        line1 (4I10,1P3E13.5): node_label coord_sys_num displacement_sys_num
   75        color x y z */
   76     static const unsigned SINGLE_PRECISION_NODES = 15;
   77  
   78     /* 3) http://www.sdrl.uc.edu/pdf/test_universal_file_formats.pdf
   79        781,2411 nodes with double precision coordinates
   80        line1 (4I10):       node_label coord_sys_num displacement_sys_num color
   81        line2 (1P3D25.16):  x y z */
   82     static const unsigned DOUBLE_PRECISION_NODES0 = 781;
   83     static const unsigned DOUBLE_PRECISION_NODES1 = 2411;
   84  
   85     /* 4) http://www.sdrl.uc.edu/universal-file-formats-for-modal-analysis-testing-1
   86        /file-format-storehouse/unv_0780.htm/
   87        71, 780, 2412 element definitions
   88        line1 (8I10): element_label fe_id phys_prop_bin_num phys_prop_num
   89        mat_prop_bin_num mat_prop_num color num_of_nodes
   90        line2 (8I10): connectivity_node_labels */
   91     static const unsigned ELEMENTS0 = 71;
   92     static const unsigned ELEMENTS1 = 780;
   93     static const unsigned ELEMENTS2 = 2412;
   94  
   95     /* Mesh elements exist inside chunks 71, 780, and 2412. Each element definition
   96        includes the finite element id that describes the element type. These are
   97        used in the OpenFOAM IDEAS reader. The canonical node ordering matches that
   98        of MBCN, as suggested by the Gmsh 2.2.3 source code.*/
   99     static const int ROD0  = 11;
  100     static const int ROD1  = 171;
  101     static const int TRI0  = 41;
  102     static const int TRI1  = 91;
  103     static const int QUAD0 = 44;
  104     static const int QUAD1 = 94;
  105     static const int TET   = 111;
  106     static const int WEDGE = 112;
  107     static const int HEX   = 115;
  108 };
  109  
  110 }  // namespace moab
  111 #endif