6 |
|
ax = arg(AC-4); { ax is roughness in Up orientation } |
7 |
|
ay = arg(AC-3); |
8 |
|
|
9 |
< |
arefl(hx,hy,hz) = exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) * |
9 |
> |
arefl(hx,hy,hz) = Exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) * |
10 |
|
(hx*hx + hy*hy + hz*hz) / |
11 |
|
(PI*ax*ay*hz*hz*hz*hz); |
12 |
|
|
13 |
< |
axmit(hx,hy,hz) = exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) / |
13 |
> |
axmit(hx,hy,hz) = Exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) / |
14 |
|
(4*PI*ax*ay); |
15 |
|
|
16 |
|
dobsdf(sdx,sdy,sdz) = if(sdz, arefl(Idx+sdx, Idy+sdy, Idz+sdz), |