In my program I have a vector "present_solution", defined as
LinearAlgebraTrilinos::MPI::Vector present_solution;
const size_t dof_numbers = dof_handler.n_dofs();
IndexSet solution_partitioning(dof_numbers), solution_relevant_partitioning(
dof_numbers);
solution_partitioning = dof_handler.locally_owned_dofs();
DoFTools::extract_locally_relevant_dofs(dof_handler,
solution_relevant_partitioning);
present_solution.reinit(solution_relevant_partitioning, MPI_COMM_WORLD);
and a ConstraintMatrix "boundary_constraints" defined as
ConstraintMatrix boundary_constraints;
boundary_constraints.clear();
boundary_constraints.reinit(solution_relevant_partitioning);
DoFTools::make_hanging_node_constraints(dof_handler, boundary_constraints);
VectorTools::interpolate_boundary_values(dof_handler, 0, BoundaryValues<dim>
(), boundary_constraints);
boundary_constraints.close();
Thus I assume that both have the same dimensions. But now when calling
boundary_constraints.distribute(present_solution);
I get the error (when running with more than one mpi thread):
An error occurred in line <1367> of file </home/roland/Downloads/dealii/
include/deal.II/lac/trilinos_vector.h> in function
dealii::IndexSet dealii::TrilinosWrappers::MPI::Vector::
locally_owned_elements() const
The violated condition was:
owned_elements.size()==size()
Additional information:
The locally owned elements have not been properly initialized! This
happens for example if this object has been initialized with exactly one
overlapping IndexSet.
In order to debug this, I would like to access the element
"owned_elements", but it is declared private. Is there another way to debug
my problem?
Thanks!
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
