Hi, I am using the gradient_list function to calculate the velocity gradient in x direction on the bottom wall of backward facing step and trying to output it using the table handler. But I am getting an exception of points not found in the coarse grid. The points vector and the gradient vector size is 288.
Exception on processing: -------------------------------------------------------- An error occurred in line <4343> of file </home/kambler/soft/dealii/dealii/source/grid/grid_tools.cc> in function std::tuple<std::vector<typename dealii::Triangulation<dim, spacedim>::active_cell_iterator, std::allocator<typename dealii::Triangulation<dim, spacedim>::active_cell_iterator> >, std::vector<std::vector<dealii::Point<dim, double>, std::allocator<dealii::Point<dim, double> > >, std::allocator<std::vector<dealii::Point<dim, double>, std::allocator<dealii::Point<dim, double> > > > >, std::vector<std::vector<unsigned int, std::allocator<unsigned int> >, std::allocator<std::vector<unsigned int, std::allocator<unsigned int> > > > > dealii::GridTools::compute_point_locations(const dealii::GridTools::Cache<dim, spacedim>&, const std::vector<dealii::Point<spacedim> >&, const typename dealii::Triangulation<dim, spacedim>::active_cell_iterator&) [with int dim = 2; int spacedim = 2; typename dealii::Triangulation<dim, spacedim>::active_cell_iterator = dealii::TriaActiveIterator<dealii::CellAccessor<2, 2> >; typename dealii::Triangulation<dim, spacedim>::active_cell_iterator = dealii::TriaActiveIterator<dealii::CellAccessor<2, 2> >] The violated condition was: std::get<3>(cqmp).size() == 0 Additional information: The point <0 0> could not be found inside any of the subcells of a coarse grid cell. -------------------------------------------------------- An error occurred in line <4343> of file </home/kambler/soft/dealii/dealii/source/grid/grid_tools.cc> in function std::tuple<std::vector<typename dealii::Triangulation<dim, spacedim>::active_cell_iterator, std::allocator<typename dealii::Triangulation<dim, spacedim>::active_cell_iterator> >, std::vector<std::vector<dealii::Point<dim, double>, std::allocator<dealii::Point<dim, double> > >, std::allocator<std::vector<dealii::Point<dim, double>, std::allocator<dealii::Point<dim, double> > > > >, std::vector<std::vector<unsigned int, std::allocator<unsigned int> >, std::allocator<std::vector<unsigned int, std::allocator<unsigned int> > > > > dealii::GridTools::compute_point_locations(const dealii::GridTools::Cache<dim, spacedim>&, const std::vector<dealii::Point<spacedim> >&, const typename dealii::Triangulation<dim, spacedim>::active_cell_iterator&) [with int dim = 2; int spacedim = 2; typename dealii::Triangulation<dim, spacedim>::active_cell_iterator = dealii::TriaActiveIterator<dealii::CellAccessor<2, 2> >; typename dealii::Triangulation<dim, spacedim>::active_cell_iterator = dealii::TriaActiveIterator<dealii::CellAccessor<2, 2> >] The violated condition was: std::get<3>(cqmp).size() == 0 Additional information: The point <17 0> could not be found inside any of the subcells of a coarse grid cell. PS: I am newer to using deal ii and I checked the discussion on the question Usage of FEFieldFunction.vector_value_list on a parallel::distributed::Triangulation posted earlier but still not able to find the mistake. -- 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/73791db1-68b0-4f49-832b-5896b7563876%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.