On 8/31/22 11:09, Lucas Myers wrote:
Thanks for the help! I printed out the number of constraints after each
of your suggestions, and after setting the hanging node constraints the
number of constraints is 0 (as it should be, I'm working with a 1D mesh
that has been globally refined some number of times), but after
interpolating the boundary values the number of constraints is only 1.
Given that dim = 1 and I'm trying to enforce Dirichlet boundaries, I
reckon there should be 2 constraints. Running the std::map version of
interpolate_boundary_values, I still only get 1 constraint.
Ah, that's a commonly asked question. In 2d and 3d, the default is for
the entire boundary to have boundary ID equal to zero. In 1d, it is zero
for the left end point and one for the right end point. So you have to
call interpolate_b_v() twice, once for each of the two boundary IDs.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/ddbc0863-92d4-b8ad-59ec-9e04773f9963%40colostate.edu.
