Alex,
I was wondering if anyone might have some insight into why this isn't
working. I've also reviewed the videos on applying non-Dirichlet
boundary conditions in deal.II, and everything seems to align with the code.
Despite that, the issue still persists!
Thank you in advance for your time and help.
On Thursday, March 20, 2025 at 9:51:50 AM UTC+1 Alex wrote:
Dear Prof. Wolfgang,
Thank you very much for your kind response.
I sincerely appreciate your insights. However, I believe this aspect
has already been addressed in step-44 through the conditions
mentioned below.
// Since the constraints are different at different Newton
iterations, we
// need to clear the constraints matrix and completely rebuild
// it. However, after the first iteration, the constraints
remain the same
// and we can simply skip the rebuilding step if we do not
clear it.
if (it_nr > 1)
return;
const bool apply_dirichlet_bc = (it_nr == 0);
To verify this, I printed the values of the constraints at each
iteration. Based on my observations, the constraints are populated
with values only at Iteration "0," and in all subsequent iterations,
they remain zero.
This is the starting point on how to debug this. If I understand you
right, you expect the constraints object to contain things also at
iterations other than the first one, but it does not. You you have
something that seems wrong to you. Now work backward to find out why the
object does not contain anything. Where is the place where it's supposed
to be filled? Why is that not happening?
Best
W.
--
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 visit
https://groups.google.com/d/msgid/dealii/4cbc0ebd-8ca4-488f-92c1-c33cdcf47440%40colostate.edu.
