gmsh_file_v2.h

 
#ifndef __1fd687af60b74bf6aad50d509ecbc4da
#define __1fd687af60b74bf6aad50d509ecbc4da
 
#include <lf/mesh/mesh.h>
 
namespace lf::io {
 
struct GMshFileV2 {
  using size_type = mesh::Mesh::size_type;
  std::string VersionNumber;
  bool IsBinary = false;
  int DoubleSize = 64;
 
  struct PhysicalEntity {
    int Dimension = 0;
 
    int Number = 0;
    std::string Name;
  };
 
  std::vector<PhysicalEntity> PhysicalEntities;
 
  std::vector<std::pair<size_type, Eigen::Vector3d>> Nodes;
 
  enum class ElementType : int {
    EDGE2 = 1,     
    TRIA3 = 2,     
    QUAD4 = 3,     
    TET4 = 4,      
    HEX8 = 5,      
    PRISM6 = 6,    
    PYRAMID5 = 7,  
    EDGE3 = 8,     
    TRIA6 = 9,   
    QUAD9 = 10,  
    TET10 = 11,  
    HEX27 = 12,  
    PRISM18 =
        13,  
    PYRAMID14 =
        14,  
    POINT = 15,    
    QUAD8 = 16,    
    HEX20 = 17,    
    PRISM15 = 18,  
    PYRAMID13 = 19,  
    TRIA9 = 20,  
    TRIA10 = 21,  
    TRIA12 = 22,  
    TRIA15 = 23,  
    TRIA15_5 = 24,  
    TRIA21 = 25,  
    EDGE4 = 26,   
    EDGE5 = 27,   
    EDGE6 = 28,   
    TET20 = 29,   
    TET35 = 30,   
    TET56 = 31,   
    HEX64 = 92,   
    HEX125 = 93   
  };
 
  static const std::vector<ElementType> AllElementTypes;
 
  struct Element {
    size_type Number = 0;
    ElementType Type = ElementType::POINT;
    int PhysicalEntityNr = 0;
    int ElementaryEntityNr = 0;
 
    std::vector<int> MeshPartitions;
 
    std::vector<size_type> NodeNumbers;
  };
 
  std::vector<Element> Elements;
 
  struct PeriodicEntity {
    int Dimension = 0;
    int ElementarySlaveNr = 0;
    int ElementaryMasterNr = 0;
    std::vector<std::pair<size_type, size_type>> NodeMapping;
  };
 
  std::vector<PeriodicEntity> Periodic;
};
 
std::ostream& operator<<(std::ostream& stream, const GMshFileV2& mf);
 
std::ostream& operator<<(std::ostream& stream, GMshFileV2::ElementType et);
 
unsigned int NumNodes(GMshFileV2::ElementType et);
 
int DimOf(GMshFileV2::ElementType et);
 
base::RefEl RefElOf(GMshFileV2::ElementType et);
 
GMshFileV2 readGmshFileV2(std::string::const_iterator begin,
                          std::string::const_iterator end,
                          const std::string& version, bool is_binary,
                          int size_t_size, int one,
                          const std::string& filename);
 
}  // namespace lf::io
 
#endif  // __1fd687af60b74bf6aad50d509ecbc4da