On 12/2/21 7:10 AM, Мария Бронзова wrote:
Basically, now I only have one additional factor in my local_matrix (apart
from it's "skeleton") within the Assembly module - density22, which is a
complex value (0,1). When I remove it, the Dirichlet BC get applied as
expected through the MatrixTools::apply_boundary_values function. But when i
leave it in the system - boundary conditions are no longer correct after
solving the system.
Would you please give a hint, what is going on here?
I *think* what might be happening is that you are running into a bug that
someone else just a few weeks ago discovered:
https://github.com/dealii/dealii/issues/12916
The patch for it is here:
https://github.com/dealii/dealii/pull/13041
I realize that that must have been very frustrating to debug :-( Can you try
and see whether the issue you see goes away if you apply that patch and
recompile deal.II with it?
Best
Wolfgang
--
------------------------------------------------------------------------
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/91ae4824-6cf7-5764-d00e-3f9dccdb2742%40colostate.edu.