From bug-octave-request at bevo dot che dot wisc dot edu Mon Nov 8 12:26:28 1999 Subject: Re: Test fails - octave-2.1.20 on irix6.5 From: Joao Cardoso To: bug-octave at bevo dot che dot wisc dot edu Date: Mon, 08 Nov 1999 18:26:14 +0000 "Billinghurst, David (RTD)" wrote: > Just built octave-2.1.20 on irix6.5 with both gcc-2.95.1 and egcs 19991102 > snapshot. > > I get a failure on test poly/residue-1.m as it misses the allowed tolerances > > by an order of magnitude. Haven't done a full error analysis, but I suspect > the tolerances may be a little tight. > > The test case poly/residue-1.m is > > b = [1, 1, 1]; > a = [1, -5, 8, -4]; > [r, p, k, e] = residue (b, a); > (abs (r - [-2; 7; 3]) < sqrt (eps) > && abs (p - [2; 2; 1]) < sqrt (eps) > && isempty (k) > && e == [1; 2; 1]) > > For octave-2.1.20 this fails as > > r - [-2; 7; 3] = > -1.0658e-14 > 1.7276e-07 > 1.1990e-14 > > while sqrt (eps) = 1.4901e-08 I also have this problem with SCO. However the same test in octave-2.1.14 runs OK. merlin [jcard-10] octave GNU Octave, version 2.1.14 (i686-pc-sco3.2v5.0.4). octave:1> b = [1, 1, 1];a = [1, -5, 8, -4];[r, p, k, e] = residue (b, a); r - [-2; 7; 3] ans = 2.6645e-15 -4.4409e-15 -3.5527e-15 Joao -- Joao Cardoso | e-mail: jcardoso at inescn dot pt INESC, R. Jose Falcao 110 | tel: + 351 2 2094322 4050 Porto, Portugal | fax: + 351 2 2008487 --------------------------------------------------------------------- Octave is freely available under the terms of the GNU GPL. To ensure that development continues, see www.che.wisc.edu/octave/giftform.html Instructions for unsubscribing: www.che.wisc.edu/octave/archive.html ---------------------------------------------------------------------