6 proofrulethickness 1pt#;
7 makegrid(0pt,0pt for i:=-5pt step 1pt until 5pt: ,i endfor)
8 (0pt,0pt for i:=-5pt step 1pt until 5pt: ,i endfor);
9 proofrulethickness .1pt#;
10 makegrid(0pt,0pt for i:=-4.8pt step .2pt until 4.8pt: ,i endfor)
11 (0pt,0pt for i:=-4.8pt step .2pt until 4.8pt: ,i endfor);
16 tracingequations := tracingonline := 1;
20 def draw_staff(expr first, last, offset)=
21 pickup pencircle scaled stafflinethickness;
22 for i:= first step 1 until last:
23 draw (- staff_space, (i + offset) * staff_space) .. (4 staff_space,( i+ offset)* staff_space);
32 def scaledabout(expr point, scale) =
33 shifted -point scaled scale shifted point
38 % make a local (restored after endgroup) copy of t_var
40 def local_copy(text type, t_var)=
51 % Urgh! Want to do parametric types
54 def del_picture_stack=
55 save save_picture_stack, picture_stack_idx;
58 % better versions of Taupin/Egler savepic cmds
61 def make_picture_stack =
62 % override previous stack.
64 picture save_picture_stack[];
65 numeric picture_stack_idx;
66 picture_stack_idx := 0;
67 def push_picture(expr p) =
68 save_picture_stack[picture_stack_idx] := p ;
69 picture_stack_idx := picture_stack_idx + 1;
71 def pop_picture = save_picture_stack[decr picture_stack_idx] enddef;
72 def top_picture = save_picture_stack[picture_stack_idx] enddef;
77 % why can't I delete individual pens?
81 numeric pen_stack_idx;
83 def push_pen(expr p) =
84 save_pen_stack[pen_stack_idx] := p ;
85 pen_stack_idx := pen_stack_idx +1;
87 def pop_pen = save_pen_stack[decr pen_stack_idx] enddef;
88 def top_pen = save_pen_stack[pen_stack_idx] enddef;
91 save save_pen_stack, pen_stack_idx;
98 def soft_penstroke text t =
99 forsuffixes e = l,r: path_.e:=t; endfor
103 fill path_.l .. tension1.5 .. reverse path_.r .. tension1.5 .. cycle
109 % make a round path segment going from P to Q. 2*A is the angle that the
113 def simple_serif(expr p,q, a)=
114 p{dir(angle(q-p) -a)} .. q{ - dir(angle(p -q) + a)}
120 % err_x: drift of y axis at top
121 % err_y: drift of x axis at right
122 def distorted_ellipse(expr a,b,err_y,err_x,super) =
123 superellipse((a,err_x),(-err_y,b),(-a,-err_x),(err_y,-b),super);
129 % draw an axis aligned block making sure that edges are on pixels.
132 def draw_rounded_block (expr bottom_left, top_right, roundness) =
134 round = floor min(roundness,xpart (top_right-bottom_left),
135 ypart (top_right-bottom_left));
138 pickup pencircle scaled round;
142 z2+(round/2,round/2) = top_right;
143 z4-(round/2,round/2) = bottom_left;
148 fill bot z1 .. rt z1 --- rt z2 .. top z2 ---
149 top z3 .. lft z3 --- lft z4 .. bot z4 --- cycle;
155 def draw_block (expr bottom_left, top_right) =
156 draw_rounded_block (bottom_left, top_right, blot_diameter);
159 def draw_gridline (expr bottom_left,top_right,thickness) =
160 draw_rounded_block (bottom_left-(thickness/2,thickness/2),
161 top_right+(thickness/2,thickness/2),
166 def draw_brush(expr a,w,b,v) =
169 penpos3(w,angle(z2-z1)+90);
170 penpos4(w,angle(z2-z1));
171 penpos5(v,angle(z1-z2)+90);
172 penpos6(v,angle(z1-z2));
176 fill z3r{z3r-z5l}..z4l..{z5r-z3l}z3l..z5r{z5r-z3l}..z6l..{z3r-z5l}z5l..cycle;
182 % make a superellipsoid segment going from FROM to TO, with SUPERNESS.
183 % Take superness = sqrt(2)/2 to get a circle segment
185 % see Knuth, p. 267 and p.126
186 def super_curvelet(expr from, to, superness, dir) =
188 (superness [xpart to, xpart from], superness [ypart from,ypart to]){to - from}
190 (superness [xpart from, xpart to], superness [ypart to,ypart from]){to - from}
196 % Bulb with smooth inside curve.
198 % alpha = start direction.
199 % beta = which side to turn to.
200 % flare = diameter of the bulb
201 % line = diameter of line attachment
202 % direction = is ink on left or right side (1 or -1)
204 def flare_path(expr pos,alpha,beta,line,flare, direction) =
207 penpos1(line,180+beta+alpha);
209 penpos2(flare,180+beta+alpha);
211 penpos3(flare,0+alpha);
212 z3l=z1r+(1/2+0.43)*flare*dir(alpha+beta) ;
215 z4=z2r- line * dir(alpha);
221 p:=z1r{dir(alpha)}..z3r{dir(180+alpha-beta)}..z2l{dir(alpha+180)}
222 ..z3l{dir(180+alpha+beta)}..tension t
223 ..z4{dir(180+alpha+beta)}..z1l{dir(alpha+180)};
235 def brush(expr a,w,b,v) =
243 % Draw a (rest) crook, starting at thickness STEM in point A,
244 % ending a ball W to the left, diameter BALLDIAM
245 % ypart of the center of the ball is BALLDIAM/4 lower than ypart A
247 def balled_crook(expr a, w, balldiam, stem) =
250 penpos1(balldiam/2,-90);
251 penpos2(balldiam/2,0);
252 penpos3(balldiam/2,90);
253 penpos4(balldiam/2,180);
254 x4r=xpart a-w; y3r=ypart a+balldiam/4;
258 x5=x4r+9/8balldiam; y5r=y1r;
260 x6l=xpart a; y6l=ypart a;
261 penstroke z1e..z2e..z3e..z4e..z1e..z5e{right}..z6e;
262 penlabels(1,2,3,4,5,6);
268 currentpicture := currentpicture yscaled -1;
269 set_char_box(charbp, charwd, charht, chardp);
274 currentpicture := currentpicture scaled -1;
275 set_char_box(charwd, charbp, charht, chardp);
280 % center_factor: typically .5, the larger, the larger the radius of the bulb
281 % radius factor: how much the bulb curves inward
283 def draw_bulb(expr turndir, zl, zr, bulb_rad, radius_factor)=
290 % don't get near infinity
291 %z0 = zr + bulb_rad * (zl-zr)/length(zr -zl);
292 z0 = zr + bulb_rad /length(zr -zl) * (zl-zr);
296 z1 = z0 + radius_factor* rad * dir(ang + turndir* 100);
297 z2 = z0 + rad * dir(ang + turndir*300);
299 fill zr{dir (ang + turndir* 90)} .. z1 .. z2 -- cycle;