changelog shortlog tags changeset files revisions annotate raw

scripts/geometry/dsearchn.m

changeset 10289: 4b124317dc38
parent:1bf0ce0930be
author: John W. Eaton <jwe@octave.org>
date: Tue Feb 09 20:58:55 2010 -0500 (70 minutes ago)
permissions: -rw-r--r--
description: base_properties::set_children: account for hidden children 
     1## Copyright (C) 2007, 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} =} dsearchn (@var{x}, @var{tri}, @var{xi})

    21## @deftypefnx {Function File} {@var{idx} =} dsearchn (@var{x}, @var{tri}, @var{xi}, @var{outval})

    22## @deftypefnx {Function File} {@var{idx} =} dsearchn (@var{x}, @var{xi})

    23## @deftypefnx {Function File} {[@var{idx}, @var{d}] =} dsearchn (@dots{})

    24## Returns the index @var{idx} or the closest point in @var{x} to the elements

    25## @var{xi}.  If @var{outval} is supplied, then the values of @var{xi} that are

    26## not contained within one of the simplicies @var{tri} are set to 

    27## @var{outval}.  Generally, @var{tri} is returned from @code{delaunayn 

    28## (@var{x})}.

    29## @seealso{dsearch, tsearch}

    30## @end deftypefn

    31

    32function [idx, d] = dsearchn (x, t, xi, outval)

    33  if (nargin < 2 || nargin > 4)

    34    print_usage ();

    35  endif

    36

    37  if (nargin == 2)

    38    [idx, d] = __dsearchn__ (x, t);

    39  else

    40    [idx, d] = __dsearchn__ (x, xi);

    41    if (nargin == 4)

    42      idx2 = isnan (tsearchn (x, t, xi));

    43      idx(idx2) = outval;

    44      d(idx2) = outval;

    45    endif

    46  endif

    47endfunction

    48

    49%!shared x, tri

    50%! x = [-1,-1;-1,1;1,-1]; 

    51%! tri = [1,2,3]; 

    52%!assert (dsearchn(x,tri,[1,1/3]), 3);

    53%!assert (dsearchn(x,tri,[1,1/3],NaN), NaN);

    54%!assert (dsearchn(x,tri,[1,1/3],NA), NA);

    55%!assert (dsearchn(x,tri,[1/3,1]), 2);

    56%!assert (dsearchn(x,tri,[1/3,1],NaN), NaN);

    57%!assert (dsearchn(x,tri,[1/3,1],NA), NA);