--- ray/src/rt/rayinit.cal	1992/10/21 11:20:26	2.5
+++ ray/src/rt/rayinit.cal	1993/12/09 16:56:28	2.8
@@ -67,6 +67,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));
+Atan2(y,x) : if(x*x+y*y, atan2(y,x), 0);
 Exp(x) : if(-x-100, 0, exp(x));
 Sqrt(x) : if(x, sqrt(x), 0);
 
@@ -138,13 +139,13 @@ un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/
 		(1+t*(1.432788+t*(.189269+t*.001308))) ;
 un1`private(p) : un2`private(sqrt(-2*log(p))) ;
 
-unif2norm(p) : if( .5-p, un1`private(p), -un1`private(1-p) ) ;
+unif2norm(p) : if( .5-p, -un1`private(p), un1`private(1-p) ) ;
 
 nrand(x) = unif2norm(rand(x));
 
 			{ Local (u,v) coordinates for planar surfaces }
 crosslen`private = Nx*Nx + Ny*Ny;
-			{ U is distance from origin in XY-plane }
+			{ U is distance from projected Z-axis }
 U = if( crosslen`private - FTINY,
 		(Py*Nx - Px*Ny)/crosslen`private,
 		Px);