Mercurial > hg > octave-lyh

comparison scripts/general/interpn.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 a938cd7869b2
children 6c69f5cda32b
comparison
equal deleted inserted replaced
7670:7a5dbd31eb76 7671:4fbaba9abec1
142 elseif (strcmp (method, "nearest")) 142 elseif (strcmp (method, "nearest"))
143 yshape = size (y{1}); 143 yshape = size (y{1});
144 yidx = cell (1, nd); 144 yidx = cell (1, nd);
145 for i = 1 : nd 145 for i = 1 : nd
146 y{i} = y{i}(:); 146 y{i} = y{i}(:);
147 yidx{i} = lookup (x{i}(2:end-1), y{i}) + 1; 147 yidx{i} = lookup (x{i}, y{i}, "lr");
148 endfor 148 endfor
149 idx = cell (1,nd); 149 idx = cell (1,nd);
150 for i = 1 : nd 150 for i = 1 : nd
151 idx{i} = yidx{i} + (y{i} - x{i}(yidx{i}) > x{i}(yidx{i} + 1) - y{i}); 151 idx{i} = yidx{i} + (y{i} - x{i}(yidx{i}) > x{i}(yidx{i} + 1) - y{i});
152 endfor 152 endfor