Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Shape derivative #882

Open
wants to merge 17 commits into
base: master
Choose a base branch
from
Prev Previous commit
Next Next commit
remove
  • Loading branch information
Connor Mallon committed Mar 2, 2023
commit 67e64d3deacf972c1b9698f2fa73a132d3dd10ed
2 changes: 0 additions & 2 deletions src/FESpaces/UnconstrainedFESpaces.jl
Original file line number Diff line number Diff line change
@@ -65,15 +65,13 @@ zero_dirichlet_values(f::UnconstrainedFESpace) = allocate_vector(f.vector_type,n
get_dirichlet_dof_tag(f::UnconstrainedFESpace) = f.dirichlet_dof_tag

function scatter_free_and_dirichlet_values(f::UnconstrainedFESpace,free_values,dirichlet_values)
#=
@check eltype(free_values) == eltype(dirichlet_values) """\n
The entries stored in free_values and dirichlet_values should be of the same type.

This error shows up e.g. when trying to build a FEFunction from a vector of integers
if the Dirichlet values of the underlying space are of type Float64, or when the
given free values are Float64 and the Dirichlet values ComplexF64.
"""
=#
cell_dof_ids = get_cell_dof_ids(f)
lazy_map(Broadcasting(PosNegReindex(free_values,dirichlet_values)),cell_dof_ids)
end