Replace the mesh in the cashocs as solver (for shape optimization) wi…

…th a custom mesh.

demos/documented/cashocs_as_solver/shape_solver/demo_shape_solver.py

@@ -25,15 +25,7 @@
 
 config = cashocs.load_config("./config.ini")
 
-meshlevel = 15
-degree = 1
-dim = 2
-mesh = UnitDiscMesh.create(MPI.comm_world, meshlevel, degree, dim)
-dx = Measure("dx", mesh)
-boundary = CompiledSubDomain("on_boundary")
-boundaries = MeshFunction("size_t", mesh, dim=1)
-boundary.mark(boundaries, 1)
-ds = Measure("ds", mesh, subdomain_data=boundaries)
+mesh, subdomains, boundaries, dx, ds, dS = cashocs.import_mesh("./mesh/mesh.xdmf")
 
 V = FunctionSpace(mesh, "CG", 1)
 u = Function(V)

demos/documented/cashocs_as_solver/shape_solver/mesh/mesh.geo

@@ -0,0 +1,10 @@
+lc = 7.5e-2;
+
+SetFactory("OpenCASCADE");
+Disk(1) = {0, 0, 0, 1, 1};
+
+Mesh.MeshSizeMax = lc;
+Mesh.MeshSizeMin = 0.0;
+
+Physical Curve(1) = {1};
+Physical Surface(1) = {1};