# HG changeset patch # User jwe # Date 1192035203 0 # Node ID cc049a392a97ae15bce619e20b45b73a05eb938c # Parent 6d0d8d6215321195af6866a9d4586569266706d8 [project @ 2007-10-10 16:51:51 by jwe] diff -r 6d0d8d621532 -r cc049a392a97 scripts/ChangeLog --- a/scripts/ChangeLog Wed Oct 10 16:41:16 2007 +0000 +++ b/scripts/ChangeLog Wed Oct 10 16:53:23 2007 +0000 @@ -1,3 +1,7 @@ +2007-10-10 Ben Abbott + + * polynomial/mpoles.m: Return indx = ordr(indx), not indx(ordr). + 2007-10-10 John W. Eaton * general/num2str.m: Always allow for sign for automatically diff -r 6d0d8d621532 -r cc049a392a97 scripts/polynomial/mpoles.m --- a/scripts/polynomial/mpoles.m Wed Oct 10 16:41:16 2007 +0000 +++ b/scripts/polynomial/mpoles.m Wed Oct 10 16:53:23 2007 +0000 @@ -107,6 +107,6 @@ n = find (multp == 0, 1); endwhile multp = multp(indx); - indx = indx(ordr); + indx = ordr(indx); endfunction diff -r 6d0d8d621532 -r cc049a392a97 scripts/polynomial/residue.m --- a/scripts/polynomial/residue.m Wed Oct 10 16:41:16 2007 +0000 +++ b/scripts/polynomial/residue.m Wed Oct 10 16:53:23 2007 +0000 @@ -302,14 +302,14 @@ %! b = [1, 1, 1]; %! a = [1, -5, 8, -4]; %! [r, p, k, e] = residue (b, a); -%! assert ((abs (r - [-2; 7; 3]) < 1e-5 -%! && abs (p - [2; 2; 1]) < 1e-7 -%! && isempty (k) -%! && e == [1; 2; 1])); +%! assert (abs (r - [-2; 7; 3]) < 1e-5 +%! && abs (p - [2; 2; 1]) < 1e-7 +%! && isempty (k) +%! && e == [1; 2; 1]); %! k = [1 0]; %! [b, a] = residue (r, p, k); %! assert ((abs (b - [1, -5, 9, -3, 1]) < 1e-12 -%! && abs (a - [1, -5, 8, -4]) < 1e-12)); +%! && abs (a - [1, -5, 8, -4]) < 1e-12)); %!test %! b = [1, 0, 1]; @@ -318,11 +318,9 @@ %! r1 = [-5i; 12; +5i; 12]/54; %! r2 = conj(r1); %! p1 = [+3i; +3i; -3i; -3i]; -%! p2 = conj(p1); -%! assert ((((abs (r - r1) < 1e-7) && (abs (p - p1) < 1e-7)) -%! || ((abs (r - r2) < 1e-7) && (abs (p - p2) < 1e-7))) -%! && isempty (k) -%! && (e == [1; 2; 1; 2])); +%! assert (abs (r - r1) < 1e-7 && abs (p - p1) < 1e-7 +%! && isempty (k) +%! && e == [1; 2; 1; 2]); %! [br, ar] = residue (r, p, k); %! assert ((abs (br - b) < 1e-12 -%! && abs (ar - a) < 1e-12)); +%! && abs (ar - a) < 1e-12));