[Trilinos-Users] trouble with Ifpack block relaxation preconditioning

[Alicia M Klinvex]

Hello,

I'm having trouble with the block relaxation preconditioners in Ifpack.  I generated a very small 5x5 diagonally dominant symmetric matrix, and I try to compute a block Jacobi preconditioner to use with pseudo-block CG.  (I started with the Trilinos example PseudoBlockPrecCGEpetraExFile.cpp and modified it to use block-Jacobi preconditioning rather than incomplete Choleski.)  If I run the code with 1 MPI process and 1 right hand side, it terminates successfully.  If I try to use more than 1 MPI process (or more than 1 right hand side) the program crashes as in the following output:

./PseudoBlockPrecCGEpetraExFile.exe --filename=test.mtx --num-rhs=2 
================================================================================
Ifpack_BlockRelaxation, Jacobi
Sweeps = 1
Damping factor = 1, using zero starting solution
Number of local blocks = 1
Global number of rows            = 5

Phase           # calls   Total Time (s)       Total MFlops     MFlops/s
-----           -------   --------------       ------------     --------
Initialize()        1       0.00310206                0                0
Compute()           1      0.000350952                0                0
ApplyInverse()      0                0                0                0
================================================================================

terminate called after throwing an instance of 'int'
Abort

Does anyone know what might be causing this problem?

The code, makefile, and matrix are available here:
https://www.dropbox.com/sh/x2yvkfb92ppf0h8/jZ-sRh5fAV

Thank you,
Alicia Klinvex
aklinvex at purdue.edu