Dear DuMuX users,
After Bernd's email, with my colleague Mustapha, we used the AMG solver
for our parallel computations.
As Mustapha told, it solved the problem but only partially since the
computations stopped without explanation as shown below
Newton iteration 2 done, maximum relative shift = 20.7159
Solve: M deltax^k =
The computation does not really crash but it seems to be frozen at this
step.
So to isolate the problem, we performed the same computation in a
sequential run with the AMG solver.
The computation was OK but 10 times longer than the same one performed
with the default solver (ILU0BiCGSTABBackend).
Have you ever met this issue? Is it normal?
For the AMG solver, are there some additional parameters to give?
I want to precise that we use a sequential approach to solve two-phase
flow with geochemical transport. We solve firstly a two phase flow with
the co2 model and then a transport reactive problem.
The problem occurs in the co2 flow step and not in the transport
reactive step.
Thank you in advance for your reply.
Best regards
Etienne
Le 29/11/2016 à 19:17, Bernd Flemisch a écrit :
Hi Mustapha,
in Dumux 2.10, you have to change the linear solver to the AMGBackend
to perform parallel runs. Add in the problem file:
#include <dumux/linear/amgbackend.hh>
...
SET_TYPE_PROP(LensProblem, LinearSolver, AMGBackend<TypeTag>);
We are setting the default to AMG soon so that you don't have this
inconvenience anymore.
Let us know if this solves your problem.
Kind regards
Bernd
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
