Isabel, can you reproduce the same problem in a 20 lines code? If yes, you could send that code to me, and I'll try to see if I can spot the problem.
Luca. -- Luca Heltai <[email protected]> http://people.sissa.it/~heltai/ Scuola Internazionale Superiore di Studi Avanzati Phone: +39 040 3787 449, Office: 255 -- There are no answers, only cross references. On 29/mar/2010, at 21.24, [email protected] wrote: > Hello Luca, > > I have compiled the new version 6.2.1 for deal.II and I use its > FEFieldFunction class in the following way: > > Functions::FEFieldFunction<DIMENSION> fe_function1 > (dof_handler_total1,vector_j1); > VectorTools::interpolate(dof_handler_total2,fe_function1,vector_j2); > > but I continue obtaining the following error message: > -------------------------------------------------------- > An error occurred in line <765> of file <source/grid/grid_tools.cc> in > function > static 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] > 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 <0.0716229 0.0728616 0.00130000> could not be found inside any > of the > subcells of a coarse grid cell. > -------------------------------------------------------- > > I have also tested: > > Functions::FEFieldFunction<DIMENSION> fe_function1 > (dof_handler_total1,vector_j1); > > ConstraintMatrix hanging_node_constraints; > hanging_node_constraints.clear (); > > DoFTools::make_hanging_node_constraints > (dof_handler_total2,hanging_node_constraints); > > QGauss3<DIMENSION> quadrature_formula; > > VectorTools::project(dof_handler_total2, > hanging_node_constraints, > quadrature_formula, > fe_function1, > vector_j2); > > and I get the following error: > -------------------------------------------------------- > An error occurred in line <765> of file <source/grid/grid_tools.cc> in > function > static 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] > 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 <0.00206228 0.0925332 0.00129324> could not be found inside any > of the subcells of a coarse grid cell. > -------------------------------------------------------- > > I really don't know where the error can be. > > Thanks in advance! > Isa > >>> >>> Hello Luca, >>> >>> thank-you for answer! >>> I use the 6.1.0 version of DealII since I have a lot of code written and >>> I didn't want to change some methods in order to use the new versions. >>> >>> Is it in that version where the bug was? Is the correct version the >>> 6.2.1 one? >> >> Yes to both questions. :) >> >> To my knowledge, upgrading from 6.1.0 to 6.2.1 (and even to the subversion >> tree!), should be mostly painless. >> >> I'd suggest you to download the svn version of deal on a different >> location, build that one, and see how your code reacts to it. If you don't >> feel like subversion is your thing, then upgrade to 6.2.1... >> >> Luca. >> >> >> > > > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
