demo04
fprintf('--------------- demo04\n')
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);
fc_oogmsh.plot_PhysicalTags(Gh)
fc_oogmsh.plot_ElementaryTags(Gh)
--------------- demo04
3D meshing: mesh and partitioned mesh using 'MshFileVersion'='2.2'
---------------
1. Building a mesh file by using :
-> gmsh version 4.2.2
-> 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.2.2
[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/4.2.2/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 4.2.2, 1 node, max. 1 thread]
Info : Started on Sun Mar 24 09:33:57 2019
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 : 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.009098 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 (Surface, MeshAdapt)
Info : Meshing surface 1052 (Surface, MeshAdapt)
Info : Meshing surface 1054 (Surface, MeshAdapt)
Info : Meshing surface 1056 (Surface, MeshAdapt)
Info : Meshing surface 1058 (Surface, MeshAdapt)
Info : Meshing surface 1060 (Surface, MeshAdapt)
Info : Done meshing 2D (0.361595 s)
Info : Meshing 3D...
Info : 3D Meshing 1 volumes with 1 connected components
Info : Tetrahedrization of 3320 points in 0.039315 seconds
Info : Reconstructing mesh ...
Info : Creating surface mesh...
Info : Identifying boundary edges...
Info : Boundary Recovery...
Info : Reconstruct time : 0.090104 sec
Info : Found region 1062
Info : 0 points created - worst tet radius 3.78846 (points removed 0 0)
Info : 500 points created - worst tet radius 1.41157 (points removed 0 0)
Info : 1000 points created - worst tet radius 1.17018 (points removed 0 0)
Info : 1500 points created - worst tet radius 1.05738 (points removed 0 0)
Info : 3D point insertion terminated (5152 points created):
Info : - 0 Delaunay cavities modified for star shapeness
Info : - 0 points could not be inserted
Info : - 21729 tetrahedra created in 0.096855 sec. (224345 tets/s)
Info : Done meshing 3D (0.26465 s)
Info : Optimizing 3D mesh...
Info : Optimizing volume 1062
Info : Optimization starts (volume = 5.27044) with worst = 0.0172298 / average = 0.739253:
Info : 0.00 < quality < 0.10 : 56 elements
Info : 0.10 < quality < 0.20 : 184 elements
Info : 0.20 < quality < 0.30 : 289 elements
Info : 0.30 < quality < 0.40 : 458 elements
Info : 0.40 < quality < 0.50 : 762 elements
Info : 0.50 < quality < 0.60 : 1813 elements
Info : 0.60 < quality < 0.70 : 3407 elements
Info : 0.70 < quality < 0.80 : 5701 elements
Info : 0.80 < quality < 0.90 : 6330 elements
Info : 0.90 < quality < 1.00 : 2729 elements
Info : 520 edge swaps, 8 node relocations (volume = 5.27044): worst = 0.281935 / average = 0.753063 (0.014201 s)
Info : 525 edge swaps, 8 node relocations (volume = 5.27044): worst = 0.300297 / average = 0.75322 (0.017745 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 : 449 elements
Info : 0.40 < quality < 0.50 : 732 elements
Info : 0.50 < quality < 0.60 : 1862 elements
Info : 0.60 < quality < 0.70 : 3380 elements
Info : 0.70 < quality < 0.80 : 5739 elements
Info : 0.80 < quality < 0.90 : 6383 elements
Info : 0.90 < quality < 1.00 : 2726 elements
Info : Done optimizing 3D mesh (0.043233 s)
Info : 5165 vertices 28530 elements
Info : Writing '<fc-oogmsh>/meshes/cylinderkey-10.msh'...
Info : Done writing '<fc-oogmsh>/meshes/cylinderkey-10.msh'
Info : Stopped on Sun Mar 24 09:33:58 2019
[fc-oogmsh] Using gmsh 4.2.2 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 4.2.2
[fc-oogmsh] Using command : gmsh -2 -saveall -part 5 -string "Mesh.MshFileVersion=2.2;" <fc-oogmsh>/meshes/cylinderkey-10.msh -o <fc-oogmsh>/meshes/cylinderkey-10-part5.msh
Be patient...
[fc-oogmsh] Using gmsh 4.2.2 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
2-node line :
physical_tag= 1075, nq= 64, nme= 64
physical_tag= 1076, nq= 64, nme= 64
physical_tag= 1077, nq= 31, nme= 31
physical_tag= 1078, nq= 31, nme= 31
physical_tag= 1079, nq= 31, nme= 31
physical_tag= 1080, nq= 31, nme= 31
physical_tag= 1081, nq= 42, nme= 42
physical_tag= 1082, nq= 42, nme= 42
physical_tag= 1083, nq= 168, nme= 160
3-node triangle :
physical_tag= 1, nq= 1562, nme= 2996
physical_tag= 10, nq= 205, nme= 352
physical_tag= 11, nq= 197, nme= 336
physical_tag= 31, nq= 707, nme= 1278
physical_tag= 1000, nq= 420, nme= 676
physical_tag= 1020, nq= 55, nme= 77
physical_tag= 1021, nq= 55, nme= 77
physical_tag= 2000, nq= 420, nme= 676
physical_tag= 2020, nq= 55, nme= 77
physical_tag= 2021, nq= 56, nme= 79
4-node tetrahedron :
physical_tag= 1, nq= 5152, nme= 21271
2-node line :
elementary_tag= 3, nq= 9, nme= 8
elementary_tag= 4, nq= 9, nme= 8
elementary_tag= 5, nq= 2, nme= 1
elementary_tag= 6, nq= 11, nme= 10
elementary_tag= 7, nq= 14, nme= 13
elementary_tag= 8, nq= 6, nme= 5
elementary_tag= 9, nq= 4, nme= 3
elementary_tag= 10, nq= 6, nme= 5
elementary_tag= 11, nq= 6, nme= 5
elementary_tag= 12, nq= 4, nme= 3
elementary_tag= 13, nq= 6, nme= 5
elementary_tag= 14, nq= 11, nme= 10
elementary_tag= 15, nq= 14, nme= 13
elementary_tag= 16, nq= 5, nme= 4
elementary_tag= 17, nq= 2, nme= 1
elementary_tag= 18, nq= 4, nme= 3
elementary_tag= 19, nq= 2, nme= 1
elementary_tag= 20, nq= 4, nme= 3
elementary_tag= 29, nq= 17, nme= 16
elementary_tag= 30, nq= 17, nme= 16
elementary_tag= 31, nq= 17, nme= 16
elementary_tag= 32, nq= 17, nme= 16
elementary_tag= 65, nq= 21, nme= 20
elementary_tag= 66, nq= 21, nme= 20
elementary_tag= 67, nq= 21, nme= 20
elementary_tag= 68, nq= 21, nme= 20
elementary_tag= 69, nq= 21, nme= 20
elementary_tag= 70, nq= 21, nme= 20
elementary_tag= 71, nq= 21, nme= 20
elementary_tag= 72, nq= 21, nme= 20
elementary_tag= 1005, nq= 4, nme= 3
elementary_tag= 1007, nq= 14, nme= 13
elementary_tag= 1008, nq= 5, nme= 4
elementary_tag= 1009, nq= 11, nme= 10
elementary_tag= 1010, nq= 4, nme= 3
elementary_tag= 1011, nq= 2, nme= 1
elementary_tag= 1012, nq= 11, nme= 10
elementary_tag= 1013, nq= 4, nme= 3
elementary_tag= 1014, nq= 2, nme= 1
elementary_tag= 1015, nq= 14, nme= 13
elementary_tag= 1016, nq= 2, nme= 1
elementary_tag= 1017, nq= 4, nme= 3
elementary_tag= 1018, nq= 4, nme= 3
elementary_tag= 1019, nq= 6, nme= 5
elementary_tag= 1020, nq= 9, nme= 8
elementary_tag= 1021, nq= 6, nme= 5
elementary_tag= 1022, nq= 4, nme= 3
elementary_tag= 1023, nq= 6, nme= 5
elementary_tag= 1024, nq= 9, nme= 8
elementary_tag= 1025, nq= 6, nme= 5
elementary_tag= 1026, nq= 17, nme= 16
elementary_tag= 1027, nq= 17, nme= 16
elementary_tag= 1028, nq= 17, nme= 16
elementary_tag= 1029, nq= 17, nme= 16
3-node triangle :
elementary_tag= 27, nq= 56, nme= 79
elementary_tag= 28, nq= 55, nme= 77
elementary_tag= 34, nq= 420, nme= 676
elementary_tag= 1006, nq= 420, nme= 676
elementary_tag= 1030, nq= 55, nme= 77
elementary_tag= 1036, nq= 55, nme= 77
elementary_tag= 1038, nq= 86, nme= 124
elementary_tag= 1040, nq= 156, nme= 260
elementary_tag= 1042, nq= 155, nme= 258
elementary_tag= 1044, nq= 83, nme= 118
elementary_tag= 1046, nq= 155, nme= 258
elementary_tag= 1048, nq= 156, nme= 260
elementary_tag= 1050, nq= 409, nme= 744
elementary_tag= 1052, nq= 414, nme= 754
elementary_tag= 1054, nq= 409, nme= 744
elementary_tag= 1056, nq= 414, nme= 754
elementary_tag= 1058, nq= 197, nme= 336
elementary_tag= 1060, nq= 205, nme= 352
4-node tetrahedron :
elementary_tag= 1062, nq= 5152, nme= 21271
