fprintf('--------------- %s\n',mfilename()) fprintf('3D meshing: mesh and partitioned mesh using ''MshFileVersion''=''4.1''\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','4.1'); fprintf(' -> msh file : %s\n',fc_oogmsh.pathrep(meshfile)); fprintf('2. Reading mesh file\n') Gh=fc_oogmsh.ooGmsh4(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','4.1'); fprintf(' -> partitioned msh file : %s\n',fc_oogmsh.pathrep(pmfile)); fprintf('4. Reading the partitioned mesh file\n') pGh=fc_oogmsh.ooGmsh4(pmfile); if fc_oogmsh.isgraphics() Gh.plotPhysicalElements() Gh.plotElementaryElements() pGh.plotPhysicalElements() pGh.plotElementaryElements() pGh.plotPartitionElements() pGh.plotInterfaceElements() end
--------------- demo14 3D meshing: mesh and partitioned mesh using 'MshFileVersion'='4.1' --------------- 1. Building a mesh file by using : -> gmsh version 4.11.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 4.11.0 [fc-oogmsh] Using command : gmsh -3 -setnumber N 10 -string "Mesh.MshFileVersion=4.1;" <fc-oogmsh>/geodir/3d/cylinderkey.geo -o <fc-oogmsh>/meshes/cylinderkey-10.msh Be patient... [fc-oogmsh] gmsh output : Info : Running '/fcopt/GMSH/64bit/4.11.0/bin/gmsh -3 -setnumber N 10 -string Mesh.MshFileVersion=4.1; <fc-oogmsh>/geodir/3d/cylinderkey.geo -o <fc-oogmsh>/meshes/cylinderkey-10.msh' [Gmsh 4.11.0, 1 node, max. 1 thread] Info : Started on Wed Jan 4 18:37:44 2023 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 : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 127: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 131: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 135: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 139: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 143: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : '<fc-oogmsh>/geodir/3d/cylinderkey.geo', line 147: 'Ruled Surface' command is deprecated: use 'Surface' instead Info : Done reading '<fc-oogmsh>/geodir/3d/cylinderkey.geo' Info : Meshing 1D... Info : [ 0%] Meshing curve 3 (Circle) Info : [ 10%] Meshing curve 4 (Circle) Info : [ 10%] Meshing curve 5 (Circle) Info : [ 10%] Meshing curve 6 (Circle) Info : [ 10%] Meshing curve 7 (Circle) Info : [ 10%] Meshing curve 8 (Line) Info : [ 20%] Meshing curve 9 (Line) Info : [ 20%] Meshing curve 10 (Line) Info : [ 20%] Meshing curve 11 (Line) Info : [ 20%] Meshing curve 12 (Line) Info : [ 20%] Meshing curve 13 (Line) Info : [ 20%] Meshing curve 14 (Circle) Info : [ 30%] Meshing curve 15 (Circle) Info : [ 30%] Meshing curve 16 (Circle) Info : [ 30%] Meshing curve 17 (Circle) Info : [ 30%] Meshing curve 18 (Circle) Info : [ 30%] Meshing curve 19 (Circle) Info : [ 30%] Meshing curve 20 (Circle) Info : [ 40%] Meshing curve 29 (Circle) Info : [ 40%] Meshing curve 30 (Circle) Info : [ 40%] Meshing curve 31 (Circle) Info : [ 40%] Meshing curve 32 (Circle) Info : [ 40%] Meshing curve 65 (Line) Info : [ 40%] Meshing curve 66 (Line) Info : [ 50%] Meshing curve 67 (Line) Info : [ 50%] Meshing curve 68 (Line) Info : [ 50%] Meshing curve 69 (Line) Info : [ 50%] Meshing curve 70 (Line) Info : [ 50%] Meshing curve 71 (Line) Info : [ 50%] Meshing curve 72 (Line) Info : [ 60%] Meshing curve 73 (Line) Info : [ 60%] Meshing curve 74 (Line) Info : [ 60%] Meshing curve 75 (Line) Info : [ 60%] Meshing curve 76 (Line) Info : [ 60%] Meshing curve 1005 (Circle) Info : [ 70%] Meshing curve 1007 (Circle) Info : [ 70%] Meshing curve 1008 (Circle) Info : [ 70%] Meshing curve 1009 (Circle) Info : [ 70%] Meshing curve 1010 (Circle) Info : [ 70%] Meshing curve 1011 (Circle) Info : [ 70%] Meshing curve 1012 (Circle) Info : [ 80%] Meshing curve 1013 (Circle) Info : [ 80%] Meshing curve 1014 (Circle) Info : [ 80%] Meshing curve 1015 (Circle) Info : [ 80%] Meshing curve 1016 (Circle) Info : [ 80%] Meshing curve 1017 (Circle) Info : [ 80%] Meshing curve 1018 (Line) Info : [ 90%] Meshing curve 1019 (Line) Info : [ 90%] Meshing curve 1020 (Circle) Info : [ 90%] Meshing curve 1021 (Line) Info : [ 90%] Meshing curve 1022 (Line) Info : [ 90%] Meshing curve 1023 (Line) Info : [ 90%] Meshing curve 1024 (Circle) Info : [100%] Meshing curve 1025 (Line) Info : [100%] Meshing curve 1026 (Circle) Info : [100%] Meshing curve 1027 (Circle) Info : [100%] Meshing curve 1028 (Circle) Info : [100%] Meshing curve 1029 (Circle) Info : Done meshing 1D (Wall 0.0082922s, CPU 0.011717s) Info : Meshing 2D... Info : [ 0%] Meshing surface 27 (Plane, Frontal-Delaunay) Info : [ 10%] Meshing surface 28 (Plane, Frontal-Delaunay) Info : [ 20%] Meshing surface 34 (Plane, Frontal-Delaunay) Info : [ 20%] Meshing surface 1006 (Plane, Frontal-Delaunay) Info : [ 30%] Meshing surface 1030 (Plane, Frontal-Delaunay) Info : [ 30%] Meshing surface 1036 (Plane, Frontal-Delaunay) Info : [ 40%] Meshing surface 1038 (Plane, Frontal-Delaunay) Info : [ 40%] Meshing surface 1040 (Plane, Frontal-Delaunay) Info : [ 50%] Meshing surface 1042 (Plane, Frontal-Delaunay) Info : [ 50%] Meshing surface 1044 (Plane, Frontal-Delaunay) Info : [ 60%] Meshing surface 1046 (Plane, Frontal-Delaunay) Info : [ 70%] Meshing surface 1048 (Plane, Frontal-Delaunay) Info : [ 70%] Meshing surface 1050 (Surface, Frontal-Delaunay) Info : [ 80%] Meshing surface 1052 (Surface, Frontal-Delaunay) Info : [ 80%] Meshing surface 1054 (Surface, Frontal-Delaunay) Info : [ 90%] Meshing surface 1056 (Surface, Frontal-Delaunay) Info : [ 90%] Meshing surface 1058 (Surface, Frontal-Delaunay) Info : [100%] Meshing surface 1060 (Surface, Frontal-Delaunay) Info : Done meshing 2D (Wall 0.182341s, CPU 0.295255s) Info : Meshing 3D... Info : 3D Meshing 1 volume with 1 connected component Info : Tetrahedrizing 3311 nodes... Info : Done tetrahedrizing 3319 nodes (Wall 0.0380898s, CPU 0.033173s) Info : Reconstructing mesh... Info : - Creating surface mesh Info : - Identifying boundary edges Info : - Recovering boundary Info : Done reconstructing mesh (Wall 0.0955693s, CPU 0.083528s) Info : Found volume 1062 Info : It. 0 - 0 nodes created - worst tet radius 3.65736 (nodes removed 0 0) Info : It. 500 - 500 nodes created - worst tet radius 1.55928 (nodes removed 0 0) Info : It. 1000 - 1000 nodes created - worst tet radius 1.29536 (nodes removed 0 0) Info : It. 1500 - 1500 nodes created - worst tet radius 1.16325 (nodes removed 0 0) Info : It. 2000 - 2000 nodes created - worst tet radius 1.07267 (nodes removed 0 0) Info : It. 2500 - 2500 nodes created - worst tet radius 1.00489 (nodes removed 0 0) Info : 3D refinement terminated (5855 nodes total): Info : - 0 Delaunay cavities modified for star shapeness Info : - 0 nodes could not be inserted Info : - 26298 tetrahedra created in 0.138849 sec. (189399 tets/s) Info : Done meshing 3D (Wall 0.321504s, CPU 0.301163s) Info : Optimizing mesh... Info : Optimizing volume 1062 Info : Optimization starts (volume = 5.2705) with worst = 0.00399 / average = 0.766187: Info : 0.00 < quality < 0.10 : 65 elements Info : 0.10 < quality < 0.20 : 166 elements Info : 0.20 < quality < 0.30 : 319 elements Info : 0.30 < quality < 0.40 : 485 elements Info : 0.40 < quality < 0.50 : 780 elements Info : 0.50 < quality < 0.60 : 1422 elements Info : 0.60 < quality < 0.70 : 3099 elements Info : 0.70 < quality < 0.80 : 6671 elements Info : 0.80 < quality < 0.90 : 9083 elements Info : 0.90 < quality < 1.00 : 4208 elements Info : 540 edge swaps, 7 node relocations (volume = 5.2705): worst = 0.212297 / average = 0.779256 (Wall 0.0141548s, CPU 0.014145s) Info : 542 edge swaps, 7 node relocations (volume = 5.2705): worst = 0.300128 / average = 0.779278 (Wall 0.0177127s, CPU 0.017705s) Info : No ill-shaped tets in the mesh :-) Info : 0.00 < quality < 0.10 : 0 elements Info : 0.10 < quality < 0.20 : 0 elements Info : 0.20 < quality < 0.30 : 0 elements Info : 0.30 < quality < 0.40 : 477 elements Info : 0.40 < quality < 0.50 : 736 elements Info : 0.50 < quality < 0.60 : 1403 elements Info : 0.60 < quality < 0.70 : 3098 elements Info : 0.70 < quality < 0.80 : 6712 elements Info : 0.80 < quality < 0.90 : 9190 elements Info : 0.90 < quality < 1.00 : 4188 elements Info : Done optimizing mesh (Wall 0.0521302s, CPU 0.051464s) Info : 5868 nodes 33061 elements Info : Writing '<fc-oogmsh>/meshes/cylinderkey-10.msh'... Info : Done writing '<fc-oogmsh>/meshes/cylinderkey-10.msh' Info : Stopped on Wed Jan 4 18:37:45 2023 (From start: Wall 0.655647s, CPU 0.729705s) [fc-oogmsh] Using gmsh 4.11.0 to write MSH file format version 4.1 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 4.11.0 [fc-oogmsh] Using command : gmsh -2 -part 5 -string "Mesh.MshFileVersion=4.1;" -saveall <fc-oogmsh>/meshes/cylinderkey-10.msh -o <fc-oogmsh>/meshes/cylinderkey-10-part5.msh Be patient... [fc-oogmsh] Using gmsh 4.11.0 to write MSH file format version 4.1 in <fc-oogmsh>/meshes/cylinderkey-10-part5.msh -> partitioned msh file : <fc-oogmsh>/meshes/cylinderkey-10-part5.msh 4. Reading the partitioned mesh file Warning: Limiting legend entries to 50. Specify a vector of graphics objects to display more than 50 entries. Warning: Limiting legend entries to 50. Specify a vector of graphics objects to display more than 50 entries. Warning: Limiting legend entries to 50. Specify a vector of graphics objects to display more than 50 entries.
