xfem 2.1.2605.dev2__tar.gz → 2.1.2605.dev3__tar.gz

{xfem-2.1.2605.dev2 → xfem-2.1.2605.dev3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: xfem
-Version: 2.1.2605.dev2
+Version: 2.1.2605.dev3
 Summary: (ngs)xfem is an Add-on library to Netgen/NGSolve for unfitted/cut FEM.
 Author-Email: Christoph Lehrenfeld <lehrenfeld@math.uni-goettingen.de>
 Classifier: License :: OSI Approved :: GNU Lesser General Public License v3 (LGPLv3)

{xfem-2.1.2605.dev2 → xfem-2.1.2605.dev3}/cutint/straightcutrule.cpp

@@ -643,7 +643,18 @@ namespace xintegration
       if (dt == IF)
       {
         auto ir_interface  = new (lh) IntegrationRule(quad_untrafo.Size(),lh);
-        if (DIM == 1) TransformQuadUntrafoToIRInterface<1>(quad_untrafo, trafo, lset, ir_interface, spacetime_mode, tval);
+        if (et == ET_SEGM && DIM == 2)
+        {
+          // BND segment in a 2D mesh: the interface (IF) is a co-dim-2 point.
+          // Use point measure: weight = 1/GetMeasure() so that mip.GetMeasure()*weight = 1.
+          for (int i = 0; i < quad_untrafo.Size(); i++)
+          {
+            MappedIntegrationPoint<1,2> mip(quad_untrafo[i], trafo);
+            (*ir_interface)[i] = IntegrationPoint(quad_untrafo[i].Point(),
+                                                  quad_untrafo[i].Weight() / mip.GetMeasure());
+          }
+        }
+        else if (DIM == 1) TransformQuadUntrafoToIRInterface<1>(quad_untrafo, trafo, lset, ir_interface, spacetime_mode, tval);
         else if (DIM == 2) TransformQuadUntrafoToIRInterface<2>(quad_untrafo, trafo, lset, ir_interface, spacetime_mode, tval);
         else TransformQuadUntrafoToIRInterface<3>(quad_untrafo, trafo, lset, ir_interface, spacetime_mode, tval);
         ir = ir_interface;

{xfem-2.1.2605.dev2 → xfem-2.1.2605.dev3}/python/_version.py

@@ -18,7 +18,7 @@ version_tuple: tuple[int | str, ...]
 commit_id: str | None
 __commit_id__: str | None
 
-__version__ = version = '2.1.2605.dev2'
-__version_tuple__ = version_tuple = (2, 1, 2605, 'dev2')
+__version__ = version = '2.1.2605.dev3'
+__version_tuple__ = version_tuple = (2, 1, 2605, 'dev3')
 
-__commit_id__ = commit_id = 'g8d6a9079d'
+__commit_id__ = commit_id = 'ge7bf5996f'

{xfem-2.1.2605.dev2 → xfem-2.1.2605.dev3}/tests/pytests/test_bndcut.py

@@ -162,6 +162,61 @@ def test_bfi_assembly_2d_partial():
     assert abs(val - cut) < 1e-2, f"Expected ~{cut}, got {val}"
 
 
+# ---------------------------------------------------------------------------
+# IF+BND (co-dimension-2) tests
+# ---------------------------------------------------------------------------
+
+@pytest.mark.parametrize("cut", [0.3, 0.5, 0.7])
+def test_integrate_if_bnd_2d_count(cut):
+    """Integrate(1 * dCut(lset, IF, vb=BND)) counts boundary–interface intersections.
+
+    With lset = x - cut the interface {lset=0} is the vertical line x=cut.
+    It intersects the boundary of the unit square in exactly 2 points:
+    (cut, 0) on the bottom edge and (cut, 1) on the top edge.
+    The point measure of 2 discrete points is 2.
+    """
+    mesh = make_2d_mesh()
+    lset = GridFunction(H1(mesh, order=1))
+    InterpolateToP1(x - cut, lset)
+
+    val = Integrate(CoefficientFunction(1) * dCut(lset, IF, vb=BND), mesh)
+    assert abs(val - 2.0) < 1e-10, f"Expected 2.0 (2 intersection points), got {val}"
+
+
+def test_integrate_if_bnd_2d_function():
+    """Integrating a function over IF+BND evaluates it at the intersection points.
+
+    lset = x - 0.5 → intersections at (0.5, 0) and (0.5, 1).
+    ∫ x dσ₀ = 0.5 + 0.5 = 1.0
+    ∫ y dσ₀ = 0.0 + 1.0 = 1.0
+    """
+    mesh = make_2d_mesh()
+    lset = GridFunction(H1(mesh, order=1))
+    InterpolateToP1(x - 0.5, lset)
+
+    val_x = Integrate(x * dCut(lset, IF, vb=BND), mesh)
+    val_y = Integrate(y * dCut(lset, IF, vb=BND), mesh)
+    assert abs(val_x - 1.0) < 1e-10, f"∫x dσ₀ expected 1.0, got {val_x}"
+    assert abs(val_y - 1.0) < 1e-10, f"∫y dσ₀ expected 1.0, got {val_y}"
+
+
+def test_integrate_if_bnd_2d_definedon():
+    """definedon restricts IF+BND to a single boundary edge."""
+    mesh = make_2d_mesh()
+    lset = GridFunction(H1(mesh, order=1))
+    InterpolateToP1(x - 0.5, lset)
+
+    # Only the bottom edge should have one intersection point at (0.5, 0)
+    val_bot = Integrate(CoefficientFunction(1) * dCut(lset, IF, vb=BND,
+                        definedon=mesh.Boundaries("bottom")), mesh)
+    assert abs(val_bot - 1.0) < 1e-10, f"Expected 1 intersection on bottom, got {val_bot}"
+
+    # Only the left edge (x=0) has no intersection with x=0.5
+    val_left = Integrate(CoefficientFunction(1) * dCut(lset, IF, vb=BND,
+                         definedon=mesh.Boundaries("left")), mesh)
+    assert abs(val_left - 0.0) < 1e-10, f"Expected 0 intersections on left, got {val_left}"
+
+
 def test_bfi_definedon_bitarray():
     """definedon can be supplied as a BitArray of boundary elements."""
     cut = 0.4