Mercurial > hg > octave-nkf

annotate scripts/general/interp3.m @ 20830:b65888ec820e draft default tip gccjit

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
dmalcom gcc jit import
author Stefan Mahr <dac922@gmx.de>
date Fri, 27 Feb 2015 16:59:36 +0100
parents 7503499a252b
children
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
20368
7503499a252b doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents: 20038
diff changeset
32 ## @code{meshgrid}, in which case the sizes of @var{x}, @var{y}, @var{z}, and
7503499a252b doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents: 20038
diff changeset
33 ## @var{v} must be equal. If @var{x}, @var{y}, @var{z} are vectors describing
7503499a252b doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents: 20038
diff changeset
34 ## a cubic grid then @code{length (@var{x}) == columns (@var{v})},
7503499a252b doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents: 20038
diff changeset
35 ## @code{length (@var{y}) == rows (@var{v})}, and
7503499a252b doc: Update docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents: 20038
diff changeset
36 ## @code{length (@var{z}) == size (@var{v}, 3)}. In either case the input
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
37 ## data must be strictly monotonic.
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
38 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
39 ## 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
40 ## data matrix @var{v}, the 3-D region
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
41 ## @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
42 ## @var{z} = 1:size (@var{v}, 3)} is assumed.
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
43 ## 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
44 ## not important.
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
45 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
46 ## 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
47 ## 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
48 ## 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
49 ## @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
50 ## 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
51 ## 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
52 ## points at [0, 1/4, 1/2, 3/4, 1].
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
53 ##
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
54 ## The interpolation @var{method} is one of:
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
55 ##
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
56 ## @table @asis
17281
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted.
Rik <rik@octave.org>
parents: 15333
diff changeset
57 ## @item @qcode{"nearest"}
9070
e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi
Rik <rdrider0-list@yahoo.com>
parents: 9051
diff changeset
58 ## Return the nearest neighbor.
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10793
diff changeset
59 ##
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
60 ## @item @qcode{"linear"} (default)
9070
e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi
Rik <rdrider0-list@yahoo.com>
parents: 9051
diff changeset
61 ## Linear interpolation from nearest neighbors.
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10793
diff changeset
62 ##
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
63 ## @item @qcode{"cubic"}
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
64 ## Piecewise cubic Hermite interpolating polynomial---shape-preserving
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
65 ## interpolation with smooth first derivative (not implemented yet).
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
66 ##
17281
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted.
Rik <rik@octave.org>
parents: 15333
diff changeset
67 ## @item @qcode{"spline"}
12175
2090995ca588 Correct en-dash,em-dash instances in docstrings.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
68 ## Cubic spline interpolation---smooth first and second derivatives
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
69 ## throughout the curve.
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
70 ## @end table
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
71 ##
20000
ca7599ae464d doc: Grammarcheck documentation ahead of 4.0 release.
Rik <rik@octave.org>
parents: 19898
diff changeset
72 ## @var{extrapval} is a scalar number. It replaces values beyond the endpoints
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 ## 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
74 ## 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
75 ## @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
76 ## @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
77 ## any other @var{method} is @qcode{"NA"}.
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
78 ## @seealso{interp1, interp2, interpn, meshgrid}
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
79 ## @end deftypefn
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
80
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
81 ## 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
82
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
83 function vi = interp3 (varargin)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
84
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
85 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
86
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
87 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
88 extrapval = [];
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
89 nargs = nargin;
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
90
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
91 if (nargin < 1 || ! isnumeric (varargin{1}))
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
92 print_usage ();
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
93 endif
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
94
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
95 if (nargs > 1 && ischar (varargin{end-1}))
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
96 if (! isnumeric (varargin{end}) || ! isscalar (varargin{end}))
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
97 error ("interp3: EXTRAPVAL must be a numeric scalar");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
98 endif
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
99 extrapval = varargin{end};
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
100 method = varargin{end-1};
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
101 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
102 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
103 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
104 nargs--;
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
105 endif
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
106
18722
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
107 if (method(1) == "*")
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
108 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
109 method(1) = [];
5cf9a02732b6 Issue warning when '*' flag used with interp2,3,n.m (bug #41459).
Rik <rik@octave.org>
parents: 17744
diff changeset
110 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
111 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
112 {"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
113
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
114 if (nargs < 3)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
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
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
120 n = varargin(2:nargs);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
121 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
122 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
123 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
124 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
125 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
126 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
127 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
128
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
129 elseif (nargs == 4 && ! isvector (varargin{1}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
130 ## Calling form interp3 (v, xi, yi, zi)
7174
ff4a4cc863a9 [project @ 2007-11-14 19:16:12 by jwe]
jwe
parents: 7016
diff changeset
131 v = varargin{1};
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
132 if (ndims (v) != 3)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
133 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
134 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
135 xi = varargin(2:4);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
136 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
137 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
138 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
139 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
140 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
141 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
142 xi{i} = permute (xi{i}, [2, 1, 3]);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
143 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
144 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
145 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
146 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
147 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
148 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
149 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
150 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
151 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
152
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
153 elseif (nargs == 7)
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
154 ## 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
155 v = varargin{4};
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
156 if (ndims (v) != 3)
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
157 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
158 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
159 x = varargin(1:3);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
160 if (any (! cellfun (@isvector, x)))
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
161 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
162 if (! size_equal (x{:}, v))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
163 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
164 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
165 for i = 1 : 3
10549
95c3e38098bf Untabify .m scripts
Rik <code@nomad.inbox5.com>
parents: 9246
diff changeset
166 x{i} = permute (x{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
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
169 xi = varargin(5:7);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
170 if (any (! cellfun (@isvector, xi)))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
171 ## Meshgridded values rather than vectors
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
172 if (! size_equal (xi{:}))
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
173 error ("interp3: XI, YI, and ZI dimensions must be equal");
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
174 endif
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
175 for i = 1 : 3
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
176 xi{i} = permute (xi{i}, [2, 1, 3]);
6722
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
177 endfor
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
178 endif
5b09d433171c [project @ 2007-06-14 10:13:14 by dbateman]
dbateman
parents: 6702
diff changeset
179 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
180 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
181 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
182 else
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
183 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
184 endif
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
185 vi = ipermute (vi, [2, 1, 3]);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
186
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
187 else
10635
d1978e7364ad Print name of function in error() string messages.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
188 error ("interp3: wrong number or incorrectly formatted input arguments");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
189 endif
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
190
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
191 endfunction
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
192
14363
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents: 14359
diff changeset
193
20038
9fc020886ae9 maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents: 20000
diff changeset
194 ## FIXME: Need some demo blocks here to show off the function like interp2.m.
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
195
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
196 %!test # basic test
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
197 %! x = y = z = -1:1; y = y + 2;
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
198 %! f = @(x,y,z) x.^2 - y - z.^2;
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
199 %! [xx, yy, zz] = meshgrid (x, y, z);
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
200 %! v = f (xx,yy,zz);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
201 %! 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
202 %! [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
203 %! 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
204 %! [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
205 %! vi2 = interpn (y, x, z, v, xxi, yyi, zzi);
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
206 %! 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
207
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
208 %!test # meshgridded xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
209 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
210 %! 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
211 %! 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
212 %! [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
213 %! [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
214 %! 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
215 %! 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
216 %! 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
217
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
218 %!test # vector xi, yi, zi
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
219 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
220 %! 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
221 %! 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
222 %! 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
223 %! 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
224 %! 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
225
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
226 %!test # vector xi+1 with extrap value
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
227 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
228 %! 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
229 %! xi = zi = .6:1.6; yi = 1;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
230 %! 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
231 %! 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
232 %! 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
233
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
234 %!test # input value matrix--no x,y,z
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
235 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
236 %! 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
237 %! 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
238 %! 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
239 %! vi2 = interpn (v, yi, xi, zi,"nearest");
6702
b2391d403ed2 [project @ 2007-06-12 21:39:26 by dbateman]
dbateman
parents:
diff changeset
240 %! assert (vi, vi2);
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
241
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
242 %!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
243 %! x = z = 1:2; y = 1:3;
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
244 %! 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
245 %! xi = zi = .6:1.6; yi = 1;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
246 %! vi = interp3 (v, xi, yi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
247 %! vi2 = interpn (v, yi, xi, zi, "nearest", 3);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
248 %! assert (vi, vi2);
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
249
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
250 %!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
251 %! 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
252 %! 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
253 %! 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
254 %! 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
255 %! 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
256 %! tol = 10 * eps;
20038
9fc020886ae9 maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents: 20000
diff changeset
257 %! x=[-0.1,0,0.1]; y=x; z=x;
19859
5a59c0e1203d Modified the "extrap" option for interp2, interp3, and interpn (bug #44002).
Kai T. Ohlhus <k.ohlhus@gmail.com>
parents: 19794
diff changeset
258 %! 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
259 %! 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
260 %! 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
261 %! 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
262
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
263 %!shared z, zout, tol
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
264 %! z = zeros (3, 3, 3);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
265 %! zout = zeros (5, 5, 5);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
266 %! z(:,:,1) = [1 3 5; 3 5 7; 5 7 9];
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
267 %! z(:,:,2) = z(:,:,1) + 2;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
268 %! z(:,:,3) = z(:,:,2) + 2;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
269 %! for n = 1:5
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
270 %! 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
271 %! 2 3 4 5 6;
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
272 %! 3 4 5 6 7;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
273 %! 4 5 6 7 8;
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
274 %! 5 6 7 8 9] + (n-1);
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
275 %! end
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
276 %! tol = 10 * eps;
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
277 %!
13151
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
278 %!assert (interp3 (z), zout, tol)
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
279 %!assert (interp3 (z, "linear"), zout, tol)
e173fda06fca Fix bug #30295.
Ben Abbott <bpabbott@mac.com>
parents: 13150
diff changeset
280 %!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
281
20038
9fc020886ae9 maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents: 20000
diff changeset
282 ## Test input validation
18754
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
283 %!error interp3 ()
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
284 %!error interp3 ({1})
0ede4dbb37f1 Overhaul interp1, interp2, interp3 functions.
Rik <rik@octave.org>
parents: 18722
diff changeset
285 %!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
286 %!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
287 %!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
288 %!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
289 %!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
290 %!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
291 %!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
292 %!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
293 %!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
294