Dear community,
I’m interested in plasticity and I have been running some tests. Among them the benchmark Spiegelman et al. 2016 which is located in benchmarks/newton_solver_benchmark_set/spiegelman_et_al_2016/. When I just run with the input.prm-file it seems that aspect tries to carry out multiple time steps. However advection of the compositional field fails at the 2nd time step. Shouldn’t the problem be instantaneous anyways? (no End time is provided in the prm file)
Thank you for your help!
*** Timestep 1: t=5.86784e+12 seconds, dt=5.86784e+12 seconds
Skipping temperature solve because RHS is zero.
Solving C_1 system ... retrying linear solve with different preconditioner...
---------------------------------------------------------
TimerOutput objects finalize timed values printed to the
screen by communicating over MPI in their destructors.
Since an exception is currently uncaught, this
synchronization (and subsequent output) will be skipped
to avoid a possible deadlock.
---------------------------------------------------------
----------------------------------------------------
Exception 'ExcMessage (exception_message.str())' on rank 0 on processing:
--------------------------------------------------------
An error occurred in line <2676> of file </cluster/home/saarom/aspect/source/utilities.cc> in function
void aspect::Utilities::linear_solver_failed(const string&, const string&, const std::vector<dealii::SolverControl>&, const std::exception&, ompi_communicator_t* const&, const string&)
The violated condition was:
false
Additional information:
The iterative advection solver in Simulator::solve_advection did not
converge.
The initial residual was: 5.654641e+02
The final residual is: 5.594030e+02
The required residual for convergence is: 1.548819e-05
See output-spiegelman_et_at_2016/solver_history.txt for the full
convergence history.
The solver reported the following error:
--------------------------------------------------------
An error occurred in line <1075> of file
</cluster/home/saarom/dealii/install/include/deal.II/lac/solver_gmres.h>
in function
void dealii::SolverGMRES<VectorType>::solve(const MatrixType&,
VectorType&, const VectorType&, const PreconditionerType&) [with
MatrixType = dealii::TrilinosWrappers::SparseMatrix;
PreconditionerType = dealii::TrilinosWrappers::PreconditionILU;
VectorType = dealii::TrilinosWrappers::MPI::Vector]
The violated condition was:
iteration_state == SolverControl::success
Additional information:
Iterative method reported convergence failure in step 1000. The
residual in the last step was 559.403.
This error message can indicate that you have simply not allowed a
sufficiently large number of iterations for your iterative solver to
converge. This often happens when you increase the size of your
problem. In such cases, the last residual will likely still be very
small, and you can make the error go away by increasing the allowed
number of iterations when setting up the SolverControl object that
determines the maximal number of iterations you allow.
The other situation where this error may occur is when your matrix is
not invertible (e.g., your matrix has a null-space), or if you try to
apply the wrong solver to a matrix (e.g., using CG for a matrix that
is not symmetric or not positive definite). In these cases, the
residual in the last iteration is likely going to be large.
--------------------------------------------------------
--------------------------------------------------------
Aborting!
----------------------------------------------------
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[62154,1],0]
Exit code: 1
--------------------------------------------------------------------------