octave: liboctave/numeric/oct-rand.cc annotate

annotate liboctave/numeric/oct-rand.cc @ 21136:7cac4e7458f2

maint: clean up code around calls to current_liboctave_error_handler. Remove statements after call to handler that are no longer reachable. Place input validation first and immediately call handler if necessary. Change if/error_handler/else to if/error_handler and re-indent code. * Array-util.cc, Array.cc, CColVector.cc, CDiagMatrix.cc, CMatrix.cc, CNDArray.cc, CRowVector.cc, CSparse.cc, DiagArray2.cc, MArray.cc, PermMatrix.cc, Sparse.cc, Sparse.h, chMatrix.cc, chNDArray.cc, dColVector.cc, dDiagMatrix.cc, dMatrix.cc, dNDArray.cc, dRowVector.cc, dSparse.cc, fCColVector.cc, fCDiagMatrix.cc, fCMatrix.cc, fCNDArray.cc, fCRowVector.cc, fColVector.cc, fDiagMatrix.cc, fMatrix.cc, fNDArray.cc, fRowVector.cc, idx-vector.cc, CmplxAEPBAL.cc, CmplxCHOL.cc, CmplxGEPBAL.cc, CmplxHESS.cc, CmplxLU.cc, CmplxQR.cc, CmplxSCHUR.cc, CmplxSVD.cc, DASPK.cc, EIG.cc, LSODE.cc, Quad.cc, SparseCmplxCHOL.cc, SparseCmplxLU.cc, SparseCmplxQR.cc, SparseQR.cc, SparsedbleCHOL.cc, SparsedbleLU.cc, base-lu.cc, bsxfun-defs.cc, dbleAEPBAL.cc, dbleCHOL.cc, dbleGEPBAL.cc, dbleHESS.cc, dbleLU.cc, dbleQR.cc, dbleSCHUR.cc, dbleSVD.cc, eigs-base.cc, fCmplxAEPBAL.cc, fCmplxCHOL.cc, fCmplxLU.cc, fCmplxQR.cc, fCmplxSCHUR.cc, fEIG.cc, floatAEPBAL.cc, floatCHOL.cc, floatGEPBAL.cc, floatHESS.cc, floatLU.cc, floatQR.cc, floatSCHUR.cc, floatSVD.cc, lo-specfun.cc, oct-fftw.cc, oct-rand.cc, oct-spparms.cc, sparse-base-chol.cc, sparse-dmsolve.cc, file-ops.cc, lo-sysdep.cc, mach-info.cc, oct-env.cc, oct-syscalls.cc, cmd-edit.cc, cmd-hist.cc, data-conv.cc, lo-ieee.cc, lo-regexp.cc, oct-base64.cc, oct-shlib.cc, pathsearch.cc, singleton-cleanup.cc, sparse-util.cc, unwind-prot.cc: Remove statements after call to handler that are no longer reachable. Place input validation first and immediately call handler if necessary. Change if/error_handler/else to if/error_handler and re-indent code.

author	Rik <rik@octave.org>
date	Sat, 23 Jan 2016 13:52:03 -0800
parents	a9574e3c6e9e
children	f7121e111991

