* Thomas Weber <[EMAIL PROTECTED]> [2008-09-06 11:09]: > On Sat, Sep 06, 2008 at 10:26:20AM +0200, Rafael Laboissiere wrote: > > * Thomas Weber <[EMAIL PROTECTED]> [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 -- To UNSUBSCRIBE, email to [EMAIL PROTECTED] with a subject of "unsubscribe". Trouble? Contact [EMAIL PROTECTED]