| 1 |
< |
{ SCCSid "$SunId$ LBL" } |
| 2 |
< |
|
| 1 |
> |
{ RCSid $Id$ } |
| 2 |
|
{ |
| 3 |
|
Initialization file for Radiance. |
| 4 |
|
|
| 17 |
|
Kx, Ky, Kz - world k unit vector |
| 18 |
|
arg(n) - real arguments, arg(0) is count |
| 19 |
|
|
| 20 |
+ |
For mesh objects, the following are available: |
| 21 |
+ |
|
| 22 |
+ |
Lu, Lv - local (u,v) coordinates |
| 23 |
+ |
|
| 24 |
|
For brdf functions, the following are also available: |
| 25 |
|
|
| 26 |
|
NxP, NyP, NzP - perturbed surface normal |
| 91 |
|
and(a,b) : if( a, b, a ); |
| 92 |
|
or(a,b) : if( a, a, b ); |
| 93 |
|
not(a) : if( a, -1, 1 ); |
| 94 |
+ |
xor(a,b) : if( a, not(b), b ); |
| 95 |
|
abs(x) : if( x, x, -x ); |
| 96 |
|
sgn(x) : if( x, 1, if(-x, -1, 0) ); |
| 97 |
|
sq(x) : x*x; |
| 103 |
|
tri(n,d) : abs( d - mod(n-d,2*d) ); |
| 104 |
|
linterp(t,p0,p1) : (1-t)*p0 + t*p1; |
| 105 |
|
|
| 106 |
< |
noop(v) = v; |
| 107 |
< |
clip(v) = bound(0,v,1); |
| 108 |
< |
noneg(v) = if(v,v,0); |
| 109 |
< |
red(r,g,b) = if(r,r,0); |
| 110 |
< |
green(r,g,b) = if(g,g,0); |
| 111 |
< |
blue(r,g,b) = if(b,b,0); |
| 112 |
< |
grey(r,g,b) = noneg(.265074126*r + .670114631*g + .064811243*b); |
| 113 |
< |
clip_r(r,g,b) = bound(0,r,1); |
| 114 |
< |
clip_g(r,g,b) = bound(0,g,1); |
| 115 |
< |
clip_b(r,g,b) = bound(0,b,1); |
| 116 |
< |
clipgrey(r,g,b) = min(grey(r,g,b),1); |
| 106 |
> |
noop(v) : v; |
| 107 |
> |
clip(v) : bound(0,v,1); |
| 108 |
> |
noneg(v) : if(v,v,0); |
| 109 |
> |
red(r,g,b) : if(r,r,0); |
| 110 |
> |
green(r,g,b) : if(g,g,0); |
| 111 |
> |
blue(r,g,b) : if(b,b,0); |
| 112 |
> |
grey(r,g,b) : noneg(.265074126*r + .670114631*g + .064811243*b); |
| 113 |
> |
clip_r(r,g,b) : bound(0,r,1); |
| 114 |
> |
clip_g(r,g,b) : bound(0,g,1); |
| 115 |
> |
clip_b(r,g,b) : bound(0,b,1); |
| 116 |
> |
clipgrey(r,g,b) : min(grey(r,g,b),1); |
| 117 |
|
|
| 118 |
|
dot(v1,v2) : v1(1)*v2(1) + v1(2)*v2(2) + v1(3)*v2(3); |
| 119 |
|
cross(i,v1,v2) : select(i, v1(2)*v2(3) - v1(3)*v2(2), |
| 120 |
|
v1(3)*v2(1) - v1(1)*v2(3), |
| 121 |
|
v1(1)*v2(2) - v1(2)*v2(1)); |
| 122 |
|
|
| 123 |
< |
fade(near_val,far_val,dist) = far_val + |
| 123 |
> |
fade(near_val,far_val,dist) : far_val + |
| 124 |
|
if (16-dist, (near_val-far_val)/(1+dist*dist), 0); |
| 125 |
|
|
| 126 |
< |
bezier(p1, p2, p3, p4, t) = p1 * (1+t*(-3+t*(3-t))) + |
| 126 |
> |
bezier(p1, p2, p3, p4, t) : p1 * (1+t*(-3+t*(3-t))) + |
| 127 |
|
p2 * 3*t*(1+t*(-2+t)) + |
| 128 |
|
p3 * 3*t*t*(1-t) + |
| 129 |
|
p4 * t*t*t ; |
| 130 |
|
|
| 131 |
< |
bspline(pp, p0, p1, pn, t) = pp * (1/6+t*(-.5+t*(.5-1/6*t))) + |
| 131 |
> |
bspline(pp, p0, p1, pn, t) : pp * (1/6+t*(-.5+t*(.5-1/6*t))) + |
| 132 |
|
p0 * (2/3+t*t*(-1+.5*t)) + |
| 133 |
|
p1 * (1/6+t*(.5+t*(.5-.5*t))) + |
| 134 |
|
pn * (1/6*t*t*t) ; |
| 135 |
|
|
| 136 |
< |
turbulence(x,y,z,s) = if( s-1.01, 0, abs(noise3(x/s,y/s,z/s)*s) + |
| 136 |
> |
turbulence(x,y,z,s) : if( s-1.01, 0, abs(noise3(x/s,y/s,z/s)*s) + |
| 137 |
|
turbulence(x,y,z,2*s) ); |
| 138 |
< |
turbulencea(x,y,z,s) = if( s-1.01, 0, |
| 138 |
> |
turbulencex(x,y,z,s) : if( s-1.01, 0, |
| 139 |
|
sgn(noise3(x/s,y/s,z/s))*noise3x(x/s,y/s,z/s) + |
| 140 |
< |
turbulencea(x,y,z,2*s) ); |
| 141 |
< |
turbulenceb(x,y,z,s) = if( s-1.01, 0, |
| 140 |
> |
turbulencex(x,y,z,2*s) ); |
| 141 |
> |
turbulencey(x,y,z,s) : if( s-1.01, 0, |
| 142 |
|
sgn(noise3(x/s,y/s,z/s))*noise3y(x/s,y/s,z/s) + |
| 143 |
< |
turbulenceb(x,y,z,2*s) ); |
| 144 |
< |
turbulencec(x,y,z,s) = if( s-1.01, 0, |
| 143 |
> |
turbulencey(x,y,z,2*s) ); |
| 144 |
> |
turbulencez(x,y,z,s) : if( s-1.01, 0, |
| 145 |
|
sgn(noise3(x/s,y/s,z/s))*noise3z(x/s,y/s,z/s) + |
| 146 |
< |
turbulencec(x,y,z,2*s) ); |
| 146 |
> |
turbulencez(x,y,z,2*s) ); |
| 147 |
|
|
| 148 |
|
{ Normal distribution from uniform range (0,1) } |
| 149 |
|
|
| 150 |
< |
un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/ |
| 150 |
> |
un2`P(t) : t - (2.515517+t*(.802853+t*.010328))/ |
| 151 |
|
(1+t*(1.432788+t*(.189269+t*.001308))) ; |
| 152 |
< |
un1`private(p) : un2`private(sqrt(-2*log(p))) ; |
| 152 |
> |
un1`P(p) : un2`P(sqrt(-2*log(p))) ; |
| 153 |
|
|
| 154 |
< |
unif2norm(p) : if( .5-p, -un1`private(p), un1`private(1-p) ) ; |
| 154 |
> |
unif2norm(p) : if( .5-p, -un1`P(p), un1`P(1-p) ) ; |
| 155 |
|
|
| 156 |
|
nrand(x) = unif2norm(rand(x)); |
| 157 |
|
|
| 158 |
|
{ Local (u,v) coordinates for planar surfaces } |
| 159 |
< |
crosslen`private = Nx*Nx + Ny*Ny; |
| 159 |
> |
crosslen`P = Nx*Nx + Ny*Ny; |
| 160 |
|
{ U is distance from projected Z-axis } |
| 161 |
< |
U = if( crosslen`private - FTINY, |
| 162 |
< |
(Py*Nx - Px*Ny)/crosslen`private, |
| 161 |
> |
U = if( crosslen`P - FTINY, |
| 162 |
> |
(Py*Nx - Px*Ny)/crosslen`P, |
| 163 |
|
Px); |
| 164 |
|
{ V is defined so that N = U x V } |
| 165 |
< |
V = if( crosslen`private - FTINY, |
| 166 |
< |
Pz - Nz*(Px*Nx + Py*Ny)/crosslen`private, |
| 165 |
> |
V = if( crosslen`P - FTINY, |
| 166 |
> |
Pz - Nz*(Px*Nx + Py*Ny)/crosslen`P, |
| 167 |
|
Py); |
