Dear Juan, My best guess is that, at some point in the loop over all active cells, you are also doing some computations involving the neighbor of each cell. TriaAccessor::neighbor() returns the neighboring cell which has, at most, the same level number as the current cell: i.e., if the neighbor cell has been refined but the current cell has not, then trying to get the degrees of freedom from the neighbor will throw exactly this exception since the cell returned by the call to neighbor() is inactive.
This is only a guess. Could you show us the part of your code that triggers this exception? That will help us give you a better answer :) Thanks, David On Fri, Dec 6, 2019 at 3:53 AM Juan Felipe Giraldo <[email protected]> wrote: > > > Hi all! > > I am currently working on the linear advection problem by using a DG > formulation, similar to the problem in step 12 > https://www.dealii.org/current/doxygen/deal.II/step_12.html, but using a > mixed formulation as it is shown in step 20 > https://www.dealii.org/current/doxygen/deal.II/step_20.html. I can solve the > problem without problems using a global refinement, however, when I implement > an adaptative refinement I have the following problem after the 1st > refinement cycle > > An error occurred in line <3792> of file > </home/arsen/Documents/instaladores/dealii/include/deal.II/dofs/dof_accessor.templates.h> > in function > void dealii::DoFCellAccessor<DoFHandlerType, > lda>::get_dof_indices(std::vector<unsigned int>&) const [with DoFHandlerType > = dealii::DoFHandler<2, 2>; bool level_dof_access = false] > The violated condition was: > this->active() > Additional information: > get_dof_indices() only works on active cells. > > > I have seen that some people had this problem a couple of years ago, but I > have not seen a proper solution to solve it. > > Anyone knows how can I solve it? > > Thank you so much, > > Regards, > > Juan Felipe Giraldo > > -- > 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/38d23dcb-3028-45eb-a539-737e72c597b0%40googlegroups.com. -- 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/CABrTbYS7_D7cz8BVjObSSy3a6zAHMSCrwn75cj79XBVnOicCpw%40mail.gmail.com.
