Skip to content
Snippets Groups Projects
Select Git revision
  • 27633c09fb73bab37e13b9e15d37828357c798f3
  • master default
  • cgnsUnstructured
  • partitioning
  • poppler
  • HighOrderBLCurving
  • gmsh_3_0_4
  • gmsh_3_0_3
  • gmsh_3_0_2
  • gmsh_3_0_1
  • gmsh_3_0_0
  • gmsh_2_16_0
  • gmsh_2_15_0
  • gmsh_2_14_1
  • gmsh_2_14_0
  • gmsh_2_13_2
  • gmsh_2_13_1
  • gmsh_2_12_0
  • gmsh_2_11_0
  • gmsh_2_10_1
  • gmsh_2_10_0
  • gmsh_2_9_3
  • gmsh_2_9_2
  • gmsh_2_9_1
  • gmsh_2_9_0
  • gmsh_2_8_6
26 results

mainOcc.cpp

Blame
  • Forked from gmsh / gmsh
    Source project has a limited visibility.
    mainOcc.cpp 2.69 KiB
    // This demonstrates how to import an OpenCASCADE model WITHOUT using
    // Gmsh's own I/O layer (useful if you create OCC objects in your own
    // code, that you want to import dynamically in Gmsh).
    //
    // If you just want to load OpenCASCADE files (.brep, .step, etc.) you
    // SHOULD NOT use this interface: simply use the standard I/O API.
    
    // g++ -I/usr/local/opencascade/inc driverOCC.cpp -lGmsh
    // -L/usr/local/opencascade/lib -lTKSTEP -lTKSTEP209 -lTKSTEPAttr
    // -lTKSTEPBase -lTKIGES -lTKXSBase -lTKOffset -lTKFeat -lTKFillet
    // -lTKBool -lTKShHealing -lTKMesh -lTKHLR -lTKBO -lTKPrim -lTKTopAlgo
    // -lTKGeomAlgo -lTKBRep -lTKGeomBase -lTKG3d -lTKG2d -lTKAdvTools
    // -lTKMath -lTKernel -lm
    
    #include <stdio.h>
    #include <gmsh/Gmsh.h>
    #include <gmsh/GmshConfig.h>
    #include <gmsh/GModel.h>
    #include <gmsh/MElement.h>
    
    #if !defined(HAVE_NO_OCC_CONFIG_H)
    #include "config.h"
    #endif
    #include <TopoDS_Shape.hxx>
    #include <BRep_Tool.hxx>
    #include <BRep_Builder.hxx>
    #include <BRepTools.hxx>
    
    class mymsg : public GmshMessage{
    private:
      GModel *_model;
    public:
      mymsg(GModel *model) : _model(model), GmshMessage() {}
      void operator()(std::string level, std::string msg)
      {
        printf("level=%s msg=%s\n", level.c_str(), msg.c_str());
    
        if(level == "Fatal" || level == "Error"){
          GEntity *e = _model->getCurrentMeshEntity();
          if(e){
            printf("error occurred while meshing entity:\n");
            printf("  tag=%d\n", e->tag());
            printf("  dimension=%d\n", e->dim());
            printf("  native pointer=%p\n", e->getNativePtr());
          }
        }
        if(level == "Fatal") throw "Fatal error in Gmsh";
      }
    };
    
    int main(int argc, char **argv)
    {
      // create an OCC shape (by loading it from a brep file)
      TopoDS_Shape shape;
      BRep_Builder builder;
      BRepTools::Read(shape, argv[1], builder);
      BRepTools::Clean(shape);
    
      // initialize gmsh
      GmshInitialize(argc, argv);
    
      // create a model and set error handler 
      GModel m;
      mymsg c(&m);
      GmshSetMessageHandler(&c);
    
      // import the shape, and mesh it
      m.importOCCShape((void*)&shape);
      try{
        m.mesh(2);
      }
      catch(...){
        printf("Unrecoverable error in gmsh--aborting mesh!\n");
      }
    
      for(GModel::fiter it = m.firstFace(); it != m.lastFace(); ++it){
        GFace *f = *it;
        printf("Surface %d contains %d elements:\n", f->tag(), f->getNumMeshElements());
        for(unsigned int i = 0; i < f->getNumMeshElements(); i++){
          MElement *e = f->getMeshElement(i);
          printf("  element %d:", e->getNum());
          for(unsigned int j = 0; j < e->getNumVertices(); j++){
            MVertex *v = e->getVertex(j);
            printf(" %d (%g,%g,%g)", v->getNum(), v->x(), v->y(), v->z());
          }
          printf("\n");
        }
        printf("\n");
      }
    
      m.writeMSH("test.msh");
    
      GmshFinalize();
    }