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.
For more options, visit https://groups.google.com/d/optout.
