Open tamiko opened 8 months ago
Fix should be straightforward: We need to assemble the normal (and boundary mass) into a temporary vector with AffineConstraints::add_local_to_global()
and compress(VectorOperation::add)
after assembly.
Then, we can collect the normal and put it into our BoundaryMap.
The logic in
OfflineData<dim, Number>::construct_boundary_map
currently reads as follows:This logic only works in 2D and will fail in 3D when we encounter hanging nodes on boundaries. We might also miss to set up boundary degrees of freedom.