Summatory of domains

mesh
1

Hi, i have saved 69 domains in a list named \textit{domain}, thats because i want to create a mesh of the new domain:

golfo = domain[0] - domain[1] - ... - domain[68]

I tried writing

golfo=domain[0]
for i in range(1,len(domain)):
    golfo = golfo-domain[i]
mesh = generate_mesh(golfo,1)
plt.figure(figsize=(26,29))
plot(mesh)

but it doesn’t work (my kernel dies). Is there a way to write the summatory in a shortish manner instead of the explicit one above? Thanks!

2

I would strongly suggest not to use mshr to generate meshes, as it is no longer maintained. I would for instance recommend gmsh, which has a large variety of tutorials, both for gui and python tutorial · master · gmsh / gmsh · GitLab