ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rambpos.cal
Revision: 2.3
Committed: Fri Apr 25 22:56:06 2014 UTC (9 years, 11 months ago) by greg
Branch: MAIN
Changes since 2.2: +2 -2 lines
Log Message:
Changed ambient value color to match 50% diffuse surface

File Contents

# User Rev Content
1 greg 2.3 { RCSid $Id: rambpos.cal,v 2.2 2014/04/25 18:13:47 greg Exp $ }
2 greg 2.1 {
3     Compute ambient position factors for genambpos.pl
4     }
5     and(a,b) : if(a,b,a);
6     { level selector }
7 greg 2.2 lvl = $10;
8     wt = $11;
9     cond = and( wt-MW, if(LV+.5, and(lvl-LV+.5,LV+.5-lvl), 6.5-lvl) );
10 greg 2.1 px = $1; py = $2; pz = $3;
11     r0 = SF*$12; r1 = SF*$13;
12     nx = $4; ny = $5; nz = $6;
13     ux = $7; uy = $8; uz = $9;
14     vx = ny*uz - nz*uy;
15     vy = nz*ux - nx*uz;
16     vz = nx*uy - ny*ux;
17 greg 2.3 avr = .5/PI*$14; avg = .5/PI*$15; avb = .5/PI*$16;
18 greg 2.1 { magnified position gradient }
19     pgx = ($17*ux + $18*vx)/SF;
20     pgy = ($17*uy + $18*vy)/SF;
21     pgz = ($17*uz + $18*vz)/SF;
22     pg = sqrt(pgx*pgx + pgy*pgy + pgz*pgz);
23     { direction gradient }
24     dgx = $19*ux + $20*vx;
25     dgy = $19*uy + $20*vy;
26     dgz = $19*uz + $20*vz;
27     val = avr + avg + avb; { ambient position color }
28 greg 2.2 agr = select(lvl+1, val, 0, 0, val, 0, val, val);
29     agg = select(lvl+1, 0, val, 0, val, val, 0, val/2);
30     agb = select(lvl+1, 0, 0, val, 0, val, val, val/2);
31 greg 2.1 psiz = r0/10;
32     { base polygon vertices }
33     px1 = px + r0*ux + r1*vx + .0001*nx;
34     py1 = py + r0*uy + r1*vy + .0001*ny;
35     pz1 = pz + r0*uz + r1*vz + .0001*nz;
36     px2 = px - r0*ux + r1*vx + .0001*nx;
37     py2 = py - r0*uy + r1*vy + .0001*ny;
38     pz2 = pz - r0*uz + r1*vz + .0001*nz;
39     px3 = px - r0*ux - r1*vx + .0001*nx;
40     py3 = py - r0*uy - r1*vy + .0001*ny;
41     pz3 = pz - r0*uz - r1*vz + .0001*nz;
42     px4 = px + r0*ux - r1*vx + .0001*nx;
43     py4 = py + r0*uy - r1*vy + .0001*ny;
44     pz4 = pz + r0*uz - r1*vz + .0001*nz;
45     alen = psiz*3; { arrow length & coordinates }
46     cx0 = px + psiz*psiz/alen*pgx/pg;
47     cy0 = py + psiz*psiz/alen*pgy/pg;
48     cz0 = pz + psiz*psiz/alen*pgz/pg;
49     cr0 = psiz*sqrt(1 - psiz*psiz/(alen*alen));
50     cx1 = px + alen*pgx/pg;
51     cy1 = py + alen*pgy/pg;
52     cz1 = pz + alen*pgz/pg;