ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cal/cal/norm.cal
(Generate patch)

Comparing ray/src/cal/cal/norm.cal (file contents):
Revision 1.2 by greg, Wed Apr 25 01:06:25 2012 UTC vs.
Revision 1.3 by greg, Wed Apr 19 22:09:12 2017 UTC

# Line 13 | Line 13
13  
14   NORMF : 1/sqrt(2*PI) ;
15  
16 < Z(u) = NORMF * exp(-u*u/2) ;
16 > Z(u) : NORMF * exp(-u*u/2) ;
17  
18 < Q(u) = if( u, Q1(u), 1-Q1(-u) ) ;
19 <
20 < Q1(u) = Z(u) * Q2(1/(1+.2316419*u)) ;
21 <
22 < Q2(t) = t*(.31938153+t*(-.356563782+t*(1.781477937+
18 > Q2(t) : t*(.31938153+t*(-.356563782+t*(1.781477937+
19                  t*(-1.821255978+t*1.330274429)))) ;
20  
21 < u(p) = if( .5-p, u1(p), -u1(1-p) ) ;
21 > Q1(u) : Z(u) * Q2(1/(1+.2316419*u)) ;
22  
23 < u1(p) = u2(sqrt(-2*log(p))) ;
23 > Q(u) : if( u, Q1(u), 1-Q1(-u) ) ;
24  
25 < u2(t) = t - (2.515517+t*(.802853+t*.010328))/
25 > u2(t) : t - (2.515517+t*(.802853+t*.010328))/
26                  (1+t*(1.432788+t*(.189269+t*.001308))) ;
27 +
28 + u1(p) : u2(sqrt(-2*log(p))) ;
29 +
30 + u(p) : if( .5-p, u1(p), -u1(1-p) ) ;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines