From 91272fe40407274e8a0c9124c3752df05489bbed Mon Sep 17 00:00:00 2001 From: Nuno Nobre Date: Tue, 19 May 2026 16:54:41 +0100 Subject: [PATCH] Fix for when find_elements returns null elements --- src/mesh/mesh_function.C | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/src/mesh/mesh_function.C b/src/mesh/mesh_function.C index 1cde22000ef..7d3f474e6c2 100644 --- a/src/mesh/mesh_function.C +++ b/src/mesh/mesh_function.C @@ -405,6 +405,15 @@ void MeshFunction::discontinuous_value (const Point & p, // empty map for (const auto & element : candidate_element) { + if (!element) + { + // We'd better be in out_of_mesh_mode if we couldn't find an + // element in the mesh + libmesh_assert (_out_of_mesh_mode); + output[element] = _out_of_mesh_value; + continue; + } + const unsigned int dim = element->dim(); // define a temporary vector to store all values