(This is Debian bug #774640: https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=774640)
I crafted a minimal example of a suspicious error with DGESDD. The example can be downloaded from there:
- Fortran source: http://people.debian.org/~sebastien/bug774640.f90
- Matrix data: http://people.debian.org/~sebastien/bug774640.data
This example crashes in LAPACK 3.5.0 with the following message:
- Code: Select all
** On entry to DLASCL parameter number 4 had an illegal value
I also tested current LAPACK version from SVN (r1594), and the problem seems to be still there.
However, the same program works well with LAPACK 3.4.1.
Cheers,