Mercurial > hg > octave-terminal
changeset 11612:204d16b1d17f release-3-0-x
[3-0-0-branch @ 2008-01-11 02:58:58 by jwe]
author | jwe |
---|---|
date | Fri, 11 Jan 2008 02:58:58 +0000 |
parents | 807cf4d67705 |
children | 96d25ed408e5 |
files | scripts/ChangeLog scripts/polynomial/mpoles.m |
diffstat | 2 files changed, 9 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/ChangeLog +++ b/scripts/ChangeLog @@ -1,9 +1,13 @@ +2008-01-10 Ben Abbott <bpabbott@mac.com> + + * polynomial/mpoles.m: Avoid cases where poles could be assigned + to more than one multiplicity group. + 2008-01-09 John W. Eaton <jwe@octave.org> * plot/drawnow.m: Fail if filename includes a directory part that does not exist. - 2008-01-07 John W. Eaton <jwe@octave.org> * miscellaneous/copyfile.m, miscellaneous/movefile.m:
--- a/scripts/polynomial/mpoles.m +++ b/scripts/polynomial/mpoles.m @@ -100,6 +100,10 @@ p0 = abs (p(n)); endif k = find (dp < tol * p0); + ## Poles can only be members of one multiplicity group. + if (numel (indx)) + k = k(! ismember (k, indx)); + endif m = 1:numel (k); multp(k) = m; indx = [indx; k];