1

Hi everyone.
I am looking to generate some simple 2D meshes (triangles filling a big rectangle with some smaller square holes). I used for years the DistMesh for Matlab. Now I am trying to use the Julia port. There are several that do not seem to be updated or to have all the functionality (e.g. Distmesh2d.jl does not output the connectivity of the elements).
I am looking for some examples of DistMesh.jl in 2D or other mesh generator.
Any help?
Kind regards.
M.

2

Are you aware that Meshes.jl can discretize geometries like Box? Check the docs and search for discretize or simplexify

1 Like
3

DelaunayTriangulation.jl doesn’t work exactly the same as DistMesh.jl, but if you wanted to triangulate for example (as you say) a rectangle with holes, you could use it like

using DelaunayTriangulation, CairoMakie
outer = [(0.0, 0.0), (1.0, 0.0), (1.0, 1.0), (0.0, 1.0), (0.0, 0.0)]
inner_1 = [(0.2, 0.2), (0.2, 0.5), (0.5, 0.5), (0.5, 0.2), (0.2, 0.2)]
inner_2 = [(0.7, 0.2), (0.7, 0.5), (0.9, 0.5), (0.9, 0.2), (0.7, 0.2)]
boundary = [[outer], [inner_1], [inner_2]]
boundary_nodes, points = convert_boundary_points_to_indices(boundary)
tri = triangulate(points; boundary_nodes)

# Unrefined 
fig = Figure()
ax = Axis(fig[1, 1], width = 400, height = 400, title = "Unrefined")
triplot!(ax, tri)

# Refined 
refine!(tri; max_area = 0.01) # you can also do custom constraints with the custom_constraint keyword 
ax = Axis(fig[1, 2], width = 400, height = 400, title = "Unrefined")
triplot!(ax, tri)
resize_to_layout!(fig)
fig

mXd8LkQ
mXd8LkQ748×392 36.1 KB

The docs give plenty of examples of working with the tri output (e.g. working with the connectivities)