Hi Simon

you can write function of setup_system like other tutorial setp, after
refine mesh, rerun function of setup_system().

Best
Chen

Simon <[email protected]> 于2024年6月10日周一 14:13写道:

> No, the program that produced the error is serial not parallel.
>
> Best,
> Simon
>
> Wolfgang Bangerth schrieb am Montag, 10. Juni 2024 um 04:30:33 UTC+2:
>
>> 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/66d6a3fd-eb61-45ab-9ca1-509ed97bdc39n%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/66d6a3fd-eb61-45ab-9ca1-509ed97bdc39n%40googlegroups.com?utm_medium=email&utm_source=footer>
> .
>

-- 
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/CADr3OdKDZqxoreYS44wst-S5Vc%3DypSLe0PNvC5gCXtR0%2B-cp4A%40mail.gmail.com.

Reply via email to