On 7/19/19 9:49 AM, 'Maxi Miller' via deal.II User Group wrote: > I can not use that function, else I will get > An error occurred in line <3543> of file > </opt/dealii/include/deal.II/grid/tria_accessor.templates.h> in function > bool dealii::CellAccessor<dim, spacedim>::is_locally_owned() const [with > int dim = 2; int spacedim = 2] > The violated condition was: > this->active() > Additional information: > is_locally_owned() can only be called on active cells!
Ah, I see. I missed that part. > In my multigrid-assembly-function I have to loop over all cells, not only the > currently active cells (as it is stated in example 50: /we don't just loop > over all active cells, but in fact all cells, active or not/). Thus I end up > in that function. Hm, that's a question for the parallel multigrid folks then. What do the parallel multigrid functions do? Could you just have your if() inside that function -- i.e., it is called on all cells, but does nothing if the cell doesn't allow you to do anything? 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/1a96dbbc-ae23-505a-9f66-6473d0ddf586%40colostate.edu.
