Mercurial > octave-nkf
changeset 7360:1e716cb4b2a3
[project @ 2008-01-11 02:59:02 by jwe]
| author | jwe |
|---|---|
| date | Fri, 11 Jan 2008 02:59:02 +0000 |
| parents | 8fb8d6985395 |
| children | a2870fd8ac58 |
| files | scripts/ChangeLog scripts/polynomial/mpoles.m |
| diffstat | 2 files changed, 9 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/ChangeLog Fri Jan 11 02:40:46 2008 +0000 +++ b/scripts/ChangeLog Fri Jan 11 02:59:02 2008 +0000 @@ -1,3 +1,8 @@ +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-10 John W. Eaton <jwe@octave.org> * plot/gnuplot_binary.m: New file.
--- a/scripts/polynomial/mpoles.m Fri Jan 11 02:40:46 2008 +0000 +++ b/scripts/polynomial/mpoles.m Fri Jan 11 02:59:02 2008 +0000 @@ -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];