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
