changelog shortlog tags changeset files revisions annotate raw

scripts/geometry/dsearch.m

changeset 10289: 4b124317dc38
parent:1bf0ce0930be
author: John W. Eaton <jwe@octave.org>
date: Tue Feb 09 20:58:55 2010 -0500 (54 minutes ago)
permissions: -rw-r--r--
description: base_properties::set_children: account for hidden children 
     1## Copyright (C) 2007, 2008, 2009 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## Returns 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 but ignored for compatibility.

    25## @seealso{dsearchn, tsearch}

    26## @end deftypefn

    27

    28function idx = dsearch (x, y, t, xi, yi, s)

    29  if (nargin < 5 || nargin > 6)

    30    print_usage ();

    31  endif

    32  idx = __dsearchn__ ([x(:), y(:)], [xi(:), yi(:)]);

    33endfunction

    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);