Mercurial > hg > octave-nkf

diff scripts/statistics/distributions/trnd.m @ 14683:57569a35765c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Simplify class handling ('single' vs. 'double') for statistical distribution random functions. * betarnd.m, cauchy_rnd.m, chi2rnd.m, exprnd.m, frnd.m, gamrnd.m, geornd.m, lognrnd.m, nbinrnd.m, normrnd.m, poissrnd.m, trnd.m, unidrnd.m, unifrnd.m, wblrnd.m: Use class argument to random number function rather than explicitly or implicitly casting output to class single.
author Rik <octave@nomad.inbox5.com>
date Thu, 24 May 2012 09:50:28 -0700
parents f3d52523cde1
children d63878346099
line wrap: on
line diff
--- a/scripts/statistics/distributions/trnd.m
+++ b/scripts/statistics/distributions/trnd.m
@@ -77,7 +77,7 @@
 
   if (isscalar (n))
     if ((n > 0) && (n < Inf))
-      rnd = randn (sz) ./ sqrt (2*randg (n/2, sz) / n);
+      rnd = randn (sz, cls) ./ sqrt (2*randg (n/2, sz, cls) / n);
     else
       rnd = NaN (sz, cls);
     endif
@@ -85,7 +85,7 @@
     rnd = NaN (sz, cls);
 
     k = (n > 0) & (n < Inf);
-    rnd(k) = randn (sum (k(:)), 1) ./ sqrt (2*randg (n(k)/2) ./ n(k))(:);
+    rnd(k) = randn (sum (k(:)), 1, cls) ./ sqrt (2*randg (n(k)/2, cls) ./ n(k))(:);
   endif
 
 endfunction