octave-nkf: scripts/general/interp3.m annotate

annotate scripts/general/interp3.m @ 19630:0e1f5a750d00

maint: Periodic merge of gui-release to default.

author	John W. Eaton <jwe@octave.org>
date	Tue, 20 Jan 2015 10:24:46 -0500
parents	0ede4dbb37f1 446c46af4b42
children	db92e7e28e1f

rev	line source
17744 d63878346099 maint: Update copyright notices for release. John W. Eaton <jwe@octave.org> parents: 17281 diff changeset	1 ## Copyright (C) 2007-2013 David Bateman
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	2 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	3 ## This file is part of Octave.
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	4 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	6 ## under the terms of the GNU General Public License as published by
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 7001 diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 7001 diff changeset	8 ## your option) any later version.
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	9 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	13 ## General Public License for more details.
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	14 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	15 ## You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 7001 diff changeset	16 ## along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 7001 diff changeset	17 ## <http://www.gnu.org/licenses/>.
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	18
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	19 ## -- texinfo --
11563 3c6e8aaa9555 Grammarcheck m-files before 3.4 release. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	20 ## @deftypefn {Function File} {@var{vi} =} interp3 (@var{x}, @var{y}, @var{z}, @var{v}, @var{xi}, @var{yi}, @var{zi})
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	21 ## @deftypefnx {Function File} {@var{vi} =} interp3 (@var{v}, @var{xi}, @var{yi}, @var{zi})
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	22 ## @deftypefnx {Function File} {@var{vi} =} interp3 (@var{v}, @var{n})
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	23 ## @deftypefnx {Function File} {@var{vi} =} interp3 (@var{v})
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	24 ## @deftypefnx {Function File} {@var{vi} =} interp3 (@dots{}, @var{method})
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	25 ## @deftypefnx {Function File} {@var{vi} =} interp3 (@dots{}, @var{method}, @var{extrapval})
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	26 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	27 ## Three-dimensional interpolation.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	28 ##
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	29 ## Interpolate reference data @var{x}, @var{y}, @var{z}, @var{v} to determine
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	30 ## @var{vi} at the coordinates @var{xi}, @var{yi}, @var{zi}. The reference
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	31 ## data @var{x}, @var{y}, @var{z} can be matrices, as returned by
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	32 ## @code{meshgrid}, in which case the sizes of
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	33 ## @var{x}, @var{y}, @var{z}, and @var{v} must be equal. If @var{x}, @var{y},
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	34 ## @var{z} are vectors describing a cubic grid then
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	35 ## @code{length (@var{x}) == columns (@var{v})},
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	36 ## @code{length (@var{y}) == rows (@var{v})},
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	37 ## and @code{length (@var{z}) == size (@var{v}, 3)}. In either case the input
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	38 ## data must be strictly monotonic.
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	39 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	40 ## If called without @var{x}, @var{y}, @var{z}, and just a single reference
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	41 ## data matrix @var{v}, the 3-D region
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	42 ## @code{@var{x} = 1:columns (@var{v}), @var{y} = 1:rows (@var{v}),
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	43 ## @var{z} = 1:size (@var{v}, 3)} is assumed.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	44 ## This saves memory if the grid is regular and the distance between points is
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	45 ## not important.
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	46 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	47 ## If called with a single reference data matrix @var{v} and a refinement
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	48 ## value @var{n}, then perform interpolation over a 3-D grid where each original
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	49 ## interval has been recursively subdivided @var{n} times. This results in
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	50 ## @code{2^@var{n}-1} additional points for every interval in the original
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	51 ## grid. If @var{n} is omitted a value of 1 is used. As an example, the
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	52 ## interval [0,1] with @code{@var{n}==2} results in a refined interval with
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	53 ## points at [0, 1/4, 1/2, 3/4, 1].
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	54 ##
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	55 ## The interpolation @var{method} is one of:
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	56 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	57 ## @table @asis
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 15333 diff changeset	58 ## @item @qcode{"nearest"}
9070 e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi Rik <rdrider0-list@yahoo.com> parents: 9051 diff changeset	59 ## Return the nearest neighbor.
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10793 diff changeset	60 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	61 ## @item @qcode{"linear"} (default)
9070 e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi Rik <rdrider0-list@yahoo.com> parents: 9051 diff changeset	62 ## Linear interpolation from nearest neighbors.
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10793 diff changeset	63 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	64 ## @item @qcode{"pchip"}
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	65 ## Piecewise cubic Hermite interpolating polynomial---shape-preserving
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	66 ## interpolation with smooth first derivative (not implemented yet).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	67 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 15333 diff changeset	68 ## @item @qcode{"cubic"}
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	69 ## Cubic interpolation (same as @qcode{"pchip"} [not implemented yet]).
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10793 diff changeset	70 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 15333 diff changeset	71 ## @item @qcode{"spline"}
12175 2090995ca588 Correct en-dash,em-dash instances in docstrings. Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	72 ## Cubic spline interpolation---smooth first and second derivatives
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	73 ## throughout the curve.
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	74 ## @end table
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	75 ##
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	76 ## If @var{extrapval} is a number, then replace values beyond the endpoints
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	77 ## with that number. When unspecified, @var{extrapval} defaults to @code{NA}.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	78 ## Note that if @var{extrapval} is used, @var{method} must be specified as well.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	79 ## @seealso{interp1, interp2, interpn, meshgrid}
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	80 ## @end deftypefn
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	81
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	82 ## FIXME: Need to validate N argument (maybe change interpn).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	83 ## FIXME: Need to add support for 'pchip' method (maybe change interpn).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	84 ## FIXME: Need to add support for "extrap" string value (maybe change interpn).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	85
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	86 function vi = interp3 (varargin)
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	87
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	88 method = "linear";
6742 ebf96cc00ee9 [project @ 2007-06-18 16:27:52 by jwe] jwe parents: 6722 diff changeset	89 extrapval = NA;
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	90 nargs = nargin;
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	91
13151 e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	92 if (nargin < 1 \|\| ! isnumeric (varargin{1}))
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	93 print_usage ();
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	94 endif
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	95
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	96 if (ischar (varargin{end}))
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	97 method = varargin{end};
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	98 nargs--;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	99 elseif (nargs > 1 && ischar (varargin{end-1}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	100 ## FIXME: No support for "extrap" string
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	101 if (! isnumeric (varargin{end}) \|\| ! isscalar (varargin{end}))
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	102 error ("interp3: EXTRAPVAL must be a numeric scalar");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	103 endif
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	104 extrapval = varargin{end};
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	105 method = varargin{end-1};
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	106 nargs -= 2;
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	107 endif
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	108
18621 5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	109 if (method(1) == "*")
5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	110 warning ("interp3: ignoring unsupported '*' flag to METHOD");
5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	111 method(1) = [];
5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	112 endif
5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	113
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	114 if (nargs < 3)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	115 ## Calling form interp3 (v) OR interp3 (v, n)
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	116 v = varargin{1};
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	117 if (ndims (v) != 3)
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	118 error ("interp3: V must be a 3-D array of values");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	119 endif
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	120 n = varargin(2:nargs);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	121 v = permute (v, [2, 1, 3]);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	122 vi = ipermute (interpn (v, n{:}, method, extrapval), [2, 1, 3]);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	123
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	124 elseif (nargs == 4 && ! isvector (varargin{1}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	125 ## Calling form interp3 (v, xi, yi, zi)
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	126 v = varargin{1};
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	127 if (ndims (v) != 3)
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	128 error ("interp3: V must be a 3-D array of values");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	129 endif
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	130 xi = varargin(2:4);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	131 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	132 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	133 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	134 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	135 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	136 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	137 xi{i} = permute (xi{i}, [2, 1, 3]);
6722 5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	138 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	139 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	140 v = permute (v, [2, 1, 3]);
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	141 vi = ipermute (interpn (v, xi{:}, method, extrapval), [2, 1, 3]);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	142
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	143 elseif (nargs == 7)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	144 ## Calling form interp3 (x, y, z, v, xi, yi, zi)
7174 ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe] jwe parents: 7016 diff changeset	145 v = varargin{4};
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	146 if (ndims (v) != 3)
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	147 error ("interp3: V must be a 3-D array of values");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	148 endif
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	149 x = varargin(1:3);
6722 5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	150 if (any (! cellfun (@isvector, x)))
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	151 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	152 if (! size_equal (x{:}, v))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	153 error ("interp3: X, Y, Z, and V dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	154 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	155 for i = 1 : 3
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9246 diff changeset	156 x{i} = permute (x{i}, [2, 1, 3]);
6722 5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	157 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	158 endif
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	159 xi = varargin(5:7);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	160 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	161 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	162 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	163 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	164 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	165 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	166 xi{i} = permute (xi{i}, [2, 1, 3]);
6722 5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	167 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	168 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman] dbateman parents: 6702 diff changeset	169 v = permute (v, [2, 1, 3]);
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	170 vi = ipermute (interpn (x{:}, v, xi{:}, method, extrapval), [2, 1, 3]);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	171
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	172 else
10635 d1978e7364ad Print name of function in error() string messages. Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	173 error ("interp3: wrong number or incorrectly formatted input arguments");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	174 endif
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	175
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	176 endfunction
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	177
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	178
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	179 %% FIXME: Need some demo blocks here to show off the function like interp2.m.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	180
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	181 %!test # basic test
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	182 %! x = y = z = -1:1; y = y + 2;
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	183 %! f = @(x,y,z) x.^2 - y - z.^2;
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	184 %! [xx, yy, zz] = meshgrid (x, y, z);
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	185 %! v = f (xx,yy,zz);
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	186 %! xi = yi = zi = -1:0.5:1; yi = yi + 2.1;
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	187 %! [xxi, yyi, zzi] = meshgrid (xi, yi, zi);
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	188 %! vi = interp3 (x, y, z, v, xxi, yyi, zzi);
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	189 %! [xxi, yyi, zzi] = ndgrid (yi, xi, zi);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	190 %! vi2 = interpn (y, x, z, v, xxi, yyi, zzi);
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	191 %! assert (vi, vi2, 10*eps);
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	192
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	193 %!test # meshgridded xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	194 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	195 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	196 %! xi = zi = .6:1.6; yi = 1;
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	197 %! [xxi3, yyi3, zzi3] = meshgrid (xi, yi, zi);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	198 %! [xxi, yyi, zzi] = ndgrid (yi, xi, zi);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	199 %! vi = interp3 (x, y, z, v, xxi3, yyi3, zzi3, "nearest");
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	200 %! vi2 = interpn (y, x, z, v, xxi, yyi, zzi, "nearest");
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	201 %! assert (vi, vi2);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	202
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	203 %!test # vector xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	204 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	205 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	206 %! xi = zi = .6:1.6; yi = 1;
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	207 %! vi = interp3 (x, y, z, v, xi, yi, zi, "nearest");
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	208 %! vi2 = interpn (y, x, z, v, yi, xi, zi,"nearest");
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	209 %! assert (vi, vi2);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	210
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	211 %!test # vector xi+1 with extrap value
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	212 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	213 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	214 %! xi = zi = .6:1.6; yi = 1;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	215 %! vi = interp3 (x, y, z, v, xi+1, yi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	216 %! vi2 = interpn (y, x, z, v, yi, xi+1, zi, "nearest", 3);
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	217 %! assert (vi, vi2);
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	218
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	219 %!test # input value matrix--no x,y,z
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	220 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	221 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	222 %! xi = zi = .6:1.6; yi = 1;
15333 7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	223 %! vi = interp3 (v, xi, yi, zi, "nearest");
7d2eb4a01798 fix problem in interp3.m when method argument is passed, added tests, and made one test more stringent Kris Thielemans <kris.f.thielemans@gmail.com> parents: 14363 diff changeset	224 %! vi2 = interpn (v, yi, xi, zi,"nearest");
6702 b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman] dbateman parents: diff changeset	225 %! assert (vi, vi2);
13151 e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	226
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	227 %!test # input value matrix--no x,y,z, with extrap value
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	228 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	229 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	230 %! xi = zi = .6:1.6; yi = 1;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	231 %! vi = interp3 (v, xi, yi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	232 %! vi2 = interpn (v, yi, xi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	233 %! assert (vi, vi2);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	234
13151 e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	235 %!shared z, zout, tol
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	236 %! z = zeros (3, 3, 3);
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	237 %! zout = zeros (5, 5, 5);
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	238 %! z(:,:,1) = [1 3 5; 3 5 7; 5 7 9];
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	239 %! z(:,:,2) = z(:,:,1) + 2;
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	240 %! z(:,:,3) = z(:,:,2) + 2;
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	241 %! for n = 1:5
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	242 %! zout(:,:,n) = [1 2 3 4 5;
19627 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	243 %! 2 3 4 5 6;
13151 e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	244 %! 3 4 5 6 7;
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	245 %! 4 5 6 7 8;
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	246 %! 5 6 7 8 9] + (n-1);
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	247 %! end
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	248 %! tol = 10 * eps;
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	249 %!
13151 e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	250 %!assert (interp3 (z), zout, tol)
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	251 %!assert (interp3 (z, "linear"), zout, tol)
e173fda06fca Fix bug #30295. Ben Abbott <bpabbott@mac.com> parents: 13150 diff changeset	252 %!assert (interp3 (z, "spline"), zout, tol)
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	253
18621 5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	254 %% Test input validation
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	255 %!error interp3 ()
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	256 %!error interp3 ({1})
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	257 %!error <EXTRAPVAL must be a numeric scalar> interp3 (1,2,3,4,1,2,3,"linear", {1})
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	258 %!error <EXTRAPVAL must be a numeric scalar> interp3 (1,2,3,4,1,2,3,"linear", ones (2,2))
18621 5cf9a02732b6 Issue warning when '' flag used with interp2,3,n.m (bug #41459). Rik <rik@octave.org>* parents: 17744 diff changeset	259 %!warning <ignoring unsupported '\' flag> interp3 (rand (3,3,3), 1, "linear");
18640 0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	260 %!error <V must be a 3-D array> interp3 (rand (2,2))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	261 %!error <V must be a 3-D array> interp3 (rand (2,2), 1,1,1)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	262 %!error <XI, YI, and ZI dimensions must be equal> interp3 (rand (2,2,2), 1,1, ones (2,2))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	263 %!error <V must be a 3-D array> interp3 (1:2, 1:2, 1:2, rand (2,2), 1,1,1)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	264 %!error <X, Y, Z, and V dimensions must be equal> interp3 (ones(1,2,2), ones(2,2,2), ones(2,2,2), rand (2,2,2), 1,1,1)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	265 %!error <XI, YI, and ZI dimensions must be equal> interp3 (1:2, 1:2, 1:2, rand (2,2,2), 1,1, ones (2,2))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions. Rik <rik@octave.org> parents: 18621 diff changeset	266

Mercurial > octave-nkf

annotate scripts/general/interp3.m @ 19630:0e1f5a750d00