[Trilinos-Users] trouble with Ifpack block relaxation preconditioning
Alicia M Klinvex
aklinvex at purdue.edu
Mon Feb 4 19:41:23 MST 2013
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
More information about the Trilinos-Users
mailing list