Dear All,
I am not able to directly localize the - what I would call - bug, but
here is the description on how far I got.
1. I test whether a point is inside my domain with a loop across all
cells and the "cell->point_inside (point)" command. Yes my point is
inside, also visual inspection and geometry: the point is inside.
2. I run GridTools::find_active_cell_around_point<3,
dealii::DoFHandler, 3>(dealii::Mapping<3, 3> const&,
dealii::DoFHandler<3, 3> const&, dealii::Point<3, double> const&)
via the FEFieldFunction and that raises the exception
------------------------------------------------------
The violated condition was:
best_cell.first.state() == IteratorState::valid
The name and call sequence of the exception was:
ExcPointNotFound<dim>(p)
Additional Information:
The point <304.767 -57.0113 254.766> could not be found inside any of
the subcells of a coarse grid cell.
-------------------------------------------------------
3. I rerun, and the exception is repeatable: every time I run it.
Running a debugger I found that the exception is thrown after (not
during) the return statement of
-------------------------------------------------------
template <typename T, typename P>
inline
T & SmartPointer<T,P>::operator * () const
{
Assert(t != 0, ExcNotInitialized());
return *t;
}
-------------------------------------------------------
Can someone help?
THANKS,
Steffen
--------------------------------------------------------
An error occurred in line <830> of file <...
/deal.II/source/grid/grid_tools.cc> in function
std::pair<typename Container<dim, spacedim>::active_cell_iterator,
dealii::Point<spacedim> >
dealii::GridTools::find_active_cell_around_point(const
dealii::Mapping<dim, spacedim>&, const Container<dim, spacedim>&, const
dealii::Point<spacedim>&) [with int dim = 3, Container =
dealii::DoFHandler, int spacedim = 3, typename Container<dim,
spacedim>::active_cell_iterator =
dealii::TriaActiveIterator<dealii::DoFCellAccessor<dealii::DoFHandler<3,
3> > >]
The violated condition was:
best_cell.first.state() == IteratorState::valid
The name and call sequence of the exception was:
ExcPointNotFound<dim>(p)
Additional Information:
The point <304.767 -57.0113 254.766> could not be found inside any of
the subcells of a coarse grid cell.
Stacktrace:
-----------
#0 deal.II/lib/libdeal_II.g.so.7.1.0: std::pair<dealii::DoFHandler<3,
3>::active_cell_iterator, dealii::Point<3, double> >
dealii::GridTools::find_active_cell_around_point<3, dealii::DoFHandler,
3>(dealii::Mapping<3, 3> const&, dealii::DoFHandler<3, 3> const&,
dealii::Point<3, double> const&)
#1 deal.II/lib/libdeal_II.g.so.7.1.0:
dealii::Functions::FEFieldFunction<3, dealii::DoFHandler<3, 3>,
dealii::Vector<double> >::vector_gradient(dealii::Point<3, double>
const&, std::vector<dealii::Tensor<1, 3, double>,
std::allocator<dealii::Tensor<1, 3, double> > >&) const
#2 exe: ares::FEM::getStress(dealii::Point<3, double>)
#3 exe: ares::FEM::getStress(Eigen::Matrix<float, 3, 1, 0, 3, 1>)
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
