[aspect-devel] Trilinos problems ( I guess )

[Jonathan Perry-Houts]

I just updated / recompiled deal.ii and aspect on my cluster and now
aspect is failing every time with a convergence error.

I attached all the outputs from running convection-box.prm to this email
(hopefully that works with the listserv). The important part is this
stuff from stderr:

----------- 8< ---- CUT HERE ------ 8< -----
The iterative Stokes solver did not converge. It reported the following
error:
....
The iterative (top left) solver in BlockSchurPreconditioner::vmult did
not converge. It reported the following error:
....
AztecOO::Iterate error code -2: numerical breakdown
----------- 8< ---- CUT HERE ------ 8< -----

This happens for all models I've tried.

I've never run into this on my desktop, or on the cluster before. I'm
not sure what could have changed to affect the Stokes solver, but I'm
not doing anything different in my own workflow so it must have been
some system library that Trilinos uses, I guess... Any ideas which one
it could be?

For reference, I'm using Trilinos 11.12.1, deal.ii development source,
aspect development source.

Thanks in advance!

Cheers,
Jonathan

-- 
Jonathan Perry-Houts
Ph.D. Candidate
Department of Geological Sciences
1272 University of Oregon
Eugene, OR 97403-1272



-------------- next part --------------
-----------------------------------------------------------------------------
-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
--     . version 1.4.pre
--     . running in OPTIMIZED mode
--     . running with 1 MPI process
--     . using Trilinos
-----------------------------------------------------------------------------

Number of active cells: 256 (on 5 levels)
Number of degrees of freedom: 3,556 (2,178+289+1,089)

*** Timestep 0:  t=0 seconds
   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 
-------------- next part --------------
0 2.14146e+13
1 2.03663e+13
2 1.69189e+13
3 1.2777e+13
4 9.24983e+12
5 5.98204e+12
6 4.48472e+12
7 3.68882e+12
8 3.34609e+12
9 3.15623e+12
10 3.04184e+12
11 3.00375e+12
12 2.99781e+12
13 2.99645e+12
14 2.99533e+12
15 2.99492e+12
16 2.99484e+12
17 2.99481e+12
18 2.9948e+12
19 2.9948e+12
20 2.99479e+12
21 2.99479e+12
22 2.99479e+12
23 2.99479e+12
24 2.99479e+12
25 2.99479e+12
26 2.99479e+12
27 2.99479e+12
28 2.99479e+12
29 2.99479e+12
30 2.99479e+12
31 2.99479e+12

33 2.99479e+12
-------------- next part --------------


----------------------------------------------------
Exception on processing: 

--------------------------------------------------------
An error occurred in line <744> of file </usr/local/packages/aspect/dev/source/simulator/solver.cc> in function
    double aspect::Simulator<dim>::solve_stokes() [with int dim = 2]
The violated condition was: 
    false
The name and call sequence of the exception was:
    ExcMessage (std::string("The iterative Stokes solver " "did not converge. It reported the following error:\n\n") + exc.what() + "\n See " + parameters.output_directory+"solver_history.txt" + " for convergence history.")
Additional Information: 
The iterative Stokes solver did not converge. It reported the following error:


--------------------------------------------------------
An error occurred in line <347> of file </usr/local/packages/aspect/dev/source/simulator/solver.cc> in function
    void aspect::internal::BlockSchurPreconditioner<PreconditionerA, PreconditionerMp>::vmult(aspect::LinearAlgebra::BlockVector&, const BlockVector&) const [with PreconditionerA = dealii::TrilinosWrappers::PreconditionAMG; PreconditionerMp = dealii::TrilinosWrappers::PreconditionILU; aspect::LinearAlgebra::BlockVector = dealii::TrilinosWrappers::MPI::BlockVector]
The violated condition was: 
    false
The name and call sequence of the exception was:
    ExcMessage (std::string("The iterative (top left) solver in BlockSchurPreconditioner::vmult " "did not converge. It reported the following error:\n\n") + exc.what())
Additional Information: 
The iterative (top left) solver in BlockSchurPreconditioner::vmult did not converge. It reported the following error:


--------------------------------------------------------
An error occurred in line <279> of file </home13/jperryh2/packages/build/deal.II/source/lac/trilinos_solver.cc> in function
    void dealii::TrilinosWrappers::SolverBase::do_solve(const dealii::TrilinosWrappers::PreconditionBase&)
The violated condition was: 
    false
The name and call sequence of the exception was:
    ExcMessage("AztecOO::Iterate error code -2: " "numerical breakdown")
Additional Information: 
AztecOO::Iterate error code -2: numerical breakdown
--------------------------------------------------------

--------------------------------------------------------

 See output/solver_history.txt for convergence history.
--------------------------------------------------------

Aborting!
----------------------------------------------------
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 538 bytes
Desc: OpenPGP digital signature
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20150909/97ead525/attachment.sig>