Hello Luca,
thank-you for answer!
I use the 6.1.0 version of DealII since I have a lot of code written and
I didn't want to change some methods in order to use the new versions.
Is it in that version where the bug was? Is the correct version the
6.2.1 one?
Thanks in advance.
My best
Isa
Luca Heltai escribió:
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
