I am having that error message when I run in parallel with 2 processor like
this
mpirun -n 2 ./darcyparl
and when I run with one MPI process with mpirun -n 1 ./darcyparl I have
this error message
An error occurred in line <474> of file
</home/franckm/darcyparl/darcyparl.cc> in function
void Step43::LinearSolvers::InverseMatrix<Matrix,
Preconditioner>::vmult(VectorType&, const VectorType&) const [with
VectorType = dealii::TrilinosWrappers::MPI::Vector; Matrix =
dealii::TrilinosWrappers::SparseMatrix; Preconditioner =
dealii::TrilinosWrappers::PreconditionIC]
The violated condition was:
false
Additional information:
--------------------------------------------------------
An error occurred in line <590> of file
</home/franckm/deal.ii-candi/deal.II-v8.5.0/include/deal.II/lac/solver_cg.h>
in function
void dealii::SolverCG<VectorType>::solve(const MatrixType&,
VectorType&, const VectorType&, const PreconditionerType&) [with MatrixType
= dealii::TrilinosWrappers::SparseMatrix; PreconditionerType =
dealii::TrilinosWrappers::PreconditionIC; VectorType =
dealii::TrilinosWrappers::MPI::Vector]
The violated condition was:
false
Additional information:
Iterative method reported convergence failure in step 81. The residual in
the last step was 49.0941.
This error message can indicate that you have simply not allowed a
sufficiently large number of iterations for your iterative solver to
converge.
how large should be the number of iterations for it to converge?
linear solver being the same as in step-43
Best
F75
On Saturday, 15 July 2017 23:12:17 UTC+2, Daniel Arndt wrote:
>
> franck75,
>
> Some hints for debugging this:
> Do you observe this error only if you run in parallel or also with only
> one MPI process?
> Use gdb to find the last line in your program that triggers the Assert.
> Which line is it?
> Try to find out where the numbers 81 and 45 come from. The best candidates
> would probably be
> the number of locally relevant or locally owned DoFs.
>
> This should give you a good hint on what is going wrong. Maybe you are
> just using a
> ghosted vector in a place where you would need a totally distributed one
> instead.
>
> Best,
> Daniel
>
--
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.
