### Think there is a bug in Zgeev

Posted:

**Fri May 21, 2010 12:44 pm**My task is to perform the iterative diagonalization of a complex-valued square matrix. Occasionally we try to diagonalize matrices that are "nearly diagonal": Most but not all of the off diagonal elements are either 0 or nearly 0.

For the nearly diagonal matrices, BLAS version of Zgeev works fine. CULA, however, gives me a "segmentation fault" error while calling culaZgeev routine.

I attached both BLAS and CULA versions of my program with a sample input matrix.

In both CULA and BLAS:

diag.c - calls BLAS or CULA routine

diagtest.c - wrapper. sets up variables and calls diag.

omegai.dat - imaginary part of input matrix

omegar.dat - real part of input matrix

Makefile - (make diagtest)

In BLAS version:

*.f - fortran dependencies

For the nearly diagonal matrices, BLAS version of Zgeev works fine. CULA, however, gives me a "segmentation fault" error while calling culaZgeev routine.

I attached both BLAS and CULA versions of my program with a sample input matrix.

In both CULA and BLAS:

diag.c - calls BLAS or CULA routine

diagtest.c - wrapper. sets up variables and calls diag.

omegai.dat - imaginary part of input matrix

omegar.dat - real part of input matrix

Makefile - (make diagtest)

In BLAS version:

*.f - fortran dependencies