test/renders/aniso.cal

{ RCSid $Id: aniso.cal,v 1.1 2018/11/17 22:09:12 greg Exp $ }
{ Ward-Geisler-Moroder-Duer anisotropic BSDF model }

{ Final 5 real arguments are alpha_x and alpha_y roughnesses and Up vector }

ax = arg(AC-4);         { ax is roughness in Up orientation }
ay = arg(AC-3);

arefl(hx,hy,hz) = exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) *
                        (hx*hx + hy*hy + hz*hz) /
                        (PI*ax*ay*hz*hz*hz*hz);

axmit(hx,hy,hz) = exp(-(hx*hx/(ax*ax) + hy*hy/(ay*ay))/(hz*hz)) /
                        (4*PI*ax*ay);

dobsdf(sdx,sdy,sdz) = if(sdz, arefl(Idx+sdx, Idy+sdy, Idz+sdz),
                        axmit(Idx+sdx, Idy+sdy, Idz-sdz) / sqrt(-sdz*Idz) );

absdf(ox,oy,oz) = dobsdf(Ldx(ox,oy,oz), Ldy(ox,oy,oz), Ldz(ox,oy,oz));
Revision:	1.2
Committed:	Sat Dec 1 02:06:04 2018 UTC (6 years, 5 months ago) by greg
Branch:	MAIN
Changes since 1.1:	+4 -25 lines
Log Message:	Finished additions for BSDF's, etc. to complete primitive testing
#	User	Rev	Content
1	greg	1.2	{ RCSid $Id: aniso.cal,v 1.1 2018/11/17 22:09:12 greg Exp $ }
2	greg	1.1	{ Ward-Geisler-Moroder-Duer anisotropic BSDF model }
3
4			{ Final 5 real arguments are alpha_x and alpha_y roughnesses and Up vector }
5
6			ax = arg(AC-4); { ax is roughness in Up orientation }
7			ay = arg(AC-3);
8
9			arefl(hx,hy,hz) = exp(-(hxhx/(axax) + hyhy/(ayay))/(hzhz))
10			(hxhx + hyhy + hz*hz) /
11			(PIaxayhzhzhzhz);
12
13			axmit(hx,hy,hz) = exp(-(hxhx/(axax) + hyhy/(ayay))/(hz*hz)) /
14			(4PIax*ay);
15
16	greg	1.2	dobsdf(sdx,sdy,sdz) = if(sdz, arefl(Idx+sdx, Idy+sdy, Idz+sdz),
17			axmit(Idx+sdx, Idy+sdy, Idz-sdz) / sqrt(-sdz*Idz) );
18	greg	1.1
19	greg	1.2	absdf(ox,oy,oz) = dobsdf(Ldx(ox,oy,oz), Ldy(ox,oy,oz), Ldz(ox,oy,oz));