demo14
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
--------------- publish
3D meshing: mesh and partitioned mesh using 'MshFileVersion'='4.1'
---------------
1. Building a mesh file by using :
-> gmsh version 4.5.1
-> 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.5.1
[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.5.1/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.5.1, 1 node, max. 1 thread]
Info : Started on Fri Mar 20 15:02:54 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 : '<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 (0.045101 s)
Info : Meshing 2D...
Info : [ 0 %] Meshing surface 27 (Plane, Frontal)
Info : [ 10 %] Meshing surface 28 (Plane, Frontal)
Info : [ 20 %] Meshing surface 34 (Plane, Frontal)
Info : [ 20 %] Meshing surface 1006 (Plane, Frontal)
Info : [ 30 %] Meshing surface 1030 (Plane, Frontal)
Info : [ 30 %] Meshing surface 1036 (Plane, Frontal)
Info : [ 40 %] Meshing surface 1038 (Plane, Frontal)
Info : [ 40 %] Meshing surface 1040 (Plane, Frontal)
Info : [ 50 %] Meshing surface 1042 (Plane, Frontal)
Info : [ 50 %] Meshing surface 1044 (Plane, Frontal)
Info : [ 60 %] Meshing surface 1046 (Plane, Frontal)
Info : [ 70 %] Meshing surface 1048 (Plane, Frontal)
Info : [ 70 %] Meshing surface 1050 (Surface, Frontal)
Info : [ 80 %] Meshing surface 1052 (Surface, Frontal)
Info : [ 80 %] Meshing surface 1054 (Surface, Frontal)
Info : [ 90 %] Meshing surface 1056 (Surface, Frontal)
Info : [ 90 %] Meshing surface 1058 (Surface, Frontal)
Info : [100 %] Meshing surface 1060 (Surface, Frontal)
Info : Done meshing 2D (0.398657 s)
Info : Meshing 3D...
Info : 3D Meshing 1 volume with 1 connected component
Info : Tetrahedrization of 3318 nodes in 0.034338 seconds
Info : Reconstructing mesh...
Info : - Creating surface mesh
Info : - Identifying boundary edges
Info : - Recovering boundary
Info : Done reconstructing mesh (0.085455 s)
Info : Found volume 1062
Info : It. 0 - 0 nodes created - worst tet radius 3.67875 (nodes removed 0 0)
Info : It. 500 - 500 nodes created - worst tet radius 1.56204 (nodes removed 0 0)
Info : It. 1000 - 1000 nodes created - worst tet radius 1.29036 (nodes removed 0 0)
Info : It. 1500 - 1500 nodes created - worst tet radius 1.16878 (nodes removed 0 0)
Info : It. 2000 - 2000 nodes created - worst tet radius 1.06982 (nodes removed 0 0)
Info : It. 2500 - 2500 nodes created - worst tet radius 1.00348 (nodes removed 0 0)
Info : 3D refinement terminated (5834 nodes total):
Info : - 0 Delaunay cavities modified for star shapeness
Info : - 0 nodes could not be inserted
Info : - 26100 tetrahedra created in 0.134153 sec. (194553 tets/s)
Info : Done meshing 3D (0.300233 s)
Info : Optimizing mesh...
Info : Optimizing volume 1062
Info : Optimization starts (volume = 5.27052) with worst = 0.00864879 / average = 0.767998:
Info : 0.00 < quality < 0.10 : 61 elements
Info : 0.10 < quality < 0.20 : 186 elements
Info : 0.20 < quality < 0.30 : 291 elements
Info : 0.30 < quality < 0.40 : 482 elements
Info : 0.40 < quality < 0.50 : 682 elements
Info : 0.50 < quality < 0.60 : 1404 elements
Info : 0.60 < quality < 0.70 : 3148 elements
Info : 0.70 < quality < 0.80 : 6441 elements
Info : 0.80 < quality < 0.90 : 9148 elements
Info : 0.90 < quality < 1.00 : 4257 elements
Info : 536 edge swaps, 5 node relocations (volume = 5.27052): worst = 0.300291 / average = 0.78108 (0.013859 s)
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 : 476 elements
Info : 0.40 < quality < 0.50 : 645 elements
Info : 0.50 < quality < 0.60 : 1391 elements
Info : 0.60 < quality < 0.70 : 3112 elements
Info : 0.70 < quality < 0.80 : 6466 elements
Info : 0.80 < quality < 0.90 : 9264 elements
Info : 0.90 < quality < 1.00 : 4251 elements
Info : Done optimizing mesh (0.045486 s)
Info : 5847 nodes 32860 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:02:54 2020
[fc-oogmsh] Using gmsh 4.5.1 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.5.1
[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.5.1 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
