Mercurial > octave-nkf
view scripts/sparse/sprandn.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 | 0e1f5a750d00 |
children | 4197fc428c7d |
line wrap: on
line source
## Copyright (C) 2004-2013 Paul Kienzle ## ## 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/>. ## ## Original version by Paul Kienzle distributed as free software in the ## public domain. ## -*- texinfo -*- ## @deftypefn {Function File} {} sprandn (@var{m}, @var{n}, @var{d}) ## @deftypefnx {Function File} {} sprandn (@var{m}, @var{n}, @var{d}, @var{rc}) ## @deftypefnx {Function File} {} sprandn (@var{s}) ## Generate a sparse matrix with normally distributed random values. ## ## The size of the matrix is @var{m}x@var{n} with a density of values @var{d}. ## @var{d} must be between 0 and 1. Values will be normally distributed with a ## mean of 0 and a variance of 1. ## ## If called with a single matrix argument, a sparse matrix is generated with ## random values wherever the matrix @var{s} is nonzero. ## ## If called with a scalar fourth argument @var{rc}, a random sparse matrix ## with reciprocal condition number @var{rc} is generated. If @var{rc} is ## a vector, then it specifies the first singular values of the generated ## matrix (@code{length (@var{rc}) <= min (@var{m}, @var{n})}). ## ## @seealso{sprand, sprandsym, randn} ## @end deftypefn ## Author: Paul Kienzle <pkienzle@users.sf.net> function s = sprandn (m, n, d, rc) if (nargin == 1 ) s = __sprand_impl__ (m, @randn); elseif ( nargin == 3) s = __sprand_impl__ (m, n, d, "sprandn", @randn); elseif (nargin == 4) s = __sprand_impl__ (m, n, d, rc, "sprandn", @randn); else print_usage (); endif endfunction %% Test 3-input calling form %!test %! s = sprandn (4, 10, 0.1); %! assert (size (s), [4, 10]); %! assert (nnz (s) / numel (s), 0.1); %% Test 4-input calling form %!test %! d = rand (); %! s1 = sprandn (100, 100, d, 0.4); %! rc = [5, 4, 3, 2, 1, 0.1]; %! s2 = sprandn (100, 100, d, rc); %! s3 = sprandn (6, 4, d, rc); %! assert (svd (s2)'(1:length (rc)), rc, sqrt (eps)); %! assert (1/cond (s1), 0.4, sqrt (eps)); %! assert (nnz (s1) / (100*100), d, 0.02); %! assert (nnz (s2) / (100*100), d, 0.02); %! assert (svd (s3)', [5 4 3 2], sqrt (eps)); %% Test 1-input calling form %!test %! s = sprandn (sparse ([1 2 3], [3 2 3], [2 2 2])); %! [i, j] = find (s); %! assert (sort (i), [1 2 3]'); %! assert (sort (j), [2 3 3]'); %% Test input validation %!error sprandn () %!error sprandn (1, 2) %!error sprandn (1, 2, 3, 4) %!error <M must be an integer greater than 0> sprandn (ones (3), 3, 0.5) %!error <M must be an integer greater than 0> sprandn (3.5, 3, 0.5) %!error <M must be an integer greater than 0> sprandn (0, 3, 0.5) %!error <N must be an integer greater than 0> sprandn (3, ones (3), 0.5) %!error <N must be an integer greater than 0> sprandn (3, 3.5, 0.5) %!error <N must be an integer greater than 0> sprandn (3, 0, 0.5) %!error <D must be between 0 and 1> sprandn (3, 3, -1) %!error <D must be between 0 and 1> sprandn (3, 3, 2) %!error <RC must be a scalar or vector> sprandn (2, 2, 0.2, ones (3,3)) %!error <RC must be between 0 and 1> sprandn (2, 2, 0.2, -1) %!error <RC must be between 0 and 1> sprandn (2, 2, 0.2, 2) %% Test very large, very low density matrix doesn't fail %!test %! s = sprandn (1e6,1e6,1e-7);