# HG changeset patch # User jwe # Date 1109119257 0 # Node ID 955ed0dcdc5fc4c0150159b7436245ed85838c81 # Parent 8ca032643f554541dbe9d1bd2efed1f18b6e3961 [project @ 2005-02-23 00:40:57 by jwe] diff -r 8ca032643f55 -r 955ed0dcdc5f scripts/ChangeLog --- a/scripts/ChangeLog Wed Feb 23 00:18:58 2005 +0000 +++ b/scripts/ChangeLog Wed Feb 23 00:40:57 2005 +0000 @@ -1,3 +1,10 @@ +2005-02-22 John W. Eaton + + * polynomial/residue.m: Force prepad to always create row vectors. + + * polynomial/poly.m: Quick return if m is 0. + From Carmen Navarrete . + 2005-02-21 David Bateman * statistics/distributions/poisson_rnd.m: fix for lambda of zero. diff -r 8ca032643f55 -r 955ed0dcdc5f scripts/polynomial/poly.m --- a/scripts/polynomial/poly.m Wed Feb 23 00:18:58 2005 +0000 +++ b/scripts/polynomial/poly.m Wed Feb 23 00:40:57 2005 +0000 @@ -40,6 +40,7 @@ n = max (size (x)); if (m == 0) y = 1; + return; elseif (m == 1) v = x; elseif (m == n) diff -r 8ca032643f55 -r 955ed0dcdc5f scripts/polynomial/residue.m --- a/scripts/polynomial/residue.m Wed Feb 23 00:18:58 2005 +0000 +++ b/scripts/polynomial/residue.m Wed Feb 23 00:40:57 2005 +0000 @@ -294,7 +294,7 @@ else cp = p (1:mpi-1); endif - rhs (1, rhi:rhi+lp-1) = prepad (poly (cp), lp); + rhs (1, rhi:rhi+lp-1) = prepad (poly (cp), lp, 0, 2); rhi = rhi + lp; endfor mpi = mpi + M (dpi); @@ -302,11 +302,11 @@ endwhile if (MM > 1) for index = 2:MM - lhs (index, :) = prepad (polyderiv (lhs (index-1, :)), lb); + lhs (index, :) = prepad (polyderiv (lhs (index-1, :)), lb, 0, 2); ind = 1; for rhi = 1:lp cp = rhs (index-1, ind:ind+lp-1); - rhs (index, ind:ind+lp-1) = prepad (polyderiv (cp), lp); + rhs (index, ind:ind+lp-1) = prepad (polyderiv (cp), lp, 0, 2); ind = ind + lp; endfor endfor