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

Ben Abbott bpabbott at mac.com
Mon Jan 14 02:49:10 CST 2008


On Jan 14, 2008, at 2:25 AM, Jason Riedy wrote:

> 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).
>

Based upon the posts below, it may be that the lapack-3.1.0 sources  
are fixed

	https://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=80&start=0&postdays=0&postorder=asc&highlight=dgelsd
	https://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=278

I'll try to verify. In the meanwhile, if someone knows what version  
they are using and is failing/working, please let me know.

Ben


More information about the Bug-octave mailing list