Skip to content

BUG: test_high_order_location_warped_interior_facet #4627

@qbisi

Description

@qbisi

Describe the bug
test_high_order_location_warped_interior_facet failed on the assertion
assert mesh.locate_cell([0.25, 0.605], tolerance=0.0001) == 3

Steps to Reproduce
Steps to reproduce the behavior:

pytest tests/firedrake/regression/test_locate_cell.py -k test_high_order_location_warped_interior_facet

Expected behavior

assert mesh.locate_cell([0.25, 0.605], tolerance=0.0001,cell_ignore=3) == 1
assert mesh.locate_cell([0.25, 0.605], tolerance=0.0001,cell_ignore=1) == None

However, the first assert success, while the second failed (assert 3 == None)

A further test show that [0.25, 0.6125](as supposed) is the edge point for cell 1

assert mesh.locate_cell([0.25, 0.6127], tolerance=0.0001, cell_ignore=3) == None
assert mesh.locate_cell([0.25, 0.6125], tolerance=0.0001, cell_ignore=3) == 1

while [0.25,0.5](remain unchnaged) is the edge point for cell 3

assert mesh.locate_cell([0.25, 0.5], tolerance=0.0001, cell_ignore=1) == 3
assert mesh.locate_cell([0.25, 0.499], tolerance=0.0001, cell_ignore=1) == None

BTW [1/6,0.6],[1/3,0.6] is the edge point for both cell 1,3 as supposed.

Error message
FAILED tests/firedrake/regression/test_locate_cell.py::test_high_order_location_warped_interior_facet - AssertionError: assert np.int32(3) == 1

Environment:

  • OS: nixos-25.05
  • Python version: python3.13

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions