lapack error (example)

[Globe Trotter]

Here is an example: please try it and see what you get.

Included are files svdd-demo.c and array.h and test.dat

To compile, use

% gcc svdd-demo.c -ansi -Wall -pedantic -lm -llapack

% a.out


reading in matrix A from test.dat:
    2.117079     2.992091 
    2.106365     1.970807 
    2.131384     3.130390 
    2.344235     2.667219 
    1.838804     3.243629 


 ** On entry to DGESDD parameter number 12 had an illegal value
STOP 0

If you change m to <=4 it works fine. So there is an error for m "substantially
larger" than n (=2 here). Works fine on SuSE and RHEL3/4.

Thanks!

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com