On 6/8/24 07:51, Simon wrote:
Before refinement, dofHandler stored 138 DoFs, after refinement 202.
I could not reproduce the error on a minimal example.
However, it is not clear to me why
make_hanging_constraints() reports an error associated with dof index 167,
if I cleared the constraints object before hand and can be certain that my
DoFHandler has 202 DoFs.
So the sequence of distributing DoFs, followed by clearing constraints and
calling make hanging node constraints should work independently of how my
program arrived at these lines of code. Is this not true?
Any general hints as for the problem source and narrowing down the problem?
Did you initialize the AffineConstraints object with index sets for the
locally owned and/or locally relevant DoFs? If so, you will have to
re-initialize it with the now larger/different number of DoFs you have after
mesh refinement.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/1350d410-7ae1-494f-9adf-50fd323746bd%40colostate.edu.
