Mercurial > hg > octave-nkf

diff scripts/polynomial/ppval.m @ 7671:4fbaba9abec1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
implement compiled binary lookup
author Jaroslav Hajek <highegg@gmail.com>
date Fri, 28 Mar 2008 15:53:09 -0400
parents a1dbe9d80eee
children eb63fbe60fab
line wrap: on
line diff
--- a/scripts/polynomial/ppval.m
+++ b/scripts/polynomial/ppval.m
@@ -39,7 +39,7 @@
     transposed = (columns (xi) == 1);
     xi = xi(:);
     xn = length (xi);
-    idx = lookup (pp.x(2:pp.n), xi) + 1;
+    idx = lookup (pp.x, xi, "lr");
     dx = (xi - pp.x(idx)).';
     dx = reshape (dx(ones(1,prod(pp.d)),:),[pp.d,xn]);
     c = reshape (pp.P(:,1), pp.n, prod (pp.d));