Hi Dmitry, looks like you are running into a deadlock. Looks like a bug in the exception handler of the Newton (I assume https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/blob/master/dumux/nonlinear/newtonsolver.hh#L548 is wrong because the exception cannot be assumed to be thrown on all processes). Thanks for reporting this. Would you open a bug report?
Best wishes Timo On 29 Jun 2022, at 19:36, Dmitry Pavlov <[email protected]<mailto:[email protected]>> wrote: Hello, In one porous flow simulation, I get this message: Newton iteration 6 done, residual = 3.1654e+13 Solve: M deltax^k = rNewton: Caught exception: "FMatrixError [luDecomposition:/home/dpavlov/DUMUX/dune-common/dune/common/densematrix.hh:939]: matrix is singular" ... and it is fine. Matrix can be singular, no objection here. After this, the Newton solver reports that it did not converge and restarts with a smaller time step, and eventually it proceeds. However, when I run the same simulation with MPI, I get this Newton iteration 8 done, residual = 7.2437e+13 Solve: M deltax^k = and the program hangs. So it seems that it misses the singularity in the MPI mode and tries to proceed with the current time step, to no avail. How this can be cured? Best regards, Dmitry _______________________________________________ DuMux mailing list [email protected]<mailto:[email protected]> https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________ DuMux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
