Mercurial > octave-nkf
view scripts/linear-algebra/bandwidth.m @ 19631:db92e7e28e1f
strip trailing whitespace from most source files
* NEWS, doc/interpreter/contributors.in, doc/interpreter/func.txi,
doc/interpreter/genpropdoc.m, doc/interpreter/octave_logo.eps,
doc/interpreter/plot.txi, doc/interpreter/stmt.txi,
examples/data/Makefile.am, libinterp/corefcn/data.cc,
libinterp/corefcn/debug.cc, libinterp/corefcn/error.cc,
libinterp/corefcn/file-io.cc, libinterp/corefcn/gl-render.cc,
libinterp/corefcn/graphics.cc, libinterp/corefcn/graphics.in.h,
libinterp/corefcn/load-path.cc, libinterp/corefcn/pr-output.cc,
libinterp/corefcn/pt-jit.cc, libinterp/corefcn/strfind.cc,
libinterp/corefcn/toplev.cc, libinterp/corefcn/toplev.h,
libinterp/corefcn/urlwrite.cc, libinterp/corefcn/variables.cc,
libinterp/octave-value/ov-classdef.cc,
libinterp/octave-value/ov-classdef.h, libinterp/octave.cc,
libinterp/parse-tree/lex.h, libinterp/parse-tree/oct-parse.in.yy,
libinterp/parse-tree/pt-classdef.h, liboctave/system/file-ops.cc,
liboctave/system/oct-env.cc, m4/acinclude.m4,
scripts/deprecated/finite.m, scripts/deprecated/fmod.m,
scripts/deprecated/fnmatch.m, scripts/deprecated/luinc.m,
scripts/deprecated/octave_tmp_file_name.m, scripts/deprecated/syl.m,
scripts/deprecated/usage.m, scripts/general/inputParser.m,
scripts/general/interp1.m, scripts/general/interp2.m,
scripts/general/interp3.m, scripts/general/isequal.m,
scripts/general/private/__isequal__.m, scripts/geometry/voronoi.m,
scripts/image/image.m, scripts/image/imshow.m,
scripts/image/ind2rgb.m, scripts/linear-algebra/bandwidth.m,
scripts/linear-algebra/isbanded.m, scripts/miscellaneous/bzip2.m,
scripts/miscellaneous/cast.m, scripts/miscellaneous/copyfile.m,
scripts/miscellaneous/delete.m, scripts/miscellaneous/fullfile.m,
scripts/miscellaneous/getappdata.m, scripts/miscellaneous/gunzip.m,
scripts/miscellaneous/isappdata.m, scripts/miscellaneous/ls.m,
scripts/miscellaneous/mex.m, scripts/miscellaneous/movefile.m,
scripts/miscellaneous/orderfields.m, scripts/miscellaneous/recycle.m,
scripts/miscellaneous/rmappdata.m, scripts/miscellaneous/setfield.m,
scripts/miscellaneous/symvar.m, scripts/miscellaneous/tar.m,
scripts/miscellaneous/tmpnam.m, scripts/miscellaneous/unpack.m,
scripts/miscellaneous/ver.m, scripts/miscellaneous/what.m,
scripts/miscellaneous/xor.m, scripts/miscellaneous/zip.m,
scripts/optimization/fminbnd.m, scripts/optimization/sqp.m,
scripts/path/private/getsavepath.m, scripts/path/savepath.m,
scripts/pkg/pkg.m, scripts/pkg/private/installed_packages.m,
scripts/plot/draw/plotyy.m, scripts/plot/draw/polar.m,
scripts/plot/draw/private/__quiver__.m,
scripts/plot/draw/private/__scatter__.m,
scripts/plot/draw/private/__stem__.m, scripts/plot/draw/surface.m,
scripts/plot/draw/surfnorm.m, scripts/plot/util/copyobj.m,
scripts/plot/util/hgload.m, scripts/plot/util/hgsave.m,
scripts/plot/util/isprop.m, scripts/plot/util/linkprop.m,
scripts/plot/util/private/__go_draw_axes__.m, scripts/set/setdiff.m,
scripts/set/union.m, scripts/signal/periodogram.m,
scripts/sparse/eigs.m, scripts/sparse/ilu.m, scripts/sparse/qmr.m,
scripts/sparse/sprand.m, scripts/sparse/sprandn.m,
scripts/specfun/beta.m, scripts/specfun/ellipke.m,
scripts/specfun/isprime.m, scripts/statistics/base/lscov.m,
scripts/testfun/__run_test_suite__.m, scripts/testfun/test.m:
Strip trailing whitespace.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 20 Jan 2015 10:29:54 -0500 |
parents | 9addb5ad9426 |
children | 4197fc428c7d |
line wrap: on
line source
## Copyright (C) 2014 Massimiliano Fasi ## ## 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{bw} =} bandwidth (@var{A}, @var{type}) ## @deftypefnx {Function File} {[@var{lower}, @var{upper}] =} bandwidth (@var{A}) ## Compute the bandwidth of @var{A}. ## ## The @var{type} argument is the string @qcode{"lower"} for the lower ## bandwidth and @qcode{"upper"} for the upper bandwidth. ## If no @var{type} is specified return both the lower and upper bandwidth ## of @var{A}. ## ## The lower/upper bandwidth of a matrix is the number of ## subdiagonals/superdiagonals with nonzero entries. ## ## @seealso{isbanded, isdiag, istril, istriu} ## @end deftypefn ## Author: Massimiliano Fasi function [lower, upper] = bandwidth (A, type) if (! ((nargin == 1 && nargout == 2) || (nargin == 2 && nargout <= 1))) print_usage (); endif if (! isnumeric (A) && ! islogical (A) || ndims (A) != 2) error ("bandwidth: A must be a 2-D numeric or logical matrix"); elseif (nargin == 2 && ! (strcmp (type, "lower") || strcmp (type, "upper"))) error ('bandwidth: TYPE must be "lower" or "upper"'); endif if (nargin == 1) [i, j] = find (A); if (isempty (i)) lower = upper = 0; else lower = max (i - j); upper = max (j - i); endif else [i, j] = find (A); if (isempty (i)) lower = 0; elseif (strcmp (type, "lower")) lower = max (i - j); else lower = max (j - i); endif endif endfunction %!test %! [a,b] = bandwidth (speye (100)); %! assert ([a,b] == [0,0]); %! assert (bandwidth (speye (100), "upper"), 0); %! assert (bandwidth (speye (100), "lower"), 0); %!test %! A = [2 3 0 0 0; 1 2 3 0 0; 0 1 2 3 0; 0 0 1 2 3; 0 0 0 1 2]; %! [a,b] = bandwidth (A); %! assert ([a,b] == [1,1]); %! assert (bandwidth (A, "lower"), 1); %! assert (bandwidth (A, "upper"), 1); %!assert (bandwidth ([], "lower"), 0) %!assert (bandwidth ([], "upper"), 0) %!assert (bandwidth (zeros (3,3), "lower"), 0) %!assert (bandwidth (zeros (3,3), "upper"), 0) %!assert (bandwidth (ones (5,5), "lower"), 4) %!assert (bandwidth (ones (5,5), "upper"), 4) %!test %! [a,b] = bandwidth ([]); %! assert ([a,b] == [0,0]); %!test %! [a,b] = bandwidth (zeros (3,3)); %! assert ([a,b] == [0,0]); %!test %! [a,b] = bandwidth (ones (5,5)); %! assert ([a,b] == [4,4]); %% Test input validation %!error bandwidth () %!error bandwidth (1,2,3) %!error [a,b,c] = bandwidth (ones (2)) %!error [a,b] = bandwidth (ones (2), "upper") %!error <A must be a 2-D numeric or logical> bandwidth ("string", "lower") %!error <A must be a 2-D numeric or logical> bandwidth (ones (3,3,3), "lower") %!error <TYPE must be "lower" or "upper"> bandwidth (ones (2), "uper") %!error <TYPE must be "lower" or "upper"> bandwidth (ones (2), "uppper")