Skip to content

Commit

Permalink
Fix in order to have an _add_constraints() function which is consistent
Browse files Browse the repository at this point in the history
with the logic of the local Gridap FESpace constructor

Before the fix, whenever we were not passing the conformity kw-arg,
we were building a H1-conforming space only for the conforming mesh
interfaces and we were not gluing at all the non-conforming ones
  • Loading branch information
amartinhuertas committed Oct 11, 2023
1 parent 2de20a1 commit 35deee4
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/FESpaces.jl
Original file line number Diff line number Diff line change
Expand Up @@ -775,7 +775,7 @@ function _add_constraints(model::GridapDistributed.DistributedDiscreteModel{Dc},
spaces_wo_constraints;
conformity=nothing,
kwargs...) where {Dc}
if (conformity!=nothing && conformity!=:L2)
if (conformity==nothing || conformity!=:L2)
ref_constraints = _build_constraint_coefficients_matrix_in_ref_space(Dc, reffe)

face_subface_ldof_to_cell_ldof = Vector{Vector{Vector{Vector{Int32}}}}(undef, Dc-1)
Expand Down

0 comments on commit 35deee4

Please sign in to comment.