Hello Luca,
I have compiled the new version 6.2.1 for deal.II and I use its
FEFieldFunction class in the following way:
Functions::FEFieldFunction<DIMENSION> fe_function1
(dof_handler_total1,vector_j1);
VectorTools::interpolate(dof_handler_total2,fe_function1,vector_j2);
but I continue obtaining the following error message:
--------------------------------------------------------
An error occurred in line <765> of file <source/grid/grid_tools.cc> in
function
static std::pair<typename Container<dim, spacedim>::active_cell_iterator,
dealii::Point<spacedim> >
dealii::GridTools::find_active_cell_around_point(const
dealii::Mapping<dim, spacedim>&, const Container<dim, spacedim>&, const
dealii::Point<spacedim>&) [with int dim = 3, Container =
dealii::DoFHandler, int
spacedim = 3]
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 also tested:
Functions::FEFieldFunction<DIMENSION> fe_function1
(dof_handler_total1,vector_j1);
ConstraintMatrix hanging_node_constraints;
hanging_node_constraints.clear ();
DoFTools::make_hanging_node_constraints
(dof_handler_total2,hanging_node_constraints);
QGauss3<DIMENSION> quadrature_formula;
VectorTools::project(dof_handler_total2,
hanging_node_constraints,
quadrature_formula,
fe_function1,
vector_j2);
and I get the following error:
--------------------------------------------------------
An error occurred in line <765> of file <source/grid/grid_tools.cc> in
function
static std::pair<typename Container<dim,
spacedim>::active_cell_iterator, dealii::Point<spacedim> >
dealii::GridTools::find_active_cell_around_point(const
dealii::Mapping<dim, spacedim>&, const Container<dim, spacedim>&,
const dealii::Point<spacedim>&) [with int dim = 3, Container =
dealii::DoFHandler, int spacedim = 3]
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.00206228 0.0925332 0.00129324> could not be found inside any
of the subcells of a coarse grid cell.
--------------------------------------------------------
I really don't know where the error can be.
Thanks in advance!
Isa
>>
>> 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?
>
> Yes to both questions. :)
>
> To my knowledge, upgrading from 6.1.0 to 6.2.1 (and even to the subversion
> tree!), should be mostly painless.
>
> I'd suggest you to download the svn version of deal on a different
> location, build that one, and see how your code reacts to it. If you don't
> feel like subversion is your thing, then upgrade to 6.2.1...
>
> Luca.
>
>
>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
