Displaying 1 result from an estimated 1 matches for "box_muller".
2003 Apr 25
4
Kinderman-Ramage (PR#2846)
...tt = 0.479727404222441-0.595507138015940*fmin2(u2,u3);
if(fmax2(u2,u3) <= 0.805577924423817)
return (u2<u3) ? tt : -tt;
+ /** <added> **/
+ if(0.053377549506886*fabs(u2-u3) <= g(tt))
+ return (u2<u3) ? tt : -tt;
+ /** </added> **/
}
case BOX_MULLER:
if(BM_norm_keep != 0.0) { /* An exact test is intentional */
----------------------------------------------------------
$ uname -a
Linux 2.4.18-3 #1 Thu Apr 18 07:37:53 EDT 2002 i686 unknown
$ R
R : Copyright 2003, The R Development Core Team
Version 1.7.0 (2003-04-16)
R is free software...