Hi Javier, Wolfgang and I managed to find and fix the bug you were most likely seeing. It is checked in in the development version of deal.II. Let me know if it does not fix your problem.
On Sat, May 26, 2012 at 5:22 PM, Timo Heister <[email protected]> wrote: > 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/ -- Timo Heister http://www.math.tamu.edu/~heister/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
