Hi Doug! On Wednesday, August 21, 2019 at 4:00:49 AM UTC+2, Doug wrote: > > 8134: void dealii::parallel::distributed::SolutionTransfer<dim, > VectorType, DoFHandlerType>::unpack_callback(const typename > dealii::parallel::distributed::Triangulation<dim, DoFHandlerType:: > space_dimension>::cell_iterator&, typename > dealii::parallel::distributed::Triangulation<dim, DoFHandlerType:: > space_dimension>::CellStatus, const > boost::iterator_range<__gnu_cxx::__normal_iterator<const char*, > std::vector<char, std::allocator<char> > > >&, std::vector<VectorType*>&) > [with int dim = 3; VectorType = > dealii::LinearAlgebra::distributed::Vector<double>; DoFHandlerType = > dealii::hp::DoFHandler<3, 3>; typename > dealii::parallel::distributed::Triangulation<dim, DoFHandlerType:: > space_dimension>::cell_iterator = > dealii::TriaIterator<dealii::CellAccessor<3, 3> >; typename > dealii::parallel::distributed::Triangulation<dim, DoFHandlerType:: > space_dimension>::CellStatus = dealii::Triangulation<3, 3>::CellStatus] >
It seems that the assertion fails only in the 3D case, which is quite interesting. Are you using the deal.II-9.1 release version, or do you work with the current master branch? I'll try to reproduce the error once I'm at my desk. It looks like either p::d::SolutionTransfer picks the wrong finite element while packing the solution on cells to be coarsened, or hp::DoFHandler chooses the wrong finite element for the parent cell in case of coarsening. I'll investigate. Marc -- 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/294c934e-5582-4c4a-a9e9-f895fa8639bb%40googlegroups.com.
