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

Reply via email to