cell_iterator = GridTools::find_active_cell_around_point(q, dof_handler, Point<3>(x, y, z) );
works fine. The segementation fault happens at this point cell_iterator->set_material_id(value); I don't understand why this should happen. Pietro DoFHandler<3> dof_handler(triangulation); MappingQ1<3> q; std::pair<DoFHandler<3>::active_cell_iterator, Point<3> > cell_iterator; cell_iterator = GridTools::find_active_cell_around_point(q, dof_handler, Point<3>(x, y, z) ); cell_iterator->set_material_id(value); --- On Fri, 12/10/10, Pietro Maximoff <[email protected]> wrote: From: Pietro Maximoff <[email protected]> Subject: Re: [deal.II] Finding a cell a point belongs to To: "Luca Heltai" <[email protected]> Cc: [email protected] Date: Friday, December 10, 2010, 9:17 AM Is this correct usage because I get segmentation fault for this test. DoFHandler<3> dof_handler(triangulation); MappingQ1<3> q; std::pair<DoFHandler<3>::active_cell_iterator, Point<3> > cell_iterator; cell_iterator = GridTools::find_active_cell_around_point(q, dof_handler, Point<3>(x, y, z) ); Pietro --- On Fri, 12/10/10, Pietro Maximoff <[email protected]> wrote: From: Pietro Maximoff <[email protected]> Subject: Re: [deal.II] Finding a cell a point belongs to To: "Luca Heltai" <[email protected]> Cc: [email protected] Date: Friday, December 10, 2010, 9:01 AM Thanks Luca. Just what I was looking for. Pietro --- On Fri, 12/10/10, Luca Heltai <[email protected]> wrote: From: Luca Heltai <[email protected]> Subject: Re: [deal.II] Finding a cell a point belongs to To: "Pietro Maximoff" <[email protected]> Cc: [email protected] Date: Friday, December 10, 2010, 7:23 AM The function GridTools::find_active_cell_around_point should do what you ask for. http://www.dealii.org/developer/doxygen/deal.II/classGridTools.html#aa4b56da0763f08676ba92f80cab491da Luca. -- Luca Heltai <[email protected]> http://people.sissa.it/~heltai/ Scuola Internazionale Superiore di Studi Avanzati Phone: +39 040 3787 449, Office: 732 -- There are no answers, only cross references On 10/dic/2010, at 16.17, Pietro Maximoff wrote: > Hi > > Is there a way to find the cell a point belongs to without having to iterate over all active cells? > > Pietro > > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii -----Inline Attachment Follows----- _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii -----Inline Attachment Follows----- _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
