ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rambpos.cal
(Generate patch)

Comparing ray/src/util/rambpos.cal (file contents):
Revision 2.5 by greg, Fri May 9 23:49:05 2014 UTC vs.
Revision 2.6 by greg, Thu May 15 17:30:53 2014 UTC

# Line 7 | Line 7 | and(a,b) : if(a,b,a);
7   lvl = $10;
8   wt = $11;
9   r0 = SF*$12; r1 = SF*$13;
10 + eps = r0*.015;
11   acond = and( r0-5e-5, and( wt-MW, if(LV+.5, and(lvl-LV+.5,LV+.5-lvl), 6.5-lvl) ) );
11 px = $1; py = $2; pz = $3;
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 + px = $1 + nx*eps;
18 + py = $2 + ny*eps;
19 + pz = $3 + nz*eps;
20   avr = .5/PI*$14; avg = .5/PI*$15; avb = .5/PI*$16;
21                                  { magnified position gradient }
22   pgx = ($17*ux + $18*vx)/SF;
# Line 25 | Line 28 | pcond = and( acond, pg - .0001 );
28   dgx = $19*ux + $20*vx;
29   dgy = $19*uy + $20*vy;
30   dgz = $19*uz + $20*vz;
31 + dg = sqrt(dgx*dgx + dgy*dgy + dgz*dgz);
32   val = avr + avg + avb;          { ambient position color }
33   agr = select(lvl+1, val, 0, 0, val, 0, val, val);
34   agg = select(lvl+1, 0, val, 0, val, val, 0, val/2);
35   agb = select(lvl+1, 0, 0, val, 0, val, val, val/2);
36   psiz = r0/10;
37 < dcond = and( acond, dgx*dgx+dgy*dgy+dgz*dgz - 1e-8 );
37 > dcond = and( acond, dg - .0001 );
38                                  { base polygon vertices }
39 < px1 = px + r0*ux + r1*vx + .0001*nx;
40 < py1 = py + r0*uy + r1*vy + .0001*ny;
41 < pz1 = pz + r0*uz + r1*vz + .0001*nz;
42 < px2 = px - r0*ux + r1*vx + .0001*nx;
43 < py2 = py - r0*uy + r1*vy + .0001*ny;
44 < pz2 = pz - r0*uz + r1*vz + .0001*nz;
45 < px3 = px - r0*ux - r1*vx + .0001*nx;
46 < py3 = py - r0*uy - r1*vy + .0001*ny;
47 < pz3 = pz - r0*uz - r1*vz + .0001*nz;
48 < px4 = px + r0*ux - r1*vx + .0001*nx;
49 < py4 = py + r0*uy - r1*vy + .0001*ny;
50 < pz4 = pz + r0*uz - r1*vz + .0001*nz;
39 > px1 = px + r0*ux + r1*vx;
40 > py1 = py + r0*uy + r1*vy;
41 > pz1 = pz + r0*uz + r1*vz;
42 > px2 = px - r0*ux + r1*vx;
43 > py2 = py - r0*uy + r1*vy;
44 > pz2 = pz - r0*uz + r1*vz;
45 > px3 = px - r0*ux - r1*vx;
46 > py3 = py - r0*uy - r1*vy;
47 > pz3 = pz - r0*uz - r1*vz;
48 > px4 = px + r0*ux - r1*vx;
49 > py4 = py + r0*uy - r1*vy;
50 > pz4 = pz + r0*uz - r1*vz;
51   alen = psiz*3;                  { arrow length & coordinates }
52   cx0 = px + psiz*psiz/alen*pgx/pg;
53   cy0 = py + psiz*psiz/alen*pgy/pg;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines