| 5 |
|
# |
| 6 |
|
set td=/usr/tmp/fc$$ |
| 7 |
|
onintr quit |
| 8 |
< |
set mult=470 |
| 8 |
> |
set mult=179 |
| 9 |
|
set label=Nits |
| 10 |
|
set scale=1000 |
| 11 |
|
set decades=0 |
| 79 |
|
end |
| 80 |
|
mkdir $td |
| 81 |
|
cat > $td/pc0.cal <<_EOF_ |
| 82 |
+ |
PI : 3.14159265358979323846 ; |
| 83 |
|
scale : $scale ; |
| 84 |
|
mult : $mult ; |
| 85 |
|
ndivs : $ndivs ; |
| 87 |
|
or(a,b) : if(a,a,b); |
| 88 |
|
EPS : 1e-7; |
| 89 |
|
neq(a,b) : if(a-b-EPS,1,b-a-EPS); |
| 90 |
< |
btwn(a,b) : if(a-x,-1,b-x); |
| 90 |
> |
btwn(a,x,b) : if(a-x,-1,b-x); |
| 91 |
|
clip(x) : if(x-1,1,if(x,x,0)); |
| 92 |
|
frac(x) : x - floor(x); |
| 93 |
|
boundary(a,b) : neq(floor(ndivs*a),floor(ndivs*b)); |
| 94 |
|
|
| 95 |
< |
isconta = or(boundary(vleft,vright),boundary(vabove,vbelow)); |
| 95 |
> |
isconta = if(btwn(0,v,1),or(boundary(vleft,vright),boundary(vabove,vbelow)),-1); |
| 96 |
|
iscontb = if(btwn(0,v,1),btwn(.4,frac(ndivs*v),.6),-1); |
| 97 |
|
|
| 98 |
|
ro = if(in,clip($redv),ra); |
