Dear Chengjiang,
Yes, you need to run reinit() again when the boundary id on certain
boundaries changes, because the boundary id gets extracted to the data
structures of matrix-free face integrals, and faces of the same ID will
get grouped together.
Best,
Martin
Am 01.08.24 um 04:56 schrieb Chengjiang Yin:
BTW, the changing boundary is not Dirichlet BC so that I do not need
to rebuild my AffineConstraints.
在2024年7月30日星期二 UTC+8 10:41:34<Chengjiang Yin> 写道：
Hi all,
I need to change my boundary id on part of the mesh boundary when
evolving the equations. And I also have a boundary operator that
acts on specific boundary_id (absorption boundary).
In that case, do I need to *reinit* my matrix free every time I
update the boundary_id? I use *tasks_parallel_scheme=none*. Does
the matrix free only cache the mapping info of boundary face or
does it pre partitioned the face it into different groups?
Regards,
Chengjiang Yin
--
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/0a311f3d-cf29-4582-9a59-30484d1be4d0n%40googlegroups.com
<https://groups.google.com/d/msgid/dealii/0a311f3d-cf29-4582-9a59-30484d1be4d0n%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/1a4dfbd1-ee85-456c-89ed-727b7e35edcc%40gmail.com.
