so for those who are interested, I found a partial solution in how
VoronoiDelaunay plots meshes. Basically, one can interrupt lines and
polygons using NaNs. For example:
points = Tuple{Float64, Float64}[]
for n = 1:size(T, 2)
p = [X[:, T[:, n]] X[:, T[1,n]]]
for m = 1:size(p, 2)
push!(points, tuple(p[1, m], p[2, m]))
end
push!(points, tuple(NaN, NaN))
end
compose(context(units=ub), polygon(points), fill(elcol), stroke(linecol),
linewidth(lwidth) )
This works well for me - for now. Though I am sure one can improve this
further.
Christoph