phiFEM 0.7.2.dev2__tar.gz → 0.8.0.dev1__tar.gz

{phifem-0.7.2.dev2 → phifem-0.8.0.dev1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.5
 Name: phiFEM
-Version: 0.7.2.dev2
+Version: 0.8.0.dev1
 Summary: phiFEM: convenience tools to use the phiFEM immersed boundary finite element method with dolfinx.
 Project-URL: Homepage, https://github.com/PhiFEM/phiFEM
 Project-URL: Issues, https://github.com/PhiFEM/phiFEM/issues

{phifem-0.7.2.dev2 → phifem-0.8.0.dev1}/demo/interface-elasticity/main.py

@@ -108,12 +108,18 @@ mesh = dfx.mesh.create_rectangle(
     MPI.COMM_WORLD, np.asarray(bbox).T, [nx, ny], cell_type
 )
 
+cell_type = mesh.topology.cell_name()
+detection_element = element(
+    "Lagrange",
+    cell_type,
+    1,
+)
+detection_space = dfx.fem.functionspace(mesh, detection_element)
+
 results = {"dof": [], "H10 relative error": [], "L2 relative error": []}
 for i in range(num_iterations):
-    x_ufl = ufl.SpatialCoordinate(mesh)
-    detection_levelset = levelset(x_ufl)
     cells_tags, facets_tags, _, d_bdry, _ = compute_tags_measures(
-        mesh, detection_levelset, detection_degree, box_mode=True
+        mesh, levelset, detection_space, box_mode=True
     )
 
     gdim = mesh.geometry.dim

{phifem-0.7.2.dev2 → phifem-0.8.0.dev1}/demo/neumann/square/main.py

@@ -59,14 +59,14 @@ detection_levelset_h.interpolate(detection_levelset)
 
 if mesh_type == "bg":
     cells_tags, facets_tags, _, ds_bdy, _ = compute_tags_measures(
-        bg_mesh, detection_levelset_h, 1, box_mode=True
+        bg_mesh, detection_levelset_h, bg_levelset_space, box_mode=True
     )
     mesh = bg_mesh
     # To get the ds measure on the bounday of Omega_h
     ds = ds_bdy(100)
 elif mesh_type == "sub":
     cells_tags, facets_tags, mesh, _, _ = compute_tags_measures(
-        bg_mesh, detection_levelset_h, 1, box_mode=False
+        bg_mesh, detection_levelset_h, bg_levelset_space, box_mode=False
     )
     ds = ufl.Measure("ds", domain=mesh)