Are you running the code in debug mode? What is the error you are getting? When does it blow up? What message is printed when it blows up?
We need a little more information to understand what is it that you are doing wrong... Most likely, the point you are looking for is not inside the domain (95% of the time, if you are in optimized mode, the code just seg faults in this case). If you run in debug mode, then an error should be printed that tells you what went wrong. If not, then the seg fault is likely to be outside of deal.II (i.e., access to non initialized data). Run it in debug mode first. Ah. I almost forgot. Run it in debug mode first. 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 18.17, Pietro Maximoff wrote: > 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 _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
