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

annotate scripts/general/interp1.m @ 17744:d63878346099

maint: Update copyright notices for release.

author	John W. Eaton <jwe@octave.org>
date	Wed, 23 Oct 2013 22:09:27 -0400
parents	4448cc742880
children	b06675ef40f2

rev	line source
17744 d63878346099 maint: Update copyright notices for release. John W. Eaton <jwe@octave.org> parents: 17307 diff changeset	1 ## Copyright (C) 2000-2013 Paul Kienzle
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	2 ## Copyright (C) 2009 VZLU Prague
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	3 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	4 ## This file is part of Octave.
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	5 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	6 ## Octave is free software; you can redistribute it and/or modify it
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	7 ## 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	8 ## 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	9 ## your option) any later version.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	10 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	11 ## Octave is distributed in the hope that it will be useful, but
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	12 ## WITHOUT ANY WARRANTY; without even the implied warranty of
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	14 ## General Public License for more details.
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	15 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	16 ## 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	17 ## 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	18 ## <http://www.gnu.org/licenses/>.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	19
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	20 ## -- texinfo --
10846 a4f482e66b65 Grammarcheck more of the documentation. Rik <octave@nomad.inbox5.com> parents: 10821 diff changeset	21 ## @deftypefn {Function File} {@var{yi} =} interp1 (@var{x}, @var{y}, @var{xi})
10820 c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	22 ## @deftypefnx {Function File} {@var{yi} =} interp1 (@var{y}, @var{xi})
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	23 ## @deftypefnx {Function File} {@var{yi} =} interp1 (@dots{}, @var{method})
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	24 ## @deftypefnx {Function File} {@var{yi} =} interp1 (@dots{}, @var{extrap})
14359 7277fe922e99 doc: Use Octave preference for double quote in docstrings in scripts/ Rik <octave@nomad.inbox5.com> parents: 14335 diff changeset	25 ## @deftypefnx {Function File} {@var{pp} =} interp1 (@dots{}, "pp")
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	26 ##
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	27 ## One-dimensional interpolation. Interpolates to determine the value of
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	28 ## @var{yi} at the points, @var{xi}. If not specified, @var{x} is taken
17087 bb3a1ec7cc4b doc: Clarify matrix and ND-array args to interpolation functions Mike Miller <mtmiller@ieee.org> parents: 16816 diff changeset	29 ## to be the indices of @var{y}. If @var{y} is a matrix or an N-dimensional
bb3a1ec7cc4b doc: Clarify matrix and ND-array args to interpolation functions Mike Miller <mtmiller@ieee.org> parents: 16816 diff changeset	30 ## array, the interpolation is performed on each column of @var{y}.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	31 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	32 ## Method is one of:
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	33 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	34 ## @table @asis
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	35 ## @item @qcode{"nearest"}
9070 e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi Rik <rdrider0-list@yahoo.com> parents: 9051 diff changeset	36 ## Return the nearest neighbor.
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10820 diff changeset	37 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	38 ## @item @qcode{"linear"}
9070 e9dc2ed2ec0f Cleanup documentation for poly.texi, interp.texi, geometry.texi Rik <rdrider0-list@yahoo.com> parents: 9051 diff changeset	39 ## Linear interpolation from nearest neighbors
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10820 diff changeset	40 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	41 ## @item @qcode{"pchip"}
11536 702dbd0c53f5 Add undocumented ppder, ppint, ppjumps functions to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	42 ## Piecewise cubic Hermite interpolating polynomial
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10820 diff changeset	43 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	44 ## @item @qcode{"cubic"}
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	45 ## Cubic interpolation (same as @code{pchip})
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 10820 diff changeset	46 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	47 ## @item @qcode{"spline"}
12175 2090995ca588 Correct en-dash,em-dash instances in docstrings. Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	48 ## Cubic spline interpolation---smooth first and second derivatives
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	49 ## throughout the curve
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	50 ## @end table
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	51 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	52 ## Appending '*' to the start of the above method forces @code{interp1}
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	53 ## to assume that @var{x} is uniformly spaced, and only @code{@var{x}(1)}
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	54 ## and @code{@var{x}(2)} are referenced. This is usually faster,
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	55 ## and is never slower. The default method is @qcode{"linear"}.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	56 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	57 ## If @var{extrap} is the string @qcode{"extrap"}, then extrapolate values
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	58 ## beyond the endpoints. If @var{extrap} is a number, replace values beyond
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	59 ## the endpoints with that number. If @var{extrap} is missing, assume NA.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	60 ##
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	61 ## If the string argument @qcode{"pp"} is specified, then @var{xi} should not be
11536 702dbd0c53f5 Add undocumented ppder, ppint, ppjumps functions to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	62 ## supplied and @code{interp1} returns the piecewise polynomial that
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	63 ## can later be used with @code{ppval} to evaluate the interpolation.
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	64 ## There is an equivalence, such that @code{ppval (interp1 (@var{x},
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	65 ## @var{y}, @var{method}, @qcode{"pp"}), @var{xi}) == interp1 (@var{x}, @var{y},
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	66 ## @var{xi}, @var{method}, @qcode{"extrap"})}.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	67 ##
10711 fbd7843974fa Periodic grammar check of documentation files to ensure common format. Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	68 ## Duplicate points in @var{x} specify a discontinuous interpolant. There
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	69 ## may be at most 2 consecutive points with the same value.
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	70 ## If @var{x} is increasing, the default discontinuous interpolant is
16816 12005245b645 doc: Periodic grammarcheck of documentation. Rik <rik@octave.org> parents: 15468 diff changeset	71 ## right-continuous. If @var{x} is decreasing, the default discontinuous
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	72 ## interpolant is left-continuous.
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	73 ## The continuity condition of the interpolant may be specified by using
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	74 ## the options, @qcode{"-left"} or @qcode{"-right"}, to select a left-continuous
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	75 ## or right-continuous interpolant, respectively.
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	76 ## Discontinuous interpolation is only allowed for @qcode{"nearest"} and
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	77 ## @qcode{"linear"} methods; in all other cases, the @var{x}-values must be
bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 17174 diff changeset	78 ## unique.
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	79 ##
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	80 ## An example of the use of @code{interp1} is
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	81 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	82 ## @example
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	83 ## @group
8507 cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	84 ## xf = [0:0.05:10];
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	85 ## yf = sin (2pixf/5);
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	86 ## xp = [0:10];
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	87 ## yp = sin (2pixp/5);
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	88 ## lin = interp1 (xp, yp, xf);
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	89 ## spl = interp1 (xp, yp, xf, "spline");
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	90 ## cub = interp1 (xp, yp, xf, "cubic");
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	91 ## near = interp1 (xp, yp, xf, "nearest");
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	92 ## plot (xf, yf, "r", xf, lin, "g", xf, spl, "b",
cadc73247d65 style fixes John W. Eaton <jwe@octave.org> parents: 7671 diff changeset	93 ## xf, cub, "c", xf, near, "m", xp, yp, "r*");
14327 4d917a6a858b doc: Use Octave coding conventions in @example blocks of docstrings. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	94 ## legend ("original", "linear", "spline", "cubic", "nearest");
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	95 ## @end group
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	96 ## @end example
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	97 ##
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	98 ## @seealso{interpft}
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	99 ## @end deftypefn
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	100
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	101 ## Author: Paul Kienzle
376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	102 ## Date: 2000-03-25
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	103 ## added 'nearest' as suggested by Kai Habel
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	104 ## 2000-07-17 Paul Kienzle
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	105 ## added '*' methods and matrix y
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	106 ## check for proper table lengths
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	107 ## 2002-01-23 Paul Kienzle
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	108 ## fixed extrapolation
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	109
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	110 function yi = interp1 (x, y, varargin)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	111
10820 c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	112 if (nargin < 2 \|\| nargin > 6)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	113 print_usage ();
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	114 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	115
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	116 method = "linear";
6742 ebf96cc00ee9 [project @ 2007-06-18 16:27:52 by jwe] jwe parents: 6721 diff changeset	117 extrap = NA;
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	118 xi = [];
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	119 ispp = false;
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	120 firstnumeric = true;
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	121 rightcontinuous = [];
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	122
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	123 if (nargin > 2)
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	124 for i = 1:length (varargin)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	125 arg = varargin{i};
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	126 if (ischar (arg))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	127 arg = tolower (arg);
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	128 if (strcmp ("extrap", arg))
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	129 extrap = "extrap";
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	130 elseif (strcmp ("pp", arg))
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	131 ispp = true;
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	132 elseif (any (strcmp ({"right", "-right"}, arg)))
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	133 rightcontinuous = true;
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	134 elseif (any (strcmp ({"left", "-left"}, arg)))
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	135 rightcontinuous = false;
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	136 else
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	137 method = arg;
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	138 endif
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	139 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	140 if (firstnumeric)
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	141 xi = arg;
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	142 firstnumeric = false;
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	143 else
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	144 extrap = arg;
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9929 diff changeset	145 endif
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	146 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	147 endfor
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	148 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	149
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	150 if (isempty (xi) && firstnumeric && ! ispp)
10820 c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	151 xi = y;
c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	152 y = x;
17088 c8586e6fc60d interp1.m: Fix interpolation on ND-arrays (bug #35860) Mike Miller <mtmiller@ieee.org> parents: 17087 diff changeset	153 if (isvector (y))
c8586e6fc60d interp1.m: Fix interpolation on ND-arrays (bug #35860) Mike Miller <mtmiller@ieee.org> parents: 17087 diff changeset	154 x = 1:numel (y);
c8586e6fc60d interp1.m: Fix interpolation on ND-arrays (bug #35860) Mike Miller <mtmiller@ieee.org> parents: 17087 diff changeset	155 else
c8586e6fc60d interp1.m: Fix interpolation on ND-arrays (bug #35860) Mike Miller <mtmiller@ieee.org> parents: 17087 diff changeset	156 x = 1:rows (y);
c8586e6fc60d interp1.m: Fix interpolation on ND-arrays (bug #35860) Mike Miller <mtmiller@ieee.org> parents: 17087 diff changeset	157 endif
10820 c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	158 endif
c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	159
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	160 ## reshape matrices for convenience
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	161 x = x(:);
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	162 nx = rows (x);
9769 9a1c4fe44af8 small interp1 simplification Jaroslav Hajek <highegg@gmail.com> parents: 9754 diff changeset	163 szx = size (xi);
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	164 if (isvector (y))
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	165 y = y(:);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	166 endif
13141 e81ddf9cacd5 maint: untabify and remove trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 12608 diff changeset	167
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	168 szy = size (y);
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	169 y = y(:,:);
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	170 [ny, nc] = size (y);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	171 xi = xi(:);
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	172
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	173 ## determine sizes
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	174 if (nx < 2 \|\| ny < 2)
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	175 error ("interp1: table too short");
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	176 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	177
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	178 ## check whether x is sorted; sort if not.
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	179 if (! issorted (x, "either"))
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	180 [x, p] = sort (x);
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	181 y = y(p,:);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	182 endif
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	183
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	184 if (isempty (rightcontinuous))
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	185 ## If not specified, set the continuity condition
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	186 if (x(end) < x(1))
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	187 rightcontinuous = false;
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	188 else
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	189 rightcontinuous = true;
17306 09543e9c8f40 Use explicit form of end (endif, endfor, etc.) in core m-files. Rik <rik@octave.org> parents: 17281 diff changeset	190 endif
15193 2d5e4d283688 Modified interp1.m file to check whether X has distinct values or not. Vivek Dogra <vivek.dogra.iitd@gmail.com> parents: 14868 diff changeset	191 endif
2d5e4d283688 Modified interp1.m file to check whether X has distinct values or not. Vivek Dogra <vivek.dogra.iitd@gmail.com> parents: 14868 diff changeset	192
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	193 if ((rightcontinuous && (x(end) < x(1)))
17307 4448cc742880 Use '!' not operator rather than '~' in core m-files. Rik <rik@octave.org> parents: 17306 diff changeset	194 \|\| (! rightcontinuous && (x(end) > x(1))))
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	195 ## Switch between left-continuous and right-continuous
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	196 x = flipud (x);
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	197 y = flipud (y);
17306 09543e9c8f40 Use explicit form of end (endif, endfor, etc.) in core m-files. Rik <rik@octave.org> parents: 17281 diff changeset	198 endif
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	199
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	200 starmethod = method(1) == "*";
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	201
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	202 if (starmethod)
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	203 dx = x(2) - x(1);
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	204 else
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	205 jumps = x(1:end-1) == x(2:end);
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	206 have_jumps = any (jumps);
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	207 if (have_jumps)
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	208 if (any (strcmp (method, {"nearest", "linear"})))
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	209 if (any (jumps(1:nx-2) & jumps(2:nx-1)))
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	210 error ("interp1: extra points in discontinuities");
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	211 endif
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	212 else
15468 6437fa7263dd use ' instead of ` in error messages, warnings and most comments John W. Eaton <jwe@octave.org> parents: 15218 diff changeset	213 error ("interp1: discontinuities not supported for method '%s'", method);
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	214 endif
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	215 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	216 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	217
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	218 ## Proceed with interpolating by all methods.
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	219 switch (method)
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	220 case "nearest"
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	221 pp = mkpp ([x(1); (x(1:nx-1)+x(2:nx))/2; x(nx)],
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	222 shiftdim (y, 1), szy(2:end));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	223 pp.orient = "first";
13141 e81ddf9cacd5 maint: untabify and remove trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 12608 diff changeset	224
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	225 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	226 yi = pp;
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	227 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	228 yi = ppval (pp, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	229 endif
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	230 case "*nearest"
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	231 pp = mkpp ([x(1), x(1)+[0.5:(nx-1)]*dx, x(nx)],
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	232 shiftdim (y, 1), szy(2:end));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	233 pp.orient = "first";
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	234 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	235 yi = pp;
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	236 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	237 yi = ppval (pp, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	238 endif
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	239 case "linear"
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	240 dy = diff (y);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	241 dx = diff (x);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	242 dx = repmat (dx, [1 size(dy)(2:end)]);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	243 coefs = [(dy./dx).'(:), y(1:nx-1, :).'(:)];
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	244 xx = x;
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	245
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	246 if (have_jumps)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	247 ## Omit zero-size intervals.
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	248 coefs(jumps, :) = [];
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	249 xx(jumps) = [];
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	250 endif
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	251
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	252 pp = mkpp (xx, coefs, szy(2:end));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	253 pp.orient = "first";
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	254
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	255 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	256 yi = pp;
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	257 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	258 yi = ppval (pp, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	259 endif
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	260
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	261 case "*linear"
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	262 dy = diff (y);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	263 coefs = [(dy/dx).'(:), y(1:nx-1, :).'(:)];
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	264 pp = mkpp (x, coefs, szy(2:end));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	265 pp.orient = "first";
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	266
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	267 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	268 yi = pp;
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	269 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	270 yi = ppval (pp, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	271 endif
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	272
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	273 case {"pchip", "pchip", "cubic", "cubic"}
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	274 if (nx == 2 \|\| starmethod)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	275 x = linspace (x(1), x(nx), ny);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	276 endif
13141 e81ddf9cacd5 maint: untabify and remove trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 12608 diff changeset	277
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	278 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	279 y = shiftdim (reshape (y, szy), 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	280 yi = pchip (x, y);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	281 yi.orient = "first";
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	282 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	283 y = shiftdim (y, 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	284 yi = pchip (x, y, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	285 if (! isvector (y))
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	286 yi = shiftdim (yi, 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	287 endif
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	288 endif
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	289 case {"spline", "*spline"}
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	290 if (nx == 2 \|\| starmethod)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	291 x = linspace (x(1), x(nx), ny);
14492 7ce925166af6 Fix inconsistent orientation of output matrix for interp1. Olaf Till <i7tiol@t-online.de> parents: 14363 diff changeset	292 endif
13141 e81ddf9cacd5 maint: untabify and remove trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 12608 diff changeset	293
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	294 if (ispp)
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	295 y = shiftdim (reshape (y, szy), 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	296 yi = spline (x, y);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	297 yi.orient = "first";
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	298 else
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	299 y = shiftdim (y, 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	300 yi = spline (x, y, reshape (xi, szx));
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	301 if (! isvector (y))
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	302 yi = shiftdim (yi, 1);
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	303 endif
14492 7ce925166af6 Fix inconsistent orientation of output matrix for interp1. Olaf Till <i7tiol@t-online.de> parents: 14363 diff changeset	304 endif
17174 c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	305 otherwise
c3c1ebfaa7dc maint: Use common indentation for switch statement. Rik <rik@octave.org> parents: 17088 diff changeset	306 error ("interp1: invalid method '%s'", method);
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	307 endswitch
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	308
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	309 if (! ispp)
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	310 if (! ischar (extrap))
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	311 ## determine which values are out of range and set them to extrap,
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	312 ## unless extrap == "extrap".
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	313 minx = min (x(1), x(nx));
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	314 maxx = max (x(1), x(nx));
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	315
4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	316 outliers = xi < minx \| ! (xi <= maxx); # this catches even NaNs
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	317 if (size_equal (outliers, yi))
59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	318 yi(outliers) = extrap;
59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	319 yi = reshape (yi, szx);
59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	320 elseif (!isvector (yi))
14492 7ce925166af6 Fix inconsistent orientation of output matrix for interp1. Olaf Till <i7tiol@t-online.de> parents: 14363 diff changeset	321 yi(outliers, :) = extrap;
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	322 else
59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	323 yi(outliers.') = extrap;
59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	324 endif
9754 4219e5cf773d improve interp1 and pchip Jaroslav Hajek <highegg@gmail.com> parents: 9070 diff changeset	325 endif
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	326 endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	327
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	328 endfunction
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	329
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	330
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	331 %!demo
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	332 %! clf;
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	333 %! xf = 0:0.05:10; yf = sin (2pixf/5);
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	334 %! xp = 0:10; yp = sin (2pixp/5);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	335 %! lin = interp1 (xp,yp,xf, "linear");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	336 %! spl = interp1 (xp,yp,xf, "spline");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	337 %! cub = interp1 (xp,yp,xf, "pchip");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	338 %! near= interp1 (xp,yp,xf, "nearest");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	339 %! plot (xf,yf,"r",xf,near,"g",xf,lin,"b",xf,cub,"c",xf,spl,"m",xp,yp,"r*");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	340 %! legend ("original", "nearest", "linear", "pchip", "spline");
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	341 %! %--------------------------------------------------------
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	342 %! % confirm that interpolated function matches the original
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	343
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	344 %!demo
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	345 %! clf;
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	346 %! xf = 0:0.05:10; yf = sin (2pixf/5);
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	347 %! xp = 0:10; yp = sin (2pixp/5);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	348 %! lin = interp1 (xp,yp,xf, "*linear");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	349 %! spl = interp1 (xp,yp,xf, "*spline");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	350 %! cub = interp1 (xp,yp,xf, "*cubic");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	351 %! near= interp1 (xp,yp,xf, "*nearest");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	352 %! plot (xf,yf,"r",xf,near,"g",xf,lin,"b",xf,cub,"c",xf,spl,"m",xp,yp,"r*");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	353 %! legend ("original", "nearest", "linear", "cubic", "*spline");
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	354 %! %--------------------------------------------------------
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	355 %! % confirm that interpolated function matches the original
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	356
6721 01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	357 %!demo
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	358 %! clf;
6721 01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	359 %! t = 0 : 0.3 : pi; dt = t(2)-t(1);
01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	360 %! n = length (t); k = 100; dti = dt*n/k;
01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	361 %! ti = t(1) + [0 : k-1]*dti;
01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	362 %! y = sin (4t + 0.3) . cos (3*t - 0.1);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	363 %! ddyc = diff (diff (interp1 (t,y,ti, "cubic")) ./dti)./dti;
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	364 %! ddys = diff (diff (interp1 (t,y,ti, "spline"))./dti)./dti;
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	365 %! ddyp = diff (diff (interp1 (t,y,ti, "pchip")) ./dti)./dti;
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	366 %! plot (ti(2:end-1),ddyc,'g+', ti(2:end-1),ddys,'b*', ti(2:end-1),ddyp,'c^');
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	367 %! legend ("cubic", "spline", "pchip");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	368 %! title ("Second derivative of interpolated 'sin (4t + 0.3) . cos (3*t - 0.1)'");
6721 01036667884a [project @ 2007-06-14 06:56:41 by dbateman] dbateman parents: 6702 diff changeset	369
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	370 %!demo
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	371 %! clf;
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	372 %! xf = 0:0.05:10; yf = sin (2pixf/5) - (xf >= 5);
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	373 %! xp = [0:.5:4.5,4.99,5:.5:10]; yp = sin (2pixp/5) - (xp >= 5);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	374 %! lin = interp1 (xp,yp,xf, "linear");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	375 %! near= interp1 (xp,yp,xf, "nearest");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	376 %! plot (xf,yf,"r", xf,near,"g", xf,lin,"b", xp,yp,"r*");
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	377 %! legend ("original", "nearest", "linear");
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	378 %! %--------------------------------------------------------
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	379 %! % confirm that interpolated function matches the original
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	380
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	381 %!demo
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	382 %! clf;
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	383 %! x = 0:0.5:3;
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	384 %! x1 = [3 2 2 1];
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	385 %! x2 = [1 2 2 3];
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	386 %! y1 = [1 1 0 0];
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	387 %! y2 = [0 0 1 1];
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	388 %! h = plot (x, interp1 (x1, y1, x), 'b', x1, y1, 'sb');
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	389 %! hold on
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	390 %! g = plot (x, interp1 (x2, y2, x), 'r', x2, y2, '*r');
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	391 %! xlim ([1 3])
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	392 %! legend ([h(1), g(1)], {'left-continous', 'right-continuous'}, ...
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	393 %! 'location', 'east')
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	394 %! legend boxoff
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	395 %! %--------------------------------------------------------
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	396 %! % red curve is left-continuos and blue is right-continuous at x = 2
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	397
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	398 ##FIXME: add test for n-d arguments here
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	399
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	400 ## For each type of interpolated test, confirm that the interpolated
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	401 ## value at the knots match the values at the knots. Points away
14359 7277fe922e99 doc: Use Octave preference for double quote in docstrings in scripts/ Rik <octave@nomad.inbox5.com> parents: 14335 diff changeset	402 ## from the knots are requested, but only "nearest" and "linear"
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	403 ## confirm they are the correct values.
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	404
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	405 %!shared xp, yp, xi, style
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	406 %! xp = 0:2:10;
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	407 %! yp = sin (2pixp/5);
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	408 %! xi = [-1, 0, 2.2, 4, 6.6, 10, 11];
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	409
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	410 ## The following BLOCK/ENDBLOCK section is repeated for each style
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	411 ## nearest, linear, cubic, spline, pchip
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	412 ## The test for ppval of cubic has looser tolerance, but otherwise
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	413 ## the tests are identical.
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	414 ## Note that the block checks style and *style; if you add more tests
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	415 ## be sure to add them to both sections of each block. One test,
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	416 ## style vs. *style, occurs only in the first section.
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	417 ## There is an ENDBLOCKTEST after the final block
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	418
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	419 %!test style = "nearest";
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	420 ## BLOCK
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	421 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	422 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	423 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	424 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	425 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	426 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	427 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	428 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	429 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	430 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	431 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	432 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	433 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	434 %!error interp1 (1,1,1, style)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	435 %!assert (interp1 (xp,[yp',yp'],xi,style),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	436 %! interp1 (xp,[yp',yp'],xi,["",style]),100eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	437 %!test style = ["*",style];
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	438 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	439 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	440 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	441 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	442 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	443 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	444 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	445 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	446 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	447 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	448 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	449 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	450 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	451 %!error interp1 (1,1,1, style)
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	452 ## ENDBLOCK
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	453
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	454 %!test style = "linear";
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	455 ## BLOCK
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	456 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	457 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	458 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	459 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	460 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	461 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	462 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	463 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	464 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	465 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	466 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	467 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	468 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	469 %!error interp1 (1,1,1, style)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	470 %!assert (interp1 (xp,[yp',yp'],xi,style),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	471 %! interp1 (xp,[yp',yp'],xi,["",style]),100eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	472 %!test style = ['*',style];
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	473 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	474 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	475 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	476 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	477 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	478 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	479 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	480 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	481 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	482 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	483 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	484 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	485 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
15218 94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	486 %!assert (interp1 ([1 2 2 3], [1 2 3 4], 2), 3);
94d512d712e3 Modified interp1.m file to check whether X has distinct values or not. Ben Abbott <bpabbott@mac.com> parents: 15193 diff changeset	487 %!assert (interp1 ([3 2 2 1], [4 3 2 1], 2), 2);
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	488 %!error interp1 (1,1,1, style)
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	489 ## ENDBLOCK
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	490
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	491 %!test style = "cubic";
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	492 ## BLOCK
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	493 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	494 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	495 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	496 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	497 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	498 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	499 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	500 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	501 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	502 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	503 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	504 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	505 %! interp1 (xp,yp,xi,style,"extrap"),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	506 %!error interp1 (1,1,1, style)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	507 %!assert (interp1 (xp,[yp',yp'],xi,style),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	508 %! interp1 (xp,[yp',yp'],xi,["",style]),100eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	509 %!test style = ["*",style];
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	510 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	511 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	512 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	513 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	514 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	515 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	516 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	517 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	518 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	519 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	520 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	521 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	522 %! interp1 (xp,yp,xi,style,"extrap"),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	523 %!error interp1 (1,1,1, style)
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	524 ## ENDBLOCK
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	525
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	526 %!test style = "pchip";
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	527 ## BLOCK
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	528 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	529 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	530 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	531 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	532 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	533 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	534 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	535 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	536 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	537 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	538 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	539 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	540 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	541 %!error interp1 (1,1,1, style)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	542 %!assert (interp1 (xp,[yp',yp'],xi,style),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	543 %! interp1 (xp,[yp',yp'],xi,["",style]),100eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	544 %!test style = ["*",style];
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	545 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	546 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	547 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	548 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	549 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	550 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	551 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	552 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	553 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	554 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	555 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	556 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	557 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	558 %!error interp1 (1,1,1, style)
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	559 ## ENDBLOCK
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	560
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	561 %!test style = "spline";
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	562 ## BLOCK
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	563 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	564 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	565 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	566 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	567 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	568 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	569 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	570 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	571 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	572 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	573 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	574 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	575 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	576 %!error interp1 (1,1,1, style)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	577 %!assert (interp1 (xp,[yp',yp'],xi,style),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	578 %! interp1 (xp,[yp',yp'],xi,["",style]),100eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	579 %!test style = ["*",style];
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	580 %!assert (interp1 (xp, yp, [min(xp)-1, max(xp)+1],style), [NA, NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	581 %!assert (interp1 (xp,yp,xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	582 %!assert (interp1 (xp,yp,xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	583 %!assert (interp1 (xp',yp',xp',style), yp', 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	584 %!assert (interp1 (xp',yp',xp,style), yp, 100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	585 %!assert (isempty (interp1 (xp',yp',[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	586 %!assert (isempty (interp1 (xp,yp,[],style)))
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	587 %!assert (interp1 (xp,[yp',yp'],xi(:),style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	588 %! [interp1(xp,yp,xi(:),style),interp1(xp,yp,xi(:),style)])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	589 %!assert (interp1 (xp,yp,xi,style),...
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	590 %! interp1 (fliplr (xp),fliplr (yp),xi,style),100*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	591 %!assert (ppval (interp1 (xp,yp,style,"pp"),xi),
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	592 %! interp1 (xp,yp,xi,style,"extrap"),10*eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	593 %!error interp1 (1,1,1, style)
6374 419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	594 ## ENDBLOCK
419017274c1e [project @ 2007-03-01 15:57:50 by jwe] jwe parents: 6366 diff changeset	595 ## ENDBLOCKTEST
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	596
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	597 %!# test linear extrapolation
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	598 %!assert (interp1 ([1:5],[3:2:11],[0,6],"linear","extrap"), [1, 13], eps)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	599 %!assert (interp1 (xp, yp, [-1, max(xp)+1],"linear",5), [5, 5])
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	600
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	601 %!assert (interp1 (1:2,1:2,1.4,"nearest"), 1)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	602 %!assert (interp1 (1:2,1:2,1.4,"linear"), 1.4)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	603 %!assert (interp1 (1:4,1:4,1.4,"cubic"), 1.4)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	604 %!assert (interp1 (1:2,1:2,1.1, "spline"), 1.1)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	605 %!assert (interp1 (1:3,1:3,1.4,"spline"), 1.4)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	606
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	607 %!assert (interp1 (1:2:4,1:2:4,1.4,"*nearest"), 1)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	608 %!assert (interp1 (1:2:4,1:2:4,[0,1,1.4,3,4],"*linear"), [NA,1,1.4,3,NA])
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	609 %!assert (interp1 (1:2:8,1:2:8,1.4,"*cubic"), 1.4)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	610 %!assert (interp1 (1:2,1:2,1.3, "*spline"), 1.3)
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	611 %!assert (interp1 (1:2:6,1:2:6,1.4,"*spline"), 1.4)
7671 4fbaba9abec1 implement compiled binary lookup Jaroslav Hajek <highegg@gmail.com> parents: 7017 diff changeset	612
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	613 %!assert (interp1 ([3,2,1],[3,2,2],2.5), 2.5)
9929 45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	614
45c08d7c2c79 allow discontinuous interpolant in interp1 Jaroslav Hajek <highegg@gmail.com> parents: 9769 diff changeset	615 %!assert (interp1 ([1,2,2,3,4],[0,1,4,2,1],[-1,1.5,2,2.5,3.5], "linear", "extrap"), [-2,0.5,4,3,1.5])
12608 59e2460acae1 make piecewise polynomial (pp) functions more compatible Kai Habel <kai.habel@gmx.de> parents: 12459 diff changeset	616 %!assert (interp1 ([4,4,3,2,0],[0,1,4,2,1],[1.5,4,4.5], "linear"), [1.75,1,NA])
10820 c44c786f87ba interp1.m: When absent set X equal to the inices of Y. Ben Abbott <bpabbott@mac.com> parents: 10793 diff changeset	617 %!assert (interp1 (0:4, 2.5), 1.5)
14363 f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	618
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	619 %!error interp1 ()
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	620 %!error interp1 (1,1,1, "linear")
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	621 %!error interp1 (1,1,1, "*nearest")
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	622 %!error interp1 (1,1,1, "*linear")
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	623 %!error interp1 (1:2,1:2,1, "bogus")
f3d52523cde1 Use Octave coding conventions in all m-file %!test blocks Rik <octave@nomad.inbox5.com> parents: 14359 diff changeset	624

Mercurial > octave-nkf

annotate scripts/general/interp1.m @ 17744:d63878346099