Controlling mesh resolution#
Coarse global resolution#
The keyword arguments default_resolution_min and default_resolution_max set the minimum and maximum characteristic lengths used by gmsh when computing the mesh.
They are used when other resolutions are not provided, and default_resolution_max effectively sets the minimum length possible, since when conflicting resolution at a point are given, the smallest one is taken.
# %matplotlib widget
import gdsfactory as gf
import meshio
from gdsfactory.generic_tech import get_generic_pdk
from gdsfactory.pdk import get_layer_stack
from gdsfactory.technology import LayerStack
from skfem.io import from_meshio
from gplugins.gmsh.get_mesh import create_physical_mesh, get_mesh
PDK = get_generic_pdk()
PDK.activate()
gf.config.rich_output()
waveguide = gf.components.straight_pin(length=10, taper=None)
waveguide_trimmed = gf.Component()
waveguide_trimmed.add_ref(
gf.geometry.trim(
component=waveguide,
domain=[[3, -4], [3, 4], [5, 4], [5, -4]],
)
)
filtered_layer_stack = LayerStack(
layers={
k: get_layer_stack().layers[k]
for k in (
"slab90",
"core",
"via_contact",
)
}
)
filename = "mesh"
def mesh_with_physicals(mesh, filename):
mesh_from_file = meshio.read(f"{filename}.msh")
return create_physical_mesh(mesh_from_file, "triangle")
With default_resolution_max set to 1 um and default_resolution_min set to 100 nm:
mesh = get_mesh(
component=waveguide_trimmed,
type="uz",
xsection_bounds=[(4, -4), (4, 4)],
layer_stack=filtered_layer_stack,
filename=f"{filename}.msh",
background_tag="oxide",
background_padding=(2.0, 2.0, 2.0, 2.0),
default_resolution_min=0.1,
default_resolution_max=1,
)
mesh = mesh_with_physicals(mesh, filename)
mesh = from_meshio(mesh)
mesh.draw().plot()
Info : [ 0%] Union
Info : [ 10%] Union - Performing intersection of shapes
Info : [ 80%] Union - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference - Performing intersection of shapes
Info : [ 80%] Difference - Building splits of containers
Info : [ 0%] Fragments
Info : [ 10%] Fragments - Performing intersection of shapes
Info : [ 80%] Fragments - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Filling splits of edges
Info : [ 80%] Fragments - Splitting faces
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference - Adding holes
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Building splits of containers
Info : [ 80%] Fragments - Splitting faces
Info : Meshing 1D...
Info : [ 0%] Meshing curve 3 (Line)
Info : [ 0%] Meshing curve 1 (Line)
Info : [ 0%] Meshing curve 4 (Line)
Info : [ 0%] Meshing curve 2 (Line)
Info : [ 10%] Meshing curve 5 (Line)
Info : [ 10%] Meshing curve 6 (Line)
Info : [ 20%] Meshing curve 7 (Line)
Info : [ 20%] Meshing curve 8 (Line)
Info : [ 20%] Meshing curve 15 (Line)
Info : [ 20%] Meshing curve 16 (Line)
Info : [ 20%] Meshing curve 17 (Line)
Info : [ 20%] Meshing curve 18 (Line)
Info : [ 30%] Meshing curve 19 (Line)
Info : [ 30%] Meshing curve 20 (Line)
Info : [ 50%] Meshing curve 22 (Line)
Info : [ 50%] Meshing curve 23 (Line)
Info : [ 50%] Meshing curve 24 (Line)
Info : [ 50%] Meshing curve 25 (Line)
Info : [ 30%] Meshing curve 21 (Line)
Info : [ 70%] Meshing curve 27 (Line)
Info : [ 70%] Meshing curve 29 (Line)
Info : [ 50%] Meshing curve 26 (Line)
Info : [ 70%] Meshing curve 28 (Line)
Info : [ 80%] Meshing curve 30 (Line)
Info : [ 80%] Meshing curve 32 (Line)
Info : [ 90%] Meshing curve 33 (Line)
Info : [ 80%] Meshing curve 31 (Line)
Info : [ 90%] Meshing curve 34 (Line)
Info : Done meshing 1D (Wall 0.00352673s, CPU 0.007651s)
Info : Meshing 2D...
Info : [ 0%] Meshing surface 1 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 3 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 4 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 5 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 6 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 2 (Plane, Frontal-Delaunay)
Info : Done meshing 2D (Wall 0.0134749s, CPU 0.036962s)
Info : 504 nodes 1093 elements
Info : Writing 'mesh.msh'...
Info : Done writing 'mesh.msh'
Info : Writing '/tmp/tmp3se6zbuv/mesh.msh'...
Info : Done writing '/tmp/tmp3se6zbuv/mesh.msh'
[]
With default_resolution_max set to 300 nm and default_resolution_max set to 50 nm:
mesh = get_mesh(
component=waveguide_trimmed,
type="uz",
xsection_bounds=[(4, -4), (4, 4)],
layer_stack=filtered_layer_stack,
filename=f"{filename}.msh",
background_tag="oxide",
background_padding=(2.0, 2.0, 2.0, 2.0),
default_resolution_min=0.05,
default_resolution_max=0.3,
)
mesh = mesh_with_physicals(mesh, filename)
mesh = from_meshio(mesh)
mesh.draw().show()
Info : [ 0%] Union
Info : [ 10%] Union - Performing intersection of shapes
Info : [ 80%] Union - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference - Performing intersection of shapes
Info : [ 80%] Difference - Building splits of containers
Info : [ 0%] Fragments
Info : [ 10%] Fragments - Performing intersection of shapes
Info : [ 80%] Fragments - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference
Info : [ 90%] Difference
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Filling splits of edges
Info : [ 80%] Fragments - Splitting faces
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference - Adding holes
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Building splits of containers
Info : [ 80%] Fragments - Splitting faces
Info : Meshing 1D...
Info : [ 0%] Meshing curve 2 (Line)
Info : [ 0%] Meshing curve 3 (Line)
Info : [ 0%] Meshing curve 1 (Line)
Info : [ 0%] Meshing curve 4 (Line)
Info : [ 0%] Meshing curve 5 (Line)
Info : [ 0%] Meshing curve 6 (Line)
Info : [ 0%] Meshing curve 7 (Line)
Info : [ 0%] Meshing curve 8 (Line)
Info : [ 0%] Meshing curve 15 (Line)
Info : [ 0%] Meshing curve 16 (Line)
Info : [ 0%] Meshing curve 17 (Line)
Info : [ 0%] Meshing curve 18 (Line)
Info : [ 0%] Meshing curve 19 (Line)
Info : [ 0%] Meshing curve 20 (Line)
Info : [ 0%] Meshing curve 22 (Line)
Info : [ 0%] Meshing curve 21 (Line)
Info : [ 0%] Meshing curve 23 (Line)
Info : [ 0%] Meshing curve 24 (Line)
Info : [ 0%] Meshing curve 25 (Line)
Info : [ 0%] Meshing curve 26 (Line)
Info : [ 0%] Meshing curve 27 (Line)
Info : [ 0%] Meshing curve 28 (Line)
Info : [ 0%] Meshing curve 29 (Line)
Info : [ 0%] Meshing curve 30 (Line)
Info : [ 0%] Meshing curve 31 (Line)
Info : [ 0%] Meshing curve 32 (Line)
Info : [ 0%] Meshing curve 33 (Line)
Info : [ 0%] Meshing curve 34 (Line)
Info : Done meshing 1D (Wall 0.0270794s, CPU 0.079419s)
Info : Meshing 2D...
Info : [ 0%] Meshing surface 1 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 3 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 2 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 5 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 6 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 4 (Plane, Frontal-Delaunay)
Info : Done meshing 2D (Wall 0.0163829s, CPU 0.042735s)
Info : 504 nodes 1093 elements
Info : Writing 'mesh.msh'...
Info : Done writing 'mesh.msh'
Info : Writing '/tmp/tmpkkm9y5lr/mesh.msh'...
Info : Done writing '/tmp/tmpkkm9y5lr/mesh.msh'
Label-wise coarse resolution control#
An advantage of finite-volume and finite-element schemes is the ability for different nodes to have different characteristics lengths.
This simply achieved to first order here by supplying a resolutions dict with keys referencing the LayerStack names, and for value a second dict with keys resolution and DistMax / SizeMax (see gmsh documentation) which control, respectively, the characteristic length within a region and the dropoff away from interfaces with this region.
For example, to refine within the core only, one could use:
resolutions = {"core": {"resolution": 0.05, "distance": 0}}
mesh = get_mesh(
component=waveguide_trimmed,
type="uz",
xsection_bounds=[(4, -4), (4, 4)],
layer_stack=filtered_layer_stack,
filename=f"{filename}.msh",
background_tag="oxide",
background_padding=(2.0, 2.0, 2.0, 2.0),
resolutions=resolutions,
)
mesh = mesh_with_physicals(mesh, filename)
mesh = from_meshio(mesh)
mesh.draw().show()
Info : [ 0%] Union
Info : [ 10%] Union - Performing intersection of shapes
Info : [ 80%] Union - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference - Performing intersection of shapes
Info : [ 80%] Difference - Building splits of containers
Info : [ 0%] Fragments
Info : [ 10%] Fragments - Performing intersection of shapes
Info : [ 80%] Fragments - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference
Info : [ 90%] Difference
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Filling splits of edges
Info : [ 80%] Fragments - Splitting faces
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference - Adding holes
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Building splits of containers
Info : [ 80%] Fragments - Splitting faces
Info : Meshing 1D...
Info : [ 0%] Meshing curve 1 (Line)
Info : [ 0%] Meshing curve 2 (Line)
Info : [ 0%] Meshing curve 3 (Line)
Info : [ 0%] Meshing curve 4 (Line)
Info : [ 0%] Meshing curve 5 (Line)
Info : [ 0%] Meshing curve 6 (Line)
Info : [ 0%] Meshing curve 7 (Line)
Info : [ 0%] Meshing curve 8 (Line)
Info : [ 0%] Meshing curve 15 (Line)
Info : [ 0%] Meshing curve 16 (Line)
Info : [ 0%] Meshing curve 17 (Line)
Info : [ 0%] Meshing curve 18 (Line)
Info : [ 0%] Meshing curve 19 (Line)
Info : [ 0%] Meshing curve 20 (Line)
Info : [ 0%] Meshing curve 21 (Line)
Info : [ 0%] Meshing curve 22 (Line)
Info : [ 50%] Meshing curve 23 (Line)
Info : [ 50%] Meshing curve 24 (Line)
Info : [ 50%] Meshing curve 25 (Line)
Info : [ 50%] Meshing curve 26 (Line)
Info : [ 50%] Meshing curve 27 (Line)
Info : [ 50%] Meshing curve 28 (Line)
Info : [ 70%] Meshing curve 29 (Line)
Info : [ 70%] Meshing curve 30 (Line)
Info : [ 70%] Meshing curve 31 (Line)
Info : [ 70%] Meshing curve 32 (Line)
Info : [ 90%] Meshing curve 33 (Line)
Info : [ 90%] Meshing curve 34 (Line)
Info : Done meshing 1D (Wall 0.00195572s, CPU 0.002442s)
Info : Meshing 2D...
Info : [ 0%] Meshing surface 4 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 1 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 2 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 3 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 5 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 6 (Plane, Frontal-Delaunay)
Info : Done meshing 2D (Wall 0.0408996s, CPU 0.097417s)
Info : 788 nodes 1695 elements
Info : Writing 'mesh.msh'...
Info : Done writing 'mesh.msh'
Info : Writing '/tmp/tmps_30z91s/mesh.msh'...
Info : Done writing '/tmp/tmps_30z91s/mesh.msh'
Adding a dropoff at the interface:
resolutions = {"core": {"resolution": 0.05, "DistMax": 1, "SizeMax": 0.2}}
mesh = get_mesh(
component=waveguide_trimmed,
type="uz",
xsection_bounds=[(4, -4), (4, 4)],
layer_stack=filtered_layer_stack,
filename=f"{filename}.msh",
background_tag="oxide",
background_padding=(2.0, 2.0, 2.0, 2.0),
resolutions=resolutions,
)
mesh = mesh_with_physicals(mesh, filename)
mesh = from_meshio(mesh)
mesh.draw().show()
Info : [ 0%] Union
Info : [ 10%] Union - Performing intersection of shapes
Info : [ 80%] Union - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference - Performing intersection of shapes
Info : [ 80%] Difference - Building splits of containers
Info : [ 0%] Fragments
Info : [ 10%] Fragments - Performing intersection of shapes
Info : [ 80%] Fragments - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference
Info : [ 90%] Difference
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Filling splits of edges
Info : [ 80%] Fragments - Splitting faces
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference - Adding holes
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Building splits of containers
Info : [ 80%] Fragments - Splitting faces
Info : Meshing 1D...
Info : [ 0%] Meshing curve 1 (Line)
Info : [ 0%] Meshing curve 3 (Line)
Info : [ 0%] Meshing curve 2 (Line)
Info : [ 0%] Meshing curve 4 (Line)
Info : [ 0%] Meshing curve 5 (Line)
Info : [ 0%] Meshing curve 6 (Line)
Info : [ 0%] Meshing curve 7 (Line)
Info : [ 0%] Meshing curve 8 (Line)
Info : [ 0%] Meshing curve 15 (Line)
Info : [ 0%] Meshing curve 16 (Line)
Info : [ 0%] Meshing curve 17 (Line)
Info : [ 0%] Meshing curve 18 (Line)
Info : [ 0%] Meshing curve 19 (Line)
Info : [ 0%] Meshing curve 20 (Line)
Info : [ 0%] Meshing curve 21 (Line)
Info : [ 0%] Meshing curve 22 (Line)
Info : [ 0%] Meshing curve 23 (Line)
Info : [ 0%] Meshing curve 24 (Line)
Info : [ 0%] Meshing curve 25 (Line)
Info : [ 0%] Meshing curve 26 (Line)
Info : [ 0%] Meshing curve 27 (Line)
Info : [ 0%] Meshing curve 28 (Line)
Info : [ 0%] Meshing curve 29 (Line)
Info : [ 0%] Meshing curve 30 (Line)
Info : [ 0%] Meshing curve 31 (Line)
Info : [ 0%] Meshing curve 32 (Line)
Info : [ 0%] Meshing curve 33 (Line)
Info : [ 0%] Meshing curve 34 (Line)
Info : Done meshing 1D (Wall 0.0342376s, CPU 0.094385s)
Info : Meshing 2D...
Info : [ 0%] Meshing surface 2 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 3 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 4 (Plane, Frontal-Delaunay)
Info : [ 20%] Meshing surface 5 (Plane, Frontal-Delaunay)
Info : [ 20%] Meshing surface 6 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 1 (Plane, Frontal-Delaunay)
Info : Done meshing 2D (Wall 0.0244079s, CPU 0.054272s)
Info : 1054 nodes 2251 elements
Info : Writing 'mesh.msh'...
Info : Done writing 'mesh.msh'
Info : Writing '/tmp/tmpaw3_16va/mesh.msh'...
Info : Done writing '/tmp/tmpaw3_16va/mesh.msh'
Refining multiple elements simultaneously:
resolutions = {
"core": {"resolution": 0.05, "DistMax": 1, "SizeMax": 0.2},
"slab90": {"resolution": 0.02, "DistMax": 1, "SizeMax": 0.2},
"via_contact": {"resolution": 0.2},
"oxide": {"resolution": 1},
}
mesh = get_mesh(
component=waveguide_trimmed,
type="uz",
xsection_bounds=[(4, -4), (4, 4)],
layer_stack=filtered_layer_stack,
filename=f"{filename}.msh",
background_tag="oxide",
background_padding=(2.0, 2.0, 2.0, 2.0),
resolutions=resolutions,
)
mesh = mesh_with_physicals(mesh, filename)
mesh = from_meshio(mesh)
mesh.draw().show()
Info : [ 0%] Union
Info : [ 10%] Union - Performing intersection of shapes
Info : [ 80%] Union - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference - Performing intersection of shapes
Info : [ 80%] Difference - Building splits of containers
Info : [ 0%] Fragments
Info : [ 10%] Fragments - Performing intersection of shapes
Info : [ 80%] Fragments - Building splits of containers
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference
Info : [ 90%] Difference
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Filling splits of edges
Info : [ 80%] Fragments - Splitting faces
Info : [ 0%] Difference
Info : [ 10%] Difference
Info : [ 20%] Difference
Info : [ 30%] Difference
Info : [ 40%] Difference
Info : [ 50%] Difference
Info : [ 60%] Difference
Info : [ 70%] Difference - Filling splits of edges
Info : [ 80%] Difference - Making faces
Info : [ 90%] Difference - Adding holes
Info : [ 0%] Fragments
Info : [ 10%] Fragments
Info : [ 20%] Fragments
Info : [ 30%] Fragments
Info : [ 40%] Fragments
Info : [ 50%] Fragments
Info : [ 60%] Fragments
Info : [ 70%] Fragments - Building splits of containers
Info : [ 80%] Fragments - Splitting faces
Info : Meshing 1D...
Info : [ 0%] Meshing curve 1 (Line)
Info : [ 0%] Meshing curve 4 (Line)
Info : [ 0%] Meshing curve 2 (Line)
Info : [ 0%] Meshing curve 5 (Line)
Info : [ 0%] Meshing curve 3 (Line)
Info : [ 0%] Meshing curve 6 (Line)
Info : [ 0%] Meshing curve 7 (Line)
Info : [ 0%] Meshing curve 8 (Line)
Info : [ 0%] Meshing curve 15 (Line)
Info : [ 0%] Meshing curve 16 (Line)
Info : [ 0%] Meshing curve 17 (Line)
Info : [ 30%] Meshing curve 18 (Line)
Info : [ 40%] Meshing curve 19 (Line)
Info : [ 40%] Meshing curve 20 (Line)
Info : [ 40%] Meshing curve 21 (Line)
Info : [ 50%] Meshing curve 22 (Line)
Info : [ 50%] Meshing curve 23 (Line)
Info : [ 60%] Meshing curve 24 (Line)
Info : [ 60%] Meshing curve 25 (Line)
Info : [ 60%] Meshing curve 26 (Line)
Info : [ 70%] Meshing curve 27 (Line)
Info : [ 70%] Meshing curve 28 (Line)
Info : [ 70%] Meshing curve 29 (Line)
Info : [ 70%] Meshing curve 30 (Line)
Info : [ 80%] Meshing curve 31 (Line)
Info : [ 80%] Meshing curve 32 (Line)
Info : [ 90%] Meshing curve 33 (Line)
Info : [ 90%] Meshing curve 34 (Line)
Info : Done meshing 1D (Wall 0.0360117s, CPU 0.121951s)
Info : Meshing 2D...
Info : [ 0%] Meshing surface 1 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 3 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 2 (Plane, Frontal-Delaunay)
Info : [ 0%] Meshing surface 4 (Plane, Frontal-Delaunay)
Info : [ 20%] Meshing surface 5 (Plane, Frontal-Delaunay)
Info : [ 20%] Meshing surface 6 (Plane, Frontal-Delaunay)
Info : Done meshing 2D (Wall 0.170126s, CPU 0.286213s)
Info : 8395 nodes 17685 elements
Info : Writing 'mesh.msh'...
Info : Done writing 'mesh.msh'
Info : Writing '/tmp/tmpbkwmg4nu/mesh.msh'...
Info : Done writing '/tmp/tmpbkwmg4nu/mesh.msh'
