On 12/08/2017 09:14 AM, Marek Čapek wrote:
--------------------------------------------------------
An error occurred in line <1099> of file
</home/mcapek/candis/candi_8_5/tmp/unpack/deal.II-v8.5.0/include/deal.II/lac/trilinos_vector_base.h>
in function
dealii::IndexSet
dealii::TrilinosWrappers::VectorBase::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.
Stacktrace:
-----------
#0
/home/mcapek/candis/candi_8_5/deal.II-v8.5.0/lib/libdeal_II.g.so.8.5.0:
dealii::TrilinosWrappers::VectorBase::locally_owned_elements() const
#1
/home/mcapek/candis/candi_8_5/deal.II-v8.5.0/lib/libdeal_II.g.so.8.5.0:
dealii::LinearAlgebra::ReadWriteVector<double>::import(dealii::TrilinosWrappers::MPI::Vector
const&, dealii::VectorOperation::values,
std::shared_ptr<dealii::LinearAlgebra::CommunicationPatternBase const>)
#2
/home/mcapek/candis/candi_8_5/deal.II-v8.5.0/lib/libdeal_II.g.so.8.5.0:
dealii::LinearAlgebra::distributed::Vector<double>::operator=(dealii::TrilinosWrappers::MPI::Vector
const&)
#3 ./main: NSSystem<3>::assemble_and_solve_system(SolutionCrate&,
SolutionCrate&, double)
#4 ./main: NSSystem<3>::compute_solution(SolutionCrate&,
SolutionCrate&, double, double)
#5 ./main: NSSystem<3>::compute_solution_get_dt(SolutionCrate&,
SolutionCrate&, double, double)
#6 ./main: Main<3>::run()
#7 ./main: main
--------------------------------
Marek -- I don't actually think that your problem happens when computing
the norm. In fact, as the backtrace above shows, it happens in your
assemble_and_solve_system() function in a location where you are
assigning a Trilinos vector to a LinearAlgebra::distributed::Vector. The
error message to me would suggest that one of the two vectors does not
have the correct size after mesh refinement.
It is often useful to carefully read the error message and stack trace.
It is also often useful to run a program in a debugger, because you can
then inspect the state of the vectors at the place where the problem
happens, and infer which variable may have been wrongly initialized.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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.
