3 |
|
{ |
4 |
|
Initialization file for Radiance. |
5 |
|
|
6 |
– |
4/14/86 |
7 |
– |
} |
8 |
– |
|
9 |
– |
{ |
6 |
|
The following are predefined: |
7 |
|
|
8 |
|
Dx, Dy, Dz - ray direction |
9 |
|
Nx, Ny, Nz - surface normal |
10 |
|
Px, Py, Pz - intersection point |
11 |
|
T - distance from start |
12 |
+ |
Ts - single ray (shadow) distance |
13 |
|
Rdot - ray dot product |
14 |
|
S - world scale |
15 |
|
Tx, Ty, Tz - world origin |
24 |
|
RdotP - perturbed ray dot product |
25 |
|
CrP, CgP, CbP - perturbed material color |
26 |
|
|
27 |
+ |
For prism1 and prism2 types, the following are available: |
28 |
+ |
|
29 |
+ |
DxA, DyA, DzA - direction to target light source |
30 |
+ |
|
31 |
|
Library functions: |
32 |
|
|
33 |
|
if(a, b, c) - if a positive, return b, else c |
71 |
|
bound(a,x,b) : if(a-x, a, if(x-b, b, x)); |
72 |
|
Acos(x) : acos(bound(-1,x,1)); |
73 |
|
Asin(x) : asin(bound(-1,x,1)); |
74 |
< |
Exp(x) : if(-x-60, 0, exp(x)); |
74 |
> |
Atan2(y,x) : if(x*x+y*y, atan2(y,x), 0); |
75 |
> |
Exp(x) : if(-x-100, 0, exp(x)); |
76 |
|
Sqrt(x) : if(x, sqrt(x), 0); |
77 |
|
|
78 |
|
{ Useful constants } |
97 |
|
|
98 |
|
noop(v) = v; |
99 |
|
clip(v) = bound(0,v,1); |
100 |
< |
noneg(v) = max(0,v); |
101 |
< |
red(r,g,b) = r; |
102 |
< |
green(r,g,b) = g; |
103 |
< |
blue(r,g,b) = b; |
104 |
< |
grey(r,g,b) = .3*r + .59*g + .11*b; |
100 |
> |
noneg(v) = if(v,v,0); |
101 |
> |
red(r,g,b) = if(r,r,0); |
102 |
> |
green(r,g,b) = if(g,g,0); |
103 |
> |
blue(r,g,b) = if(b,b,0); |
104 |
> |
grey(r,g,b) = noneg(.265074126*r + .670114631*g + .064811243*b); |
105 |
|
clip_r(r,g,b) = bound(0,r,1); |
106 |
|
clip_g(r,g,b) = bound(0,g,1); |
107 |
|
clip_b(r,g,b) = bound(0,b,1); |
108 |
< |
clipgrey(r,g,b) = bound(0,grey(r,g,b),1); |
108 |
> |
clipgrey(r,g,b) = min(grey(r,g,b),1); |
109 |
|
|
110 |
|
dot(v1,v2) : v1(1)*v2(1) + v1(2)*v2(2) + v1(3)*v2(3); |
111 |
|
cross(i,v1,v2) : select(i, v1(2)*v2(3) - v1(3)*v2(2), |
136 |
|
turbulencec(x,y,z,s) = if( s-1.01, 0, |
137 |
|
sgn(noise3(x/s,y/s,z/s))*noise3c(x/s,y/s,z/s) + |
138 |
|
turbulencec(x,y,z,2*s) ); |
139 |
+ |
|
140 |
+ |
{ Normal distribution from uniform range (0,1) } |
141 |
+ |
|
142 |
+ |
un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/ |
143 |
+ |
(1+t*(1.432788+t*(.189269+t*.001308))) ; |
144 |
+ |
un1`private(p) : un2`private(sqrt(-2*log(p))) ; |
145 |
+ |
|
146 |
+ |
unif2norm(p) : if( .5-p, -un1`private(p), un1`private(1-p) ) ; |
147 |
+ |
|
148 |
+ |
nrand(x) = unif2norm(rand(x)); |
149 |
+ |
|
150 |
+ |
{ Local (u,v) coordinates for planar surfaces } |
151 |
+ |
crosslen`private = Nx*Nx + Ny*Ny; |
152 |
+ |
{ U is distance from projected Z-axis } |
153 |
+ |
U = if( crosslen`private - FTINY, |
154 |
+ |
(Py*Nx - Px*Ny)/crosslen`private, |
155 |
+ |
Px); |
156 |
+ |
{ V is defined so that N = U x V } |
157 |
+ |
V = if( crosslen`private - FTINY, |
158 |
+ |
Pz - Nz*(Px*Nx + Py*Ny)/crosslen`private, |
159 |
+ |
Py); |