#1035 Parallel AMG in non-overlapping ISTLBackend ISTLBackend_NOVLP_CG_AMG_SSOR crashes
Metadata
| Property | Value |
|---|---|
| Reported by | Eike Mueller (em459@bath.ac.uk) |
| Reported at | Feb 2, 2012 14:49 |
| Type | Bug Report |
| Version | Git (pre2.4) [autotools] |
| Operating System | Linux |
| Last edited by | Eike Mueller (em459@bath.ac.uk) |
| Last edited at | Feb 3, 2012 13:48 |
| Closed by | Eike Mueller (em459@bath.ac.uk) |
| Closed at | Feb 3, 2012 13:48 |
| Closed in version | Unknown |
| Resolution | Not a bug |
| Comment |
Description
The parallel AMG preconditioner crashes with following error message for the attached test code, which solves a Poisson problem in a cube (with YaspGrid) or on a spherical shell (implemented with ALUGrid) (see source dune_alugridnovlpamg.cc):
Dune reported error: MathError [mgc:/u/n/em459/Library/Dune_trunk_30Jan2012/include/dune/istl/paamg/amg.hh:670]: Coarse solver did not converge mpirun: killing job...
This happened both for the YaspGrid and the ALUGrid implementation when run on 8 processors. For YaspGrid it also crashes on two or four processes, but runs fine on just one process. However, on very rare occasions the 8 core run on YaspGrid did complete.
For YaspGrid, the code works if I increase the number of grid points in each direction from 16 to 32.
Use Dune revisions (trunk): dune-common: r6559 dune-grid: r7861 dune-geometry: r38 dune-pdelab: r1705 dune-istl: r1522 dune-localfunctions: r1037
I also attach the config.log file from the configure of the test code and a full output log from an 8 core run with YaspGrid (run.log). The compiler I used is g++44 (GCC) 4.4.4 20100726 (Red Hat 4.4.4-13).