octave_packages/m/geometry/dsearch.m

   1 ## Copyright (C) 2007-2012 David Bateman
   2 ##
   3 ## This file is part of Octave.
   4 ##
   5 ## Octave is free software; you can redistribute it and/or modify it
   6 ## under the terms of the GNU General Public License as published by
   7 ## the Free Software Foundation; either version 3 of the License, or (at
   8 ## your option) any later version.
   9 ##
  10 ## Octave is distributed in the hope that it will be useful, but
  11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
  12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13 ## General Public License for more details.
  14 ##
  15 ## You should have received a copy of the GNU General Public License
  16 ## along with Octave; see the file COPYING.  If not, see
  17 ## <http://www.gnu.org/licenses/>.
  18
  19 ## -*- texinfo -*-
  20 ## @deftypefn  {Function File} {@var{idx} =} dsearch (@var{x}, @var{y}, @var{tri}, @var{xi}, @var{yi})
  21 ## @deftypefnx {Function File} {@var{idx} =} dsearch (@var{x}, @var{y}, @var{tri}, @var{xi}, @var{yi}, @var{s})
  22 ## Return the index @var{idx} or the closest point in @code{@var{x}, @var{y}}
  23 ## to the elements @code{[@var{xi}(:), @var{yi}(:)]}.  The variable @var{s} is
  24 ## accepted for compatibility but is ignored.
  25 ## @seealso{dsearchn, tsearch}
  26 ## @end deftypefn
  27
  28 function idx = dsearch (x, y, tri, xi, yi, s)
  29   if (nargin < 5 || nargin > 6)
  30     print_usage ();
  31   endif
  32   idx = __dsearchn__ ([x(:), y(:)], [xi(:), yi(:)]);
  33 endfunction
  34
  35 %!shared x, y, tri
  36 %! x = [-1;-1;1];
  37 %! y = [-1;1;-1];
  38 %! tri = [1,2,3];
  39 %!assert (dsearch(x,y,tri,1,1/3), 3);
  40 %!assert (dsearch(x,y,tri,1/3,1), 2);