Hello Prof. Bangerth,
Thanks for your reply.
"you show the backtrace, but what is the actual error message you observe?"
Did you mean what is the thrown error message in the command line?
input: mpirun -np 2 ./v1.0
output :
corrupted double-linked list
Decompressing section '.debug_rnglists'...inflated from 6492533 to 22203369
bytes
it stopped at this line
this->local_assemble_volume(cell, scratch, copy_data);
in the worker function
auto worker =
[this](const typename DoFHandler<dim>::active_cell_iterator &cell,
ScratchData & scratch,
PerTaskData & copy_data) {
this->local_assemble_volume(cell, scratch, copy_data);
};
Did I understand your question clear?
Lex
On Wednesday, October 25, 2023 at 7:35:05 AM UTC-7 Wolfgang Bangerth wrote:
>
> Lex,
> you show the backtrace, but what is the actual error message you observe?
> Best
> WB
>
> On 10/24/23 23:08, Lex Lee wrote:
> > *** Caution: EXTERNAL Sender ***
> >
> > Hello deal.ii develpoers,
> >
> >
> > In my FSI problem, the parallization code solved 2 sets of equations
> using the
> > "workstream" class interfacing with PETSc. Equation 1 will get known
> values of
> > variables in equation 2, and so does equation 2.
> >
> > The solution for equation 2 is initiallized in this way.
> >
> > hp_index_set = dof_handler.locally_owned_dofs();
> > hp_relevant_set = DoFTools::extract_locally_relevant_dofs(dof_handler);
> > solution.reinit(hp_index_set,
> > hp_relevant_set,
> > mpi_communicator);
> >
> >
> > In the "assemble_local" part, I used this line to get parts of the
> solution
> > values.
> > fe_values[extractor_displacement].get_function_gradients(solution,
> grad_u);
> >
> >
> > In the "workstream" part, there are worker, copier, and workstream
> function.
> >
> > PerTaskData cp;
> >
> > auto worker =
> > [this](const typename DoFHandler<dim>::active_cell_iterator &cell,
> > ScratchData & scratch,
> > PerTaskData &
> copy_data) {
> > this->local_assemble_volume(cell, scratch, copy_data);
> > };
> > auto copier = [this](const PerTaskData ©_data) {
> > this->copy_local_to_global_volume(copy_data);
> > };
> >
> > WorkStream::run(CellFilter(IteratorFilters::LocallyOwnedCell(),
> > volume_dof_handler.begin_active()),
> > CellFilter(IteratorFilters::LocallyOwnedCell(),
> > volume_dof_handler.end()),
> > worker,
> > copier,
> > sd,
> > cp);
> > volume_system_matrix.compress(VectorOperation::add);
> > volume_system_rhs.compress(VectorOperation::add);
> >
> >
> > Note that I used a cell filter to make the code run on locally owned
> cells,
> > instead of ghosted cells.However, an error will occur in the "get
> function
> > values" line. It failed to call back the function
> >
> > get_dof_values(values, local_values.begin(), local_values.end());
> >
> > in the source file. Refer to the picture below.
> >
> > The reason is deal.ii somehow tries to get values on ghosted cells,
> ignoring
> > the cell filter.
> >
> > Screenshot from 2023-10-24 21-34-45.png
> >
> > Screenshot from 2023-10-24 21-35-05.png
> > This puzzled me a lot for a long time. Any suggestion?
> >
> > Regards,
> > Lex
> >
> > --
> > The deal.II project is located at http://www.dealii.org/
> > <http://www.dealii.org/>
> > For mailing list/forum options, see
> > https://groups.google.com/d/forum/dealii?hl=en
> > <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]
> > <mailto:[email protected]>.
> > To view this discussion on the web visit
> >
> https://groups.google.com/d/msgid/dealii/d408514c-0f75-4448-8279-ac550830c6c1n%40googlegroups.com
>
> <
> https://groups.google.com/d/msgid/dealii/d408514c-0f75-4448-8279-ac550830c6c1n%40googlegroups.com?utm_medium=email&utm_source=footer
> >.
>
> --
> ------------------------------------------------------------------------
> 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/0870dd82-ef4c-461f-9872-c39df9469da7n%40googlegroups.com.
