Strengthen/bug-fix compute_facet_owners

NEWS.md

@@ -5,8 +5,12 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
-## [Unreleased]
+## [0.20.6] - 2026-04-29
 
+### Fixed
+
+- Strengthen compute_facet_owners(...) so that the owner of a facet does not depend on how the cells shared by a facet are listed in the mesh data structure, but on the (local or global) cell identifiers of the cells sharing the facet. Since PR[#1291](https://github.com/gridap/Gridap.jl/pull/1291).
+
 ## [0.20.5] - 2026-04-28
 
 ### Fixed

Project.toml

@@ -1,7 +1,7 @@
 name = "Gridap"
 uuid = "56d4f2e9-7ea1-5844-9cf6-b9c51ca7ce8e"
 authors = ["Santiago Badia <santiago.badia@monash.edu>", "Francesc Verdugo <f.verdugo.rojano@vu.nl>", "Alberto F. Martin <alberto.f.martin@anu.edu.au>"]
-version = "0.20.5"
+version = "0.20.6"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"

src/FESpaces/Pullbacks.jl

@@ -188,17 +188,25 @@ function evaluate!(cache,k::NormalSignMap,reffe,facet_own_dofs,cell)
   return Diagonal(dof_sign)
 end
 
-function compute_facet_owners(model::DiscreteModel{Dc}) where {Dc}
+function compute_facet_owners(model::DiscreteModel{Dc}, local_to_global::Function=identity) where {Dc}
   topo = get_grid_topology(model)
   facet_to_cell = get_faces(topo, Dc-1, Dc)
 
   nfacets = num_faces(topo, Dc-1)
   owners = Vector{Int32}(undef, nfacets)
   for facet in 1:nfacets
     facet_cells = view(facet_to_cell, facet)
-    owners[facet] = first(facet_cells)
+    max_gid=-1
+    max_lid=-1
+    for cell in facet_cells
+        gid = local_to_global(cell)
+        if gid > max_gid
+          max_gid = gid
+          max_lid = cell
+        end
+    end
+    owners[facet] = max_lid
   end
-
   return owners
 end