Then I get
An error occurred in line <754> of file <~/Downloads/dealii/include/deal.II/
lac/constraint_matrix.templates.h> in function
void dealii::internal::{anonymous}::import_vector_with_ghost_elements(
const dealii::TrilinosWrappers::MPI::Vector&, const dealii::IndexSet&, const
dealii::IndexSet&, dealii::TrilinosWrappers::MPI::Vector&, dealii::internal
::bool2type<false>)
The violated condition was:
!vec.has_ghost_elements()
Additional information:
You are trying an operation on a vector that is only allowed if the
vector has no ghost elements, but the vector you are operating on does have
ghost elements. Specifically, vectors with ghost elements are read-only and
cannot appear in operations that write into the
se vectors.
See the glossary entry on 'Ghosted vectors' for more information.
Stacktrace:
-----------
#0 /opt/dealII/lib/libdeal_II.g.so.9.0.0-pre:
#1 /opt/dealII/lib/libdeal_II.g.so.9.0.0-pre: void
dealii::ConstraintMatrix::distribute<dealii::TrilinosWrappers::MPI::Vector>(dealii::TrilinosWrappers::MPI::Vector&)
const
#2 main: Step15::MinimalSurfaceProblem<2>::run()
#3 main: main
--------------------------------------------------------
Calling MPI_Abort now.
To break execution in a GDB session, execute 'break MPI_Abort' before
running. You can also put the following into your ~/.gdbinit:
set breakpoint pending on
break MPI_Abort
set breakpoint pending auto
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 255.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
Calling MPI_Abort now.
To break execution in a GDB session, execute 'break MPI_Abort' before
running. You can also put the following into your ~/.gdbinit:
set breakpoint pending on
break MPI_Abort
set breakpoint pending auto
[linux-lb8c:17949] 1 more process has sent help message help-mpi-api.txt /
mpi-abort
[linux-lb8c:17949] Set MCA parameter "orte_base_help_aggregate" to 0 to see
all help / error messages
at the same line (.distribute())
Am Sonntag, 10. September 2017 23:25:07 UTC+2 schrieb Timo Heister:
>
> Maxi,
>
> can you try replacing
> > present_solution.reinit(solution_relevant_partitioning, MPI_COMM_WORLD);
> by
> > present_solution.reinit(solution_partitioning,
> solution_relevant_partitioning, MPI_COMM_WORLD);
> ?
>
> On Sun, Sep 10, 2017 at 3:46 PM, 'Maxi Miller' via deal.II User Group
> <dea...@googlegroups.com <javascript:>> wrote:
> > 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
> https://urldefense.proofpoint.com/v2/url?u=http-3A__www.dealii.org_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=AhRdZ_UVqz83aYxbDsBIq5YyD9M9NAy6-hnCaGUmA20&s=D1aNzlOq6SC-rbTJzL5NI7WYom6R609iP-sWpsqUzas&e=
>
> > For mailing list/forum options, see
> >
> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=AhRdZ_UVqz83aYxbDsBIq5YyD9M9NAy6-hnCaGUmA20&s=sMVWJCGxV_WAg6GGoIWl5y3R4pyGQzxivam0NjzmAe8&e=
>
> > ---
> > 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 dealii+un...@googlegroups.com <javascript:>.
> > For more options, visit
> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=AhRdZ_UVqz83aYxbDsBIq5YyD9M9NAy6-hnCaGUmA20&s=c0d5V5EPM8_gh4Lh5_uhq962B54bq2TPyTqk1Vf-GgM&e=
>
> .
>
>
>
> --
> Timo Heister
> http://www.math.clemson.edu/~heister/
>
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
