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]

Reply via email to