steve
/
factor
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

random: gauss-random-float is just normal-random-float.

db4
John Benediktsson 2012-03-30 17:55:18 -07:00
parent 1a8e52379a
commit 1e7382b893
1 changed files with 7 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
basis/random/random.factor
View File

@ -111,11 +111,6 @@ ERROR: too-many-samples seq n ;
: normal-random-float ( mean sigma -- n ) : normal-random-float ( mean sigma -- n )
(cos-random-float) (log-sqrt-random-float) * * + ; (cos-random-float) (log-sqrt-random-float) * * + ;
{
{ [ os windows? ] [ "random.windows" require ] }
{ [ os unix? ] [ "random.unix" require ] }
} cond
: lognormal-random-float ( mean sigma -- n ) : lognormal-random-float ( mean sigma -- n )
normal-random-float exp ; normal-random-float exp ;
@ -128,12 +123,13 @@ ERROR: too-many-samples seq n ;
: pareto-random-float ( alpha -- n ) : pareto-random-float ( alpha -- n )
[ 0. 1. uniform-random-float 1 swap - ] dip [ 1. swap / ] bi@ ^ ; [ 0. 1. uniform-random-float 1 swap - ] dip [ 1. swap / ] bi@ ^ ;
: gauss-random-float ( mean sigma -- n )
0. 1. uniform-random-float 1 swap - log -2 * sqrt
(cos-random-float) * * + ;
: beta-random-float ( alpha beta -- n ) : beta-random-float ( alpha beta -- n )
[ 1. gauss-random-float ] dip over zero? [ 1. normal-random-float ] dip over zero?
[ 2drop 0 ] [ 1. gauss-random-float dupd + / ] if ; [ 2drop 0 ] [ 1. normal-random-float dupd + / ] if ;
{
{ [ os windows? ] [ "random.windows" require ] }
{ [ os unix? ] [ "random.unix" require ] }
} cond
"random.mersenne-twister" require "random.mersenne-twister" require