blas error and crash when doing ldivide/rdivide of a scalar by a vector

Jason Riedy ejr at cs.berkeley.edu
Mon Jan 14 01:25:51 CST 2008


And John W. Eaton writes:
>
> How?  We call DGELSD with lwork == -1 and ask it to tell us how much
> work space to use.

Yeah, it's more painful than I thought, especially since some
people (well, me) switch LAPACK libraries at run time.

And worse, the documentation in {C,Z}GELSD is completely busted
so I can't claim that some obscure formula *must* work.  Ugh.
The problems stretch way back; I don't know how they lurked
unnoticed.  There may be no sane Octave-level workaround, but
I'll get it fixed in LAPACK before the next release (whenever
that is, no control over that part).

Jason


More information about the Bug-octave mailing list