rev	line source
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	1 /*
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	2
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 18573 diff changeset	3 Copyright (C) 2003-2015 John W. Eaton
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	4
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	5 This file is part of Octave.
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	6
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	8 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6435 diff changeset	9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6435 diff changeset	10 option) any later version.
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	11
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	15 for more details.
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	16
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	17 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: 6435 diff changeset	18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6435 diff changeset	19 <http://www.gnu.org/licenses/>.
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	20
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	21 */
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	22
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	23 #ifdef HAVE_CONFIG_H
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	24 #include <config.h>
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	25 #endif
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	26
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	27 #include <map>
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	28 #include <vector>
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	29
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	30 #include <stdint.h>
6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	31
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	32 #include "data-conv.h"
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	33 #include "f77-fcn.h"
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	34 #include "lo-error.h"
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	35 #include "lo-ieee.h"
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	36 #include "lo-mappers.h"
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	37 #include "mach-info.h"
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	38 #include "oct-locbuf.h"
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	39 #include "oct-rand.h"
4415 ed8c4aaa8648 [project @ 2003-05-16 21:20:33 by jwe] jwe parents: 4308 diff changeset	40 #include "oct-time.h"
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	41 #include "randgamma.h"
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	42 #include "randmtzig.h"
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	43 #include "randpoisson.h"
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	44 #include "singleton-cleanup.h"
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	45
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	46 extern "C"
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	47 {
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	48 F77_RET_T
6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	49 F77_FUNC (dgennor, DGENNOR) (const double&, const double&, double&);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	50
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	51 F77_RET_T
6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	52 F77_FUNC (dgenunf, DGENUNF) (const double&, const double&, double&);
6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	53
6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	54 F77_RET_T
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	55 F77_FUNC (dgenexp, DGENEXP) (const double&, double&);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	56
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	57 F77_RET_T
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	58 F77_FUNC (dignpoi, DIGNPOI) (const double&, double&);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	59
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	60 F77_RET_T
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	61 F77_FUNC (dgengam, DGENGAM) (const double&, const double&, double&);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	62
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	63 F77_RET_T
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	64 F77_FUNC (setall, SETALL) (const int32_t&, const int32_t&);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	65
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	66 F77_RET_T
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	67 F77_FUNC (getsd, GETSD) (int32_t&, int32_t&);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	68
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	69 F77_RET_T
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	70 F77_FUNC (setsd, SETSD) (const int32_t&, const int32_t&);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	71
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4543 diff changeset	72 F77_RET_T
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	73 F77_FUNC (setcgn, SETCGN) (const int32_t&);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	74 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	75
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	76 octave_rand *octave_rand::instance = 0;
6326 129ef2ae319f [project @ 2007-02-18 14:01:27 by dbateman] dbateman parents: 5828 diff changeset	77
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	78 octave_rand::octave_rand (void)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	79 : current_distribution (uniform_dist), use_old_generators (false),
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	80 rand_states ()
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	81 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	82 initialize_ranlib_generators ();
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	83
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	84 initialize_mersenne_twister ();
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	85 }
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	86
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	87 bool
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	88 octave_rand::instance_ok (void)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	89 {
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	90 bool retval = true;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	91
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	92 if (! instance)
13983 7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	93 {
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	94 instance = new octave_rand ();
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	95
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	96 if (instance)
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	97 singleton_cleanup_list::add (cleanup_instance);
7dd7cccf0757 clean up memory allocated for singletons before exit John W. Eaton <jwe@octave.org> parents: 11586 diff changeset	98 }
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	99
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	100 if (! instance)
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 20232 diff changeset	101 (*current_liboctave_error_handler)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 20232 diff changeset	102 ("unable to create octave_rand object!");
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	103
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	104 return retval;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	105 }
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	106
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	107 double
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	108 octave_rand::do_seed (void)
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	109 {
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	110 union d2i { double d; int32_t i[2]; };
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	111 union d2i u;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	112
6435 71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	113 oct_mach_info::float_format ff = oct_mach_info::native_float_format ();
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	114
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	115 switch (ff)
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	116 {
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	117 case oct_mach_info::flt_fmt_ieee_big_endian:
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	118 F77_FUNC (getsd, GETSD) (u.i[1], u.i[0]);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	119 break;
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	120 default:
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	121 F77_FUNC (getsd, GETSD) (u.i[0], u.i[1]);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	122 break;
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	123 }
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	124
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	125 return u.d;
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	126 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	127
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	128 static int32_t
6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	129 force_to_fit_range (int32_t i, int32_t lo, int32_t hi)
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	130 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	131 assert (hi > lo && lo >= 0 && hi > lo);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	132
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	133 i = i > 0 ? i : -i;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	134
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	135 if (i < lo)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	136 i = lo;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	137 else if (i > hi)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	138 i = i % hi;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	139
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	140 return i;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	141 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	142
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	143 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	144 octave_rand::do_seed (double s)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	145 {
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	146 use_old_generators = true;
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	147
6435 71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	148 int i0, i1;
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	149 union d2i { double d; int32_t i[2]; };
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	150 union d2i u;
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	151 u.d = s;
6435 71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	152
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	153 oct_mach_info::float_format ff = oct_mach_info::native_float_format ();
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	154
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	155 switch (ff)
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	156 {
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	157 case oct_mach_info::flt_fmt_ieee_big_endian:
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	158 i1 = force_to_fit_range (u.i[0], 1, 2147483563);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	159 i0 = force_to_fit_range (u.i[1], 1, 2147483399);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	160 break;
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	161 default:
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	162 i0 = force_to_fit_range (u.i[0], 1, 2147483563);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	163 i1 = force_to_fit_range (u.i[1], 1, 2147483399);
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	164 break;
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	165 }
71888a0ab040 [project @ 2007-03-23 14:17:02 by dbateman] dbateman parents: 6326 diff changeset	166
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	167 F77_FUNC (setsd, SETSD) (i0, i1);
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	168 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	169
10709 92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	170 void
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	171 octave_rand::do_reset (void)
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	172 {
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	173 use_old_generators = true;
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	174 initialize_ranlib_generators ();
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	175 }
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	176
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	177 ColumnVector
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	178 octave_rand::do_state (const std::string& d)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	179 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	180 return rand_states[d.empty () ? current_distribution : get_dist_id (d)];
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	181 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	182
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	183 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	184 octave_rand::do_state (const ColumnVector& s, const std::string& d)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	185 {
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	186 use_old_generators = false;
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	187
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	188 int old_dist = current_distribution;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	189
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	190 int new_dist = d.empty () ? current_distribution : get_dist_id (d);
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	191
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	192 ColumnVector saved_state;
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	193
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	194 if (old_dist != new_dist)
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	195 saved_state = get_internal_state ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	196
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	197 set_internal_state (s);
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	198
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	199 rand_states[new_dist] = get_internal_state ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	200
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	201 if (old_dist != new_dist)
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	202 rand_states[old_dist] = saved_state;
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	203 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	204
10709 92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	205 void
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	206 octave_rand::do_reset (const std::string& d)
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	207 {
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	208 use_old_generators = false;
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	209
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	210 int old_dist = current_distribution;
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	211
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	212 int new_dist = d.empty () ? current_distribution : get_dist_id (d);
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	213
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	214 ColumnVector saved_state;
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	215
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	216 if (old_dist != new_dist)
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	217 saved_state = get_internal_state ();
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	218
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	219 oct_init_by_entropy ();
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	220 rand_states[new_dist] = get_internal_state ();
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	221
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	222 if (old_dist != new_dist)
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	223 rand_states[old_dist] = saved_state;
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	224 }
92a85ed5b86e Don't special case color_property type when emitting factory default (bug #30118) David Bateman <dbateman@free.fr> parents: 10350 diff changeset	225
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	226 std::string
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	227 octave_rand::do_distribution (void)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	228 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	229 std::string retval;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	230
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	231 switch (current_distribution)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	232 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	233 case uniform_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	234 retval = "uniform";
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	235 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	236
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	237 case normal_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	238 retval = "normal";
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	239 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	240
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	241 case expon_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	242 retval = "exponential";
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	243 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	244
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	245 case poisson_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	246 retval = "poisson";
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	247 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	248
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	249 case gamma_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	250 retval = "gamma";
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	251 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	252
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	253 default:
7536 4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	254 (*current_liboctave_error_handler)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	255 ("rand: invalid distribution ID = %d", current_distribution);
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	256 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	257 }
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	258
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	259 return retval;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	260 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	261
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	262 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	263 octave_rand::do_distribution (const std::string& d)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	264 {
7536 4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	265 int id = get_dist_id (d);
4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	266
4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	267 switch (id)
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	268 {
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	269 case uniform_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	270 octave_rand::uniform_distribution ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	271 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	272
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	273 case normal_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	274 octave_rand::normal_distribution ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	275 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	276
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	277 case expon_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	278 octave_rand::exponential_distribution ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	279 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	280
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	281 case poisson_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	282 octave_rand::poisson_distribution ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	283 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	284
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	285 case gamma_dist:
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	286 octave_rand::gamma_distribution ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	287 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	288
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	289 default:
7536 4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	290 (*current_liboctave_error_handler)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	291 ("rand: invalid distribution ID = %d", id);
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	292 break;
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	293 }
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	294 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	295
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	296 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	297 octave_rand::do_uniform_distribution (void)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	298 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	299 switch_to_generator (uniform_dist);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	300
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	301 F77_FUNC (setcgn, SETCGN) (uniform_dist);
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	302 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	303
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	304 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	305 octave_rand::do_normal_distribution (void)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	306 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	307 switch_to_generator (normal_dist);
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	308
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	309 F77_FUNC (setcgn, SETCGN) (normal_dist);
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	310 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	311
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	312 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	313 octave_rand::do_exponential_distribution (void)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	314 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	315 switch_to_generator (expon_dist);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	316
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	317 F77_FUNC (setcgn, SETCGN) (expon_dist);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	318 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	319
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	320 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	321 octave_rand::do_poisson_distribution (void)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	322 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	323 switch_to_generator (poisson_dist);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	324
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	325 F77_FUNC (setcgn, SETCGN) (poisson_dist);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	326 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	327
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	328 void
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	329 octave_rand::do_gamma_distribution (void)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	330 {
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	331 switch_to_generator (gamma_dist);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	332
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	333 F77_FUNC (setcgn, SETCGN) (gamma_dist);
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	334 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	335
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	336
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	337 double
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	338 octave_rand::do_scalar (double a)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	339 {
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	340 double retval = 0.0;
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	341
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	342 if (use_old_generators)
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	343 {
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	344 switch (current_distribution)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	345 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	346 case uniform_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	347 F77_FUNC (dgenunf, DGENUNF) (0.0, 1.0, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	348 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	349
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	350 case normal_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	351 F77_FUNC (dgennor, DGENNOR) (0.0, 1.0, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	352 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	353
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	354 case expon_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	355 F77_FUNC (dgenexp, DGENEXP) (1.0, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	356 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	357
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	358 case poisson_dist:
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	359 if (a < 0.0 \|\| ! xfinite (a))
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	360 retval = octave_NaN;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	361 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	362 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	363 // workaround bug in ignpoi, by calling with different Mu
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	364 F77_FUNC (dignpoi, DIGNPOI) (a + 1, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	365 F77_FUNC (dignpoi, DIGNPOI) (a, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	366 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	367 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	368
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	369 case gamma_dist:
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	370 if (a <= 0.0 \|\| ! xfinite (a))
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	371 retval = octave_NaN;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	372 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	373 F77_FUNC (dgengam, DGENGAM) (1.0, a, retval);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	374 break;
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	375
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	376 default:
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	377 (*current_liboctave_error_handler)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	378 ("rand: invalid distribution ID = %d", current_distribution);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	379 break;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	380 }
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	381 }
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	382 else
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	383 {
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	384 switch (current_distribution)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	385 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	386 case uniform_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	387 retval = oct_randu ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	388 break;
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	389
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	390 case normal_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	391 retval = oct_randn ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	392 break;
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	393
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	394 case expon_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	395 retval = oct_rande ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	396 break;
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	397
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	398 case poisson_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	399 retval = oct_randp (a);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	400 break;
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	401
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	402 case gamma_dist:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	403 retval = oct_randg (a);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	404 break;
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	405
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	406 default:
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	407 (*current_liboctave_error_handler)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	408 ("rand: invalid distribution ID = %d", current_distribution);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	409 break;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	410 }
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	411
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	412 save_state ();
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	413 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	414
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	415 return retval;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	416 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4415 diff changeset	417
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	418 float
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	419 octave_rand::do_float_scalar (float a)
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	420 {
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	421 float retval = 0.0;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	422
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	423 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	424 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	425 double da = a;
15866 c1c6502fe52b oct-rand.cc: Silence compiler warning about possible uninitialized variable. Rik <rik@octave.org> parents: 15467 diff changeset	426 double dretval = 0.0;
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	427 switch (current_distribution)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	428 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	429 case uniform_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	430 F77_FUNC (dgenunf, DGENUNF) (0.0, 1.0, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	431 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	432
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	433 case normal_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	434 F77_FUNC (dgennor, DGENNOR) (0.0, 1.0, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	435 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	436
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	437 case expon_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	438 F77_FUNC (dgenexp, DGENEXP) (1.0, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	439 break;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	440
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	441 case poisson_dist:
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	442 if (da < 0.0 \|\| ! xfinite (a))
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	443 dretval = octave_NaN;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	444 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	445 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	446 // workaround bug in ignpoi, by calling with different Mu
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	447 F77_FUNC (dignpoi, DIGNPOI) (da + 1, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	448 F77_FUNC (dignpoi, DIGNPOI) (da, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	449 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	450 break;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	451
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	452 case gamma_dist:
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	453 if (da <= 0.0 \|\| ! xfinite (a))
18573 53ba52143af7 Fix dead assignment in oct-rand.cc do_float_scalar (bug #41827). Branden Archer <b.m.archer4@gmail.com> parents: 17769 diff changeset	454 dretval = octave_NaN;
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	455 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	456 F77_FUNC (dgengam, DGENGAM) (1.0, da, dretval);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	457 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	458
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	459 default:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	460 (*current_liboctave_error_handler)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	461 ("rand: invalid distribution ID = %d", current_distribution);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	462 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	463 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	464 retval = dretval;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	465 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	466 else
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	467 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	468 switch (current_distribution)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	469 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	470 case uniform_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	471 retval = oct_float_randu ();
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	472 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	473
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	474 case normal_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	475 retval = oct_float_randn ();
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	476 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	477
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	478 case expon_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	479 retval = oct_float_rande ();
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	480 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	481
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	482 case poisson_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	483 // Keep poisson distribution in double precision for accuracy
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	484 retval = oct_randp (a);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	485 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	486
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	487 case gamma_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	488 retval = oct_float_randg (a);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	489 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	490
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	491 default:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	492 (*current_liboctave_error_handler)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	493 ("rand: invalid distribution ID = %d", current_distribution);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	494 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	495 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	496
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	497 save_state ();
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	498 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	499
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	500 return retval;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	501 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	502
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	503 Array<double>
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	504 octave_rand::do_vector (octave_idx_type n, double a)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	505 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	506 Array<double> retval;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	507
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	508 if (n > 0)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	509 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	510 retval.clear (n, 1);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	511
20229 5dfaaaae784f Deprecate Array::capacity() and Sparse::capacity() for numel() and nzmax(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	512 fill (retval.numel (), retval.fortran_vec (), a);
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	513 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	514 else if (n < 0)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	515 (*current_liboctave_error_handler) ("rand: invalid negative argument");
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	516
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	517 return retval;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	518 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	519
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	520 Array<float>
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	521 octave_rand::do_float_vector (octave_idx_type n, float a)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	522 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	523 Array<float> retval;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	524
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	525 if (n > 0)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	526 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	527 retval.clear (n, 1);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	528
20229 5dfaaaae784f Deprecate Array::capacity() and Sparse::capacity() for numel() and nzmax(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	529 fill (retval.numel (), retval.fortran_vec (), a);
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	530 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	531 else if (n < 0)
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	532 (*current_liboctave_error_handler) ("rand: invalid negative argument");
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	533
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	534 return retval;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	535 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	536
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	537 NDArray
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	538 octave_rand::do_nd_array (const dim_vector& dims, double a)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	539 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	540 NDArray retval;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	541
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	542 if (! dims.all_zero ())
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	543 {
9647 54f45f883a53 optimize & extend randperm Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	544 retval.clear (dims);
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	545
20229 5dfaaaae784f Deprecate Array::capacity() and Sparse::capacity() for numel() and nzmax(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	546 fill (retval.numel (), retval.fortran_vec (), a);
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	547 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	548
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	549 return retval;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	550 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	551
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	552 FloatNDArray
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	553 octave_rand::do_float_nd_array (const dim_vector& dims, float a)
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	554 {
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	555 FloatNDArray retval;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	556
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	557 if (! dims.all_zero ())
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	558 {
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	559 retval.clear (dims);
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	560
20229 5dfaaaae784f Deprecate Array::capacity() and Sparse::capacity() for numel() and nzmax(). Carnë Draug <carandraug@octave.org> parents: 19697 diff changeset	561 fill (retval.numel (), retval.fortran_vec (), a);
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	562 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	563
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	564 return retval;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	565 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	566
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	567 // Make the random number generator give us a different sequence every
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	568 // time we start octave unless we specifically set the seed. The
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	569 // technique used below will cycle monthly, but it it does seem to
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	570 // work ok to give fairly different seeds each time Octave starts.
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	571
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	572 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	573 octave_rand::initialize_ranlib_generators (void)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	574 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	575 octave_localtime tm;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	576 int stored_distribution = current_distribution;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	577 F77_FUNC (setcgn, SETCGN) (uniform_dist);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	578
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14655 diff changeset	579 int hour = tm.hour () + 1;
460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14655 diff changeset	580 int minute = tm.min () + 1;
460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14655 diff changeset	581 int second = tm.sec () + 1;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	582
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14655 diff changeset	583 int32_t s0 = tm.mday () * hour * minute * second;
10331 6a50aa2cdf00 never compile ranlib with 8-byte integers John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	584 int32_t s1 = hour * minute * second;
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	585
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	586 s0 = force_to_fit_range (s0, 1, 2147483563);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	587 s1 = force_to_fit_range (s1, 1, 2147483399);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	588
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	589 F77_FUNC (setall, SETALL) (s0, s1);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	590 F77_FUNC (setcgn, SETCGN) (stored_distribution);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	591 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	592
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	593 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	594 octave_rand::initialize_mersenne_twister (void)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	595 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	596 oct_init_by_entropy ();
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	597
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	598 ColumnVector s = get_internal_state ();
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	599
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	600 rand_states[uniform_dist] = s;
8124 d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	601
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	602 oct_init_by_entropy ();
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	603 s = get_internal_state ();
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	604 rand_states[normal_dist] = s;
8124 d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	605
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	606 oct_init_by_entropy ();
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	607 s = get_internal_state ();
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	608 rand_states[expon_dist] = s;
8124 d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	609
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	610 oct_init_by_entropy ();
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	611 s = get_internal_state ();
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	612 rand_states[poisson_dist] = s;
8124 d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	613
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	614 oct_init_by_entropy ();
d227d096d49e oct-rand.cc (initialize_mersenne_twister): use separate initializations for each generator Brian Gough parents: 7537 diff changeset	615 s = get_internal_state ();
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	616 rand_states[gamma_dist] = s;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	617 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	618
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	619 ColumnVector
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	620 octave_rand::get_internal_state (void)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	621 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	622 ColumnVector s (MT_N + 1);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	623
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	624 OCTAVE_LOCAL_BUFFER (uint32_t, tmp, MT_N + 1);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	625
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	626 oct_get_state (tmp);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	627
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	628 for (octave_idx_type i = 0; i <= MT_N; i++)
15018 3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	629 s.elem (i) = static_cast<double> (tmp[i]);
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	630
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	631 return s;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	632 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	633
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	634 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	635 octave_rand::save_state (void)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	636 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	637 rand_states[current_distribution] = get_internal_state ();;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	638 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	639
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	640 int
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	641 octave_rand::get_dist_id (const std::string& d)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	642 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	643 int retval = unknown_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	644
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	645 if (d == "uniform" \|\| d == "rand")
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	646 retval = uniform_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	647 else if (d == "normal" \|\| d == "randn")
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	648 retval = normal_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	649 else if (d == "exponential" \|\| d == "rande")
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	650 retval = expon_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	651 else if (d == "poisson" \|\| d == "randp")
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	652 retval = poisson_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	653 else if (d == "gamma" \|\| d == "randg")
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	654 retval = gamma_dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	655 else
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	656 (*current_liboctave_error_handler)
15466 d174210ce1ec use ' instead of ` in error messages, warnings and most comments John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	657 ("rand: invalid distribution '%s'", d.c_str ());
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	658
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	659 return retval;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	660 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	661
17258 213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	662 // Guarantee reproducible conversion of negative initialization values to
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	663 // random number algorithm. Note that Matlab employs slightly different rules.
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	664 // 1) Seed saturates at 2^32-1 for any value larger than that.
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	665 // 2) NaN, Inf are translated to 2^32-1.
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	666 // 3) -Inf is translated to 0.
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	667 static uint32_t
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	668 double2uint32 (double d)
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	669 {
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	670 uint32_t u;
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	671 static const double TWOUP32 = std::numeric_limits<uint32_t>::max() + 1.0;
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	672
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	673 if (! xfinite (d))
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	674 u = 0;
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	675 else
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	676 {
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	677 d = fmod (d, TWOUP32);
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	678 if (d < 0)
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	679 d += TWOUP32;
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	680 u = static_cast<uint32_t> (d);
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	681 }
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	682
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	683 return u;
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	684 }
213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	685
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	686 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	687 octave_rand::set_internal_state (const ColumnVector& s)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	688 {
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20229 diff changeset	689 octave_idx_type len = s.numel ();
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	690 octave_idx_type n = len < MT_N + 1 ? len : MT_N + 1;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	691
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	692 OCTAVE_LOCAL_BUFFER (uint32_t, tmp, MT_N + 1);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	693
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	694 for (octave_idx_type i = 0; i < n; i++)
17258 213ee68b59da Handle out-of-range values consistently when initializing random number generator. Philipp Kutin <philipp.kutin@gmail.com> parents: 16971 diff changeset	695 tmp[i] = double2uint32 (s.elem (i));
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	696
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	697 if (len == MT_N + 1 && tmp[MT_N] <= MT_N && tmp[MT_N] > 0)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	698 oct_set_state (tmp);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	699 else
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	700 oct_init_by_array (tmp, len);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	701 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	702
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	703 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	704 octave_rand::switch_to_generator (int dist)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	705 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	706 if (dist != current_distribution)
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	707 {
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	708 current_distribution = dist;
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	709
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	710 set_internal_state (rand_states[dist]);
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	711 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	712 }
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	713
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	714 #define MAKE_RAND(len) \
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	715 do \
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	716 { \
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	717 double val; \
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	718 for (volatile octave_idx_type i = 0; i < len; i++) \
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	719 { \
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	720 octave_quit (); \
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	721 RAND_FUNC (val); \
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	722 v[i] = val; \
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	723 } \
4308 b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	724 } \
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	725 while (0)
b738d1a02adb [project @ 2003-01-24 19:37:12 by jwe] jwe parents: diff changeset	726
7537 a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	727 void
a2950622f070 make octave_rand a proper singleton class John W. Eaton <jwe@octave.org> parents: 7536 diff changeset	728 octave_rand::fill (octave_idx_type len, double *v, double a)
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	729 {
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	730 if (len < 1)
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	731 return;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	732
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	733 switch (current_distribution)
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	734 {
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	735 case uniform_dist:
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	736 if (use_old_generators)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	737 {
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	738 #define RAND_FUNC(x) F77_FUNC (dgenunf, DGENUNF) (0.0, 1.0, x)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	739 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	740 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	741 }
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	742 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	743 oct_fill_randu (len, v);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	744 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	745
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	746 case normal_dist:
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	747 if (use_old_generators)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	748 {
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	749 #define RAND_FUNC(x) F77_FUNC (dgennor, DGENNOR) (0.0, 1.0, x)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	750 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	751 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	752 }
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	753 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	754 oct_fill_randn (len, v);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	755 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	756
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	757 case expon_dist:
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	758 if (use_old_generators)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	759 {
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	760 #define RAND_FUNC(x) F77_FUNC (dgenexp, DGENEXP) (1.0, x)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	761 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	762 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	763 }
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	764 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	765 oct_fill_rande (len, v);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	766 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	767
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	768 case poisson_dist:
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	769 if (use_old_generators)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	770 {
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	771 if (a < 0.0 \|\| ! xfinite (a))
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	772 #define RAND_FUNC(x) x = octave_NaN;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	773 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	774 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	775 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	776 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	777 // workaround bug in ignpoi, by calling with different Mu
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	778 double tmp;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	779 F77_FUNC (dignpoi, DIGNPOI) (a + 1, tmp);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	780 #define RAND_FUNC(x) F77_FUNC (dignpoi, DIGNPOI) (a, x)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	781 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	782 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	783 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	784 }
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	785 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	786 oct_fill_randp (a, len, v);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	787 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	788
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	789 case gamma_dist:
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	790 if (use_old_generators)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	791 {
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	792 if (a <= 0.0 \|\| ! xfinite (a))
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	793 #define RAND_FUNC(x) x = octave_NaN;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	794 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	795 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	796 else
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	797 #define RAND_FUNC(x) F77_FUNC (dgengam, DGENGAM) (1.0, a, x)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	798 MAKE_RAND (len);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	799 #undef RAND_FUNC
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	800 }
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	801 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	802 oct_fill_randg (a, len, v);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	803 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	804
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	805 default:
7536 4dda6fbc8ba6 oct-rand.cc: fix typo John W. Eaton <jwe@octave.org> parents: 7535 diff changeset	806 (*current_liboctave_error_handler)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	807 ("rand: invalid distribution ID = %d", current_distribution);
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	808 break;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	809 }
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	810
7533 ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	811 save_state ();
ff52243af934 save state separately for each MT random number generator John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	812
5730 109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	813 return;
109fdf7b3dcb [project @ 2006-04-03 19:18:26 by jwe] jwe parents: 5307 diff changeset	814 }
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	815
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	816 void
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	817 octave_rand::fill (octave_idx_type len, float *v, float a)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	818 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	819 if (len < 1)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	820 return;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	821
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	822 switch (current_distribution)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	823 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	824 case uniform_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	825 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	826 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	827 #define RAND_FUNC(x) F77_FUNC (dgenunf, DGENUNF) (0.0, 1.0, x)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	828 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	829 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	830 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	831 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	832 oct_fill_float_randu (len, v);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	833 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	834
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	835 case normal_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	836 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	837 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	838 #define RAND_FUNC(x) F77_FUNC (dgennor, DGENNOR) (0.0, 1.0, x)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	839 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	840 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	841 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	842 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	843 oct_fill_float_randn (len, v);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	844 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	845
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	846 case expon_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	847 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	848 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	849 #define RAND_FUNC(x) F77_FUNC (dgenexp, DGENEXP) (1.0, x)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	850 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	851 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	852 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	853 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	854 oct_fill_float_rande (len, v);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	855 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	856
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	857 case poisson_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	858 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	859 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	860 double da = a;
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	861 if (da < 0.0 \|\| ! xfinite (a))
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	862 #define RAND_FUNC(x) x = octave_NaN;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	863 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	864 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	865 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	866 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	867 // workaround bug in ignpoi, by calling with different Mu
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	868 double tmp;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	869 F77_FUNC (dignpoi, DIGNPOI) (da + 1, tmp);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	870 #define RAND_FUNC(x) F77_FUNC (dignpoi, DIGNPOI) (da, x)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	871 MAKE_RAND (len);
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	872 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	873 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	874 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	875 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	876 oct_fill_float_randp (a, len, v);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	877 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	878
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	879 case gamma_dist:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	880 if (use_old_generators)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	881 {
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	882 double da = a;
16971 259c1f295a1e Use xfinite to replace some (isinf \|\| isnan) instances in C++ code. Rik <rik@octave.org> parents: 15866 diff changeset	883 if (da <= 0.0 \|\| ! xfinite (a))
14655 43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	884 #define RAND_FUNC(x) x = octave_NaN;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	885 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	886 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	887 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	888 #define RAND_FUNC(x) F77_FUNC (dgengam, DGENGAM) (1.0, da, x)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	889 MAKE_RAND (len);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	890 #undef RAND_FUNC
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	891 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	892 else
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	893 oct_fill_float_randg (a, len, v);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	894 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	895
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	896 default:
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	897 (*current_liboctave_error_handler)
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	898 ("rand: invalid distribution ID = %d", current_distribution);
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	899 break;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	900 }
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	901
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	902 save_state ();
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	903
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	904 return;
43db83eff9db Implement single precision rand, randn, rande, randg and randp generators (bug #34351, #36293) David Bateman <dbateman@free.fr> parents: 14138 diff changeset	905 }

Mercurial > octave

annotate liboctave/numeric/oct-rand.cc @ 21136:7cac4e7458f2