Hi Javier, yes, you are not the only one with this problem. We are working on improving GridTools::find_active_cell_around_point, that should fix your problem.
Best, Timo On Wed, May 23, 2012 at 6:20 PM, Javier Muñoz <[email protected]> wrote: > Dear all > > > I'm not sure if this question has been answered before. I'm using > deal.ii 7.1.0 and I'm working with an adaptively refined mesh in 3d. For > the first step, before refinement, I extract a point value via > > VectorTools::point_value(const DoFHandler< dim, spacedim > & dof,const > InVector &fe_function,const Point<spacedim>& point) > > and there is no problem. But, after refinement, when I tried to access > again to this point, the program prompt me with an error: > > The point <-0.0715000 1.13240 0.00000> could not be found inside any of > the subcells of a coarse grid cell. > > Stacktrace: > ----------- > #0 /home/zerpiko/Documents/Simulation/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 /home/zerpiko/Documents/Simulation/deal.II/lib/libdeal_II.g.so.7.1.0: > double dealii::VectorTools::point_value<3, dealii::Vector<double>, > 3>(dealii::Mapping<3, 3> const&, dealii::DoFHandler<3, 3> const&, > dealii::Vector<double> const&, dealii::Point<3, double> const&) > #2 /home/zerpiko/Documents/Simulation/deal.II/lib/libdeal_II.g.so.7.1.0: > double dealii::VectorTools::point_value<3, dealii::Vector<double>, > 3>(dealii::DoFHandler<3, 3> const&, dealii::Vector<double> const&, > dealii::Point<3, double> const&) > #3 ./heat_pipe_with_advection: Heat_Pipe::assemble_system_3d() > #4 ./heat_pipe_with_advection: Heat_Pipe::run() > #5 ./heat_pipe_with_advection: main > > > This happens with some points and no with others. I had read in previous > emails that this could be related with the new hanging nodes generated > by the refinement. Am I doing the wrong thing or is there another way to > access the point values? > > Thanks in advance > Regards > Javier Munoz > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii -- Timo Heister http://www.math.tamu.edu/~heister/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
