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?
Thomas
>
> --
> Rafael
>
>
>
> _______________________________________________
> Pkg-octave-devel mailing list
> [EMAIL PROTECTED]
> http://lists.alioth.debian.org/mailman/listinfo/pkg-octave-devel
--
To UNSUBSCRIBE, email to [EMAIL PROTECTED]
with a subject of "unsubscribe". Trouble? Contact [EMAIL PROTECTED]
