Rajeshwari, We need a little bit more information to be able to help you. - Which version of the library are you using? - What does the mesh you are using look like? - Do you see this problem also when running in serial or only in parallel?
It would also help to reduce your test case to a size as small as possible if you want to share it with us. Best, Daniel Am Mo., 15. Juli 2019 um 13:06 Uhr schrieb Rajeshwari Kamble < [email protected]>: > 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 [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/73791db1-68b0-4f49-832b-5896b7563876%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/73791db1-68b0-4f49-832b-5896b7563876%40googlegroups.com?utm_medium=email&utm_source=footer> > . > For more options, visit https://groups.google.com/d/optout. > -- 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/CAOYDWb%2BSUCbhb012wPZ04z_%3Dir5pTgi8bCjRAzKkCbT1h-ShJQ%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
