GMSHInputModule.h

Go to the documentation of this file.

 
#pragma once
 
#include <filesystem>
#include <gmsh.h>
 
#include <jem/base/Array.h>
#include <jem/base/CString.h>
#include <jem/base/Class.h>
#include <jem/base/System.h>
#include <jem/util/ArrayBuffer.h>
#include <jem/util/HashDictionary.h>
#include <jem/util/Properties.h>
 
#include <jive/app/Module.h>
#include <jive/app/ModuleFactory.h>
#include <jive/app/Names.h>
#include <jive/fem/ElementGroup.h>
#include <jive/fem/XElementSet.h>
#include <jive/fem/XNodeSet.h>
#include <jive/model/Actions.h>
#include <jive/model/Model.h>
#include <jive/model/StateVector.h>
#include <jive/util/Assignable.h>
#include <jive/util/DenseTable.h>
#include <jive/util/DofSpace.h>
#include <jive/util/FuncUtils.h>
#include <jive/util/Globdat.h>
#include <jive/util/XTable.h>
 
using jem::Array;
using jem::max;
using jem::newInstance;
 
using jive::ALL;
using jive::idx_t;
using jive::IdxMatrix;
using jive::IdxVector;
using jive::Matrix;
using jive::Properties;
using jive::Ref;
using jive::String;
using jive::StringVector;
using jive::Vector;
 
using jive::app::Module;
using jive::fem::ElementGroup;
using jive::fem::XElementSet;
using jive::fem::XNodeSet;
using jive::model::Model;
using jive::model::StateVector;
using jive::util::Assignable;
using jive::util::DofSpace;
using jive::util::Function;
using jive::util::Globdat;
 
typedef jem::util::ArrayBuffer<idx_t> IdxBuffer;
 
static_assert(GMSH_API_VERSION_MAJOR == 4, "Wrong GMSH API version, please use 4.9.5");
static_assert(GMSH_API_VERSION_MINOR == 9, "Wrong GMSH API version, please use 4.9.5");
static_assert(GMSH_API_VERSION_PATCH == 5, "Wrong GMSH API version, please use 4.9.5");
 
class GMSHInputModule : public Module
{
public:
  JEM_DECLARE_CLASS(GMSHInputModule, Module);
 
  static const char *TYPE_NAME;       
  static const char *ORDER;           
  static const char *MESH_DIM;        
  static const char *SAVE_MSH;        
  static const char *STORE_TANGENTS;  
  static const char *ENTITY_NAMES[4]; 
  static const char *ONELAB_PROPS;    
  static const char *VERBOSE;         
  static const char *OUT_FILE;        
  static const char *OUT_TABLES;      
 
  explicit GMSHInputModule(const String &name = "GMSHInput");
 
  Status init(const Properties &conf,
              const Properties &props,
              const Properties &globdat) override;
 
  Status run(const Properties &globdat) override;
 
  void shutdown(const Properties &globdat) override;
 
  static Ref<Module> makeNew(const String &name,
                             const Properties &conf,
                             const Properties &props,
                             const Properties &globdat);
 
  static void declare();
 
protected:
  ~GMSHInputModule();
 
  void openMesh_(const String &geoFile, const idx_t order);
 
  void prepareOnelab_(const Properties &onelabProps);
 
  void createNodes_(const idx_t dim);
 
  void createElems_(const Properties &globdat);
 
  void storeTangents_(const Properties &globdat);
 
  void writeOutFile_(const Properties &globdat) const;
 
private:
  Assignable<XNodeSet> nodes_;                               
  std::unordered_map<std::size_t, idx_t> gmshToJiveNodeMap_; 
  Assignable<XElementSet> elements_;                         
  std::unordered_map<std::size_t, idx_t> gmshToJiveElemMap_; 
  IdxMatrix entities_;                                       
 
  bool verbose_;             
  bool writeOutput_;         
  Ref<Function> sampleCond_; 
  String outFile_;           
  String outExt_;            
  StringVector outTables_;   
  idx_t nodeView_;           
  idx_t elemView_;           
};