<> <> <> <> DIRECTORY Cubic USING [Bezier, Coeffs], Vector USING [Add, Div, In, Max, Min, Mul, Sub, VEC]; CubicImpl: CEDAR PROGRAM IMPORTS Vector EXPORTS Cubic = { OPEN Cubic; VEC: TYPE ~ Vector.VEC; CoeffsToBezier: PUBLIC PROCEDURE[c: Coeffs] RETURNS[b: Bezier] = { OPEN Vector,c,b; <> <> <> <> <> <> <> b0_c0; b1_Add[c0,Div[c1,3]]; b2_Add[b1,Div[Add[c1,c2],3]]; b3_Add[Add[Add[c0,c1],c2],c3]; }; BezierToCoeffs: PUBLIC PROCEDURE[b: Bezier] RETURNS[c: Coeffs] = { OPEN Vector,b,c; temp: VEC; <> <> <> <> <> c0 _ b0; c1 _ Mul[Sub[b1,b0],3]; c2 _ Sub[(temp_Mul[Sub[b2,b1],3]),c1]; c3 _ Sub[b3,Add[b0,temp]]; }; BezierPolygon: PUBLIC PROCEDURE[b: Bezier, epsilon: REAL,Proc: PROCEDURE[Vector.VEC]] = { lvlmax: CARDINAL=10; SubDivide: PROCEDURE[b: Bezier, lvl: CARDINAL] = { IF lvl>=lvlmax OR IsFlat[b,epsilon] THEN Proc[b.b3] ELSE { OPEN Vector; q1,q2,q3,qp1,qp2,q: VEC; bb: Bezier; q1_Mul[Add[b.b0,b.b1],0.5]; q2_Mul[Add[b.b1,b.b2],0.5]; q3_Mul[Add[b.b2,b.b3],0.5]; qp1_Mul[Add[q1,q2],0.5]; qp2_Mul[Add[q2,q3],0.5]; q_Mul[Add[qp1,qp2],0.5]; bb_[b0: b.b0, b1: q1, b2: qp1, b3: q]; SubDivide[bb,lvl+1]; bb_[b0: q, b1: qp2, b2: q3, b3: b.b3]; SubDivide[bb,lvl+1]; }; }; SubDivide[b,1]; }; IsFlat: PROCEDURE[b: Bezier, eps: REAL] RETURNS[BOOLEAN] = INLINE { OPEN Vector; dx,dy: REAL; d1,d2,d,bl,bh: VEC; oh: VEC=[0.5,0.5]; bh_Add[Max[b.b0,b.b3],oh]; bl_Sub[Min[b.b0,b.b3],oh]; IF ~In[b.b1,bl,bh] OR ~In[b.b2,bl,bh] THEN RETURN[FALSE]; d1_Sub[b.b1,b.b0]; d2_Sub[b.b2,b.b0]; d_Sub[b.b3,b.b0]; dx_ABS[d.x]; dy_ABS[d.y]; IF dx+dy < 1 THEN RETURN[TRUE]; IF dy < dx THEN RETURN[ ABS[d2.y-((d.y*d2.x)/d.x)]