Skip to content
Snippets Groups Projects
Select Git revision
  • 7db8eac073206d58061cfa3f5ea651c285644937
  • master default protected
  • dof-renumbering
  • test-dof-hash
  • gdemesy-master-patch-30528
  • eval-space-time
  • oscillating_multiharm
  • MH_movement
  • axisqu
  • write_vtu_and_ensight_formats
  • movingband
  • CP_1972_add_vtu_file_writing
  • mortar
  • fast_freq_sweep_Resolution
  • applyresolvent_again
  • marteaua-master-patch-54323
  • patch-1
  • binde-master-patch-08072
  • binde-master-patch-52461
  • BCGSL
  • resolvent
  • getdp_3_5_0
  • getdp_3_4_0
  • getdp_3_3_0
  • getdp_3_2_0
  • getdp_3_1_0
  • getdp_3_0_4
  • getdp_3_0_3
  • getdp_3_0_2
  • getdp_3_0_1
  • getdp_3_0_0
  • onelab_mobile_2.1.0
  • getdp_2_11_3 protected
  • getdp_2_11_2 protected
  • getdp_2_11_1 protected
  • getdp_2_11_0 protected
  • getdp_2_10_0 protected
  • getdp_2_9_2 protected
  • getdp_2_9_1 protected
  • getdp_2_9_0 protected
  • getdp_2_8_0 protected
41 results

Cal_Quantity.cpp

Blame
  • driverOCC.cpp 1.50 KiB
    // 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
    
    #if !defined(WIN32) || defined(__CYGWIN__)
    #include "config.h"
    #endif
    #include "TopoDS_Shape.hxx"
    #include "BRep_Tool.hxx"
    #include "BRep_Builder.hxx"
    #include "BRepTools.hxx"
    
    #include <stdio.h>
    #include <gmsh/Gmsh.h>
    #include <gmsh/GModel.h>
    #include <gmsh/MElement.h>
    
    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);
    
      // import the shape in gmsh and mesh it
      GmshInitialize(argc, argv);
    
      GModel m;
      m.importOCCShape((void*)&shape, 0);
      m.mesh(2);
    
      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();
    }