[aspect-devel] Problems with Mapping::transform_real_to_unit_cell

Wolfgang Bangerth bangerth at math.tamu.edu
Sun Jul 1 07:14:01 PDT 2012

> I've attached the parameter file that causes the problem.  It still throws the exception below, but I think this is a slightly different cause.
> --------------------------------------------------------
> An error occurred in line <851> of file </home/emheien/trunk/deal.II/source/grid/grid_tools.cc> in function
>      std::pair<typename Container<dim, spacedim>::active_cell_iterator, dealii::Point<dim, double> > 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<spacedim>(p)
> Additional Information:
> The point <891669. 896973. 3.23958e+06> could not be found inside any of the subcells of a coarse grid cell.

This point has a distance from the origin equal to 3477716, but in your 
parameter file you describe the inner radius as 3481000. Do you know how 
that point got there?


Wolfgang Bangerth               email:            bangerth at math.tamu.edu
                                 www: http://www.math.tamu.edu/~bangerth/

More information about the Aspect-devel mailing list