--- ray/src/rt/rayinit.cal 1991/11/12 17:10:31 2.1 +++ ray/src/rt/rayinit.cal 1991/11/25 09:51:57 2.2 @@ -66,7 +66,7 @@ noise3d(i,x,y,z) = select(i, noise3a(x,y,z), noise3b(x bound(a,x,b) : if(a-x, a, if(x-b, b, x)); Acos(x) : acos(bound(-1,x,1)); Asin(x) : asin(bound(-1,x,1)); -Exp(x) : if(-x-60, 0, exp(x)); +Exp(x) : if(-x-100, 0, exp(x)); Sqrt(x) : if(x, sqrt(x), 0); { Useful constants } @@ -135,7 +135,7 @@ turbulencec(x,y,z,s) = if( s-1.01, 0, un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/ (1+t*(1.432788+t*(.189269+t*.001308))) ; -un1`private(p) : un2`private(sqrt(log(1/p/p))) ; +un1`private(p) : un2`private(sqrt(-2*log(p))) ; unif2norm(p) : if( .5-p, un1`private(p), -un1`private(1-p) ) ;