Mercurial > hg > octave-nkf
view scripts/image/summer.m @ 20771:87b557ee8e5d
clean up and vectorize code for dense output in ode45
* scripts/ode/private/ode_rk_interpolate.m: new file
* scripts/ode/private/ode_rk_interpolate.m(hermite_quartic_interpolation):
move to internal function, use vectorization and broadcasting.
* scripts/ode/private/hermite_quartic_interpolation.m: remove file
* scripts/ode/module.mk: list added and removed files
* scripts/ode/private/integrate_adaptive.m: use new interpolation code.
author | Carlo de Falco <carlo.defalco@polimi.it> |
---|---|
date | Tue, 06 Oct 2015 19:28:59 +0200 |
parents | 7503499a252b |
children |
line wrap: on
line source
## Copyright (C) 1999-2015 Kai Habel ## ## This file is part of Octave. ## ## Octave is free software; you can redistribute it and/or modify it ## under the terms of the GNU General Public License as published by ## the Free Software Foundation; either version 3 of the License, or (at ## your option) any later version. ## ## Octave is distributed in the hope that it will be useful, but ## WITHOUT ANY WARRANTY; without even the implied warranty of ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU ## General Public License for more details. ## ## You should have received a copy of the GNU General Public License ## along with Octave; see the file COPYING. If not, see ## <http://www.gnu.org/licenses/>. ## -*- texinfo -*- ## @deftypefn {Function File} {@var{map} =} summer () ## @deftypefnx {Function File} {@var{map} =} summer (@var{n}) ## Create color colormap. This colormap varies from green to yellow. ## ## The argument @var{n} must be a scalar. ## If unspecified, the length of the current colormap, or 64, is used. ## @seealso{colormap} ## @end deftypefn ## Author: Kai Habel <kai.habel@gmx.de> ## Date: 06/03/2000 ## PKG_ADD: colormap ("register", "summer"); ## PKG_DEL: colormap ("unregister", "summer"); function map = summer (n = rows (colormap ())) if (nargin > 1) print_usage (); elseif (! isscalar (n)) error ("summer: N must be a scalar"); endif n = double (n); if (n == 1) map = [0, 0.5, 0.4]; elseif (n > 1) r = [0:(n-1)]' / (n - 1); g = 0.5 + r / 2; b = 0.4 * ones (n, 1); map = [r, g, b]; else map = zeros (0, 3); endif endfunction %!demo %! ## Show the 'summer' colormap as an image %! image (1:64, linspace (0, 1, 64), repmat ((1:64)', 1, 64)); %! axis ([1, 64, 0, 1], "ticy", "xy"); %! colormap (summer (64));