ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cal/cal/patch3w.cal
Revision: 1.1
Committed: Sat Feb 22 02:07:21 2003 UTC (21 years, 1 month ago) by greg
Branch: MAIN
CVS Tags: rad5R2, rad4R2P2, rad5R0, rad5R1, rad3R7P2, rad3R7P1, rad4R2, rad4R1, rad4R0, rad3R5, rad3R6, rad3R6P1, rad3R8, rad3R9, rad4R2P1
Log Message:
Changes and check-in for 3.5 release
Includes new source files and modifications not recorded for many years
See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release

File Contents

# Content
1 {
2 Uniform Rational Patch
3
4 02Mar90
5
6 Define Px(i,j), Py(i,j), Pz(i,j) and W(i,j), nrows, ncols
7 }
8 { default basis functions }
9 br(a,b,c,d,s) = b3(a,b,c,d,s);
10 bc(a,b,c,d,t) = b3(a,b,c,d,t);
11 b3(a,b,c,d,p) = bezier(a,b,c,d,p);
12 { default patch interpolation }
13 rint = interpolate;
14 cint = interpolate;
15 interpolate = 1;
16
17 x(s,t) = patch(s,t,Px);
18 y(s,t) = patch(s,t,Py);
19 z(s,t) = patch(s,t,Pz);
20
21 patch(s,t,p) = bc(P2(s,c0(t)+1,p), P2(s,c0(t)+2,p),
22 P2(s,c0(t)+3,p), P2(s,c0(t)+4,p), cf(t)) /
23 bc(W2(s,c0(t)+1), W2(s,c0(t)+2),
24 W2(s,c0(t)+3), W2(s,c0(t)+4), cf(t)) ;
25
26 P2(s,j,p) = br(W(1+r0(s),j)*p(1+r0(s),j), W(2+r0(s),j)*p(2+r0(s),j),
27 W(3+r0(s),j)*p(3+r0(s),j), W(4+r0(s),j)*p(4+r0(s),j), rf(s)) ;
28
29 W2(s,j) = br(W(1+r0(s),j), W(2+r0(s),j), W(3+r0(s),j), W(4+r0(s),j), rf(s)) ;
30
31 r0(s) = rmult*floor(clamp(s)*npr);
32 rf(s) = s*npr - floor(clamp(s)*npr);
33 c0(t) = cmult*floor(clamp(t)*npc);
34 cf(t) = t*npc - floor(clamp(t)*npc);
35
36 rmult = if(rint, 4, 1);
37 cmult = if(cint, 4, 1);
38 npr = if(rint, nrows/4, nrows-3);
39 npc = if(cint, ncols/4, ncols-3);
40
41 clamp(p) = if(p, if(1-p, p, .9999999), .0000001);