[Wien] Bug in seclr4.F (probably dsyxev4.f or dstebz2.f) for MKL_NUM_THREADS>1

Laurence Marks L-marks at northwestern.edu
Fri Apr 11 01:25:06 CEST 2008


Location of error appears to be line 258 of dsytrd2.f
            CALL DLATRD( UPLO, N-I+1, NB, A( I, I ), LDA, E( I ), &
                         TAU( I ), WORK, LDWORK )

Putting a "return" before this and the code works (albeit nonsense),
after this and it crashes.

On Thu, Apr 10, 2008 at 5:26 PM, Laurence Marks
<L-marks at northwestern.edu> wrote:
> The array IWORK is getting corrupted, and the deallocation on line 367
>  is failing with a SIGSEV.
>
>  Compiler: ifort 10.1.015, mkl 10.0.2.018
>  libpthread-2.5.so   glibc-2.5.s0 libgcc_s-4.1.1-20070105
>
>  --
>  Laurence Marks
>  Department of Materials Science and Engineering
>  MSE Rm 2036 Cook Hall
>  2220 N Campus Drive
>  Northwestern University
>  Evanston, IL 60208, USA
>  Tel: (847) 491-3996 Fax: (847) 491-7820
>  email: L-marks at northwestern dot edu
>  Web: www.numis.northwestern.edu
>  Commission on Electron Diffraction of IUCR
>  www.numis.northwestern.edu/IUCR_CED
>



-- 
Laurence Marks
Department of Materials Science and Engineering
MSE Rm 2036 Cook Hall
2220 N Campus Drive
Northwestern University
Evanston, IL 60208, USA
Tel: (847) 491-3996 Fax: (847) 491-7820
email: L-marks at northwestern dot edu
Web: www.numis.northwestern.edu
Commission on Electron Diffraction of IUCR
www.numis.northwestern.edu/IUCR_CED


More information about the Wien mailing list