Think there is a bug in Zgeev
3 posts
• Page 1 of 1
Think there is a bug in Zgeev
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
- wk-cof
- CULA Academic
- Posts: 2
- Joined: Wed Mar 17, 2010 7:47 am
Re:Think there is a bug in Zgeev
It doesn't look like the files were successfully attached. If you'd like, you can e-mail them to me: spagnoli@emphotonics.com
- kyle
- Administrator
- Posts: 301
- Joined: Fri Jun 12, 2009 7:47 pm
Re:Think there is a bug in Zgeev
I can confirm that this is a reproducible bug. The segmentation fault has since been fixed and will most likely be available in CULA 2.0 (final).
- kyle
- Administrator
- Posts: 301
- Joined: Fri Jun 12, 2009 7:47 pm
3 posts
• Page 1 of 1
Who is online
Users browsing this forum: No registered users and 0 guests