<> <> DIRECTORY Rope, SetBasics; LichenIntBasics: CEDAR DEFINITIONS IMPORTS SetBasics = { LNAT: TYPE ~ INT--[0 .. INT.LAST], but we can't write that--; ROPE: TYPE ~ Rope.ROPE; Int2: TYPE ~ PACKED ARRAY Dim2 OF Int; Dim2: TYPE ~ {X, Y}; Int: TYPE ~ INT16; Range: TYPE = RECORD [min, maxPlusOne: Int]; Range2: TYPE = ARRAY Dim2 OF Range; OtherDim2: ARRAY Dim2 OF Dim2 ~ [X: Y, Y: X]; ConsInt2: PROC [d1: Dim2, x1, x2: Int] RETURNS [x: Int2] ~ INLINE {x[d1] _ x1; x[OtherDim2[d1]] _ x2}; Neg: PROC [a: Int2] RETURNS [Int2] ~ INLINE {RETURN [[X: -a[X], Y: -a[Y]]]}; Add: PROC [a, b: Int2] RETURNS [Int2] ~ INLINE {RETURN [[X: a[X]+b[X], Y: a[Y]+b[Y]]]}; Sub: PROC [a, b: Int2] RETURNS [Int2] ~ INLINE {RETURN [[X: a[X]-b[X], Y: a[Y]-b[Y]]]}; InRange: PROC [i: Int2, r: Range2] RETURNS [BOOL] ~ INLINE {RETURN [i[X] IN [r[X].min .. r[X].maxPlusOne) AND i[Y] IN [r[Y].min .. r[Y].maxPlusOne)]}; Tweak: PROC [i: Int2, d: Dim2, ~ INLINE {i[d] _ i[d] + Mul: PROC [i: Int2, ~ INLINE {RETURN [[X: i[X]* Scale: PROC [a: Int2, b: Int] RETURNS [Int2] ~ INLINE {RETURN [[X: a[X]*b, Y: a[Y]*b]]}; Dot: PROC [a, b: Int2] RETURNS [Int] ~ INLINE {RETURN [a[X]*b[X]+a[Y]*b[Y]]}; Cross: PROC [a, b: Int2] RETURNS [INT] ~ INLINE {RETURN [INT[a[X]]*b[Y]-INT[a[Y]]*b[X]]}; Mod: PROC [a, mod: Int2] RETURNS [Int2] ~ INLINE {RETURN [[ X: ((a[X] MOD mod[X])+mod[X]) MOD mod[X], Y: ((a[Y] MOD mod[Y])+mod[Y]) MOD mod[Y]]]}; Div: PROC [a, b: Int2] RETURNS [Int2] ~ INLINE {RETURN [[X: a[X]/b[X], Y: a[Y]/b[Y]]]}; AddMod: PROC [a: Int2, b: Int2, mod: Int2] RETURNS [Int2] ~ INLINE {RETURN [[ X: (a[X]+b[X]+mod[X]) MOD mod[X], Y: (a[Y]+b[Y]+mod[Y]) MOD mod[Y]]]}; Area: PROC [x: Int2] RETURNS [LNAT] ~ INLINE {RETURN [x[X] * x[Y]]}; Int2Hash: PROC [x: Int2] RETURNS [CARDINAL] ~ INLINE {RETURN [17*LOOPHOLE[x[X], CARDINAL] + 257*LOOPHOLE[x[Y], CARDINAL]]}; Int2Compare: PROC [a, b: Int2] RETURNS [SetBasics.TotalComparison] ~ INLINE {RETURN SetBasics.CompareIntI[LOOPHOLE[a], LOOPHOLE[b]]}; RangeOff: PROC [r: Range, ~ INLINE {RETURN[[min: r.min+ RangeOffClip: PROC [r: Range, ~ INLINE {RETURN[[min: MAX[r.min+ ShaveRange2Top1: PROC [r: Range2, d: Dim2] RETURNS [Range2] ~ INLINE {r[d].min _ MIN[r[d].min, r[d].maxPlusOne _ r[d].maxPlusOne - 1]; RETURN [r]}; ConsRange2: PROC [d1: Dim2, x1, x2: Range] RETURNS [x: Range2] ~ INLINE {x[d1] _ x1; x[OtherDim2[d1]] _ x2}; Range2Empty: PROC [r: Range2] RETURNS [BOOL] ~ INLINE {RETURN [r[X].maxPlusOne<=r[X].min OR r[Y].maxPlusOne<=r[Y].min]}; Range2IsSingleton: PROC [r: Range2] RETURNS [BOOL] ~ INLINE {RETURN [r[X].maxPlusOne=r[X].min+1 AND r[Y].maxPlusOne=r[Y].min+1]}; Range2Min: PROC [r2: Range2] RETURNS [Int2] ~ INLINE {RETURN[[X: r2[X].min, Y: r2[Y].min]]}; Range2Off: PROC [r: Range2, ~ INLINE {RETURN[[X: RangeOff[r[X], Range2OffClip: PROC [r: Range2, ~ INLINE {RETURN[[X: RangeOffClip[r[X], Range2Included: PROC [sub, in: Range2] RETURNS [BOOL] ~ INLINE {RETURN [RangeIncluded[sub[X], in[X]] AND RangeIncluded[sub[Y], in[Y]]]}; RangeIncluded: PROC [sub, in: Range] RETURNS [BOOL] ~ INLINE {RETURN [sub.min>=in.min AND sub.maxPlusOne<=in.maxPlusOne]}; Range2Intersection: PROC [a, b: Range2] RETURNS [Range2] ~ INLINE {RETURN [[ X: [ min: MAX[a[X].min, b[X].min], maxPlusOne: MIN[a[X].maxPlusOne, b[X].maxPlusOne]], Y: [ min: MAX[a[Y].min, b[Y].min], maxPlusOne: MIN[a[Y].maxPlusOne, b[Y].maxPlusOne]]]]}; RangeLength: PROC [r: Range] RETURNS [NATURAL] ~ INLINE {RETURN [r.maxPlusOne - r.min]}; RangeArea: PROC [r: Range2] RETURNS [area: LNAT] ~ INLINE {area _ INT[RangeLength[r[X]]] * RangeLength[r[Y]]}; RangeShape: PROC [r: Range2] RETURNS [Int2] ~ INLINE {RETURN [[RangeLength[r[X]], RangeLength[r[Y]]]]}; SizeRange: PROC [size: Int2] RETURNS [Range2] ~ INLINE {RETURN [[[0, size[X]], [0, size[Y]]]]}; Int2sRange: PROC [a, b: Int2] RETURNS [r: Range2] ~ INLINE {RETURN [[ X: [MIN[a[X], b[X]], MAX[a[X], b[X]]+1], Y: [MIN[a[Y], b[Y]], MAX[a[Y], b[Y]]+1]]]}; Range2Mbb: PROC [a, b: Range2] RETURNS [Range2] ~ INLINE {RETURN [[X: RangeMbb[a[X], b[X]], Y: RangeMbb[a[Y], b[Y]]]]}; RangeMbb: PROC [a, b: Range] RETURNS [Range] ~ INLINE {RETURN [[min: MIN[a.min, b.min], maxPlusOne: MAX[a.maxPlusOne, b.maxPlusOne]]]}; Range2sIntersect: PROC [r1, r2: Range2] RETURNS [BOOL] = INLINE {RETURN [RangesIntersect[r1[X], r2[X]] AND RangesIntersect[r1[Y], r2[Y]]]}; RangesIntersect: PROC [r1, r2: Range] RETURNS [BOOL] = INLINE {RETURN [ (r1.min IN [r2.min .. r2.maxPlusOne) AND r1.maxPlusOne > r1.min) OR (r2.min IN [r1.min .. r1.maxPlusOne) AND r2.maxPlusOne > r2.min)]}; Range2Div: PROC [r: Range2, Range1Div: PROC [r: Range, <> Range2RoundA: PROC [r: Range2, Range2MulA: PROC [r: Range2, Range1MulA: PROC [r: Range, <<.maxPlusOne-1 has given phase.>> Range2RoundB: PROC [r: Range2, Range2MulB: PROC [r: Range2, Range1MulB: PROC [r: Range, <<.maxPlusOne has given phase.>> FmtRange: PROC [r: Range] RETURNS [ROPE]; }.