Hi, I am using gridap to solve helmholtz PDE. I have some passive elements which has fixed values of solution there. I have defined this elements as Dirichlet BC as below in two ways once fixing those elements surface and once by adding the outer boundary of those elements to the previous surface and i get different results. Is there any one who knows that I am doing right or the method that I am using is totally wrong? thank you

```
p_reffe = ReferenceFE(lagrangian, Float64, 0)
Q = TestFESpace(Ω, p_reffe, vector_type = Vector{Float64},
dirichlet_tags = ["NonDesignLoad","LoadLine"],
dirichlet_masks =[(true),(true)])
P1(x)=1.0
P2(x)=1.0
P = TrialFESpace(Q ,[P1,P2])
np = num_free_dofs(P) # Number of cells in design region (number of design parameters)
pf_reffe = ReferenceFE(lagrangian, Float64, 1)
Qf = TestFESpace(Ω, pf_reffe, vector_type = Vector{Float64},
dirichlet_tags = ["NonDesignLoad","LoadLine"],
dirichlet_masks =[(true),(true)])
sr = num_free_dofs(Qf)
Pf1(x)=1.0
Pf2(x)=1.0
Pf = TrialFESpace(Qf ,[Pf1,Pf2])
```

```
p_reffe = ReferenceFE(lagrangian, Float64, 0)
Q = TestFESpace(Ω, p_reffe, vector_type = Vector{Float64},
dirichlet_tags = ["NonDesignLoad"],
dirichlet_masks =[(true)])
P1(x)=1.0
P = TrialFESpace(Q ,[P1])
np = num_free_dofs(P) # Number of cells in design region (number of design parameters)
pf_reffe = ReferenceFE(lagrangian, Float64, 1)
Qf = TestFESpace(Ω, pf_reffe, vector_type = Vector{Float64},
dirichlet_tags = ["NonDesignLoad"],
dirichlet_masks =[(true)])
sr = num_free_dofs(Qf)
Pf1(x)=1.0
Pf = TrialFESpace(Qf ,[Pf1])
```