Skip to content
Merged
Show file tree
Hide file tree
Changes from 5 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions NEWS.md
Original file line number Diff line number Diff line change
Expand Up @@ -7,6 +7,10 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0

## [Unreleased]

### Fixed

- Make `compute_facet_owners(...)` more flexible, allowing the user to provide a function to select the owner from neighboring cells. Since PR[#1291](https://github.com/gridap/Gridap.jl/pull/1291).

## [0.20.5] - 2026-04-28

Comment thread
amartinhuertas marked this conversation as resolved.
### Fixed
Expand Down
2 changes: 1 addition & 1 deletion Project.toml
Original file line number Diff line number Diff line change
@@ -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"
Expand Down
7 changes: 3 additions & 4 deletions src/FESpaces/Pullbacks.jl
Original file line number Diff line number Diff line change
Expand Up @@ -188,21 +188,20 @@ 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}, select_nbor=first) where {Dc}
Comment thread
amartinhuertas marked this conversation as resolved.
Outdated
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)
@check !isempty(facet_cells) "Facet $facet has no adjacent cells"
owners[facet] = select_nbor(facet_cells)
Comment thread
amartinhuertas marked this conversation as resolved.
Outdated
end

return owners
end


#################
# DOFScalingMap #
#################
Expand Down
Loading