On 7/12/21 4:18 AM, Sylvain Mathonnière wrote:
_*Question:*_
How can I modify the cells worker loop (on dof1) initialisation to have a
second iterator linked with a second dof_handler looping at the same time ? or
how can I add a concurrentely running for loop to the cell_worker with a
different iterator ?
I checked in step 16 and 47 as well as in the doc for MeshWorkers and I did
not find a solution to this problem. I hope the question is clear enough, do
not hesitate to inquires more details if I am missing something important.
Sylvain,
it's not going to be possible to hack this into MeshWorker::mesh_loop() in an
easy way. But what you can do is this: Inside your assembler, whenever you are
given a cell iterator pointing to a cell in dof_handler_1, you need to find
the corresponding iterator for dof_handler_2. You can do this in the following
way (up to typos and misspelled function names):
DoFHandler::active_cell_iterator cell_2 (cell_1->get_triangulation(),
cell_1->level(),
cell_1->index(),
cell_1->get_dof_handler());
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/fa539e45-c082-29eb-b95e-61c602c95000%40colostate.edu.