fprintf('--------------- %s\n',mfilename()) fprintf('3D meshing: mesh and partitioned mesh using ''MshFileVersion''=''2.2''\n') fprintf('---------------\n') geofile='cylinderkey';np=5; fprintf('1. Building a mesh file by using :\n') fprintf(' -> gmsh version %s\n',fc_oogmsh.gmsh.version()) fprintf(' -> geo file : %s\n',geofile); meshfile=fc_oogmsh.gmsh.buildmesh3d(geofile,10,'verbose',3,'force',true,'MshFileVersion','2.2'); fprintf(' -> msh file : %s\n',fc_oogmsh.pathrep(meshfile)); fprintf('2. Reading mesh file\n') Gh=fc_oogmsh.ooGmsh2(meshfile); fprintf('3. Building partitioned mesh file\n') fprintf(' -> number of partitions : %d\n',np); pmfile=fc_oogmsh.gmsh.buildpartmesh(meshfile,np,'verbose',2,'force',true,'MshFileVersion','2.2'); fprintf(' -> partitioned msh file : %s\n',fc_oogmsh.pathrep(pmfile)); fprintf('4. Reading the partitioned mesh file\n') pGh=fc_oogmsh.ooGmsh2(pmfile); if fc_oogmsh.isgraphics() Gh.plotPhysicalElements() Gh.plotElementaryElements() pGh.plotPhysicalElements() pGh.plotElementaryElements() pGh.plotPartitionElements() pGh.plotInterfaceElements() end
--------------- demo04 3D meshing: mesh and partitioned mesh using 'MshFileVersion'='2.2' --------------- 1. Building a mesh file by using : -> gmsh version 2.16.0 -> geo file : cylinderkey [fc-oogmsh] Input file : <fc-oogmsh>/geodir/3d/cylinderkey.geo [fc-oogmsh] Overwritting mesh file <fc-oogmsh>/meshes/cylinderkey-10.msh [fc-oogmsh] Starting building mesh <fc-oogmsh>/meshes/cylinderkey-10.msh with gmsh 2.16.0 [fc-oogmsh] Using command : gmsh -3 -setnumber N 10 -string "Mesh.MshFileVersion=2.2;" <fc-oogmsh>/geodir/3d/cylinderkey.geo -o <fc-oogmsh>/meshes/cylinderkey-10.msh Be patient... [fc-oogmsh] gmsh output : Info : Running '/fcopt/GMSH/64bit/2.16.0/bin/gmsh -3 -setnumber N 10 -string Mesh.MshFileVersion=2.2; <fc-oogmsh>/geodir/3d/cylinderkey.geo -o <fc-oogmsh>/meshes/cylinderkey-10.msh' [Gmsh 2.16.0, 1 node, max. 1 thread] Info : Started on Fri Mar 20 15:07:38 2020 Info : Reading '<fc-oogmsh>/geodir/3d/cylinderkey.geo'... Info : Reading '<fc-oogmsh>/geodir/3d/options01_data.geo'... Info : Done reading '<fc-oogmsh>/geodir/3d/options01_data.geo' Info : Done reading '<fc-oogmsh>/geodir/3d/cylinderkey.geo' Info : Finalized high order topology of periodic connections Info : Meshing 1D... Info : Meshing curve 3 (Circle) Info : Meshing curve 4 (Circle) Info : Meshing curve 5 (Circle) Info : Meshing curve 6 (Circle) Info : Meshing curve 7 (Circle) Info : Meshing curve 8 (Line) Info : Meshing curve 9 (Line) Info : Meshing curve 10 (Line) Info : Meshing curve 11 (Line) Info : Meshing curve 12 (Line) Info : Meshing curve 13 (Line) Info : Meshing curve 14 (Circle) Info : Meshing curve 15 (Circle) Info : Meshing curve 16 (Circle) Info : Meshing curve 17 (Circle) Info : Meshing curve 18 (Circle) Info : Meshing curve 19 (Circle) Info : Meshing curve 20 (Circle) Info : Meshing curve 29 (Circle) Info : Meshing curve 30 (Circle) Info : Meshing curve 31 (Circle) Info : Meshing curve 32 (Circle) Info : Meshing curve 65 (Line) Info : Meshing curve 66 (Line) Info : Meshing curve 67 (Line) Info : Meshing curve 68 (Line) Info : Meshing curve 69 (Line) Info : Meshing curve 70 (Line) Info : Meshing curve 71 (Line) Info : Meshing curve 72 (Line) Info : Meshing curve 73 (Line) Info : Meshing curve 74 (Line) Info : Meshing curve 75 (Line) Info : Meshing curve 76 (Line) Info : Meshing curve 1005 (Circle) Info : Meshing curve 1007 (Circle) Info : Meshing curve 1008 (Circle) Info : Meshing curve 1009 (Circle) Info : Meshing curve 1010 (Circle) Info : Meshing curve 1011 (Circle) Info : Meshing curve 1012 (Circle) Info : Meshing curve 1013 (Circle) Info : Meshing curve 1014 (Circle) Info : Meshing curve 1015 (Circle) Info : Meshing curve 1016 (Circle) Info : Meshing curve 1017 (Circle) Info : Meshing curve 1018 (Line) Info : Meshing curve 1019 (Line) Info : Meshing curve 1020 (Circle) Info : Meshing curve 1021 (Line) Info : Meshing curve 1022 (Line) Info : Meshing curve 1023 (Line) Info : Meshing curve 1024 (Circle) Info : Meshing curve 1025 (Line) Info : Meshing curve 1026 (Circle) Info : Meshing curve 1027 (Circle) Info : Meshing curve 1028 (Circle) Info : Meshing curve 1029 (Circle) Info : Done meshing 1D (0.033809 s) Info : Meshing 2D... Info : Meshing surface 27 (Plane, Delaunay) Info : Meshing surface 28 (Plane, Delaunay) Info : Meshing surface 34 (Plane, Delaunay) Info : Meshing surface 1006 (Plane, Delaunay) Info : Meshing surface 1030 (Plane, Delaunay) Info : Meshing surface 1036 (Plane, Delaunay) Info : Meshing surface 1038 (Plane, Delaunay) Info : Meshing surface 1040 (Plane, Delaunay) Info : Meshing surface 1042 (Plane, Delaunay) Info : Meshing surface 1044 (Plane, Delaunay) Info : Meshing surface 1046 (Plane, Delaunay) Info : Meshing surface 1048 (Plane, Delaunay) Info : Meshing surface 1050 (Ruled surface, MeshAdapt) Info : Meshing surface 1052 (Ruled surface, MeshAdapt) Info : Meshing surface 1054 (Ruled surface, MeshAdapt) Info : Meshing surface 1056 (Ruled surface, MeshAdapt) Info : Meshing surface 1058 (Ruled surface, MeshAdapt) Info : Meshing surface 1060 (Ruled surface, MeshAdapt) Info : Done meshing 2D (0.438618 s) Info : Meshing 3D... Info : Delaunay Meshing 1 volumes with 1 connected components Info : Meshing volume 1062 (Delaunay) Delaunizing vertices... Delaunay seconds: 0.031235 Creating surface mesh ... Surface mesh seconds: 0.007388 Recovering boundaries... Boundary recovery seconds: 0.013792 Removing exterior tetrahedra ... Exterior tets removal seconds: 0.005188 Recovering Delaunayness... Delaunay recovery seconds: 0.007394 Optimizing mesh... Optimization seconds: 0.002693 Writing nodes. Writing elements. Writing faces. Writing edges. Output seconds: 0.004626 Total running seconds: 0.072627 Statistics: Input points: 3256 Input facets: 6512 Input segments: 9768 Input holes: 0 Input regions: 0 Mesh points: 3256 Mesh tetrahedra: 9963 Mesh faces: 23182 Mesh edges: 16475 Mesh faces on exterior boundary: 6512 Mesh edges on exterior boundary: 9768 Mesh faces on input facets: 6512 Mesh edges on input segments: 9768 Info : 3256 points 16475 edges and 6512 faces in the initial mesh Info : Found region 1062 Info : 0 points created - Worst tet radius is 3.89251 (PTS removed 0 0) Info : 3D point insertion terminated (5129 points created): Info : - 1 Delaunay cavities modified for star shapeness Info : - 0 points could not be inserted Info : - 21679 tetrahedra created in 0.116062 sec. (186788 tets/sec.) Info : Done meshing 3D (0.269016 s) Info : 5142 vertices 28826 elements Info : Writing '<fc-oogmsh>/meshes/cylinderkey-10.msh'... Info : Done writing '<fc-oogmsh>/meshes/cylinderkey-10.msh' Info : Stopped on Fri Mar 20 15:07:39 2020 [fc-oogmsh] Using gmsh 2.16.0 to write MSH file format version 2.2 in <fc-oogmsh>/meshes/cylinderkey-10.msh -> msh file : <fc-oogmsh>/meshes/cylinderkey-10.msh 2. Reading mesh file 3. Building partitioned mesh file -> number of partitions : 5 [fc-oogmsh] Input file : <fc-oogmsh>/meshes/cylinderkey-10.msh [fc-oogmsh] Overwritting mesh file <fc-oogmsh>/meshes/cylinderkey-10-part5.msh [fc-oogmsh] Starting building mesh <fc-oogmsh>/meshes/cylinderkey-10-part5.msh with gmsh 2.16.0 [fc-oogmsh] Using command : gmsh -2 -part 5 -string "Mesh.MshFileVersion=2.2;" -saveall <fc-oogmsh>/meshes/cylinderkey-10.msh -o <fc-oogmsh>/meshes/cylinderkey-10-part5.msh Be patient... [fc-oogmsh] Using gmsh 2.16.0 to write MSH file format version 2.2 in <fc-oogmsh>/meshes/cylinderkey-10-part5.msh -> partitioned msh file : <fc-oogmsh>/meshes/cylinderkey-10-part5.msh 4. Reading the partitioned mesh file