As another follow-up, when R says > svd(h9,0,0) Error in La.svd(x, nu, nv) : BLAS/LAPACK routine 'DGEBRD' gave error code -10 >
then this really comes from DGEBRD in dlapack2.f where ELSE IF( LWORK.LT.MAX( 1, M, N ) .AND. .NOT.LQUERY ) THEN INFO = -10 lwork is insufficiently sized, going back to my initial bug report about DGESDD giving the faulty size estimate. This really is the same bug. And as Ben showed in a follow-up Lapack 3.1 from experimental seems to cure it. Dirk -- Three out of two people have difficulties with fractions. -- To UNSUBSCRIBE, email to [EMAIL PROTECTED] with a subject of "unsubscribe". Trouble? Contact [EMAIL PROTECTED]