Hi Daniel, Below I have attached a text file which contains the backtrace message. It looks like the error originates from dealii::KellyErrorEstimator<2,2>::estimate and could involve dealii::hp::FESubfaceValues and QProjector. Could this be related to the way the triangulated mesh is handled in the estimator?
Do you have any suggestions on how to debug this further? Thanks, Devansh. On Monday, 17 March 2025 at 21:37:02 UTC+5:30 d.arnd...@gmail.com wrote: > Devanash, > > What does a backtrace look like when running in a debugger? > > Best, > Daniel > > On Mon, Mar 17, 2025 at 9:42 AM Devansh Sonigra > <deva...@tifrbng.res.in> wrote: > > > > Hello everyone, > > > > I am working with deal.II and encountering an issue with > KellyErrorEstimator. I started with a code that solves the Laplace equation > in a non-convex domain using adaptive refinement. The original code was > designed for rectangular mesh, but I updated it to work with a triangular > mesh. While KellyErrorEstimator runs fine on the initial mesh and first > iteration of refinement, the code breaks on the second iteration, and the > error seems to originate from KellyErrorEstimator. > > > > Has anyone faced a similar issue or have any suggestions on what might > be causing this? Any insights on debugging would be greatly appreciated. > > > > Thanks in advance! > > > > Best, > > Devansh > > > > Output: > > > > Number of active cells: 2400 > > > > Total number of cells: 2400 > > > > Number of degrees of freedom: 1281 > > > > 35 CG iterations needed to obtain convergence. > > > > Computing error norms > > > > Refining grid based on error estimator > > > > Debug message - Before running KellyErrorEstimator > > > > Debug message - After running KellyErrorEstimator > > > > Number of active cells: 2634 > > > > Total number of cells: 2712 > > > > Number of degrees of freedom: 1414 > > > > 53 CG iterations needed to obtain convergence. > > > > Computing error norms > > > > Refining grid based on error estimator > > > > Debug message - Before running KellyErrorEstimator > > > > Segmentation fault: 11 > > > > > > -- > > 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+un...@googlegroups.com. > > To view this discussion visit > https://groups.google.com/d/msgid/dealii/b4386a30-773c-4bc3-bd5d-68fb843b11cen%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 dealii+unsubscr...@googlegroups.com. To view this discussion visit https://groups.google.com/d/msgid/dealii/5121b014-5a78-4151-81cc-11662fe1632cn%40googlegroups.com.
* thread #1, queue = 'com.apple.main-thread', stop reason = signal SIGKILL * frame #0: 0x00007fff6fead929 libsystem_platform.dylib`_platform_memmove$VARIANT$Haswell + 41 frame #1: 0x00007fff6cfba353 libc++.1.dylib`std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >::__grow_by_and_replace(unsigned long, unsigned long, unsigned long, unsigned long, unsigned long, unsigned long, char const*) + 163 frame #2: 0x00007fff6cfba261 libc++.1.dylib`std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >::append(char const*, unsigned long) + 93 frame #3: 0x0000000105f4dc20 libdeal_II.g.9.6.0.dylib`std::__1::basic_string<char, std::__1::char_traits<char>, std::__1::allocator<char> >::operator+=(char const*) + 16 frame #4: 0x0000000105f4d9b3 libdeal_II.g.9.6.0.dylib`dealii::ExceptionBase::print_stack_trace(std::__1::basic_ostream<char, std::__1::char_traits<char> >&) const + 483 frame #5: 0x0000000105f4d4ad libdeal_II.g.9.6.0.dylib`dealii::ExceptionBase::generate_message() const + 349 frame #6: 0x0000000105f4d315 libdeal_II.g.9.6.0.dylib`dealii::ExceptionBase::what() const + 37 frame #7: 0x0000000105f4de3a libdeal_II.g.9.6.0.dylib`dealii::deal_II_exceptions::internals::abort(dealii::ExceptionBase const&) + 10 frame #8: 0x0000000100516d39 libdeal_II.g.9.6.0.dylib`void dealii::deal_II_exceptions::internals::issue_error_noreturn<dealii::StandardExceptions::ExcNotImplemented>(dealii::deal_II_exceptions::internals::ExceptionHandling, char const*, int, char const*, char const*, char const*, dealii::StandardExceptions::ExcNotImplemented) + 201 frame #9: 0x00000001061f76fe libdeal_II.g.9.6.0.dylib`dealii::QProjector<2>::DataSetDescriptor::subface(dealii::ReferenceCell const&, unsigned int, unsigned int, unsigned char, unsigned int, dealii::internal::SubfaceCase<2>) + 174 frame #10: 0x000000010206531c libdeal_II.g.9.6.0.dylib`dealii::MappingFE<2, 2>::fill_fe_subface_values(dealii::TriaIterator<dealii::CellAccessor<2, 2> > const&, unsigned int, unsigned int, dealii::Quadrature<1> const&, dealii::Mapping<2, 2>::InternalDataBase const&, dealii::internal::FEValuesImplementation::MappingRelatedData<2, 2>&) const + 348 frame #11: 0x00000001022e30cf libdeal_II.g.9.6.0.dylib`dealii::FESubfaceValues<2, 2>::do_reinit(unsigned int, unsigned int) + 287 frame #12: 0x00000001022eda0e libdeal_II.g.9.6.0.dylib`void dealii::FESubfaceValues<2, 2>::reinit<false>(dealii::TriaIterator<dealii::DoFCellAccessor<2, 2, false> > const&, unsigned int, unsigned int) + 382 frame #13: 0x0000000106ee4ab0 libdeal_II.g.9.6.0.dylib`void dealii::hp::FESubfaceValues<2, 2>::reinit<false>(dealii::TriaIterator<dealii::DoFCellAccessor<2, 2, false> > const&, unsigned int, unsigned int, unsigned int, unsigned int, unsigned int) + 256 frame #14: 0x0000000100c6b720 libdeal_II.g.9.6.0.dylib`void dealii::internal::integrate_over_irregular_face<double, 2, 2>(dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::DoFHandler<2, 2>::face_iterator, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::DoFHandler<2, 2>::face_iterator>, std::__1::allocator<std::__1::pair<dealii::DoFHandler<2, 2>::face_iterator const, std::__1::vector<double, std::__1::allocator<double> > > > >&, dealii::DoFHandler<2, 2>::active_cell_iterator const&, unsigned int, dealii::hp::FEFaceValues<2, 2>&, dealii::hp::FESubfaceValues<2, 2>&, dealii::KellyErrorEstimator<2, 2>::Strategy) + 496 frame #15: 0x0000000100c6a9da libdeal_II.g.9.6.0.dylib`void dealii::internal::estimate_one_cell<double, 2, 2>(dealii::DoFHandler<2, 2>::active_cell_iterator const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::DoFHandler<2, 2>::face_iterator, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::DoFHandler<2, 2>::face_iterator>, std::__1::allocator<std::__1::pair<dealii::DoFHandler<2, 2>::face_iterator const, std::__1::vector<double, std::__1::allocator<double> > > > >&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::KellyErrorEstimator<2, 2>::Strategy) + 266 frame #16: 0x0000000100c6a8cc libdeal_II.g.9.6.0.dylib`void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)::operator()(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&) const + 28 frame #17: 0x0000000100c6a88e libdeal_II.g.9.6.0.dylib`decltype(std::__1::forward<double>(fp)(std::__1::forward<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)&>(fp0)...)) std::__1::__invoke<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)&, dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&>(double&&, void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)&...) + 78 frame #18: 0x0000000100c6a81e libdeal_II.g.9.6.0.dylib`void std::__1::__invoke_void_return_wrapper<void>::__call<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)&, dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&>(double&&...) + 78 frame #19: 0x0000000100c6a79e libdeal_II.g.9.6.0.dylib`std::__1::__function::__alloc_func<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), std::__1::allocator<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>, void (dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>::operator()(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&) + 78 frame #20: 0x0000000100c69de7 libdeal_II.g.9.6.0.dylib`std::__1::__function::__func<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), std::__1::allocator<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>, void (dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>::operator()(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&) + 71 frame #21: 0x0000000100c714fb libdeal_II.g.9.6.0.dylib`std::__1::__function::__value_func<void (dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>::operator()(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&) const + 75 frame #22: 0x0000000100c6f103 libdeal_II.g.9.6.0.dylib`std::__1::function<void (dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&)>::operator()(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&) const + 67 frame #23: 0x0000000100c6edc8 libdeal_II.g.9.6.0.dylib`void dealii::WorkStream::internal::tbb_no_coloring::run<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::std_cxx20::type_identity<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > >::type const&, double, void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::internal::ParallelData<2, 2, double> const&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&, unsigned int, unsigned int)::'lambda'(dealii::WorkStream::internal::tbb_no_coloring::IteratorRangeToItemStream<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >::ItemType*)::operator()(dealii::WorkStream::internal::tbb_no_coloring::IteratorRangeToItemStream<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >::ItemType*) const + 328 frame #24: 0x0000000100c6ec50 libdeal_II.g.9.6.0.dylib`tbb::interface6::internal::concrete_filter<dealii::WorkStream::internal::tbb_no_coloring::IteratorRangeToItemStream<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >::ItemType*, dealii::WorkStream::internal::tbb_no_coloring::IteratorRangeToItemStream<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >::ItemType*, void dealii::WorkStream::internal::tbb_no_coloring::run<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::std_cxx20::type_identity<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > >::type const&, double, void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::internal::ParallelData<2, 2, double> const&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&, unsigned int, unsigned int)::'lambda'(dealii::WorkStream::internal::tbb_no_coloring::IteratorRangeToItemStream<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >::ItemType*)>::operator()(void*) + 48 frame #25: 0x0000000100482a03 libdeal_II.g.9.6.0.dylib`tbb::internal::stage_task::execute() + 115 frame #26: 0x000000010048a756 libdeal_II.g.9.6.0.dylib`tbb::internal::custom_scheduler<tbb::internal::IntelSchedulerTraits>::local_wait_for_all(tbb::task&, tbb::task*) + 598 frame #27: 0x000000010048790c libdeal_II.g.9.6.0.dylib`tbb::internal::generic_scheduler::local_spawn_root_and_wait(tbb::task*, tbb::task*&) + 268 frame #28: 0x00000001004879d0 libdeal_II.g.9.6.0.dylib`tbb::internal::generic_scheduler::spawn_root_and_wait(tbb::task&, tbb::task*&) + 32 frame #29: 0x00000001004836af libdeal_II.g.9.6.0.dylib`tbb::task::spawn_root_and_wait(tbb::task&) + 47 frame #30: 0x00000001004835e1 libdeal_II.g.9.6.0.dylib`tbb::pipeline::run(unsigned long, tbb::task_group_context&) + 145 frame #31: 0x00000001005409a5 libdeal_II.g.9.6.0.dylib`tbb::interface6::parallel_pipeline(unsigned long, tbb::interface6::filter_t<void, void> const&, tbb::task_group_context&) + 53 frame #32: 0x000000010053ad1b libdeal_II.g.9.6.0.dylib`tbb::interface6::parallel_pipeline(unsigned long, tbb::interface6::filter_t<void, void> const&) + 75 frame #33: 0x0000000100c65dc2 libdeal_II.g.9.6.0.dylib`void dealii::WorkStream::internal::tbb_no_coloring::run<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::std_cxx20::type_identity<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > >::type const&, double, void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::internal::ParallelData<2, 2, double> const&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&, unsigned int, unsigned int) + 306 frame #34: 0x0000000100c55e7e libdeal_II.g.9.6.0.dylib`void dealii::WorkStream::run<void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::internal::ParallelData<2, 2, double>&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > >&), void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> >, dealii::internal::ParallelData<2, 2, double>, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > >(dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > const&, dealii::std_cxx20::type_identity<dealii::TriaActiveIterator<dealii::DoFCellAccessor<2, 2, false> > >::type const&, double, void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy)::'lambda'(std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&), dealii::internal::ParallelData<2, 2, double> const&, std::__1::map<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> >, std::__1::vector<double, std::__1::allocator<double> >, std::__1::less<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > >, std::__1::allocator<std::__1::pair<dealii::TriaIterator<dealii::DoFAccessor<1, 2, 2, false> > const, std::__1::vector<double, std::__1::allocator<double> > > > > const&, unsigned int, unsigned int) + 190 frame #35: 0x0000000100c551b5 libdeal_II.g.9.6.0.dylib`void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::hp::QCollection<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy) + 853 frame #36: 0x0000000100c54d94 libdeal_II.g.9.6.0.dylib`void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::Quadrature<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ArrayView<dealii::ReadVector<double> const*, dealii::MemorySpace::Host> const&, dealii::ArrayView<dealii::Vector<float>*, dealii::MemorySpace::Host>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy) + 116 frame #37: 0x0000000100c54ccf libdeal_II.g.9.6.0.dylib`void dealii::KellyErrorEstimator<2, 2>::estimate<double>(dealii::Mapping<2, 2> const&, dealii::DoFHandler<2, 2> const&, dealii::Quadrature<1> const&, std::__1::map<unsigned int, dealii::Function<2, double> const*, std::__1::less<unsigned int>, std::__1::allocator<std::__1::pair<unsigned int const, dealii::Function<2, double> const*> > > const&, dealii::ReadVector<double> const&, dealii::Vector<float>&, dealii::ComponentMask const&, dealii::Function<2, double> const*, unsigned int, unsigned int, unsigned int, dealii::KellyErrorEstimator<2, 2>::Strategy) + 175 frame #38: 0x0000000100013f39 demo`LaplaceProblem<2>::refine_grid(this=<unavailable>) at demo.cc:280:4 frame #39: 0x0000000100003fb8 demo`LaplaceProblem<2>::run(this=0x00007ffeefbfe310, ncell=size=8, ndofs=size=8, L2_error=size=8, H1_error=size=8) at demo.cc:319:10 frame #40: 0x000000010000390a demo`main at demo.cc:342:12 frame #41: 0x00007fff6fcb7cc9 libdyld.dylib`start + 1