Isabel, I had a similar problem some time ago for a bug in a previous version of deal.II (which I fixed then). The bug was manifesting when the point was close to a hanging node. If this is the case, can I ask you what version of Deal you are using?
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 13.26, Isabel Gil wrote: > To the email placed below, I give more information in case anyone could help > me. > > ***************************************************************************** > Both meshes describe the same domain. > > Tringulation 1, dof_hander1 have: > > Total number of dofs: 39664. > 29748 dofs with a material value of 1. > 9916 dofs with a material value of 2. > > 9916 dofs with boundary_value of 4. > 9916 dofs with boundary_value of 3. > 1648 dofs with boundary_value of 1. > > On another hand, Triangulation 2, dof_hander2 have: > > Total number of dofs: 39660. > 29745 dofs with a material value of 1. > 9915 dofs with a material value of 2. > > 9915 dofs with boundary_value of 3. > 6644 dofs with boundary_value of 4. > 4855 dofs with boundary_value of 1. > > I would like to remind you that what I want to get is the following > interpolation: > > Functions::FEFieldFunction<DIMENSION> fe_function1 > (dof_handler_total1,vector1); > VectorTools::interpolate(dof_handler_total2,fe_function1,vector2); > > ***************************************************************************** > > Hi! > > I have run the debugger version of my program and the information about the > error I obtain is: > > -------------------------------------------------------- > An error occurred in line <759> of file <source/grid/grid_tools.cc> in > function > static std::pair<typename Container<dim>::active_cell_iterator, > dealii::Point<dim> > dealii::GridTools::find_active_cell_around_point(const > dealii::Mapping<dim>&, const Container<dim>&, const dealii::Point<dim>&) > [with int dim = 3, Container = dealii::DoFHandler] > 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 revised my meshes and both of them work with exactly the same domain > and, according to FEFieldFunction documentation, the only requirement was > "the points at which this function can be evaluated must be inside the domain > of the dof_handler". > > I have revised my meshes and the point to which the error message makes > reference doesn't exist but there are several points around it, so, I think > it could be interpolated, couldn't it? > > Any advice will be welcome! > > Best > Isa > > Wolfgang Bangerth escribió: >>> ****** >>> p0_29592: p4_error: interrupt SIGSEGV: 11 >>> ****** >>> >>> >> >> You are accessing an invalid memory location. You need to run the program >> in a debugger to see at which location in your program this happens. >> >> W. >> >> ------------------------------------------------------------------------- >> Wolfgang Bangerth email: >> [email protected] >> >> www: >> http://www.math.tamu.edu/~bangerth/ >> >> >> >> > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
