[Pkg-octave-devel] Bug#348480: Bug#348480: Bug#348480: sample data

[Rafael Laboissiere]

* Thomas Weber <thomas.weber.mail at gmail.com> [2008-09-06 11:09]:

> On Sat, Sep 06, 2008 at 10:26:20AM +0200, Rafael Laboissiere wrote:
> > * Thomas Weber <thomas.weber.mail at gmail.com> [2008-05-17 23:22]:
> > 
> > > Upstream changed the polyfit algorithm in the 3.1 branch to a QR based
> > > algorithm. Your sample is included as a test case, so once we come to
> > > that version, this bug should be fixed. 
> > 
> > It does not seem to be the case:
> > 
> > $ octave3.1 -q
> > octave3.1:1> inputData = load( "test.dat" );
> > octave3.1:2> x = inputData( : , 1 );
> > octave3.1:3> y = inputData( : , 2 );
> > octave3.1:4> P = polyfit(x,y,10);
> > warning: dgelsd: rank deficient 13x11 matrix, rank = 3
> > octave3.1:5> P
> > P =
> > 
> >    1.9998e-24
> >    5.2878e-21
> >    3.5391e-18
> >   -8.9277e-21
> >    1.5014e-23
> >   -2.1042e-26
> >    2.6541e-29
> >   -3.1245e-32
> >    3.5032e-35
> >   -3.7875e-38
> >    3.9812e-41
> > 
> > This is with version 3.1.51.  Thomas, do you know whether it is fixed in hg?
> 
> It should be in 3.1.51, the sample is one test in polyfit.m. I've just
> copied the polyfit from the ss-3-1-51 tag in my 3.0 installation and it
> seems to work fine; however, my P is a row vector?

You are right, the problem is fixed in 3.1.51.  I had a wrong polyfit.m in
my path when I did the test above.  Sorry.

-- 
Rafael