ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/falsecolor.csh
(Generate patch)

Comparing ray/src/px/falsecolor.csh (file contents):
Revision 1.2 by greg, Thu May 23 14:22:00 1991 UTC vs.
Revision 1.3 by greg, Thu May 23 16:41:38 1991 UTC

# Line 8 | Line 8 | onintr quit
8   set mult=470
9   set label=Nits
10   set scale=1000
11 + set decades=0
12   set redv='2*v-1'
13   set grnv='if(v-.5,2-2*v,2*v)'
14   set bluv='1-2*v'
# Line 28 | Line 29 | while ($#argv > 0)
29                  shift argv
30                  set label="$argv[1]"
31                  breaksw
32 +        case -log:
33 +                shift argv
34 +                set decades=$argv[1]
35 +                breaksw
36          case -r:
37                  shift argv
38                  set redv="$argv[1]"
# Line 65 | Line 70 | while ($#argv > 0)
70                  set ndivs="$argv[1]"
71                  breaksw
72          default:
73 <                echo bad option "'$argv[1]'"
73 >                echo bad option "'$argv[1]'" >/dev/tty
74                  exit 1
75          endsw
76          shift argv
# Line 83 | Line 88 | btwn(a,b) : if(a-x,-1,b-x);
88   frac(x) : x - floor(x);
89   boundary(a,b) : neq(floor(ndivs*a),floor(ndivs*b));
90  
91 < red=$redv;
92 < grn=$grnv;
93 < blu=$bluv;
91 > map(x) = x;
92 > red = $redv;
93 > grn = $grnv;
94 > blu = $bluv;
95  
96 < v = li(1)*(mult/scale);
97 < vleft = li(1,-1,0)*(mult/scale);
98 < vright = li(1,1,0)*(mult/scale);
99 < vabove = li(1,0,1)*(mult/scale);
100 < vbelow = li(1,0,-1)*(mult/scale);
96 > v = map(li(1)*(mult/scale));
97 > vleft = map(li(1,-1,0)*(mult/scale));
98 > vright = map(li(1,1,0)*(mult/scale));
99 > vabove = map(li(1,0,1)*(mult/scale));
100 > vbelow = map(li(1,0,-1)*(mult/scale));
101   isconta = or(boundary(vleft,vright),boundary(vabove,vbelow));
102   iscontb = if(btwn(0,v,1),btwn(.4,frac(ndivs*v),.6),0);
103  
# Line 114 | Line 120 | if ("$cpict" == "") then
120   else if ("$cpict" == "$picture") then
121          set cpict=
122   endif
123 + if ($decades > 0) then
124 +        set pcargs=($pcargs -e "map(x)=log10(x)/$decades+1")
125 +        set imap="imap(y)=10^((y-1)*$decades)"
126 + else
127 +        set imap="imap(y)=y"
128 + endif
129   pcomb $pcargs -e 'v=(y+.5)/200;vleft=v;vright=v' \
130                  -e 'vbelow=(y-.5)/200;vabove=(y+1.5)/200' \
131                  -e 'ra=0;ga=0;ba=0' -x 100 -y 200 \
132                  > $tempdir/scol.pic
133 < (echo $label; cnt $ndivs |rcalc -e '$1='"($scale)/$ndivs*($ndivs"'-.5-$1)') \
133 > ( echo $label; cnt $ndivs \
134 >                | rcalc -e '$1='"($scale)*imap(($ndivs-.5-"'$1'")/$ndivs)" \
135 >                -e "$imap" ) \
136          | psign -cf 1 1 1 -cb 0 0 0 -h `ev "floor(2*200/$ndivs+.5)"` \
137          | pfilt -1 -x /2 -y /2 > $tempdir/slab.pic
138   pcomb $pcargs $picture $cpict \

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines