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
