DIRECTORY
Atom, Cubic2, CubicPaths, CubicPathsExtras, CubicSplines, GGBasicTypes, GGBoundBox, GGCubic2, GGError, GGInterfaceTypes, GGModelTypes, GGParseIn, GGParseOut, GGSegment, GGSegmentTypes, GGStatistics, GGTransform, GGVector, Imager, ImagerBackdoor, ImagerPath, ImagerTransformation, IO, RealFns, Rope;

GGSegmentImplA: CEDAR PROGRAM
IMPORTS Atom, CubicSplines, CubicPaths, CubicPathsExtras, GGBoundBox, GGCubic2, GGError, GGParseIn, GGParseOut, GGSegment, GGStatistics, GGTransform, GGVector, IO, Imager, ImagerPath, ImagerTransformation, RealFns, Rope
EXPORTS GGSegment = BEGIN

VEC: TYPE = Imager.VEC;
X: NAT = CubicSplines.X;
Y: NAT = CubicSplines.Y;

Bezier: TYPE = Cubic2.Bezier;
BezierRef: TYPE = REF Cubic2.Bezier;
BitVector: TYPE = GGModelTypes.BitVector;
BoundBox: TYPE = GGModelTypes.BoundBox;
BuildPathProc: TYPE = GGSegmentTypes.BuildPathProc;
BuildPathTransformProc: TYPE = GGSegmentTypes.BuildPathTransformProc;
ClassDef: TYPE = REF ClassDefRec;
ClassDefRec: TYPE = GGSegment.ClassDefRec;
Circle: TYPE = GGBasicTypes.Circle;
Edge: TYPE = GGBasicTypes.Edge;
GargoyleData: TYPE = GGInterfaceTypes.GargoyleData;
Line: TYPE = GGBasicTypes.Line;
Point: TYPE = GGBasicTypes.Point;
Segment: TYPE = GGSegmentTypes.Segment;
SegmentClass: TYPE = REF SegmentClassObj;
SegmentClassObj: TYPE = GGSegmentTypes.SegmentClassObj;
SegmentObj: TYPE = GGSegmentTypes.SegmentObj;
SelectedObjectData: TYPE = GGSegmentTypes.SelectedObjectData;
SelectionClass: TYPE = GGSegmentTypes.SelectionClass;
Vector: TYPE = GGBasicTypes.Vector;

defaultStrokeWidth: REAL _ 2.0;
segmentClasses: LIST OF ClassDef;

Problem: ERROR [msg: Rope.ROPE] = GGError.Problem;


RegisterSegmentClass: PUBLIC PROC [classDef: ClassDef] = {
segmentClasses _ CONS[classDef, segmentClasses];
};

FetchSegmentClass: PUBLIC PROC [type: ATOM] RETURNS [class: SegmentClass] = {
FOR l: LIST OF ClassDef _ segmentClasses, l.rest UNTIL l=NIL DO
IF l.first.type=type THEN RETURN[l.first.class];
ENDLOOP;
ERROR Problem[msg: IO.PutFR["SegmentClass %g is unknown", [rope[Atom.GetPName[type]]] ]];
};

CopySegment: PUBLIC PROC [seg: Segment] RETURNS [copy: Segment] = {
copyData: REF ANY;
copyData _ seg.class.copyData[seg];
copy _ NEW[SegmentObj _ [
class: seg.class,
looks: seg.looks,
strokeWidth: seg.strokeWidth,
color: seg.color,
lo: seg.lo, hi: seg.hi,
TselectedInFull: seg.TselectedInFull,
bBox: GGBoundBox.CopyBoundBox[seg.bBox],
data: copyData,
props: seg.props
]];
};

ReverseSegment: PUBLIC PROC [seg: Segment] = {
temp: Point;
temp _ seg.lo;
seg.lo _ seg.hi;
seg.hi _ temp;
seg.class.reverse[seg];
};

OpenUpSegment: PUBLIC PROC [seg: Segment]  = {
IF seg.class.type=$CubicSpline THEN {
data: CubicSplineData _ NARROW[seg.data];
IF data.type=cyclicAL THEN {
data.type _ naturalAL;
data.path _ CubicPaths.PathFromCubic[CubicSplines.MakeSpline[data.cps, data.type]]; -- remake path data
};
};
};

TransformSegment: PUBLIC PROC [seg: Segment, transform: ImagerTransformation.Transformation] =  {
seg.lo _ GGTransform.Transform[transform, seg.lo];
seg.hi _ GGTransform.Transform[transform, seg.hi];
seg.class.transform[seg, transform];
};

TranslateSegment: PUBLIC PROC [seg: Segment, vector: Vector] =  {
transform: ImagerTransformation.Transformation;
transform _ ImagerTransformation.Translate[[vector.x, vector.y]];
TransformSegment[seg, transform];
};

MoveEndPointSegment: PUBLIC PROC [seg: Segment, lo: BOOL, newPoint: Point] = {
IF lo THEN seg.lo _ newPoint
ELSE seg.hi _ newPoint;
seg.class.endPointMoved[seg, lo, newPoint];
};

ConicData: TYPE = REF ConicDataRec;
ConicDataRec: TYPE = RECORD [
p1: VEC,
s: REAL,
p1Selected: SelectedObjectData,
path: CubicPaths.Path
];

BuildConicClass: PROC [] RETURNS [class: SegmentClass] = {
OPEN GGSegment;
class _ NEW[SegmentClassObj _ [
type: $Conic,
copyData: ConicCopyData,
reverse: ConicReverse,
boundBox: ConicBoundBox,
tightBox: ConicTightBox,
transform: ConicTransform,
endPointMoved: ConicEndPointMoved,
controlPointMoved: ConicControlPointMoved,
buildPath: ConicBuildPath,
buildPathTransform: ConicBuildPathTransform,
controlPointGet: ConicControlPointGet,
controlPointCount: ConicControlPointCount,
controlPointFieldSet: ConicFieldSet,
controlPointFieldGet: ConicFieldGet,
closestPoint: CurveClosestPoint,
closestControlPoint: ConicClosestControlPoint,
closestPointAndTangent: NoOpClosestPointAndTangent,
lineIntersection: NoOpLineIntersection,
circleIntersection: NoOpCircleIntersection,
asSimpleCurve: NoOpAsSimpleCurve,
addJoint: GGSegment.NoOpAddJoint,
describe: ConicDescribe,
fileOut: ConicFileOut,
fileIn: ConicFileIn,
setStrokeWidth: ConicSetStrokeWidth
]];
};

ConicSetStrokeWidth: PROC [seg: Segment, strokeWidth: REAL] = {
seg.strokeWidth _ strokeWidth;
UpdateConicBoundBox[seg];
};

MakeConic: PUBLIC PROC [p0, p1, p2: Point, r: REAL, props: LIST OF REF ANY] RETURNS [seg: Segment] = {
data: ConicData _ NEW[ConicDataRec _ [
p1: [p1.x,p1.y], s: r, p1Selected: [FALSE, FALSE, FALSE],
path: PathFromConic[p0, p1, p2, r]
]];
seg _ NEW[SegmentObj _ [
class: GGSegment.FetchSegmentClass[$Conic],
looks: NIL,
strokeWidth: defaultStrokeWidth,
color: Imager.black,
lo: p0, hi: p2,
bBox: GGBoundBox.CreateBoundBox[0,0,0,0], -- gets updated one line down.
data: data,
props: props
]];
UpdateConicBoundBox[seg];
};

ConicBoundBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ seg.bBox;
};

ConicTightBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ GGBoundBox.NullBoundBox[];
UpdateConicTightBox[seg, bBox];
};

UpdateConicTightBox: PROC [seg: Segment, bBox: BoundBox] = {
data: ConicData _ NARROW[seg.data];
CurveBoundBox[seg, bBox];
};

UpdateConicBoundBox: PROC [seg: Segment] = {
data: ConicData _ NARROW[seg.data];
cpHalf: REAL _ GGModelTypes.halfJointSize + 1;
strokeWidth: REAL _ seg.strokeWidth;
pad: REAL _ MAX[cpHalf, strokeWidth];
UpdateConicTightBox[seg, seg.bBox];
GGBoundBox.EnlargeByPoint[seg.bBox, [data.p1.x, data.p1.y] ]; -- CP is far from curve
GGBoundBox.EnlargeByOffset[seg.bBox, pad];
};

ConicCopyData: GGSegmentTypes.CopyDataProc = {
conicData: ConicData _ NARROW[seg.data];
new: ConicData _ NEW[ConicDataRec];
new.p1 _ conicData.p1;
new.s _ conicData.s;
new.p1Selected _ conicData.p1Selected;
new.path _ CubicPathsExtras.CopyPath[conicData.path];
RETURN[new];
};

ConicReverse: PROC [seg: Segment] = {
data: ConicData _ NARROW[seg.data];
CubicPathsExtras.ReversePath[data.path];
};

ConicBuildPath: PROC [seg: Segment, lineTo: ImagerPath.LineToProc, curveTo: ImagerPath.CurveToProc, conicTo: ImagerPath.ConicToProc, arcTo: ImagerPath.ArcToProc] = {
data: ConicData _ NARROW[seg.data];
conicTo[data.p1, seg.hi, data.s];
};

ConicBuildPathTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation, entire, lo, hi: BOOL, controlPoints: BitVector, lineTo: ImagerPath.LineToProc, curveTo: ImagerPath.CurveToProc, conicTo: ImagerPath.ConicToProc, arcTo: ImagerPath.ArcToProc] = {
data: ConicData _ NARROW[seg.data];
p0,p1,p2: VEC;
IF entire THEN {
[p0, p2] _ TransformEndPoints[seg.lo, seg.hi, TRUE, TRUE, transform];
p1 _ ImagerTransformation.Transform[transform, data.p1];
}
ELSE {
[p0, p2] _ TransformEndPoints[seg.lo, seg.hi, lo, hi, transform];
p1 _ IF controlPoints[0] THEN ImagerTransformation.Transform[transform, data.p1]
ELSE data.p1;
};
conicTo[p1,p2,data.s];
};
ConicTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation] = {
data: ConicData _ NARROW[seg.data];
data.p1 _ ImagerTransformation.Transform[transform, data.p1];
CubicPaths.TransformPath[data.path, transform];
UpdateConicBoundBox[seg];
};

ConicEndPointMoved: PROC [seg: Segment, lo: BOOL, newPoint: Point] = {
data: ConicData _ NARROW[seg.data];
p0: Point _ IF lo THEN newPoint ELSE seg.lo;
p1: Point _ [data.p1.x, data.p1.y];
p2: Point _ IF lo THEN seg.hi ELSE newPoint;
data.path _ PathFromConic[p0, p1, p2, data.s];
UpdateConicBoundBox[seg];
};

ConicControlPointMoved: PROC [seg: Segment, transform: ImagerTransformation.Transformation, controlPointNum: NAT] = {
p1Vec: VEC;
p1: Point;
data: ConicData _ NARROW[seg.data];
IF controlPointNum#0 THEN ERROR;
p1Vec _ ImagerTransformation.Transform[transform, data.p1];
p1 _ [p1Vec.x, p1Vec.y];
data.path _ PathFromConic[seg.lo, p1, seg.hi, data.s];
data.p1 _ p1Vec;
UpdateConicBoundBox[seg];
};

TransformEndPoints: PROC [loPt, hiPt: Point, lo, hi: BOOL, transform: ImagerTransformation.Transformation] RETURNS [newLo, newHi: VEC] = {
IF lo THEN newLo _ ImagerTransformation.Transform[transform, [loPt.x,loPt.y]]
ELSE newLo _ [loPt.x,loPt.y];
IF hi THEN newHi _ ImagerTransformation.Transform[transform, [hiPt.x,hiPt.y]]
ELSE newHi _ [hiPt.x,hiPt.y];
};

ConicDescribe: PROC [seg: Segment, self, lo, hi: BOOL, cps: BitVector] RETURNS [rope: Rope.ROPE] = {
rope _ "Parabolic Conic";
};

ConicFileOut: PROC [seg: Segment, f: IO.STREAM] = {
data: ConicData _ NARROW[seg.data];
p1: Point _ [data.p1.x, data.p1.y];
GGParseOut.WritePoint[f, p1];
f.PutF[" %g", [real[data.s]]];
};

ConicFileIn: PROC [f: IO.STREAM, loPoint, hiPoint: Point, version: REAL] RETURNS [seg: Segment] = {
s: REAL;
p1: Point;
p1 _ GGParseIn.ReadPoint[f];
s _ GGParseIn.ReadBlankAndReal[f];
seg _ MakeConic[loPoint, p1, hiPoint, s, NIL];
};

ConicFieldSet: GGSegmentTypes.ControlPointFieldSetProc = {
data: ConicData _ NARROW[seg.data];
IF controlPointNum # 0 THEN ERROR;
SELECT selectClass FROM
normal => data.p1Selected.normal _ selected;
hot => data.p1Selected.hot _ selected;
active => data.p1Selected.active _ selected;
ENDCASE => ERROR;
};

ConicFieldGet: GGSegmentTypes.ControlPointFieldGetProc = {
data: ConicData _ NARROW[seg.data];
IF controlPointNum # 0 THEN ERROR;
SELECT selectClass FROM
normal => selected _ data.p1Selected.normal;
hot => selected _ data.p1Selected.hot;
active => selected _ data.p1Selected.active;
ENDCASE => ERROR;
};

ConicControlPointGet: PROC [seg: Segment, controlPointNum: NAT] RETURNS [point: Point] = {
data: ConicData _ NARROW[seg.data];
IF controlPointNum#0 THEN ERROR;
RETURN[ [data.p1.x, data.p1.y] ];
};

ConicControlPointCount: PROC [seg: Segment] RETURNS [controlPointCount: NAT] = {
data: ConicData _ NARROW[seg.data]; -- just for NarrowRefFault if there is a problem
RETURN[1];
};
ConicClosestControlPoint: PROC [seg: Segment, testPoint: Point, tolerance: REAL] RETURNS [point: Point, controlPointNum: NAT, success: BOOL] = {
data: ConicData _ NARROW[seg.data];
success _ TRUE;
controlPointNum _ 0;
point _ [data.p1.x, data.p1.y];
};

BezierData: TYPE = REF BezierDataRec;
BezierDataRec: TYPE = RECORD [
b1Selected, b2Selected: SelectedObjectData,
path: CubicPaths.Path];

BuildBezierClass: PROC [] RETURNS [class: SegmentClass] = {
OPEN GGSegment;
class _ NEW[SegmentClassObj _ [
type: $Bezier,
copyData: BZCopyData,
reverse: BZReverse,
boundBox: BZBoundBox,
tightBox: BZTightBox,
transform: BZTransform,
endPointMoved: BZEndPointMoved,
controlPointMoved: BZControlPointMoved,
buildPath: CurveBuildPath,
buildPathTransform: BZBuildPathTransform,
controlPointGet: BZControlPointGet,
controlPointCount: BZControlPointCount,
controlPointFieldSet: BZFieldSet,
controlPointFieldGet: BZFieldGet,
closestPoint: CurveClosestPoint,
closestControlPoint: BZClosestControlPoint,
closestPointAndTangent: GGSegment.NoOpClosestPointAndTangent,
lineIntersection: GGSegment.NoOpLineIntersection,
circleIntersection: GGSegment.NoOpCircleIntersection,
asSimpleCurve: NoOpAsSimpleCurve,
addJoint: BZAddJoint,
describe: BZDescribe,
fileOut: BZFileOut,
fileIn: BZFileIn,
setStrokeWidth: BZSetStrokeWidth
]];
};

BZSetStrokeWidth: PROC [seg: Segment, strokeWidth: REAL] = {
seg.strokeWidth _ strokeWidth;
UpdateBZBoundBox[seg];
};

MakeBezier: PUBLIC PROC [p0, p1, p2, p3: Point, props: LIST OF REF ANY] RETURNS [seg: Segment] = {
data: BezierData _ NEW[BezierDataRec _ [
path: NEW[CubicPaths.PathRec]
]];
data.path.cubics _ NEW[CubicPaths.PathSequence[1]];
data.path.cubics[0] _ [p0, p1, p2, p3];
CubicPathsExtras.UpdateBounds[data.path];
seg _ NEW[SegmentObj _ [
class: FetchSegmentClass[$Bezier],
looks: NIL,
strokeWidth: defaultStrokeWidth,
color: Imager.black,
lo: p0, hi: p3,
bBox: GGBoundBox.CreateBoundBox[0,0,0,0], -- gets updated one line down.
data: data,
props: props
]];
UpdateBZBoundBox[seg];
};

BZBoundBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ seg.bBox;
};

BZTightBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ GGBoundBox.NullBoundBox[];
UpdateBZTightBox[seg, bBox];
};

UpdateBZBoundBox: PROC [seg: Segment] = {
cpHalf: REAL _ GGModelTypes.halfJointSize + 1;
strokeWidth: REAL _ seg.strokeWidth;
pad: REAL _ MAX[cpHalf, strokeWidth];
data: BezierData _ NARROW[seg.data];
cubics: REF CubicPaths.PathSequence _ data.path.cubics;
seg.bBox^ _ [
loX: MIN[cubics[0].b0.x, cubics[0].b1.x, cubics[0].b2.x, cubics[0].b3.x],
loY: MIN[cubics[0].b0.y, cubics[0].b1.y, cubics[0].b2.y, cubics[0].b3.y],
hiX: MAX[cubics[0].b0.x, cubics[0].b1.x, cubics[0].b2.x, cubics[0].b3.x],
hiY: MAX[cubics[0].b0.y, cubics[0].b1.y, cubics[0].b2.y, cubics[0].b3.y],
null: FALSE, infinite: FALSE];
GGBoundBox.EnlargeByOffset[seg.bBox, pad];
};

UpdateBZTightBox: PROC [seg: Segment, bBox: BoundBox] = {
CurveBoundBox[seg, bBox];
};

BZCopyData: GGSegmentTypes.CopyDataProc = {
bzData: BezierData _ NARROW[seg.data];
new: BezierData _ NEW[BezierDataRec];
new.b1Selected _ bzData.b1Selected;
new.b2Selected _ bzData.b2Selected;
new.path _ CubicPathsExtras.CopyPath[bzData.path];
RETURN[new];
};

BZReverse: PROC [seg: Segment] = {
data: BezierData _ NARROW[seg.data];
CubicPathsExtras.ReversePath[data.path];
};
BZTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation] = {
data: BezierData _ NARROW[seg.data];
CubicPaths.TransformPath[data.path, transform];
UpdateBZBoundBox[seg];
};

BZEndPointMoved: PROC [seg: Segment, lo: BOOL, newPoint: Point] = {
data: BezierData _ NARROW[seg.data];
IF lo THEN data.path.cubics[0].b0 _ [newPoint.x, newPoint.y]
ELSE data.path.cubics[0].b3 _ [newPoint.x, newPoint.y];
UpdateBZBoundBox[seg];
};

BZControlPointMoved: PROC [seg: Segment, transform: ImagerTransformation.Transformation, controlPointNum: NAT] = {
data: BezierData _ NARROW[seg.data];
SELECT controlPointNum FROM
0 => data.path.cubics[0].b1 _ ImagerTransformation.Transform[transform, data.path.cubics[0].b1];
1 => data.path.cubics[0].b2 _ ImagerTransformation.Transform[transform, data.path.cubics[0].b2];
ENDCASE => ERROR;
UpdateBZBoundBox[seg];
};

BZBuildPathTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation, entire, lo, hi: BOOL, controlPoints: BitVector, lineTo: ImagerPath.LineToProc, curveTo: ImagerPath.CurveToProc, conicTo: ImagerPath.ConicToProc, arcTo: ImagerPath.ArcToProc] = {
path: CubicPaths.Path _ GetPath[seg];
moveTo: ImagerPath.MoveToProc _ {};
IF entire THEN {
CubicPaths.TransformPath[path, transform];
CubicPaths.EnumeratePath[path, moveTo, curveTo];
CubicPaths.TransformPath[path, ImagerTransformation.Invert[transform]];
}
ELSE {
oldLo, oldHi, oldC1, oldC2: VEC;
[oldLo, oldHi, oldC1, oldC2] _ BZTransformPath[path, transform, lo, hi, controlPoints];
CubicPaths.EnumeratePath[path, moveTo, curveTo];
BZRestorePath[path, oldLo, oldHi, oldC1, oldC2];
};
};
BZTransformPath: PROC [path: CubicPaths.Path, transform: ImagerTransformation.Transformation, lo, hi: BOOL, controlPoints: BitVector] RETURNS [oldLo, oldHi, oldC1, oldC2: VEC] ~ {
oldLo _ path.cubics[0].b0;
oldHi _ path.cubics[0].b3;
oldC1 _ path.cubics[0].b1;
oldC2 _ path.cubics[0].b2;
IF lo THEN path.cubics[0].b0 _ ImagerTransformation.Transform[transform, oldLo];
IF hi THEN path.cubics[0].b3 _ ImagerTransformation.Transform[transform, oldHi];
IF controlPoints[0] THEN path.cubics[0].b1 _ ImagerTransformation.Transform[transform, path.cubics[0].b1];
IF controlPoints[1] THEN path.cubics[0].b2 _ ImagerTransformation.Transform[transform, path.cubics[0].b2];
};

BZRestorePath: PROC [path: CubicPaths.Path, oldLo, oldHi, oldC1, oldC2: VEC] ~ {
path.cubics[0].b0 _ oldLo;
path.cubics[0].b3 _ oldHi;
path.cubics[0].b1 _ oldC1;
path.cubics[0].b2 _ oldC2;
};

BZControlPointGet: PROC [seg: Segment, controlPointNum: NAT] RETURNS [point: Point] = {
data: BezierData _ NARROW[seg.data];
SELECT controlPointNum FROM
0 => RETURN [data.path.cubics[0].b1];
1 => RETURN [data.path.cubics[0].b2];
ENDCASE => ERROR;
};

BZControlPointCount: PROC [seg: Segment] RETURNS [controlPointCount: NAT] = {
data: BezierData _ NARROW[seg.data]; -- just for NarrowRefFault if there is a problem
RETURN[2];
};

BZFieldSet: GGSegmentTypes.ControlPointFieldSetProc = {
data: BezierData _ NARROW[seg.data];
SELECT controlPointNum FROM
0 => {
SELECT selectClass FROM
normal => data.b1Selected.normal _ selected;
hot => data.b1Selected.hot _ selected;
active => data.b1Selected.active _ selected;
ENDCASE => ERROR;
};
1 => {
SELECT selectClass FROM
normal => data.b2Selected.normal _ selected;
hot => data.b2Selected.hot _ selected;
active => data.b2Selected.active _ selected;
ENDCASE => ERROR;
};
ENDCASE => ERROR;
};

BZFieldGet: GGSegmentTypes.ControlPointFieldGetProc = {
data: BezierData _ NARROW[seg.data];
SELECT controlPointNum FROM
0 => {
SELECT selectClass FROM
normal => selected _ data.b1Selected.normal;
hot => selected _ data.b1Selected.hot;
active => selected _ data.b1Selected.active;
ENDCASE => ERROR;
};
1 => {
SELECT selectClass FROM
normal => selected _ data.b2Selected.normal;
hot => selected _ data.b2Selected.hot;
active => selected _ data.b2Selected.active;
ENDCASE => ERROR;
};
ENDCASE => ERROR;
};
BZClosestControlPoint: PROC [seg: Segment, testPoint: Point, tolerance: REAL] RETURNS [point: Point, controlPointNum: NAT, success: BOOL] = {
data: BezierData _ NARROW[seg.data];
d0: REAL _ GGVector.DistanceSquared[testPoint, data.path.cubics[0].b1];
d1: REAL _ GGVector.DistanceSquared[testPoint, data.path.cubics[0].b2];
IF d0 < d1 THEN {
d0 _ RealFns.SqRt[d0];
point _ data.path.cubics[0].b1;
controlPointNum _ 0;
success _ (d0 <= tolerance);
}
ELSE {
d1 _ RealFns.SqRt[d1];
point _ data.path.cubics[0].b2;
controlPointNum _ 1;
success _ (d1 <= tolerance);
};
};

BZAddJoint: PUBLIC PROC [seg: Segment, pos: Point] RETURNS [seg1, seg2: Segment] = {
data: BezierData _ NARROW[seg.data];
bezier, bezier1, bezier2: Cubic2.Bezier;
alpha: REAL;
bezier _ data.path.cubics[0];
alpha _ GGCubic2.GetParam[bezier, pos];
[bezier1, bezier2] _ GGCubic2.AlphaSplit[bezier, alpha];
seg1 _ MakeBezier[bezier1.b0, bezier1.b1, bezier1.b2, bezier1.b3, NIL];
seg2 _ MakeBezier[bezier2.b0, bezier2.b1, bezier2.b2, bezier2.b3, NIL];
};

BZDescribe: PROC [seg: Segment, self, lo, hi: BOOL, cps: BitVector] RETURNS [rope: Rope.ROPE] = {
rope _ "Bezier";
};

BZFileOut: PROC [seg: Segment, f: IO.STREAM] = {
data: BezierData _ NARROW[seg.data];
GGParseOut.WritePoint[f, data.path.cubics[0].b1];
f.PutChar[IO.SP];
GGParseOut.WritePoint[f, data.path.cubics[0].b2];
};

BZFileIn: PROC [f: IO.STREAM, loPoint, hiPoint: Point, version: REAL] RETURNS [seg: Segment] = {
p1, p2: Point;
p1 _ GGParseIn.ReadPoint[f];
p2 _ GGParseIn.ReadPoint[f];
seg _ MakeBezier[loPoint, p1, p2, hiPoint, NIL];
};


CubicSplineData: TYPE = REF CubicSplineDataRec;
CubicSplineDataRec: TYPE = RECORD [
cps: CubicSplines.KnotSequence,
type: CubicSplines.SplineType,
path: CubicPaths.Path,
cpsSelected: SelectionSequence];

SelectionSequence: TYPE = REF SelectionSequenceObj;
SelectionSequenceObj: TYPE = RECORD [
seq: SEQUENCE len: NAT OF SelectedObjectData
];

BuildCubicSplineClass: PROC [] RETURNS [class: SegmentClass] = {
OPEN GGSegment;
class _ NEW[SegmentClassObj _ [
type: $CubicSpline,
copyData: CSCopyData,
reverse: CSReverse,
boundBox: CSBoundBox,
tightBox: CSTightBox,
transform: CSTransform,
endPointMoved: CSEndPointMoved,
controlPointMoved: CSControlPointMoved,
buildPath: CurveBuildPath,
buildPathTransform: CSBuildPathTransform,
controlPointGet: CSControlPointGet,
controlPointCount: CSControlPointCount,
controlPointFieldSet: CSFieldSet,
controlPointFieldGet: CSFieldGet,
closestPoint: CurveClosestPoint,
closestControlPoint: CSClosestControlPoint,
closestPointAndTangent: NoOpClosestPointAndTangent,
lineIntersection: NoOpLineIntersection,
circleIntersection: NoOpCircleIntersection,
asSimpleCurve: NoOpAsSimpleCurve,
describe: CSDescribe,
fileOut: CSFileOut,
fileIn: CSFileIn,
addJoint: CSAddJoint,
setStrokeWidth: CSSetStrokeWidth
]];
};

CSSetStrokeWidth: PROC [seg: Segment, strokeWidth: REAL] = {
seg.strokeWidth _ strokeWidth;
UpdateCSBoundBox[seg];
};

MakeCubicSpline: PUBLIC PROC [cps: CubicSplines.KnotSequence, type: CubicSplines.SplineType, props: LIST OF REF ANY] RETURNS[seg: Segment] = {
coeffs: CubicSplines.CoeffsSequence;
data: CubicSplineData;
newCPs: CubicSplines.KnotSequence;
last: NAT _ cps.length-1;
SELECT type FROM -- determine correct number of control points
naturalAL, cyclicAL => newCPs _ NEW[CubicSplines.KnotSequenceRec[cps.length]];
bspline => newCPs _ NEW[CubicSplines.KnotSequenceRec[cps.length+4]];
ENDCASE => ERROR;
data _ NEW[CubicSplineDataRec _ [
cps: newCPs, type: type,
cpsSelected: NEW[SelectionSequenceObj[cps.length]]
]];
seg _ NEW[SegmentObj _ [
class: FetchSegmentClass[$CubicSpline],
looks: NIL,
strokeWidth: defaultStrokeWidth,
color: Imager.black,
lo: [cps[0][X], cps[0][Y]],
hi: [cps[last][X], cps[last][Y]],
bBox: GGBoundBox.CreateBoundBox[0,0,0,0], -- gets updated one line down.
data: data,
props: props
]];
CSEndPointSet[seg: seg, lo: TRUE, point: FPToPoint[cps[0]]];
CSEndPointSet[seg: seg, lo: FALSE, point: FPToPoint[cps[last]]];
FOR i:INT IN [0..seg.class.controlPointCount[seg]) DO
CSControlPointSet[seg, i, FPToPoint[cps[i+1]]];
ENDLOOP;
coeffs _ CubicSplines.MakeSpline[data.cps, type];
data.path _ CubicPaths.PathFromCubic[coeffs]; -- build a path of beziers and compute bboxes.
UpdateCSBoundBox[seg: seg, boundsOK: TRUE];
};

CSBoundBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ seg.bBox;
};

CSTightBox: PROC [seg: Segment] RETURNS [bBox: BoundBox] = {
bBox _ GGBoundBox.NullBoundBox[];
CurveBoundBox[seg, bBox];
};

UpdateCSBoundBox: PROC [seg: Segment, boundsOK: BOOL _ FALSE] = {
data: CubicSplineData _ NARROW[seg.data];
cpHalf: REAL _ GGModelTypes.halfJointSize + 1;
strokeWidth: REAL _ seg.strokeWidth;
pad: REAL _ MAX[cpHalf, strokeWidth];
IF data.type = bspline THEN KnotBoundBox[seg]
ELSE {
CurveBoundBox[seg, seg.bBox, boundsOK];
GGBoundBox.EnlargeByOffset[seg.bBox, pad];
};
};

CSCopyData: GGSegmentTypes.CopyDataProc = {
csData: CubicSplineData _ NARROW[seg.data];
new: CubicSplineData _ NEW[CubicSplineDataRec];
cpCount: NAT _ seg.class.controlPointCount[seg];
new.cps _ NEW[CubicSplines.KnotSequenceRec[csData.cps.length]];
new.cpsSelected _ NEW[SelectionSequenceObj[cpCount]];
FOR i: NAT IN [0..new.cps.length) DO
new.cps[i] _ csData.cps[i];
ENDLOOP;
FOR i: NAT IN [0..cpCount) DO
new.cpsSelected[i] _ csData.cpsSelected[i];
ENDLOOP;
new.type _ csData.type;
new.path _ CubicPathsExtras.CopyPath[csData.path];
RETURN[new];
};

CSReverse: PROC [seg: Segment] = {
data: CubicSplineData _ NARROW[seg.data];
temp: CubicSplines.FPCoords;
last: NAT _ data.cps.length-1;
FOR i: NAT IN [0..data.cps.length/2) DO
temp _ data.cps[i];
data.cps[i] _ data.cps[last-i];
data.cps[last-i] _ temp;
ENDLOOP;
CubicPathsExtras.ReversePath[data.path];
};

CSBuildPathTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation, entire, lo, hi: BOOL, controlPoints: BitVector, lineTo: ImagerPath.LineToProc, curveTo: ImagerPath.CurveToProc, conicTo: ImagerPath.ConicToProc, arcTo: ImagerPath.ArcToProc] = {
data: CubicSplineData _ NARROW[seg.data];
moveTo: ImagerPath.MoveToProc _ {};
IF entire THEN {
CubicPaths.TransformPath[data.path, transform];
CubicPaths.EnumeratePath[data.path, moveTo, curveTo];
CubicPaths.TransformPath[data.path, ImagerTransformation.Invert[transform]];
}
ELSE {
coeffs: CubicSplines.CoeffsSequence;
cpsCopy: CubicSplines.KnotSequence _ NEW[CubicSplines.KnotSequenceRec[data.cps.length]];	-- save copy of cps
FOR i:INT IN [0..data.cps.length) DO
cpsCopy[i] _ data.cps[i];
ENDLOOP;
FOR i:INT IN [0..seg.class.controlPointCount[seg]) DO
IF controlPoints[i] THEN CSControlPointSet[seg, i, ImagerTransformation.Transform[transform, CSControlPointGet[seg, i]]];
ENDLOOP;
IF lo THEN CSEndPointSet[seg: seg, lo: TRUE, point: ImagerTransformation.Transform[transform, seg.lo]];
IF hi THEN CSEndPointSet[seg: seg, lo: FALSE, point: ImagerTransformation.Transform[transform, seg.hi]];
coeffs _ CubicSplines.MakeSpline[data.cps, data.type];
data.path _ CubicPaths.PathFromCubic[coeffs];
CubicPaths.EnumeratePath[data.path, moveTo, curveTo];
data.cps _ cpsCopy;
coeffs _ CubicSplines.MakeSpline[data.cps, data.type];    -- restore old spline data
data.path _ CubicPaths.PathFromCubic[coeffs];
};
};

CSTransform: PROC [seg: Segment, transform: ImagerTransformation.Transformation] = {
data: CubicSplineData _ NARROW[seg.data];
FOR i: NAT IN [0..data.cps.length) DO
vec: VEC _ [data.cps[i][X], data.cps[i][Y]];
vec _ ImagerTransformation.Transform[transform, vec];
data.cps[i][X] _ vec.x;
data.cps[i][Y] _ vec.y;
ENDLOOP;
CubicPaths.TransformPath[data.path, transform];
UpdateCSBoundBox[seg];
};

CSEndPointMoved: PROC [seg: Segment, lo: BOOL, newPoint: Point] = {
data: CubicSplineData _ NARROW[seg.data];
coeffs: CubicSplines.CoeffsSequence;
CSEndPointSet[seg, lo, newPoint];
coeffs _ CubicSplines.MakeSpline[data.cps, data.type];  -- must recompute coeffs & path
data.path _ CubicPaths.PathFromCubic[coeffs];
UpdateCSBoundBox[seg];
};

CSControlPointMoved: PROC [seg: Segment, transform: ImagerTransformation.Transformation, controlPointNum: NAT] = {
data: CubicSplineData _ NARROW[seg.data];
coeffs: CubicSplines.CoeffsSequence;
CSControlPointSet[seg, controlPointNum, ImagerTransformation.Transform[transform, CSControlPointGet[seg, controlPointNum]]];
coeffs _ CubicSplines.MakeSpline[data.cps, data.type];
data.path _ CubicPaths.PathFromCubic[coeffs];
UpdateCSBoundBox[seg];
};

CSDescribe: PROC [seg: Segment, self, lo, hi: BOOL, cps: BitVector] RETURNS [rope: Rope.ROPE] = {
data: CubicSplineData _ NARROW[seg.data];
SELECT data.type FROM
naturalAL => rope _ "Natural";
cyclicAL => rope _ "Cyclic Natural";
bspline => rope _ "BSpline";
ENDCASE => ERROR;
};

CSFileOut: PROC [seg: Segment, f: IO.STREAM] = {
data: CubicSplineData _ NARROW[seg.data];
cpCount: NAT _ seg.class.controlPointCount[seg];
SELECT data.type FROM
naturalAL => f.PutF["Type: Natural "];
cyclicAL => f.PutF["Type: CyclicNatural "];
bspline => f.PutF["Type: BSpline "];
ENDCASE => ERROR;
f.PutF["%g ", [integer[cpCount]] ];
FOR i: INT IN [0..cpCount) DO
GGParseOut.WritePoint[f, CSControlPointGet[seg, i]];
f.PutChar[' ];
ENDLOOP;
};

CSFileIn: PROC [f: IO.STREAM, loPoint, hiPoint: Point, version: REAL] RETURNS [seg: Segment] = {
cps: CubicSplines.KnotSequence;
type: CubicSplines.SplineType;
typeWord: Rope.ROPE;
cyclic, ok: BOOL;
IF version >= 8608.12 THEN {
GGParseIn.ReadBlankAndRope[f, "Type:"];
typeWord _ GGParseIn.ReadBlankAndWord[f];
IF Rope.Equal[typeWord, "Natural"] THEN type _ naturalAL
ELSE IF Rope.Equal[typeWord, "CyclicNatural"] THEN type _ cyclicAL
ELSE IF Rope.Equal[typeWord, "BSpline"] THEN type _ bspline
ELSE ERROR;
}
ELSE IF version >= 8607.22 THEN {
GGParseIn.ReadBlankAndRope[f, "Cyclic:"];
[truth: cyclic, good: ok] _ GGParseIn.ReadBOOL[f, version];
IF NOT ok THEN ERROR; -- Bad format in input file.  Expecting BOOL.
IF cyclic THEN type _ cyclicAL ELSE type _ naturalAL;
}
ELSE type _ naturalAL;
cps _ NEW[CubicSplines.KnotSequenceRec[GGParseIn.ReadBlankAndNAT[f]+2]];
FOR i: INT IN [1..cps.length-1) DO
cps[i] _ PointToFP[GGParseIn.ReadPoint[f]];
ENDLOOP;
IF type = cyclicAL AND (loPoint # hiPoint) THEN ERROR; -- both joints of cyclic spline should have the same position
cps[0] _ PointToFP[loPoint];
cps[cps.length-1] _ PointToFP[hiPoint];
seg _ MakeCubicSpline[cps, type, NIL];
};

CSFieldSet: GGSegmentTypes.ControlPointFieldSetProc = {
data: CubicSplineData _ NARROW[seg.data];
SELECT selectClass FROM
normal => data.cpsSelected[controlPointNum].normal _ selected;
hot => data.cpsSelected[controlPointNum].hot _ selected;
active => data.cpsSelected[controlPointNum].active _ selected;
ENDCASE => ERROR;
};

CSFieldGet: GGSegmentTypes.ControlPointFieldGetProc = {
data: CubicSplineData _ NARROW[seg.data];
SELECT selectClass FROM
normal => selected _ data.cpsSelected[controlPointNum].normal;
hot => selected _ data.cpsSelected[controlPointNum].hot;
active => selected _ data.cpsSelected[controlPointNum].active;
ENDCASE => ERROR;
};

CSControlPointGet: PROC [seg: Segment, controlPointNum: NAT] RETURNS [point: Point] = {
data: CubicSplineData _ NARROW[seg.data];
SELECT data.type FROM
naturalAL, cyclicAL => {
IF controlPointNum >= data.cps.length-2 THEN ERROR;
RETURN [FPToPoint[data.cps[controlPointNum+1]]];
};
bspline => {
IF controlPointNum >= data.cps.length-4 THEN ERROR;
RETURN [FPToPoint[data.cps[controlPointNum+3]]];
};
ENDCASE => ERROR;
};

CSControlPointCount: PROC [seg: Segment] RETURNS [controlPointCount: NAT] = {
data: CubicSplineData _ NARROW[seg.data]; -- just for NarrowRefFault if there is a problem
SELECT data.type FROM
naturalAL, cyclicAL => RETURN[data.cps.length-2]; -- 2 true control points are joints
bspline => RETURN[data.cps.length-6] -- 6 true control points are joints
ENDCASE => ERROR;
};
CSClosestControlPoint: PROC [seg: Segment, testPoint: Point, tolerance: REAL] RETURNS [point: Point, controlPointNum: NAT, success: BOOL] = {
data: CubicSplineData _ NARROW[seg.data];
nextDist: REAL;
minDist: REAL;
IF seg.class.controlPointCount[seg] = 0 THEN RETURN[[0,0], 999, FALSE];
minDist _ GGVector.DistanceSquared[testPoint, CSControlPointGet[seg, 0]];
controlPointNum _ 0;
FOR i:INT IN [1..seg.class.controlPointCount[seg]) DO
IF (nextDist _ GGVector.DistanceSquared[testPoint, CSControlPointGet[seg, i]]) < minDist THEN {
minDist _ nextDist;
controlPointNum _ i; 
};
ENDLOOP;
minDist _ RealFns.SqRt[minDist];
success _ (minDist <= tolerance);
point _ CSControlPointGet[seg, controlPointNum];
};

CSAddJoint: PROC [seg: Segment, pos: Point] RETURNS [seg1, seg2: Segment] = {
data: CubicSplineData _ NARROW[seg.data];
data1, data2: CubicSplineData;
oldCPCount: INT _ seg.class.controlPointCount[seg];
nearest: INT _ FindNearestPath[data.path, pos];
cpsInSeg1: INT _ PathNumToCPIndex[seg, nearest]; -- The number of cps that will be in seg1
cpsInSeg2: INT _ oldCPCount - cpsInSeg1;
coeffs: CubicSplines.CoeffsSequence;
seg1 _ CopySegment[seg];
seg2 _ CopySegment[seg];
data1 _ NARROW[seg1.data];
data2 _ NARROW[seg2.data];
IF data.type=cyclicAL THEN data1.type _ data2.type _ naturalAL; --cyclic splines are changed
data1.cps _ NEW[CubicSplines.KnotSequenceRec[data.cps.length-cpsInSeg2]];
data1.cpsSelected _ NEW[SelectionSequenceObj[cpsInSeg1]];
CSEndPointSet[seg1, TRUE, seg.lo];
CSEndPointSet[seg1, FALSE, pos];
seg1.hi _ pos;
FOR i:INT IN [0..cpsInSeg1) DO
CSControlPointSet[seg1, i, CSControlPointGet[seg, i]];
data1.cpsSelected[i] _ data.cpsSelected[i];
ENDLOOP;
coeffs _ CubicSplines.MakeSpline[data1.cps, data1.type];
data1.path _ CubicPaths.PathFromCubic[coeffs];
data2.cps _ NEW[CubicSplines.KnotSequenceRec[data.cps.length-cpsInSeg1]];
data2.cpsSelected _ NEW[SelectionSequenceObj[cpsInSeg2]];
CSEndPointSet[seg2, TRUE, pos];
CSEndPointSet[seg2, FALSE, seg.hi];
seg2.lo _ pos;
FOR i:INT IN [cpsInSeg1..oldCPCount) DO
CSControlPointSet[seg2, i-cpsInSeg1, CSControlPointGet[seg, i]];
data2.cpsSelected[i-cpsInSeg1] _ data.cpsSelected[i];
ENDLOOP;
coeffs _ CubicSplines.MakeSpline[data2.cps, data2.type];
data2.path _ CubicPaths.PathFromCubic[coeffs];
UpdateCSBoundBox[seg1]; -- added by KAP. December 22, 1986
UpdateCSBoundBox[seg2]; -- added by KAP. December 22, 1986
};

CSControlPointAdd: PUBLIC PROC [seg: Segment, pos: Point, gargoyleData: GargoyleData] RETURNS [Segment] = {
data: CubicSplineData _ NARROW[seg.data];
oldCPCount: INT _ seg.class.controlPointCount[seg];
nearest: INT _ FindNearestPath[data.path, pos];
newCPIndex: INT _ PathNumToCPIndex[seg, nearest]; -- newCPIndex is the index of the CP in the old segment that the new CP will precede
newSeg: Segment _ CopySegment[seg];
newSegData: CubicSplineData _ NARROW[newSeg.data];
coeffs: CubicSplines.CoeffsSequence;
newSegData.cps _ NEW[CubicSplines.KnotSequenceRec[data.cps.length+1]];
newSegData.cpsSelected _ NEW[SelectionSequenceObj[oldCPCount+1]];
CSEndPointSet[newSeg, TRUE, seg.lo];
CSEndPointSet[newSeg, FALSE, seg.hi];
FOR i:INT IN [0..newCPIndex) DO
CSControlPointSet[newSeg, i, CSControlPointGet[seg, i]];
newSegData.cpsSelected[i] _ data.cpsSelected[i];
ENDLOOP;
CSControlPointSet[newSeg, newCPIndex, pos];
CSFieldSet[newSeg, newCPIndex, TRUE, normal]; --newly created cp will be selected
CSFieldSet[newSeg, newCPIndex, FALSE, hot];
CSFieldSet[newSeg, newCPIndex, FALSE, active];
FOR i:INT IN [newCPIndex..oldCPCount) DO
CSControlPointSet[newSeg, i+1, CSControlPointGet[seg, i]];
newSegData.cpsSelected[i+1] _ data.cpsSelected[i];
ENDLOOP;
coeffs _ CubicSplines.MakeSpline[newSegData.cps, newSegData.type];
newSegData.path _ CubicPaths.PathFromCubic[coeffs];
UpdateCSBoundBox[newSeg]; -- added by KAP. December 22, 1986
RETURN [newSeg];
};

CSControlPointDelete: PUBLIC PROC [seg: Segment, cpVec: BitVector, gargoyleData: GargoyleData] RETURNS [Segment] = {
CountSelected: PROC [bitVec: BitVector] RETURNS [selectedCount: INT _ 0] = {
FOR i:INT IN [0..bitVec.len) DO
IF bitVec[i] THEN selectedCount _ selectedCount + 1;
ENDLOOP;
};
data: CubicSplineData _ NARROW[seg.data];
currentCPIndex: INT _ 0;
oldCPCount: INT _ seg.class.controlPointCount[seg];
cpSelectedCount: INT _ CountSelected[cpVec];
newSeg: Segment _ CopySegment[seg];
coeffs: CubicSplines.CoeffsSequence;
newSegData: CubicSplineData _ NARROW[newSeg.data];
newSegData.cps _ NEW[CubicSplines.KnotSequenceRec[data.cps.length-cpSelectedCount]];
newSegData.cpsSelected _ NEW[SelectionSequenceObj[data.cps.length-cpSelectedCount]];
CSEndPointSet[newSeg, TRUE, seg.lo];
CSEndPointSet[newSeg, FALSE, seg.hi];
FOR i:INT IN [0..oldCPCount) DO
IF NOT cpVec[i] THEN {
CSControlPointSet[newSeg, currentCPIndex, CSControlPointGet[seg, i]];
newSegData.cpsSelected[currentCPIndex] _ data.cpsSelected[i];
currentCPIndex _ currentCPIndex + 1;
};
ENDLOOP;
coeffs _ CubicSplines.MakeSpline[newSegData.cps, newSegData.type];
newSegData.path _ CubicPaths.PathFromCubic[coeffs];
UpdateCSBoundBox[newSeg]; -- added by KAP. December 22, 1986
RETURN [newSeg];
};
CSControlPointSet: PROC [seg: Segment, controlPointNum: NAT, point: Point] = {
data: CubicSplineData _ NARROW[seg.data];
SELECT data.type FROM
naturalAL, cyclicAL => data.cps[controlPointNum+1] _ PointToFP[point];
bspline => data.cps[controlPointNum+3] _ PointToFP[point];
ENDCASE => ERROR;
};

CSEndPointSet: PROC [seg: Segment, lo: BOOL, point: Point] = {
data: CubicSplineData _ NARROW[seg.data];
SELECT data.type FROM
naturalAL => {
IF lo THEN data.cps[0] _ PointToFP[point]
ELSE data.cps[data.cps.length-1] _ PointToFP[point];
};
cyclicAL => data.cps[0] _ data.cps[data.cps.length-1] _ PointToFP[point];
bspline => {
IF lo THEN data.cps[0] _ data.cps[1] _ data.cps[2] _ PointToFP[point]
ELSE data.cps[data.cps.length-1] _ data.cps[data.cps.length-2] _ data.cps[data.cps.length-3] _ PointToFP[point];
};
ENDCASE => ERROR;
};

CurveBoundBox: PROC [seg: Segment, bBox: BoundBox, boundsOK: BOOL _ FALSE] = {
path: CubicPaths.Path _ GetPath[seg];
r: Imager.Rectangle;
IF NOT boundsOK THEN CubicPathsExtras.UpdateBounds[path]; -- the bbox of the curve itself
r _ path.bounds;
GGBoundBox.UpdateBoundBox[bBox, r.x, r.y, r.x+r.w, r.y+r.h];
};
CurveMaskStroke: PROC [seg: Segment, dc: Imager.Context] = {
path: CubicPaths.Path _ GetPath[seg];
pathProc: ImagerPath.PathProc = {CubicPaths.EnumeratePath[path, moveTo, curveTo]};
Imager.SetStrokeEnd[dc, round];
Imager.MaskStroke[dc, pathProc];
};
CurveBuildPath: PROC [seg: Segment, lineTo: ImagerPath.LineToProc, curveTo: ImagerPath.CurveToProc, conicTo: ImagerPath.ConicToProc, arcTo: ImagerPath.ArcToProc] = {
path: CubicPaths.Path _ GetPath[seg];
moveTo: ImagerPath.MoveToProc _ {};
CubicPaths.EnumeratePath[path, moveTo, curveTo];
};

useBBox: BOOL _ TRUE;
CurveClosestPoint: PROC [seg: Segment, testPoint: Point, tolerance: REAL] RETURNS [point: Point, success: BOOL] = {
path: CubicPaths.Path;
IF useBBox THEN {
bigBox: GGBasicTypes.BoundBoxObj _ [seg.bBox.loX-tolerance, seg.bBox.loY-tolerance, seg.bBox.hiX+tolerance, seg.bBox.hiY+tolerance, FALSE, FALSE];
IF NOT PointIsInBox[testPoint, bigBox] THEN RETURN[[0,0], FALSE];
};
GGStatistics.StartInterval[$CurveClosestPoint, GGStatistics.GlobalTable[]];
path _ GetPath[seg];
[point, success] _ GGCubic2.ClosestPointAnalytic[testPoint, path, tolerance];
GGStatistics.StopInterval[$CurveClosestPoint, GGStatistics.GlobalTable[]];
};

PointIsInBox: PROC [test: Point, box: GGBasicTypes.BoundBoxObj] RETURNS [BOOL] = {
RETURN[ NOT (test.x < box.loX OR test.x > box.hiX OR test.y < box.loY OR test.y > box.hiY) ];
};

KnotBoundBox: PROC [seg: Segment] = {
data: CubicSplineData _ NARROW[seg.data];
cps: CubicSplines.KnotSequence _ data.cps;
minX, maxX: REAL _ cps[0][X];
minY, maxY: REAL _ cps[0][Y];
FOR i:INT IN [1..cps.length) DO
IF maxX < cps[i][X] THEN maxX _ cps[i][X];
IF minX > cps[i][X] THEN minX _ cps[i][X];
IF maxY < cps[i][Y] THEN maxY _ cps[i][Y];
IF minY > cps[i][Y] THEN minY _ cps[i][Y];
ENDLOOP;
seg.bBox^ _ [loX: minX, loY: minY, hiX: maxX, hiY: maxY, null: FALSE, infinite: FALSE];
GGBoundBox.AllowForJoints[seg.bBox];
};

FPToPoint: PROC [FP: CubicSplines.FPCoords] RETURNS [point: Point] = {
RETURN [[FP[X], FP[Y]]];
};

PointToFP: PROC [point: Point] RETURNS [CubicSplines.FPCoords] = {
RETURN [[point.x, point.y]];
};

FindNearestPath: PROC [path: CubicPaths.Path, pos: Point] RETURNS [INT] = {
miniPath: CubicPaths.Path _ NEW[CubicPaths.PathRec];
miniPath.cubics _ NEW[CubicPaths.PathSequence[1]];
FOR i:INT IN [0..path.cubics.length) DO
miniPath.cubics[0] _ path.cubics[i];
CubicPathsExtras.UpdateBounds[miniPath];
IF CubicPaths.PointOnPath[pos, miniPath, 0.5] THEN RETURN [i];
ENDLOOP;
RETURN [-1];
};

PathNumToCPIndex: PROC [seg: Segment, pathNum: INT] RETURNS [cpIndex: INT] = {
data: CubicSplineData _ NARROW[seg.data];
cpCount: INT _ seg.class.controlPointCount[seg];
SELECT data.type FROM
cyclicAL, naturalAL => cpIndex _ pathNum;
bspline => cpIndex _ pathNum-1;
ENDCASE => ERROR;
IF cpIndex < 0 THEN cpIndex _ 0
ELSE IF cpIndex > cpCount THEN cpIndex _ cpCount;
RETURN;
};

GetPath: PROC [seg: Segment] RETURNS [path: CubicPaths.Path] ~ {
WITH seg.data SELECT FROM
type: BezierData => path _ type.path;
type: ConicData => path _ type.path;
type: CubicSplineData => path _ type.path;
ENDCASE => ERROR;
RETURN[path];
};

PathFromConic: PROC [p0, p1, p2: Point, r: REAL] RETURNS [path: CubicPaths.Path] ~ {
bList: LIST OF BezierRef;
joint: VEC _ [p0.x,p0.y];
bproc: ImagerPath.CurveToProc = {
new: BezierRef _ NEW[Cubic2.Bezier _ [joint, p1, p2, p3]];
bList _ CONS[new, bList];
joint _ p3;
};
ImagerPath.ConicToCurves[[p0.x,p0.y], [p1.x,p1.y], [p2.x,p2.y], r, bproc];
path _ PathFromList[bList];
RETURN[path];
};

PathFromList: PROC [list: LIST OF BezierRef] RETURNS [CubicPaths.Path] ~ {
count: NAT _ 0;
path: CubicPaths.Path _ NEW[CubicPaths.PathRec];
FOR l: LIST OF BezierRef _ list, l.rest UNTIL l=NIL DO
count _ count+1;
ENDLOOP;
path.cubics _ NEW[CubicPaths.PathSequence[count]];
FOR l: LIST OF BezierRef _ list, l.rest UNTIL l=NIL DO
count _ count-1;	--count starts at number of entries, not last entry
path.cubics[count] _ l.first^;
ENDLOOP;
CubicPathsExtras.UpdateBounds[path];
RETURN[path];
};

Init: PROC [] = {
classDef: ClassDef _ NEW[ClassDefRec _ [type: $Bezier, class: BuildBezierClass[]]];
RegisterSegmentClass[classDef];
classDef _ NEW[ClassDefRec _ [type: $CubicSpline, class: BuildCubicSplineClass[]]];
RegisterSegmentClass[classDef];
classDef _ NEW[ClassDefRec _ [type: $Conic, class: BuildConicClass[]]];
RegisterSegmentClass[classDef];
};

InitStats: PROC [] = {
interval: GGStatistics.Interval;
interval _ GGStatistics.CreateInterval[$CurveClosestPoint];
GGStatistics.AddInterval[interval, GGStatistics.GlobalTable[]];
};

Init[];
InitStats[];

END.
μGGSegmentImplA.mesa
Copyright c 1985, 1986 by Xerox Corporation.  All rights reserved.
Last edited by Pier on December 22, 1986 2:36:13 pm PST
Contents:  Procedures which implement the Gargoyle segment classes.
Bier, January 26, 1987 4:00:36 pm PST
David Kurlander August 28, 1986 8:02:58 pm PDT

Generic Segment Class Procedures
Changes all of the cyclic segments in the traj to non-cyclic segs
A convenience routine which does a TransformSegment, where transform is a simple translation.

The Conic Class
Transforming
Drawing
Control Point Access
Hit Testing
Editing
Textual Description
Editing
Let m be the midpoint of [p0, p2].
Let p be the point on [m, p1] such that Length[m, p]/Length[m, p1] = r.
The curve starts at p0, passes through p, and ends at p2.
It is a line if r=0; an ellipse if 0<r<1/2; a parabola if r=1/2; a hyperbola if 1/2<r<1.
The curve is bounded by the triangle [p0, p1, p2].
GGSegmentTypes.BoundBoxProc
GGSegmentTypes.BoundBoxProc
Updates bound box which bounds this segment (including stroke width and control points).
GGSegmentTypes.BoundBoxProc
GGSegmentTypes.CopyDataProc
Copies the data in the "data" field of seg.  This data is class-dependent.
GGSegmentTypes.ReverseProc
The "lo" end and "hi" end have switched roles.  Update data structures as necessary.
Conic Drawing

GGSegmentTypes.BuildPathProc
GGSegmentTypes.BuildPathTransformProc
Like BuildPathProc but you should assume that transform has been applied to your lower joint iff lo = TRUE and to your higher joint iff hi = TRUE.  This is for rubberbanding filled areas.

Conic Transforming

GGSegmentTypes.TransformProc
Apply the given transformation to all internal data of the segment.  It is now in a new position, orientation, scaling, or skewing.
GGSegmentTypes.EndPointMovedProc
GGSegmentTypes.ControlPointMovedProc

Conic Description
GGSegmentTypes.FileOutProc
GGSegmentTypes.FileInProc
Conic Parts

[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]
[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]
GGSegmentTypes.ControlPointGetProc
GGSegmentTypes.ControlPointCountProc

Conic Hit Testing

GGSegmentTypes.ClosestControlPointProc

The Bezier Class
Transformations
Drawing
Control Point Access
Hit Testing
Editing
Textual Description
Editing
Create a Bezier segment which passed thru p0 and p3, guided by p1 and p2.
GGSegmentTypes.BoundBoxProc.
Find the bounding hull of the control points.  This will bound the curve as well.
Creates a new box which bounds this Bezier, and does NOT allow for stroke width and control point size.
GGSegmentTypes.CopyDataProc
Copies the data in the "data" field of seg.  This data is class-dependent.
GGSegmentTypes.ReverseProc
The "lo" end and "hi" end have switched roles.  Update data structures as necessary.

Bezier Transforming

GGSegmentTypes.TransformProc
Apply the given transformation to all internal data of the segment.  It is now in a new position, orientation, scaling, or skewing.
GGSegmentTypes.EndPointMovedProc
GGSegmentTypes.ControlPointMovedProc
Bezier Drawing

GGSegmentTypes.BuildPathTransformProc
Like BuildPathProc but you should assume that transform has been applied to your lower joint iff lo = TRUE and to your higher joint iff hi = TRUE.  This is for rubberbanding filled areas.

Bezier Control Point Access

GGSegmentTypes.ControlPointGetProc
GGSegmentTypes.ControlPointCountProc
[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]
[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]

Beier Hit Testing

GGSegmentTypes.ClosestControlPointProc
-- Consider only the middle two control points -- other two are treated as joints
Beier Editing

GGSegmentTypes.AddJointProc
Beier Description

GGSegmentTypes.FileOutProc
GGSegmentTypes.FileInProc
The Cubic Spline Class
Transformations
Drawing
Control Point Access
Hit Testing
Textual Description
Editing
Makes a segment of the types defined in CubicSplines.  cps should be a sequence of knots (including endpoints).  We'll tweek the data to be in the correct form.
copy in new cp data
Now copy in new cp data. Must do this last so that the set routines work.
Fit bounding box
The bound box which bounds this segment, NOT including stroke width and control points.
GGSegmentTypes.BoundBoxProc
GGSegmentTypes.CopyDataProc
Copies the data in the "data" field of seg.  This data is class-dependent.
GGSegmentTypes.ReverseProc
The "lo" end and "hi" end have switched roles.  Update data structures as necessary.
CS Drawing

GGSegmentTypes.BuildPathTransformProc
Like BuildPathProc but you should assume that transform has been applied to your lower joint iff lo = TRUE and to your higher joint iff hi = TRUE.  This is for rubberbanding filled areas.
GGSegmentTypes.TransformProc
Apply the given transformation to all internal data of the segment.  It is now in a new position, orientation, scaling, or skewing.
GGSegmentTypes.EndPointMovedProc
GGSegmentTypes.ControlPointMovedProc
CS Description

GGSegmentTypes.FileInProc
GGSegmentTypes.FileInProc
CS Parts

[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]
[seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]
CS Part Generators

GGSegmentTypes.ControlPointGetProc
GGSegmentTypes.ControlPointCountProc

CS Hit Testing

GGSegmentTypes.ClosestControlPointProc
-- Consider only the inner control points (edgepoints are treated as joints)
CS Editing

Returns the CubicSpline segments which are formed by placing a joint in the specified position.  Note that this position must lie on the spline, since we use this to determine cp ordering
--create and fill in the new CP array of seg1
--create and fill in the new CP array of seg2
Returns a CubicSpline segment with a control point added.  Note that the specified position must lie on the spline, since control point ordering is determined from this.
create and fill in the new CP array
Returns a CubicSpline Segment with selected control points deleted.

Curve Auxiliary Procedures

GGSegmentTypes.BoundBoxProc
Creates a new box which bounds this segment (NOT including stroke width).

GGSegmentTypes.MaskStrokeProc
Draw yourself into dc as a stroke.  Assume that strokeWidth, and color are already set.

GGSegmentTypes.BuildPathProc
Assume that the Imager's current point is at your lower joint.  Call the procedures given to draw yourself (moveTo is left out -- you shouldn't need it).

GGSegmentTypes.ClosestPointProc
Used for hit testing.  Find the nearest point on seg to testPoint.
GGSegmentTypes.BoundBoxProc
Creates a boundbox of the segment's knots.  Segment must be a cubic spline.
Auxiliary Path Procs

Finds which on which subpath a point.  Returns -1 if none.
--Given a subpath number, figures out which control point "divides" this subpath from the next.
Put the cubics in the sequence.  CONS put them on the list "backwards"
Init Procs

Κ0(˜
codeš
œ™Kš
œB™BKš
œ7™7KšΟnœ;™CKš
œ%™%Kšœ+Οk™.—K™š
ž	˜	Kšœ˜žœ˜ͺ—K˜šœžœ
ž˜Kšžœ™
žœ9˜Ϋ
Kšžœ
ž˜—˜Kšžœžœ
žœ
˜Kš
œžœ˜Kš
œžœ˜K˜Kšœžœ˜Kšœžœžœ˜$Kšœžœ˜)Kšœ
žœ˜'Kšœžœ ˜3Kšœžœ)˜EKšœ
žœžœ
˜!Kšœ
žœ˜*Kšœžœ˜#Kšœžœ˜Kšœžœ!˜3Kšœžœ˜Kšœžœ˜!Kšœ	žœ˜'Kšœžœžœ˜)Kšœžœ"˜7Kšœžœ˜-Kšœžœ%˜=Kšœžœ!˜5Kšœžœ˜#—˜Kšœžœ˜Kšœžœ
žœ
˜!—K˜Kšœžœžœ˜2K˜š
œ ™ K˜—šœžœ
žœ˜:Kšœžœ˜0Kš
œ˜K˜—š
œžœ
žœžœžœ˜Mšžœžœ
žœ#žœžœ
ž˜?Kšžœžœ
žœ˜0Kšžœ
˜—KšžœžœD˜YK˜K˜—šœžœ
žœžœ˜CKšœ
žœ
žœ
˜Kš
œ#˜#šœžœ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ%˜%Kš
œ(˜(Kš
œ˜Kš
œ˜Kš
œ˜—K˜K˜—šœžœ
žœ˜.K˜K˜K˜K˜K˜K˜K˜—š
œžœ
žœ˜.K™Ašžœžœ˜%Jšœžœ˜)šžœžœ˜Jš
œ˜JšœTΟc˜gJ˜—J˜—J˜J˜—šœžœ
žœD˜aKš
œ2˜2Kš
œ2˜2Kš
œ$˜$K˜K˜—šœžœ
žœ$˜AKš
œ]™]K˜/Kš
œA˜AKš
œ!˜!K˜K˜—šœžœ
žœžœ˜NKšžœžœ˜Kšžœ˜Kš
œ+˜+K˜—K™š
Οb™K˜—Kšœžœžœ˜#šœžœžœ˜Kšœžœ
˜Kšœžœ
˜Kš
œ˜Kš
œ˜Kš
œ˜K˜—šœžœžœ˜:KšΠbk ˜šœžœ˜K˜
K˜K˜Kš
œ˜˜Kš
 ™—K˜K˜"š
œ*˜*Kš
 ™—Kš
œ˜š
œ,˜,Kš
 ™—Kš
œ&˜&Kš
œ*˜*K˜$˜$Kš
 ™—K˜ Kš
œ.˜.Kš
œ3˜3Jš
œ'˜'Jš
œ+˜+š
œ!˜!Kš
 ™—˜!Kš
 ™—K˜Kš
œ˜˜K™—Kš
œ#˜#K˜—K˜K˜—šœžœžœ˜?Kš
œ˜Kš
œ˜K˜K˜—š	œžœ
žœžœ	žœ
žœ
žœ
žœžœ˜fKš
œ"™"Kš
œG™GKš
œ9™9Kš
œX™XKš
œ2™2šœžœ˜&Kšœ$žœžœžœ˜9Kš
œ"˜"K˜—šœžœ˜K˜+Kšœžœ
˜Kš
œ ˜ K˜Kš
œ˜Kšœ*Ÿ˜HKš
œ˜Kš
œ˜Kš
œ˜—Kš
œ˜K˜K˜—š
œžœžœ˜?Kš
œ™Kš
œ˜Kš
œ˜K˜—š
œžœžœ˜?Kš
œ™Kš
œ!˜!Kš
œ˜Kš
œ˜K˜—šœžœ#˜<K™XKšœžœ˜#Kš
œ˜K˜K˜—šœžœ˜,Kš
œ™Kšœžœ˜#Kšœžœ"˜.Kšœ
žœ˜$Kšœžœžœ˜%Kš
œ#˜#Kšœ>Ÿ˜UKš
œ*˜*Kš
œ˜K˜—š
œ!˜.Kš
œ™K™JKšœžœ˜(Kšœžœ˜#K˜K˜Kš
œ&˜&K˜5Kšžœ˜K˜K˜—šœžœ˜%Kš
œ™K™TKšœžœ˜#K˜(K˜K˜—š
œ
™
K™—šœžœ‘
˜₯
Kš
œ™Kšœžœ˜#Kš
œ!˜!K˜K˜—šœžœPžœ
˜ŽKš
œ%™%Kšœfžœ#žœ*™»
Kšœžœ˜#Kšœ
žœ
˜šžœžœ˜Kšœ.žœžœ
˜EKš
œ8˜8K˜
—šžœ˜Kš
œA˜AKšœžœžœ3˜PKšžœ	˜
K˜—K˜K˜K™—š
œ™K™—šœžœC˜WKš
œ™K™ƒ
Kšœžœ˜#Kš
œ=˜=Kš
œ/˜/Kš
œ˜K˜K˜—šœžœžœ˜FKš
œ ™ Kšœžœ˜#Kšœžœžœ
žœ˜,Kš
œ#˜#Kšœžœžœžœ
˜,Kš
œ.˜.Kš
œ˜Kš
œ˜K˜—šœžœQžœ˜uKš
œ$™$Kšœžœ
˜Kš
œ
˜
Kšœžœ˜#Kšžœžœ
žœ
˜ Kš
œ;˜;Kš
œ˜Kš
œ6˜6Kš
œ˜Kš
œ˜Kš
œ˜K˜—š
œžœžœ2žœžœ˜Š
KšžœžœC˜MKšžœ˜KšžœžœC˜MKšžœ˜K˜K˜K™—Kš
œ™š

œžœžœžœ
žœ˜dKš
œ˜K˜K˜—šœžœžœ
žœ˜3Kš
œ™Kšœžœ˜#Kš
œ#˜#Kš
œ˜Kš
œ˜K˜K˜—•StartOfExpansionm -- [f: STREAM, loPoint: GGBasicTypes.Point, hiPoint: GGBasicTypes.Point] RETURNS [seg: GGModelTypes.Segment]šœžœžœ
žœ$žœžœ˜cKš
œ™Kšœžœ
˜Kš
œ
˜
Kš
œ˜K˜"Kšœ)žœ˜.K˜K˜—š
œ™K™—–s -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]š
œ-˜:Kšœ/žœžœ-™oKšœžœ˜#Kšžœžœ
žœ
˜"šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜K˜—–| -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]š
œ-˜:Kšœ/žœ.žœžœ
™xKšœžœ˜#Kšžœžœ
žœ
˜"šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜K˜—šœžœ!žœžœ˜ZKš
œ"™"Kšœžœ˜#Kšžœžœ
žœ
˜ Kšžœ˜!Kš
œ˜K˜—šœžœžœžœ˜PJš
œ$™$KšœžœŸ0˜TKšžœ˜
Kš
œ˜K™—š
œ™K™—šœžœ-žœžœ!žœžœ˜
Kš
œ&™&Kšœžœ˜#Kšœ
žœ
˜Kš
œ˜Kš
œ˜K˜—K™š
 ™K˜—Kšœžœžœ˜%šœžœžœ˜K˜+Kš
œ˜K˜—šœžœžœ˜;Kš‘ ˜šœžœ˜K˜K˜K˜Kš
œ˜˜Kš
 ™—K˜K˜š
œ'˜'Kš
 ™—Kš
œ˜š
œ)˜)Kš
 ™—Kš
œ#˜#Kš
œ'˜'Kš
œ!˜!˜!Kš
 ™—Kš
œ ˜ Kš
œ+˜+Kš
œ=˜=Jš
œ1˜1Jš
œ5˜5š
œ!˜!Kš
 ™—˜Kš
 ™—Kš
œ˜K˜˜Kš
 ™—Kš
œ ˜ K˜—K˜K˜—šœžœžœ˜<Kš
œ˜Kš
œ˜K˜K˜—š
œžœ
žœ žœ
žœ
žœ
žœžœ˜bK™Išœžœ˜(Kšœžœ˜K˜—Kšœžœ˜3Kš
œ'˜'Kš
œ)˜)šœžœ˜K˜"Kšœžœ
˜Kš
œ ˜ K˜Kš
œ˜Kšœ*Ÿ˜HKš
œ˜Kš
œ˜Kš
œ˜—K˜K˜K˜—š
œžœžœ˜<Kš
œ˜K˜K˜—š
œžœžœ˜<Kš
œ!˜!Kš
œ˜K˜K˜—šœžœ˜)Kš
œ™K™QKšœžœ"˜.Kšœ
žœ˜$Kšœžœžœ˜%Kšœžœ˜$Kšœžœ,˜7š
œ
˜
KšœžœA˜IKšœžœA˜IKšœžœA˜IKšœžœA˜IKšœžœžœ˜—Kš
œ*˜*Kš
œ˜K˜—šœžœ#˜9Kšœ5žœ/™gKš
œ˜K˜K˜—š
œ!˜+Kš
œ™K™JKšœžœ˜&Kšœžœ˜%Kš
œ#˜#Kš
œ#˜#Kš
œ2˜2Kšžœ˜K˜K˜—š	œžœ˜"Kš
œ™K™TKšœžœ˜$K˜(K˜K™—š
œ™K™—šœžœC˜TKš
œ™K™ƒ
Kšœžœ˜$Kš
œ/˜/Kš
œ˜K˜K˜—šœžœžœ˜CKš
œ ™ Kšœžœ˜$Kšžœžœ2˜<Kšžœ3˜7Kš
œ˜K˜K˜—šœžœQžœ˜rKš
œ$™$Kšœžœ˜$šžœž˜Kš
œ`˜`Kš
œ`˜`Kšžœžœ
˜—Kš
œ˜Kš
œ˜K˜—š
œ™K™—šœžœPžœ
˜‹Kš
œ%™%Kšœfžœ#žœ*™»
Kš
œ%˜%Kš
œ#˜#šžœžœ˜Kš
œ*˜*K˜0Kš
œG˜GK˜
—šžœ˜Kšœžœ
˜ Kš
œW˜WKš
œ0˜0Kš
œ0˜0Kš
œ˜—Kš
œ˜K™—š
œžœQžœžœžœ˜³
Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜KšžœžœF˜PKšžœžœF˜PKšžœžœR˜jKšžœžœR˜jKš
œ˜K˜—š
œžœ5žœ˜PKš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜K˜K˜—š
œ™K™—šœžœ!žœžœ˜Wš
œ"™"Kšœžœ˜$šžœž˜Kšœžœ˜%Kšœžœ˜%Kšžœžœ
˜—Kš
œ˜—K˜—šœžœžœžœ˜MKš
œ$™$KšœžœŸ0˜UKšžœ˜
Kš
œ˜K˜—–s -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]š
œ-˜7Kšœ/žœžœ-™oKšœžœ˜$Kšžœž˜˜šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜—˜šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜—Kšžœžœ
˜K˜K˜—–| -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]š
œ-˜7Kšœ/žœ.žœžœ
™xKšœžœ˜$Kšžœž˜˜šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜—˜šžœ
ž˜Kš
œ,˜,Kš
œ&˜&Kš
œ,˜,Kšžœžœ
˜—K˜—Kšžœžœ
˜K˜K™—š
œ™K™—šœžœ-žœžœ!žœžœ˜
Kš
œ&™&Kšœžœ˜$JšŸ1œ ™QKšœžœ?˜GKšœžœ?˜Gšžœ	žœ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ
˜
—šžœ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ˜—Kš
œ˜K˜—š
œ
™
K™—š
œžœ
žœžœ˜TKš
œ™Kšœžœ˜$K˜(Kšœžœ
˜Kš
œ˜Kš
œ'˜'Kš
œ8˜8KšœBžœ˜GKšœBžœ˜GK˜K˜—š
œ™K™—š

œžœžœžœ
žœ˜aKš
œ˜K˜K˜—š	œžœžœ
žœ˜0Kš
œ™Kšœžœ˜$Kš
œ1˜1Kšœ
žœ
žœ˜Kš
œ1˜1K˜K˜—šœžœžœ
žœ$žœžœ˜`Kš
œ™K˜K˜K˜Kšœ+žœ˜0K˜—K˜š
 ™K˜—Kšœžœžœ˜/šœžœžœ˜#Kš
œ˜Kš
œ˜Kš
œ˜Kš
œ ˜ K˜—K–s -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selected: BOOL, selectClass: GGSegmentTypes.SelectionClass]šœžœžœ˜3šœžœžœ˜%Kšœžœžœ
žœ˜,K˜K˜—šœžœžœ˜@Kš‘ ˜šœžœ˜K˜K˜K˜K˜˜Kš
 ™—K˜K˜š
œ'˜'Kš
 ™—Kš
œ˜š
œ)˜)Kš
 ™—Kš
œ#˜#Kš
œ'˜'Kš
œ!˜!š
œ!˜!Kš
 ™—K˜ Kš
œ+˜+Kš
œ3˜3Jš
œ'˜'Jš
œ+˜+š
œ!˜!Kš
 ™—K˜Kš
œ˜˜Kš
 ™—K˜Kš
œ ˜ K˜—K˜K˜—šœžœžœ˜<Kš
œ˜Kš
œ˜K˜K˜—šœžœ
žœHžœ
žœ
žœ
žœžœ˜Ž
K™ 
Kš
œ$˜$Kš
œ˜Kš
œ"˜"Kšœžœ˜šžœžœ
Ÿ-˜>Kšœ žœ+˜NKšœžœ-˜DKšžœžœ
˜—K™šœžœ˜!Kš
œ˜Kšœ
žœ"˜2K˜—šœžœ˜K˜'Kšœžœ
˜Kš
œ ˜ K˜K˜K˜!Kšœ*Ÿ˜HKš
œ˜Kš
œ˜Kš
œ˜—K™IKšœžœ˜<Kšœžœ˜@šžœžœ
žœ'ž˜5K˜/Kšžœ
˜—K˜1Kšœ.Ÿ.˜\K™Kšœ%žœ˜+K˜K˜—š
œžœžœ˜<Kš
œ˜K˜K˜—š
œžœžœ˜<K™WKš
œ!˜!Kš
œ˜K˜K˜—šœžœžœžœ˜AKš
œ™Kšœžœ˜)Kšœžœ"˜.Kšœ
žœ˜$Kšœžœžœ˜%Kšžœžœ˜-šžœ˜Kš
œ'˜'Kš
œ*˜*K˜—K˜K˜—š
œ!˜+Kš
œ™K™JKšœžœ˜+Kšœžœ˜/Kšœ	žœ$˜0Kšœ
žœ2˜?Kšœžœ ˜5šžœžœ
žœž˜$Kš
œ˜Kšžœ
˜—šžœžœ
žœž˜Kš
œ+˜+Kšžœ
˜—Kš
œ˜Kš
œ2˜2Kšžœ˜K˜K˜—š	œžœ˜"Kš
œ™K™TKšœžœ˜)K˜Kšœžœ˜šžœžœ
žœž˜'K˜K˜K˜Kšžœ
˜—K˜(K˜K˜—š
œ
™
K™—šœžœPžœ
˜‹Kš
œ%™%Kšœfžœ#žœ*™»
Kšœžœ˜)Kš
œ#˜#šžœžœ˜Kš
œ/˜/K˜5Kš
œL˜LK˜
—šžœ˜Kš
œ$˜$Kšœ%žœ1Ÿ˜lšžœžœ
žœž˜$K˜Kšžœ
˜—šžœžœ
žœ'ž˜5Kšžœžœa˜y—Kšžœ
˜Kšžœžœžœ<˜gKšžœžœžœ<˜hJ˜6J˜-Kš
œ5˜5K˜Kšœ:Ÿ˜TJ˜-Kš
œ˜—Kš
œ˜K˜—šœžœC˜TKš
œ™K™ƒ
Kšœžœ˜)šžœžœ
žœž˜%Kšœžœ$˜,Kš
œ5˜5Kš
œ˜Kš
œ˜Kšžœ
˜—Kš
œ/˜/K˜K˜K˜—šœžœžœ˜CKš
œ ™ Kšœžœ˜)K˜$J˜!Kšœ8Ÿ˜WKš
œ-˜-K˜K˜K˜—šœžœQžœ˜rKš
œ$™$Kšœžœ˜)K˜$K˜|K˜6K˜-K˜Kš
œ˜K˜—š
œ™K™—–m -- [f: STREAM, loPoint: GGBasicTypes.Point, hiPoint: GGBasicTypes.Point] RETURNS [seg: GGModelTypes.Segment]š

œžœžœžœ
žœ˜aKšœžœ˜)šžœž˜Kš
œ˜Kš
œ$˜$Kš
œ˜Kšžœžœ
˜—K˜K˜—š	œžœžœ
žœ˜0Kš
œ™Kšœžœ˜)Kšœ	žœ$˜0šžœž˜K˜&K˜+K˜$Kšžœžœ
˜—Kš
œ#˜#šžœžœ
žœž˜K˜4K˜Kšžœ
˜—K˜K˜—–m -- [f: STREAM, loPoint: GGBasicTypes.Point, hiPoint: GGBasicTypes.Point] RETURNS [seg: GGModelTypes.Segment]šœžœžœ
žœ$žœžœ˜`Kš
œ™Kš
œ˜K˜Kšœžœ
˜Kšœžœ
˜šžœžœ˜K˜'K˜)Kšžœ!žœ˜8Kšžœ
žœ'žœ˜BKšžœ
žœ!žœ˜;Kšžœ
žœ
˜K˜
—šžœ
žœžœ˜!K˜)K˜;Kš	žœ
žœžœ
žœŸ-˜CKšžœžœžœ˜5K˜
—Kšžœ˜Kšœžœ?˜Hšžœžœ
žœž˜"Kš
œ+˜+Kšžœ
˜—Kš	žœžœžœ
žœŸ=˜tKš
œ˜Kš
œ'˜'Kšœ!žœ˜&K˜K˜—š
œ™K™—š
œ-˜7Kšœ/žœžœ-™oKšœžœ˜)šžœ
ž˜Kš
œ>˜>Kš
œ8˜8Kš
œ>˜>Kšžœžœ
˜—K˜K˜—–| -- [seg: GGSegmentTypes.Segment, controlPointNum: NAT, selectClass: GGSegmentTypes.SelectionClass] RETURNS [selected: BOOL]š
œ-˜7Kšœ/žœ.žœžœ
™xKšœžœ˜)šžœ
ž˜Kš
œ>˜>Kš
œ8˜8Kš
œ>˜>Kšžœžœ
˜—K˜K˜—š
œ™K™—šœžœ!žœžœ˜WKš
œ"™"Kšœžœ˜)šžœž˜˜Kšžœ&žœ
žœ
˜3Kšžœ*˜0K˜—˜Kšžœ&žœ
žœ
˜3Kšžœ*˜0K˜—Kšžœžœ
˜Kš
œ˜K˜——šœžœžœžœ˜MKš
œ$™$KšœžœŸ0˜Zšžœž˜KšœžœŸ#˜UKšœžœŸ#˜HKšžœžœ
˜—Kš
œ˜K™—š
œ™K™—šœžœ-žœžœ!žœžœ˜
Kš
œ&™&Kšœžœ˜)Jš
ŸL™LJšœ
žœ
˜Kšœ	žœ
˜Kšžœ&žœ
žœ
žœ˜GKš
œI˜IK˜šžœžœ
žœ'ž˜5šžœWžœ˜_K˜K˜K˜—Kšžœ
˜—K˜ K˜!K˜0Kš
œ˜K˜—š
œ
™
K™—š
œžœžœ˜MKš
œ»
™»
Jšœžœ˜)J˜Jšœžœ$˜3Jšœ	žœ#˜/Jšœžœ#Ÿ)˜ZJšœžœ˜(Kš
œ$˜$J˜J˜Jšœžœ˜Jšœžœ˜Jšžœžœ&Ÿ˜\Jš
Πbc-™-Jšœžœ:˜IJšœžœ"˜9Jšœžœ
˜"Jšœžœ˜ J˜šžœžœ
žœž˜J˜6J˜+Jšžœ
˜—J˜8J˜.Jš
’-™-Jšœžœ:˜IJšœžœ"˜9Jšœžœ˜Jšœžœ
˜#J˜šžœžœ
žœž˜'J˜@J˜5Jšžœ
˜—J˜8J˜.KšœŸ"˜:KšœŸ"˜:J˜J˜—šœžœ
žœ8žœ˜kKš
ϩ
™©
Jšœžœ˜)Jšœžœ$˜3Jšœ	žœ#˜/Jšœžœ#ŸT˜†
J˜#Jšœžœ˜2Kš
œ$˜$Jšœžœ™#Jšœžœ2˜FJšœžœ%˜AJšœžœ
˜$Jšœžœ
˜%šžœžœ
žœž˜J˜8J˜0Jšžœ
˜—J˜+JšœžœŸ#˜QJšœžœ˜+Jšœžœ
˜.šžœžœ
žœž˜(J˜:J˜2Jšžœ
˜—J˜BJ˜3KšœŸ"˜<Jšžœ
˜J˜J˜—šœžœ
žœ>žœ˜tK™Cš
œžœžœžœ	˜Lšžœžœ
žœž˜Kšžœžœ#˜4Kšžœ
˜—K˜—Kšœžœ˜)Jšœžœ˜Kšœžœ$˜3Kšœžœ˜,K˜#Kš
œ$˜$Kšœžœ˜2Kšœžœ@˜TKšœžœ8˜TJšœžœ
˜$Jšœžœ
˜%šžœžœ
žœž˜šžœ
žœ
žœ˜J˜EJ˜=J˜$J˜Jšžœ
˜——J˜BJ˜3KšœŸ"˜<Jšžœ
˜Jš
œ˜K™—šœžœ!žœ˜NKšœžœ˜)šžœž˜K˜FK˜:Kšžœžœ
˜—K˜K˜—š
œžœžœ˜>Kšœžœ˜)šžœž˜˜Kšžœžœ˜)Kšžœ0˜4K˜—K˜I˜Kšžœžœ;˜EKšžœl˜pK˜—Kšžœžœ
˜—K˜—K˜š
œ™K™—š
œžœ*žœžœ˜NKš
œ™Kšœ-žœ™IKš
œ%˜%K˜Kšžœ
žœ
žœ&Ÿ˜YKš
œ˜Kš
œ<˜<K˜K™—šœžœ'˜<Kš
œ™K™WKš
œ%˜%K˜RKš
œ˜Kš
œ ˜ K˜K™—šœžœ‘
˜₯
Kš
œ™KšœŸ™™
Kš
œ%˜%Kš
œ#˜#K˜0K˜—K˜šœ	žœžœ
˜K™—š
œžœ-žœžœžœ˜sKš
œ™Kš
œB™BKš
œ˜–[bBox: 
GGBasicTypes.BoundBox]šžœ	žœ˜K–[bBox: 
GGBasicTypes.BoundBox]šœ„
žœžœ˜’
Kš
žœ
žœ!žœ
žœžœ˜AK˜—Kš
œK˜KKš
œ˜Kšœ œ˜MKš
œJ˜JK˜K˜—šœžœ.žœžœ˜RKš
žœžœžœžœžœ˜]K˜K˜—šœžœ˜%Kš
œ™K™KKšœžœ˜)K˜*Kšœžœ
˜Kšœžœ
˜šžœžœ
žœž˜Kšžœžœ˜*Kšžœžœ˜*Kšžœžœ˜*Kšžœžœ˜*Kšžœ
˜—Kšœ?žœžœ˜WKš
œ$˜$K˜—K˜š
 ™K™—š	œžœžœžœ˜FKšžœžœžœ˜K˜K˜—š	œžœžœ˜BKšžœ˜K˜—šœžœ%žœžœ˜KJš
œ:™:Jšœžœ˜4Kšœžœ˜2šžœžœ
žœž˜'J˜$J˜(Jšžœ,žœ
žœ˜>Jšžœ
˜—Jšžœ˜J˜J˜—š
œžœžœžœžœ˜NJš
Ÿ_™_Jšœžœ˜)Jšœ	žœ$˜0šžœž˜J˜)J˜Jšžœžœ
˜—Jšžœ
žœ˜Jšžœ
žœžœ˜1Jšžœ
˜J˜J˜—šœžœžœ˜@šžœ
žœ
ž˜Kš
œ%˜%Kš
œ$˜$Kš
œ*˜*Kšžœžœ
˜—Kšžœ˜
K˜K˜—š
œžœžœžœ˜TKšœžœ
žœ˜Kšœžœ˜˜!Kšœžœ&˜:Kšœžœ
˜Kš
œ˜K˜—Kš
œJ˜JK˜Kšžœ˜
K˜K˜—š
œžœžœ
žœžœ˜JKšœžœ˜Kšœžœ˜0šžœžœ
žœžœžœ
ž˜6K˜Kšžœ
˜—Kšœžœ!˜2Kšœ!žœ!™Fšžœžœ
žœžœžœ
ž˜6KšœŸ3˜DKš
œ˜Kšžœ
˜—Kš
œ$˜$Kšžœ˜
K˜—K˜š
 
™
K™—šœžœ˜Kšœžœ;˜SKš
œ˜KšœžœE˜SKš
œ˜Kšœžœ9˜GKš
œ˜K˜K˜—š	œžœ˜Kš
œ ˜ Kš
œ;˜;Kš
œ?˜?K˜—K˜K˜Kš
œ˜K˜Kšžœ
˜—…—›Φηκ