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

Sergei Steshenko sergstesh at yahoo.com
Mon Jan 14 18:18:07 CST 2008


--- Ben Abbott <bpabbott at mac.com> wrote:

> 
> 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
> _______________________________________________
> Bug-octave mailing list
> Bug-octave at octave.org
> https://www.cae.wisc.edu/mailman/listinfo/bug-octave
> 

I am using lapack-3.1.1, and I do see the bug in 'octave'.

--Sergei.


Applications From Scratch: http://appsfromscratch.berlios.de/


      ____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ 



More information about the Bug-octave mailing list