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

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.5
 Name: phiFEM
-Version: 0.7.2.dev1
+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
@@ -11,8 +11,8 @@ License-File: COPYING
 License-File: COPYING.lesser
 Classifier: Operating System :: OS Independent
 Classifier: Programming Language :: Python :: 3
-Requires-Python: >=3.11
-Requires-Dist: fenics-dolfinx>=0.9.0
+Requires-Python: >=3.9
+Requires-Dist: fenics-dolfinx==0.9.0
 Description-Content-Type: text/markdown
 
 <p align="center" width="100%">
@@ -118,7 +118,7 @@ To run the demos first clone the repository.
 The author Raphaël Bulle guarantees that his contributions to this project are generative-ai-free.
 
 <p align="center" width="50%">
-    <img width="30%" src="https://github.com/PhiFEM/phiFEM/raw/main/logo/sw_wbg.png">
+    <img width="30%" src="https://github.com/PhiFEM/phiFEM/raw/main/logo/image-572.png">
 </p>
 
 For more information, see [here (in french)](https://doi.org/10.58079/1591s).
@@ -133,4 +133,4 @@ You should have received a copy of the GNU Lesser General Public License along w
 
 ## Logo credit
 
-Frédérique Lecourtier ([https://mimesis.inria.fr/members/frederique-lecourtier/](https://mimesis.inria.fr/members/frederique-lecourtier/))
+Frédérique Lecourtier ([https://mimesis.inria.fr/members/frederique-lecourtier/](https://mimesis.inria.fr/members/frederique-lecourtier/))

{phifem-0.7.2.dev1 → phifem-0.8.0.dev1}/README.md

@@ -101,7 +101,7 @@ To run the demos first clone the repository.
 The author Raphaël Bulle guarantees that his contributions to this project are generative-ai-free.
 
 <p align="center" width="50%">
-    <img width="30%" src="https://github.com/PhiFEM/phiFEM/raw/main/logo/sw_wbg.png">
+    <img width="30%" src="https://github.com/PhiFEM/phiFEM/raw/main/logo/image-572.png">
 </p>
 
 For more information, see [here (in french)](https://doi.org/10.58079/1591s).
@@ -116,4 +116,4 @@ You should have received a copy of the GNU Lesser General Public License along w
 
 ## Logo credit
 
-Frédérique Lecourtier ([https://mimesis.inria.fr/members/frederique-lecourtier/](https://mimesis.inria.fr/members/frederique-lecourtier/))
+Frédérique Lecourtier ([https://mimesis.inria.fr/members/frederique-lecourtier/](https://mimesis.inria.fr/members/frederique-lecourtier/))

{phifem-0.7.2.dev1 → 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.dev1 → phifem-0.8.0.dev1}/demo/neumann/square/main.py

@@ -58,13 +58,15 @@ detection_levelset_h = dfx.fem.Function(bg_levelset_space)
 detection_levelset_h.interpolate(detection_levelset)
 
 if mesh_type == "bg":
-    cells_tags, facets_tags, _, ds, _ = compute_tags_measures(
-        bg_mesh, detection_levelset_h, 1, box_mode=True
+    cells_tags, facets_tags, _, ds_bdy, _ = compute_tags_measures(
+        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
+    cells_tags, facets_tags, mesh, _, _ = compute_tags_measures(
+        bg_mesh, detection_levelset_h, bg_levelset_space, box_mode=False
     )
     ds = ufl.Measure("ds", domain=mesh)
 

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

@@ -65,12 +65,14 @@ detection_levelset_h = dfx.fem.Function(bg_levelset_space)
 detection_levelset_h.interpolate(detection_levelset)
 
 if mesh_type == "bg":
-    cells_tags, facets_tags, _, ds, _ = compute_tags_measures(
+    cells_tags, facets_tags, _, ds_bdy, _ = compute_tags_measures(
         bg_mesh, detection_levelset_h, 1, 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(
+    cells_tags, facets_tags, mesh, _, _ = compute_tags_measures(
         bg_mesh, detection_levelset_h, 1, box_mode=False
     )
     ds = ufl.Measure("ds", domain=mesh)

{phifem-0.7.2.dev1 → phifem-0.8.0.dev1}/demo/strong-dirichlet/flower/main.py

@@ -56,12 +56,14 @@ detection_levelset_h = dfx.fem.Function(bg_levelset_space)
 detection_levelset_h.interpolate(detection_levelset)
 
 if mesh_type == "bg":
-    cells_tags, facets_tags, _, ds, _ = compute_tags_measures(
+    cells_tags, facets_tags, _, ds_bdy, _ = compute_tags_measures(
         bg_mesh, detection_levelset_h, 1, 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(
+    cells_tags, facets_tags, mesh, _, _ = compute_tags_measures(
         bg_mesh, detection_levelset_h, 1, box_mode=False
     )
     ds = ufl.Measure("ds", domain=mesh)
@@ -140,11 +142,10 @@ pc = ksp.getPC()
 pc.setType("lu")
 
 # When solving on the background mesh, we need mumps to handle the null space of the matrix
-if mesh_type == "bg":
-    pc.setFactorSolverType("mumps")
-    pc.setFactorSetUpSolverType()
-    pc.getFactorMatrix().setMumpsIcntl(icntl=24, ival=1)
-    pc.getFactorMatrix().setMumpsIcntl(icntl=25, ival=0)
+pc.setFactorSolverType("mumps")
+pc.setFactorSetUpSolverType()
+pc.getFactorMatrix().setMumpsIcntl(icntl=24, ival=1)
+pc.getFactorMatrix().setMumpsIcntl(icntl=25, ival=0)
 
 """
 ===============================

{phifem-0.7.2.dev1 → phifem-0.8.0.dev1}/demo/weak-dirichlet/flower/main.py

@@ -53,7 +53,7 @@ detection_levelset_h = dfx.fem.Function(bg_levelset_space)
 detection_levelset_h.interpolate(detection_levelset)
 
 if mesh_type == "bg":
-    cells_tags, facets_tags, _, ds, _ = compute_tags_measures(
+    cells_tags, facets_tags, _, ds_bdy, _ = compute_tags_measures(
         bg_mesh,
         detection_levelset_h,
         1,
@@ -61,8 +61,10 @@ if mesh_type == "bg":
         single_layer_cut=True,
     )
     mesh = bg_mesh
+    # To get the ds measure on the boundary of Omega_h
+    ds = ds_bdy(100)
 elif mesh_type == "sub":
-    cells_tags, facets_tags, mesh, _, _, _ = compute_tags_measures(
+    cells_tags, facets_tags, mesh, _, _ = compute_tags_measures(
         bg_mesh,
         detection_levelset_h,
         1,

phifem-0.8.0.dev1/phifem-env.yml

@@ -0,0 +1,9 @@
+name: phifem
+channels:
+  - conda-forge
+dependencies:
+  - python>=3.9.0
+  - numpy>=1.21.0
+  - fenics-dolfinx==0.9.0
+  - pip:
+    - phifem==0.7.0