On 7/19/19 9:15 AM, 'Maxi Miller' via deal.II User Group wrote:
> I am trying to port example 15 to a MPI-supported version, using the
> geometric
> multigrid method as preconditioner (to get more familiar with the
> implementation). I followed example 50 while doing that, but currently I get
> the error
> An error occurred in line <3866> of file
> <~/Downloads/git-files/dealii/include/deal.II/dofs/dof_accessor.templates.h>
> in function
> void dealii::DoFCellAccessor<DoFHandlerType, lda>::get_dof_values(const
> InputVector&, ForwardIterator, ForwardIterator) const [with InputVector =
> dealii::TrilinosWrappers::MPI::Vector; ForwardIterator = double*;
> DoFHandlerType = dealii::DoFHandler<2, 2>; bool level_dof_access = true]
> The violated condition was:
> this->is_artificial() == false
> Additional information:
> Can't ask for DoF indices on artificial cells.
> as soon as I try to run my program using several MPI threads. I tried to
> guard
> that part using if (cell->level_subdomain_id() ==
> triangulation.locally_owned_subdomain()) as suggested in step-50, but that
> did
> not help. Why am I still trying to request data from artificial cells, even
> though I guarded that code part (at least I assume it guards it accordingly).
Are you sure you guarded the right place? Have you checked with a debugger
where the problem actually happens?
Best
WB
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c12c4887-a6ee-a3ad-d19e-5ca53d448afc%40colostate.edu.
