fprintf('--------------- %s\n',mfilename()) fprintf('2D meshing with gmsh.buildpartrectangle function using ''MshFileVersion''=''4.1''\n') fprintf('---------------\n') fprintf('1. Using gmsh.buildpartrectangle function :\n') meshfile=fc_oogmsh.gmsh.buildpartrectangle(2,1,3,2,10,'verbose',4,'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) if fc_oogmsh.isgraphics() Gh.plotPhysicalElements() Gh.plotElementaryElements() end
--------------- demo13 2D meshing with gmsh.buildpartrectangle function using 'MshFileVersion'='4.1' --------------- 1. Using gmsh.buildpartrectangle function : [fc-oogmsh] Input file : <fc-oogmsh>/geodir/2d/rectanglepart.geo [fc-oogmsh] Overwritting mesh file <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh [fc-oogmsh] Starting building mesh <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh with gmsh 4.5.1 [fc-oogmsh] Using command : '/fcopt/GMSH/64bit/4.5.1/bin/gmsh' -2 -string "Mesh.MshFileVersion=4.1;" -setnumber N 10 -setnumber NX 3 -setnumber NY 2 -setnumber LX 2 -setnumber LY 1 '/home/cuvelier/Travail/Recherch/Matlab/fc-oogmsh/geodir/2d/rectanglepart.geo' -o '/home/cuvelier/Travail/Recherch/Matlab/fc-oogmsh/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh' Be patient... [fc-oogmsh] Starting building mesh <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh with gmsh 4.5.1 [fc-oogmsh] Using command : gmsh -2 -string "Mesh.MshFileVersion=4.1;" -setnumber N 10 -setnumber NX 3 -setnumber NY 2 -setnumber LX 2 -setnumber LY 1 <fc-oogmsh>/geodir/2d/rectanglepart.geo -o <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh Be patient... [fc-oogmsh] gmsh output : Info : Running '/fcopt/GMSH/64bit/4.5.1/bin/gmsh -2 -string Mesh.MshFileVersion=4.1; -setnumber N 10 -setnumber NX 3 -setnumber NY 2 -setnumber LX 2 -setnumber LY 1 <fc-oogmsh>/geodir/2d/rectanglepart.geo -o <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh' [Gmsh 4.5.1, 1 node, max. 1 thread] Info : Started on Fri Mar 20 15:08:37 2020 Info : Reading '<fc-oogmsh>/geodir/2d/rectanglepart.geo'... Info : Reading '<fc-oogmsh>/geodir/2d/partitions01_data.geo'... Info : Done reading '<fc-oogmsh>/geodir/2d/partitions01_data.geo' Info : Reading '<fc-oogmsh>/geodir/2d/partitions_shape.geo'... Info : Done reading '<fc-oogmsh>/geodir/2d/partitions_shape.geo' Info : Done reading '<fc-oogmsh>/geodir/2d/rectanglepart.geo' Info : Meshing 1D... Info : [ 0 %] Meshing curve 1 (Line) Info : [ 10 %] Meshing curve 2 (Line) Info : [ 20 %] Meshing curve 3 (Line) Info : [ 20 %] Meshing curve 4 (Line) Info : [ 30 %] Meshing curve 5 (Line) Info : [ 30 %] Meshing curve 6 (Line) Info : [ 40 %] Meshing curve 7 (Line) Info : [ 50 %] Meshing curve 8 (Line) Info : [ 50 %] Meshing curve 9 (Line) Info : [ 60 %] Meshing curve 10 (Line) Info : [ 60 %] Meshing curve 11 (Line) Info : [ 70 %] Meshing curve 12 (Line) Info : [ 80 %] Meshing curve 13 (Line) Info : [ 80 %] Meshing curve 14 (Line) Info : [ 90 %] Meshing curve 15 (Line) Info : [ 90 %] Meshing curve 16 (Line) Info : [100 %] Meshing curve 17 (Line) Info : Done meshing 1D (0.005939 s) Info : Meshing 2D... Info : [ 0 %] Meshing surface 19 (Plane, Frontal) Info : [ 20 %] Meshing surface 21 (Plane, Frontal) Info : [ 40 %] Meshing surface 23 (Plane, Frontal) Info : [ 50 %] Meshing surface 25 (Plane, Frontal) Info : [ 70 %] Meshing surface 27 (Plane, Frontal) Info : [ 90 %] Meshing surface 29 (Plane, Frontal) Info : Done meshing 2D (0.065868 s) Info : 308 nodes 667 elements Info : Writing '<fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh'... Info : Done writing '<fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh' Info : Stopped on Fri Mar 20 15:08:37 2020 [fc-oogmsh] Using gmsh 4.5.1 to write MSH file format version 4.1 in <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh -> msh file : <fc-oogmsh>/meshes/rectanglepart-Lx2.000-Ly1.000-Nx3-Ny2-N10.msh 2. Reading mesh file Gh = fc_oogmsh.ooGmsh4 with properties: q: (2x308 double) nq: 308 double dim: 2 double d: 2 double toGlobal: (1x308 double) meshfile: (1x101 char) partitionnedfile: 0 logical MeshFormat: (1x1 struct) PhysicalNames: [] Entities: (1x1 struct) PartitionedEntities: [] Nodes: (1x1 struct) Elements: (1x1 struct) PeriodicLinks: [] Info: (1x1 struct) orders: 1 double