Mercurial > hg > octave-nkf

annotate scripts/general/interp3.m @ 19898:4197fc428c7d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Update copyright notices for 2015.
author John W. Eaton <jwe@octave.org>
date Wed, 11 Feb 2015 14:19:08 -0500
parents 5a59c0e1203d
children ca7599ae464d
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
19898
4197fc428c7d maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents: 19859
diff changeset
1 ## Copyright (C) 2007-2015 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})
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
27 ## Three-dimensional interpolation.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
28 ##
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
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: 18722
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: 18722
diff changeset
32 ## @code{meshgrid}, in which case the sizes of
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
diff changeset
34 ## @var{z} are vectors describing a cubic grid then
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
35 ## @code{length (@var{x}) == columns (@var{v})},
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
36 ## @code{length (@var{y}) == rows (@var{v})},
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
diff changeset
38 ## data must be strictly monotonic.
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
39 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
diff changeset
41 ## data matrix @var{v}, the 3-D region
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
diff changeset
43 ## @var{z} = 1:size (@var{v}, 3)} is assumed.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
diff changeset
45 ## not important.
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
46 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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: 18722
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: 18722
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: 18722
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: 18722
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: 18722
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: 18722
diff changeset
53 ## points at [0, 1/4, 1/2, 3/4, 1].
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
54 ##
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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 ##
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
64 ## @item @qcode{"cubic"}
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
65 ## Piecewise cubic Hermite interpolating polynomial---shape-preserving
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
66 ## interpolation with smooth first derivative (not implemented yet).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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{"spline"}
12175
2090995ca588 Correct en-dash,em-dash instances in docstrings.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
69 ## Cubic spline interpolation---smooth first and second derivatives
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
70 ## throughout the curve.
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
71 ## @end table
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
72 ##
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
73 ## @var{extrapval} is a scalar number. It replaces values beyond the endpoints
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
74 ## with @var{extrapval}. Note that if @var{extrapval} is used, @var{method}
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
75 ## must be specified as well. If @var{extrapval} is omitted and the
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
76 ## @var{method} is @qcode{"spline"}, then the extrapolated values of the
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
77 ## @qcode{"spline"} are used. Otherwise the default @var{extrapval} value for
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
78 ## any other @var{method} is @qcode{"NA"}.
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
82 ## FIXME: Need to add support for 'cubic' method (maybe change interpn).
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
83
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
84 function vi = interp3 (varargin)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
85
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
86 narginchk (1,9);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
87
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
88 method = "linear";
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
89 extrapval = [];
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
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
96 if (nargs > 1 && ischar (varargin{end-1}))
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
97 if (! isnumeric (varargin{end}) || ! isscalar (varargin{end}))
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
98 error ("interp3: EXTRAPVAL must be a numeric scalar");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
99 endif
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
100 extrapval = varargin{end};
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
101 method = varargin{end-1};
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
102 nargs -= 2;
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
103 elseif (ischar (varargin{end}))
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
104 method = varargin{end};
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
105 nargs--;
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
106 endif
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
107
18722
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
108 if (method(1) == "*")
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
109 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
110 method(1) = [];
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
111 endif
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
112 method = validatestring (method, ...
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
113 {"nearest", "linear", "cubic", "spline"});
18722
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
114
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
115 if (nargs < 3)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
116 ## Calling form interp3 (v) OR interp3 (v, n)
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
117 v = varargin{1};
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
118 if (ndims (v) != 3)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
119 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
120 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
121 n = varargin(2:nargs);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
122 v = permute (v, [2, 1, 3]);
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
123 if (isempty (extrapval))
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
124 vi = interpn (v, n{:}, method);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
125 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
126 vi = interpn (v, n{:}, method, extrapval);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
127 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
128 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
129
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
130 elseif (nargs == 4 && ! isvector (varargin{1}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
131 ## Calling form interp3 (v, xi, yi, zi)
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
132 v = varargin{1};
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
133 if (ndims (v) != 3)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
134 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
135 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
136 xi = varargin(2:4);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
137 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
138 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
139 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
140 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
141 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
142 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
143 xi{i} = permute (xi{i}, [2, 1, 3]);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
144 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
145 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
146 v = permute (v, [2, 1, 3]);
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
147 if (isempty (extrapval))
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
148 vi = interpn (v, xi{:}, method);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
149 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
150 vi = interpn (v, xi{:}, method, extrapval);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
151 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
152 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
153
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
154 elseif (nargs == 7)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
155 ## 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
156 v = varargin{4};
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
157 if (ndims (v) != 3)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
158 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
159 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
160 x = varargin(1:3);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
161 if (any (! cellfun (@isvector, x)))
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
162 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
163 if (! size_equal (x{:}, v))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
164 error ("interp3: X, Y, Z, and V dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
165 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
166 for i = 1 : 3
10549
95c3e38098bf Untabify .m scripts
Rik <code@nomad.inbox5.com>
parents: 9246
diff changeset
167 x{i} = permute (x{i}, [2, 1, 3]);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
168 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
169 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
170 xi = varargin(5:7);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
171 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
172 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
173 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
174 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
175 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
176 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
177 xi{i} = permute (xi{i}, [2, 1, 3]);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
178 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
179 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
180 v = permute (v, [2, 1, 3]);
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
181 if (isempty (extrapval))
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
182 vi = interpn (x{:}, v, xi{:}, method);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
183 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
184 vi = interpn (x{:}, v, xi{:}, method, extrapval);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
185 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
186 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
187
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
188 else
10635
d1978e7364ad Print name of function in error() string messages.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
189 error ("interp3: wrong number or incorrectly formatted input arguments");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
190 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
191
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
192 endfunction
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
193
14363
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents: 14359
diff changeset
194
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
195 %% 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: 18722
diff changeset
196
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
197 %!test # basic test
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
198 %! x = y = z = -1:1; y = y + 2;
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
199 %! f = @(x,y,z) x.^2 - y - z.^2;
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
200 %! [xx, yy, zz] = meshgrid (x, y, z);
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
201 %! v = f (xx,yy,zz);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
202 %! 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
203 %! [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
204 %! 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
205 %! [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
206 %! vi2 = interpn (y, x, z, v, xxi, yyi, zzi);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
207 %! 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
208
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
209 %!test # meshgridded xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
210 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
211 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
19794
db92e7e28e1f strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 19793
diff changeset
212 %! 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
213 %! [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
214 %! [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
215 %! vi = interp3 (x, y, z, v, xxi3, yyi3, zzi3, "nearest");
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
216 %! 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
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
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
219 %!test # vector xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
220 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
221 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
19794
db92e7e28e1f strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 19793
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 (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
224 %! 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
225 %! 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
226
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
227 %!test # vector xi+1 with extrap value
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
228 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
229 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
19794
db92e7e28e1f strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 19793
diff changeset
230 %! xi = zi = .6:1.6; yi = 1;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
231 %! vi = interp3 (x, y, z, v, xi+1, yi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
232 %! 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
233 %! 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
234
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
235 %!test # input value matrix--no x,y,z
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
236 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
237 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
19794
db92e7e28e1f strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 19793
diff changeset
238 %! 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
239 %! 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
240 %! vi2 = interpn (v, yi, xi, zi,"nearest");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
241 %! assert (vi, vi2);
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
242
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
243 %!test # input value matrix--no x,y,z, with extrap value
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
244 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
245 %! v = ones ([3,2,2]); v(:,2,1) = [7;5;4]; v(:,1,2) = [2;3;5];
19794
db92e7e28e1f strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 19793
diff changeset
246 %! xi = zi = .6:1.6; yi = 1;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
247 %! vi = interp3 (v, xi, yi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
248 %! vi2 = interpn (v, yi, xi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
249 %! assert (vi, vi2);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
250
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
251 %!test # extrapolation
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
252 %! X=[0,0.5,1]; Y=X; Z=X;
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
253 %! V = zeros (3,3,3);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
254 %! V(:,:,1) = [1 3 5; 3 5 7; 5 7 9];
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
255 %! V(:,:,2) = V(:,:,1) + 2;
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
256 %! V(:,:,3) = V(:,:,2) + 2;
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
257 %! tol = 10 * eps;
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
258 %! x=[-0.1,0,0.1]; y=x; z=x;
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
259 %! assert(interp3(X,Y,Z,V,x,y,z,"spline"), [-0.2, 1.0, 2.2]',tol);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
260 %! assert(interp3(X,Y,Z,V,x,y,z,"linear"), [NA, 1.0, 2.2]',tol);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
261 %! assert(interp3(X,Y,Z,V,x,y,z,"spline", 0), [0, 1.0, 2.2]',tol);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
262 %! assert(interp3(X,Y,Z,V,x,y,z,"linear", 0), [0, 1.0, 2.2]',tol);
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
263
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
264 %!shared z, zout, tol
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
265 %! z = zeros (3, 3, 3);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
266 %! zout = zeros (5, 5, 5);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
267 %! z(:,:,1) = [1 3 5; 3 5 7; 5 7 9];
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
268 %! z(:,:,2) = z(:,:,1) + 2;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
269 %! z(:,:,3) = z(:,:,2) + 2;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
270 %! for n = 1:5
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
271 %! zout(:,:,n) = [1 2 3 4 5;
19790
446c46af4b42 strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents: 17744
diff changeset
272 %! 2 3 4 5 6;
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
273 %! 3 4 5 6 7;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
274 %! 4 5 6 7 8;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
275 %! 5 6 7 8 9] + (n-1);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
276 %! end
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
277 %! tol = 10 * eps;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
278 %!
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
279 %!assert (interp3 (z), zout, tol)
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
280 %!assert (interp3 (z, "linear"), zout, tol)
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
281 %!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
282
18722
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
283 %% Test input validation
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
284 %!error interp3 ()
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
285 %!error interp3 ({1})
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
286 %!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: 18722
diff changeset
287 %!error <EXTRAPVAL must be a numeric scalar> interp3 (1,2,3,4,1,2,3,"linear", ones (2,2))
18722
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
288 %!warning <ignoring unsupported '\*' flag> interp3 (rand (3,3,3), 1, "*linear");
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
289 %!error <V must be a 3-D array> interp3 (rand (2,2))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
290 %!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: 18722
diff changeset
291 %!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: 18722
diff changeset
292 %!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: 18722
diff changeset
293 %!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: 18722
diff changeset
294 %!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: 18722
diff changeset
295