[Dune] Error in implementing PAAMG using dune-istl

[Markus Blatt]

On Thu, Dec 14, 2017 at 01:46:04PM +0100, Kumar, Paras wrote:
> ***********************************************************************
> [i10hpc2:52204] *** Process received signal ***
> [i10hpc2:52204] Signal: Segmentation fault (11)
> [i10hpc2:52204] Signal code: Address not mapped (1)
> [i10hpc2:52204] Failing at address: 0x10
> [i10hpc2:52204] [ 0]
> /lib/x86_64-linux-gnu/libpthread.so.0(+0x11390)[0x7f932a47b390]
> [i10hpc2:52204] [ 1] ./my_Ex-05(_ZNK4Dune3Imp31compressed_base_array_unmanagedINS_11FieldMatrixIdLi1ELi1EEESaIS3_EEixEm+0x30)[0x63328a]
> [i10hpc2:52204] [ 2] ./my_Ex-05(_ZN4Dune3Amg19SymmetricDependencyINS_10BCRSMatrixINS_11FieldMatrixIdLi1ELi1EEESaIS4_EEENS0_13FirstDiagonalEE7initRowERKNS_3Imp27CompressedBlockVectorWindowIS4_S5_EEi+0xb3)[0x62238b]
> [i10hpc2:52204] [ 3]
> ./my_Ex-05(_ZN4Dune3Amg15buildDependencyINS0_15PropertiesGraphINS0_8SubGraphINS0_11MatrixGraphIKNS_10BCRSMatrixINS_11FieldMatrixIdLi1ELi1EEESaIS7_EEEEESt6vectorIbSaIbEEEENS0_16VertexPropertiesENS0_14EdgePropertiesENS_11IdentityMapENSF_12EdgeIndexMapEEENS0_16CoarsenCriterionINS0_18SymmetricCriterionIS9_NS0_13FirstDiagonalEEEEEEEvRT_RKNT0_6MatrixESS_b+0xd7)[0x625d75]

Well, it seems like you are using SymmetricCriterion and I the
sparsity pattern of the matrix is not symmetric. (Weiredly, you are
not using SymmetricCriterion for the sequential solve.

Please try whether switching to

typedef Dune::Amg::CoarsenCriterion<Dune::Amg::UnSymmetricCriterion<DuneMatrixType,Dune::Amg::FirstDiagonal> > Criterion;


fixes the issue.

Markus

-- 
Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de
Pedettistr. 38, 85072 Eichstätt, Germany
Tel.: +49 (0) 160 97590858