From 35deee477e403365b8f21fbf3a9efdb452f23b16 Mon Sep 17 00:00:00 2001 From: "Alberto F. Martin" Date: Thu, 12 Oct 2023 10:23:40 +1100 Subject: [PATCH] Fix in order to have an _add_constraints() function which is consistent 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 --- src/FESpaces.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/FESpaces.jl b/src/FESpaces.jl index f605120..dc37615 100644 --- a/src/FESpaces.jl +++ b/src/FESpaces.jl @@ -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)