@jax-data-science/sci-image-visualizer 0.8.2 → 0.9.0

package/esm2022/lib/geometry/raster.mjs

@@ -0,0 +1,98 @@
+import { rasterizePolygon as rasterizeRing } from 'napari-js/geometry';
+/**
+ * The largest bbox (in pixels) {@link rasterizePolygon} allocates before it
+ * falls back to clipping at the window. Also the cap the wand/brush use to
+ * decide whether a region can be edited at the current zoom.
+ */
+export const MAX_RASTER_PIXELS = 4096 * 4096;
+/** `[x0, y0, x1, y1, …]` from parallel arrays, the ring layout napari-js takes. */
+function flatRing(xs, ys) {
+    const out = new Array(xs.length * 2);
+    for (let i = 0; i < xs.length; i++) {
+        out[i * 2] = xs[i];
+        out[i * 2 + 1] = ys[i];
+    }
+    return out;
+}
+/**
+ * Rasterize a closed polygon (minus its holes) into a bbox-relative mask, e.g.
+ * as the starting point for a wand/brush stroke accumulator.
+ *
+ * Pixel (px, py) covers `[px, px + 1) × [py, py + 1)` and is set when its centre
+ * is inside — the rule under which the pixel-edge outlines `geometry/contour`
+ * traces come back as exactly the pixels they were traced from.
+ *
+ * Returns the bbox origin plus the filled mask, or null if the polygon is
+ * degenerate (covers no pixel centre). The polygon's full extent is kept even
+ * where it is negative or exceeds `imageWidth`×`imageHeight` (jit-ui#102); only
+ * a bbox larger than {@link MAX_RASTER_PIXELS} is clipped to that window (rare:
+ * a region enormous at the current zoom; accepts the clip).
+ */
+export function rasterizePolygon(xpoints, ypoints, imageWidth, imageHeight, holes) {
+    if (xpoints.length < 3)
+        return null;
+    const outer = flatRing(xpoints, ypoints);
+    const inner = holes?.map((ring) => flatRing(ring.map((p) => p[0]), ring.map((p) => p[1])));
+    const bounds = { x: 0, y: 0, width: imageWidth, height: imageHeight };
+    let m;
+    try {
+        m = rasterizeRing(outer, inner, { bounds, maxPixels: MAX_RASTER_PIXELS });
+    }
+    catch (err) {
+        // Still over the cap once clipped to the window (a large export mask): allocate the clip.
+        if (!(err instanceof RangeError))
+            throw err;
+        m = rasterizeRing(outer, inner, { bounds, clip: true, maxPixels: Infinity });
+    }
+    return m && { bx: m.x, by: m.y, bw: m.width, bh: m.height, mask: m.data };
+}
+/** Returns true if the two bbox-relative masks share at least one set pixel. */
+export function masksOverlap(a, b) {
+    const ix0 = Math.max(a.bx, b.bx);
+    const iy0 = Math.max(a.by, b.by);
+    const ix1 = Math.min(a.bx + a.bw, b.bx + b.bw);
+    const iy1 = Math.min(a.by + a.bh, b.by + b.bh);
+    if (ix0 >= ix1 || iy0 >= iy1)
+        return false;
+    for (let y = iy0; y < iy1; y++) {
+        const ar = (y - a.by) * a.bw;
+        const br = (y - b.by) * b.bw;
+        for (let x = ix0; x < ix1; x++) {
+            if (a.mask[ar + (x - a.bx)] && b.mask[br + (x - b.bx)])
+                return true;
+        }
+    }
+    return false;
+}
+/**
+ * Allocate a fresh accumulator covering both masks and OR them in. The
+ * resulting mask's bbox is the smallest rectangle containing both inputs.
+ */
+export function unionMasks(a, b) {
+    const ux0 = Math.min(a.bx, b.bx);
+    const uy0 = Math.min(a.by, b.by);
+    const ux1 = Math.max(a.bx + a.bw, b.bx + b.bw);
+    const uy1 = Math.max(a.by + a.bh, b.by + b.bh);
+    const ubw = ux1 - ux0;
+    const ubh = uy1 - uy0;
+    const umask = new Uint8Array(ubw * ubh);
+    const dax = a.bx - ux0;
+    const day = a.by - uy0;
+    for (let row = 0; row < a.bh; row++) {
+        const src = row * a.bw;
+        const dst = (row + day) * ubw + dax;
+        umask.set(a.mask.subarray(src, src + a.bw), dst);
+    }
+    const dbx = b.bx - ux0;
+    const dby = b.by - uy0;
+    for (let row = 0; row < b.bh; row++) {
+        const srcRow = row * b.bw;
+        const dstRow = (row + dby) * ubw;
+        for (let col = 0; col < b.bw; col++) {
+            if (b.mask[srcRow + col])
+                umask[dstRow + (dbx + col)] = 1;
+        }
+    }
+    return { bx: ux0, by: uy0, bw: ubw, bh: ubh, mask: umask };
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoicmFzdGVyLmpzIiwic291cmNlUm9vdCI6IiIsInNvdXJjZXMiOlsiLi4vLi4vLi4vLi4vc3JjL2xpYi9nZW9tZXRyeS9yYXN0ZXIudHMiXSwibmFtZXMiOltdLCJtYXBwaW5ncyI6IkFBQUEsT0FBTyxFQUFFLGdCQUFnQixJQUFJLGFBQWEsRUFBRSxNQUFNLG9CQUFvQixDQUFDO0FBa0J2RTs7OztHQUlHO0FBQ0gsTUFBTSxDQUFDLE1BQU0saUJBQWlCLEdBQUcsSUFBSSxHQUFHLElBQUksQ0FBQztBQUU3QyxtRkFBbUY7QUFDbkYsU0FBUyxRQUFRLENBQUMsRUFBcUIsRUFBRSxFQUFxQjtJQUM1RCxNQUFNLEdBQUcsR0FBRyxJQUFJLEtBQUssQ0FBUyxFQUFFLENBQUMsTUFBTSxHQUFHLENBQUMsQ0FBQyxDQUFDO0lBQzdDLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxFQUFFLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7UUFDbkMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDbkIsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ3pCLENBQUM7SUFDRCxPQUFPLEdBQUcsQ0FBQztBQUNiLENBQUM7QUFFRDs7Ozs7Ozs7Ozs7OztHQWFHO0FBQ0gsTUFBTSxVQUFVLGdCQUFnQixDQUM5QixPQUFpQixFQUNqQixPQUFpQixFQUNqQixVQUFrQixFQUNsQixXQUFtQixFQUNuQixLQUFvQjtJQUVwQixJQUFJLE9BQU8sQ0FBQyxNQUFNLEdBQUcsQ0FBQztRQUFFLE9BQU8sSUFBSSxDQUFDO0lBQ3BDLE1BQU0sS0FBSyxHQUFHLFFBQVEsQ0FBQyxPQUFPLEVBQUUsT0FBTyxDQUFDLENBQUM7SUFDekMsTUFBTSxLQUFLLEdBQUcsS0FBSyxFQUFFLEdBQUcsQ0FBQyxDQUFDLElBQUksRUFBRSxFQUFFLENBQ2hDLFFBQVEsQ0FDTixJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsRUFDckIsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQ3RCLENBQ0YsQ0FBQztJQUNGLE1BQU0sTUFBTSxHQUFHLEVBQUUsQ0FBQyxFQUFFLENBQUMsRUFBRSxDQUFDLEVBQUUsQ0FBQyxFQUFFLEtBQUssRUFBRSxVQUFVLEVBQUUsTUFBTSxFQUFFLFdBQVcsRUFBRSxDQUFDO0lBQ3RFLElBQUksQ0FBbUMsQ0FBQztJQUN4QyxJQUFJLENBQUM7UUFDSCxDQUFDLEdBQUcsYUFBYSxDQUFDLEtBQUssRUFBRSxLQUFLLEVBQUUsRUFBRSxNQUFNLEVBQUUsU0FBUyxFQUFFLGlCQUFpQixFQUFFLENBQUMsQ0FBQztJQUM1RSxDQUFDO0lBQUMsT0FBTyxHQUFHLEVBQUUsQ0FBQztRQUNiLDBGQUEwRjtRQUMxRixJQUFJLENBQUMsQ0FBQyxHQUFHLFlBQVksVUFBVSxDQUFDO1lBQUUsTUFBTSxHQUFHLENBQUM7UUFDNUMsQ0FBQyxHQUFHLGFBQWEsQ0FBQyxLQUFLLEVBQUUsS0FBSyxFQUFFLEVBQUUsTUFBTSxFQUFFLElBQUksRUFBRSxJQUFJLEVBQUUsU0FBUyxFQUFFLFFBQVEsRUFBRSxDQUFDLENBQUM7SUFDL0UsQ0FBQztJQUNELE9BQU8sQ0FBQyxJQUFJLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxLQUFLLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxNQUFNLEVBQUUsSUFBSSxFQUFFLENBQUMsQ0FBQyxJQUFJLEVBQUUsQ0FBQztBQUM1RSxDQUFDO0FBRUQsZ0ZBQWdGO0FBQ2hGLE1BQU0sVUFBVSxZQUFZLENBQUMsQ0FBVyxFQUFFLENBQVc7SUFDbkQsTUFBTSxHQUFHLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUNqQyxNQUFNLEdBQUcsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDO0lBQ2pDLE1BQU0sR0FBRyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEVBQUUsR0FBRyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDO0lBQy9DLE1BQU0sR0FBRyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEVBQUUsR0FBRyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDO0lBQy9DLElBQUksR0FBRyxJQUFJLEdBQUcsSUFBSSxHQUFHLElBQUksR0FBRztRQUFFLE9BQU8sS0FBSyxDQUFDO0lBQzNDLEtBQUssSUFBSSxDQUFDLEdBQUcsR0FBRyxFQUFFLENBQUMsR0FBRyxHQUFHLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQztRQUMvQixNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQztRQUM3QixNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQztRQUM3QixLQUFLLElBQUksQ0FBQyxHQUFHLEdBQUcsRUFBRSxDQUFDLEdBQUcsR0FBRyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7WUFDL0IsSUFBSSxDQUFDLENBQUMsSUFBSSxDQUFDLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUMsSUFBSSxDQUFDLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLENBQUM7Z0JBQUUsT0FBTyxJQUFJLENBQUM7UUFDdEUsQ0FBQztJQUNILENBQUM7SUFDRCxPQUFPLEtBQUssQ0FBQztBQUNmLENBQUM7QUFFRDs7O0dBR0c7QUFDSCxNQUFNLFVBQVUsVUFBVSxDQUFDLENBQVcsRUFBRSxDQUFXO0lBQ2pELE1BQU0sR0FBRyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsRUFBRSxDQUFDLENBQUM7SUFDakMsTUFBTSxHQUFHLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUNqQyxNQUFNLEdBQUcsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsRUFBRSxHQUFHLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUMvQyxNQUFNLEdBQUcsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsRUFBRSxHQUFHLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUMvQyxNQUFNLEdBQUcsR0FBRyxHQUFHLEdBQUcsR0FBRyxDQUFDO0lBQ3RCLE1BQU0sR0FBRyxHQUFHLEdBQUcsR0FBRyxHQUFHLENBQUM7SUFDdEIsTUFBTSxLQUFLLEdBQUcsSUFBSSxVQUFVLENBQUMsR0FBRyxHQUFHLEdBQUcsQ0FBQyxDQUFDO0lBRXhDLE1BQU0sR0FBRyxHQUFHLENBQUMsQ0FBQyxFQUFFLEdBQUcsR0FBRyxDQUFDO0lBQ3ZCLE1BQU0sR0FBRyxHQUFHLENBQUMsQ0FBQyxFQUFFLEdBQUcsR0FBRyxDQUFDO0lBQ3ZCLEtBQUssSUFBSSxHQUFHLEdBQUcsQ0FBQyxFQUFFLEdBQUcsR0FBRyxDQUFDLENBQUMsRUFBRSxFQUFFLEdBQUcsRUFBRSxFQUFFLENBQUM7UUFDcEMsTUFBTSxHQUFHLEdBQUcsR0FBRyxHQUFHLENBQUMsQ0FBQyxFQUFFLENBQUM7UUFDdkIsTUFBTSxHQUFHLEdBQUcsQ0FBQyxHQUFHLEdBQUcsR0FBRyxDQUFDLEdBQUcsR0FBRyxHQUFHLEdBQUcsQ0FBQztRQUNwQyxLQUFLLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxJQUFJLENBQUMsUUFBUSxDQUFDLEdBQUcsRUFBRSxHQUFHLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDO0lBQ25ELENBQUM7SUFDRCxNQUFNLEdBQUcsR0FBRyxDQUFDLENBQUMsRUFBRSxHQUFHLEdBQUcsQ0FBQztJQUN2QixNQUFNLEdBQUcsR0FBRyxDQUFDLENBQUMsRUFBRSxHQUFHLEdBQUcsQ0FBQztJQUN2QixLQUFLLElBQUksR0FBRyxHQUFHLENBQUMsRUFBRSxHQUFHLEdBQUcsQ0FBQyxDQUFDLEVBQUUsRUFBRSxHQUFHLEVBQUUsRUFBRSxDQUFDO1FBQ3BDLE1BQU0sTUFBTSxHQUFHLEdBQUcsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDO1FBQzFCLE1BQU0sTUFBTSxHQUFHLENBQUMsR0FBRyxHQUFHLEdBQUcsQ0FBQyxHQUFHLEdBQUcsQ0FBQztRQUNqQyxLQUFLLElBQUksR0FBRyxHQUFHLENBQUMsRUFBRSxHQUFHLEdBQUcsQ0FBQyxDQUFDLEVBQUUsRUFBRSxHQUFHLEVBQUUsRUFBRSxDQUFDO1lBQ3BDLElBQUksQ0FBQyxDQUFDLElBQUksQ0FBQyxNQUFNLEdBQUcsR0FBRyxDQUFDO2dCQUFFLEtBQUssQ0FBQyxNQUFNLEdBQUcsQ0FBQyxHQUFHLEdBQUcsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUM7UUFDNUQsQ0FBQztJQUNILENBQUM7SUFDRCxPQUFPLEVBQUUsRUFBRSxFQUFFLEdBQUcsRUFBRSxFQUFFLEVBQUUsR0FBRyxFQUFFLEVBQUUsRUFBRSxHQUFHLEVBQUUsRUFBRSxFQUFFLEdBQUcsRUFBRSxJQUFJLEVBQUUsS0FBSyxFQUFFLENBQUM7QUFDN0QsQ0FBQyIsInNvdXJjZXNDb250ZW50IjpbImltcG9ydCB7IHJhc3Rlcml6ZVBvbHlnb24gYXMgcmFzdGVyaXplUmluZyB9IGZyb20gJ25hcGFyaS1qcy9nZW9tZXRyeSc7XG5cbi8qKlxuICogUHVyZSBiaW5hcnktbWFzayBoZWxwZXJzOiBwb2x5Z29uIHJhc3Rlcml6YXRpb24gKG5hcGFyaS1qcydzIGByYXN0ZXJpemVQb2x5Z29uYClcbiAqIGFuZCBiYm94LXJlbGF0aXZlIG1hc2sgc2V0IG9wZXJhdGlvbnMuIE5vIERPTSwgbm8gQW5ndWxhciDigJQgdGhlIG1hc2stZXhwb3J0XG4gKiB3b3JrZXIgaW1wb3J0cyB0aGlzLlxuICovXG5cbi8qKiBBIGJib3gtcmVsYXRpdmUgYmluYXJ5IG1hc2sgKDEgPSBmaWxsZWQsIDAgPSBlbXB0eSkgd2hvc2UgcGl4ZWwgKDAsMCkgc2l0c1xuICogIGF0IChieCwgYnkpIGluIHRoZSBmcmFtZSBpdCB3YXMgcmFzdGVyaXplZCBpbi4gKi9cbmV4cG9ydCBpbnRlcmZhY2UgQkJveE1hc2sge1xuICBieDogbnVtYmVyO1xuICBieTogbnVtYmVyO1xuICBidzogbnVtYmVyO1xuICBiaDogbnVtYmVyO1xuICBtYXNrOiBVaW50OEFycmF5O1xufVxuXG4vKipcbiAqIFRoZSBsYXJnZXN0IGJib3ggKGluIHBpeGVscykge0BsaW5rIHJhc3Rlcml6ZVBvbHlnb259IGFsbG9jYXRlcyBiZWZvcmUgaXRcbiAqIGZhbGxzIGJhY2sgdG8gY2xpcHBpbmcgYXQgdGhlIHdpbmRvdy4gQWxzbyB0aGUgY2FwIHRoZSB3YW5kL2JydXNoIHVzZSB0b1xuICogZGVjaWRlIHdoZXRoZXIgYSByZWdpb24gY2FuIGJlIGVkaXRlZCBhdCB0aGUgY3VycmVudCB6b29tLlxuICovXG5leHBvcnQgY29uc3QgTUFYX1JBU1RFUl9QSVhFTFMgPSA0MDk2ICogNDA5NjtcblxuLyoqIGBbeDAsIHkwLCB4MSwgeTEsIOKApl1gIGZyb20gcGFyYWxsZWwgYXJyYXlzLCB0aGUgcmluZyBsYXlvdXQgbmFwYXJpLWpzIHRha2VzLiAqL1xuZnVuY3Rpb24gZmxhdFJpbmcoeHM6IEFycmF5TGlrZTxudW1iZXI+LCB5czogQXJyYXlMaWtlPG51bWJlcj4pOiBudW1iZXJbXSB7XG4gIGNvbnN0IG91dCA9IG5ldyBBcnJheTxudW1iZXI+KHhzLmxlbmd0aCAqIDIpO1xuICBmb3IgKGxldCBpID0gMDsgaSA8IHhzLmxlbmd0aDsgaSsrKSB7XG4gICAgb3V0W2kgKiAyXSA9IHhzW2ldO1xuICAgIG91dFtpICogMiArIDFdID0geXNbaV07XG4gIH1cbiAgcmV0dXJuIG91dDtcbn1cblxuLyoqXG4gKiBSYXN0ZXJpemUgYSBjbG9zZWQgcG9seWdvbiAobWludXMgaXRzIGhvbGVzKSBpbnRvIGEgYmJveC1yZWxhdGl2ZSBtYXNrLCBlLmcuXG4gKiBhcyB0aGUgc3RhcnRpbmcgcG9pbnQgZm9yIGEgd2FuZC9icnVzaCBzdHJva2UgYWNjdW11bGF0b3IuXG4gKlxuICogUGl4ZWwgKHB4LCBweSkgY292ZXJzIGBbcHgsIHB4ICsgMSkgw5cgW3B5LCBweSArIDEpYCBhbmQgaXMgc2V0IHdoZW4gaXRzIGNlbnRyZVxuICogaXMgaW5zaWRlIOKAlCB0aGUgcnVsZSB1bmRlciB3aGljaCB0aGUgcGl4ZWwtZWRnZSBvdXRsaW5lcyBgZ2VvbWV0cnkvY29udG91cmBcbiAqIHRyYWNlcyBjb21lIGJhY2sgYXMgZXhhY3RseSB0aGUgcGl4ZWxzIHRoZXkgd2VyZSB0cmFjZWQgZnJvbS5cbiAqXG4gKiBSZXR1cm5zIHRoZSBiYm94IG9yaWdpbiBwbHVzIHRoZSBmaWxsZWQgbWFzaywgb3IgbnVsbCBpZiB0aGUgcG9seWdvbiBpc1xuICogZGVnZW5lcmF0ZSAoY292ZXJzIG5vIHBpeGVsIGNlbnRyZSkuIFRoZSBwb2x5Z29uJ3MgZnVsbCBleHRlbnQgaXMga2VwdCBldmVuXG4gKiB3aGVyZSBpdCBpcyBuZWdhdGl2ZSBvciBleGNlZWRzIGBpbWFnZVdpZHRoYMOXYGltYWdlSGVpZ2h0YCAoaml0LXVpIzEwMik7IG9ubHlcbiAqIGEgYmJveCBsYXJnZXIgdGhhbiB7QGxpbmsgTUFYX1JBU1RFUl9QSVhFTFN9IGlzIGNsaXBwZWQgdG8gdGhhdCB3aW5kb3cgKHJhcmU6XG4gKiBhIHJlZ2lvbiBlbm9ybW91cyBhdCB0aGUgY3VycmVudCB6b29tOyBhY2NlcHRzIHRoZSBjbGlwKS5cbiAqL1xuZXhwb3J0IGZ1bmN0aW9uIHJhc3Rlcml6ZVBvbHlnb24oXG4gIHhwb2ludHM6IG51bWJlcltdLFxuICB5cG9pbnRzOiBudW1iZXJbXSxcbiAgaW1hZ2VXaWR0aDogbnVtYmVyLFxuICBpbWFnZUhlaWdodDogbnVtYmVyLFxuICBob2xlcz86IG51bWJlcltdW11bXSxcbik6IEJCb3hNYXNrIHwgbnVsbCB7XG4gIGlmICh4cG9pbnRzLmxlbmd0aCA8IDMpIHJldHVybiBudWxsO1xuICBjb25zdCBvdXRlciA9IGZsYXRSaW5nKHhwb2ludHMsIHlwb2ludHMpO1xuICBjb25zdCBpbm5lciA9IGhvbGVzPy5tYXAoKHJpbmcpID0+XG4gICAgZmxhdFJpbmcoXG4gICAgICByaW5nLm1hcCgocCkgPT4gcFswXSksXG4gICAgICByaW5nLm1hcCgocCkgPT4gcFsxXSksXG4gICAgKSxcbiAgKTtcbiAgY29uc3QgYm91bmRzID0geyB4OiAwLCB5OiAwLCB3aWR0aDogaW1hZ2VXaWR0aCwgaGVpZ2h0OiBpbWFnZUhlaWdodCB9O1xuICBsZXQgbTogUmV0dXJuVHlwZTx0eXBlb2YgcmFzdGVyaXplUmluZz47XG4gIHRyeSB7XG4gICAgbSA9IHJhc3Rlcml6ZVJpbmcob3V0ZXIsIGlubmVyLCB7IGJvdW5kcywgbWF4UGl4ZWxzOiBNQVhfUkFTVEVSX1BJWEVMUyB9KTtcbiAgfSBjYXRjaCAoZXJyKSB7XG4gICAgLy8gU3RpbGwgb3ZlciB0aGUgY2FwIG9uY2UgY2xpcHBlZCB0byB0aGUgd2luZG93IChhIGxhcmdlIGV4cG9ydCBtYXNrKTogYWxsb2NhdGUgdGhlIGNsaXAuXG4gICAgaWYgKCEoZXJyIGluc3RhbmNlb2YgUmFuZ2VFcnJvcikpIHRocm93IGVycjtcbiAgICBtID0gcmFzdGVyaXplUmluZyhvdXRlciwgaW5uZXIsIHsgYm91bmRzLCBjbGlwOiB0cnVlLCBtYXhQaXhlbHM6IEluZmluaXR5IH0pO1xuICB9XG4gIHJldHVybiBtICYmIHsgYng6IG0ueCwgYnk6IG0ueSwgYnc6IG0ud2lkdGgsIGJoOiBtLmhlaWdodCwgbWFzazogbS5kYXRhIH07XG59XG5cbi8qKiBSZXR1cm5zIHRydWUgaWYgdGhlIHR3byBiYm94LXJlbGF0aXZlIG1hc2tzIHNoYXJlIGF0IGxlYXN0IG9uZSBzZXQgcGl4ZWwuICovXG5leHBvcnQgZnVuY3Rpb24gbWFza3NPdmVybGFwKGE6IEJCb3hNYXNrLCBiOiBCQm94TWFzayk6IGJvb2xlYW4ge1xuICBjb25zdCBpeDAgPSBNYXRoLm1heChhLmJ4LCBiLmJ4KTtcbiAgY29uc3QgaXkwID0gTWF0aC5tYXgoYS5ieSwgYi5ieSk7XG4gIGNvbnN0IGl4MSA9IE1hdGgubWluKGEuYnggKyBhLmJ3LCBiLmJ4ICsgYi5idyk7XG4gIGNvbnN0IGl5MSA9IE1hdGgubWluKGEuYnkgKyBhLmJoLCBiLmJ5ICsgYi5iaCk7XG4gIGlmIChpeDAgPj0gaXgxIHx8IGl5MCA+PSBpeTEpIHJldHVybiBmYWxzZTtcbiAgZm9yIChsZXQgeSA9IGl5MDsgeSA8IGl5MTsgeSsrKSB7XG4gICAgY29uc3QgYXIgPSAoeSAtIGEuYnkpICogYS5idztcbiAgICBjb25zdCBiciA9ICh5IC0gYi5ieSkgKiBiLmJ3O1xuICAgIGZvciAobGV0IHggPSBpeDA7IHggPCBpeDE7IHgrKykge1xuICAgICAgaWYgKGEubWFza1thciArICh4IC0gYS5ieCldICYmIGIubWFza1ticiArICh4IC0gYi5ieCldKSByZXR1cm4gdHJ1ZTtcbiAgICB9XG4gIH1cbiAgcmV0dXJuIGZhbHNlO1xufVxuXG4vKipcbiAqIEFsbG9jYXRlIGEgZnJlc2ggYWNjdW11bGF0b3IgY292ZXJpbmcgYm90aCBtYXNrcyBhbmQgT1IgdGhlbSBpbi4gVGhlXG4gKiByZXN1bHRpbmcgbWFzaydzIGJib3ggaXMgdGhlIHNtYWxsZXN0IHJlY3RhbmdsZSBjb250YWluaW5nIGJvdGggaW5wdXRzLlxuICovXG5leHBvcnQgZnVuY3Rpb24gdW5pb25NYXNrcyhhOiBCQm94TWFzaywgYjogQkJveE1hc2spOiBCQm94TWFzayB7XG4gIGNvbnN0IHV4MCA9IE1hdGgubWluKGEuYngsIGIuYngpO1xuICBjb25zdCB1eTAgPSBNYXRoLm1pbihhLmJ5LCBiLmJ5KTtcbiAgY29uc3QgdXgxID0gTWF0aC5tYXgoYS5ieCArIGEuYncsIGIuYnggKyBiLmJ3KTtcbiAgY29uc3QgdXkxID0gTWF0aC5tYXgoYS5ieSArIGEuYmgsIGIuYnkgKyBiLmJoKTtcbiAgY29uc3QgdWJ3ID0gdXgxIC0gdXgwO1xuICBjb25zdCB1YmggPSB1eTEgLSB1eTA7XG4gIGNvbnN0IHVtYXNrID0gbmV3IFVpbnQ4QXJyYXkodWJ3ICogdWJoKTtcblxuICBjb25zdCBkYXggPSBhLmJ4IC0gdXgwO1xuICBjb25zdCBkYXkgPSBhLmJ5IC0gdXkwO1xuICBmb3IgKGxldCByb3cgPSAwOyByb3cgPCBhLmJoOyByb3crKykge1xuICAgIGNvbnN0IHNyYyA9IHJvdyAqIGEuYnc7XG4gICAgY29uc3QgZHN0ID0gKHJvdyArIGRheSkgKiB1YncgKyBkYXg7XG4gICAgdW1hc2suc2V0KGEubWFzay5zdWJhcnJheShzcmMsIHNyYyArIGEuYncpLCBkc3QpO1xuICB9XG4gIGNvbnN0IGRieCA9IGIuYnggLSB1eDA7XG4gIGNvbnN0IGRieSA9IGIuYnkgLSB1eTA7XG4gIGZvciAobGV0IHJvdyA9IDA7IHJvdyA8IGIuYmg7IHJvdysrKSB7XG4gICAgY29uc3Qgc3JjUm93ID0gcm93ICogYi5idztcbiAgICBjb25zdCBkc3RSb3cgPSAocm93ICsgZGJ5KSAqIHVidztcbiAgICBmb3IgKGxldCBjb2wgPSAwOyBjb2wgPCBiLmJ3OyBjb2wrKykge1xuICAgICAgaWYgKGIubWFza1tzcmNSb3cgKyBjb2xdKSB1bWFza1tkc3RSb3cgKyAoZGJ4ICsgY29sKV0gPSAxO1xuICAgIH1cbiAgfVxuICByZXR1cm4geyBieDogdXgwLCBieTogdXkwLCBidzogdWJ3LCBiaDogdWJoLCBtYXNrOiB1bWFzayB9O1xufVxuIl19

package/esm2022/lib/geometry/ring.mjs

@@ -0,0 +1,201 @@
+/**
+ * Pure ring (closed polyline) helpers: bounds, containment, vertex dropping and
+ * simplification. No DOM, no Angular, no `this` — usable from workers.
+ *
+ * A ring is a pair of parallel `xs`/`ys` arrays without a repeated closing
+ * vertex; hole rings use the `[[x, y], …]` convention of {@link Polygon.holes}.
+ */
+/** The axis-aligned bounds of a ring's vertices. */
+export function ringBounds(xs, ys) {
+    let minX = Infinity, maxX = -Infinity, minY = Infinity, maxY = -Infinity;
+    for (let i = 0; i < xs.length; i++) {
+        if (xs[i] < minX)
+            minX = xs[i];
+        if (xs[i] > maxX)
+            maxX = xs[i];
+        if (ys[i] < minY)
+            minY = ys[i];
+        if (ys[i] > maxY)
+            maxY = ys[i];
+    }
+    return { minX, minY, maxX, maxY };
+}
+/**
+ * Standard ray-cast (even-odd) point-in-ring test. (px, py) and the ring
+ * vertices must be in the same coordinate system. Rings with fewer than three
+ * vertices contain nothing.
+ */
+export function pointInRing(px, py, xs, ys) {
+    const n = xs.length;
+    if (n < 3)
+        return false;
+    let inside = false;
+    for (let i = 0, j = n - 1; i < n; j = i++) {
+        const yi = ys[i], yj = ys[j];
+        const xi = xs[i], xj = xs[j];
+        const intersect = yi > py !== yj > py && px < ((xj - xi) * (py - yi)) / (yj - yi) + xi;
+        if (intersect)
+            inside = !inside;
+    }
+    return inside;
+}
+/**
+ * Point-in-region test that honours interior rings (holes): inside the
+ * exterior AND outside every hole. Coordinates must share the polygon's frame.
+ * With no holes this is exactly {@link pointInRing}.
+ */
+export function pointInPolygonWithHoles(px, py, xs, ys, holes) {
+    if (!pointInRing(px, py, xs, ys))
+        return false;
+    if (holes) {
+        for (const ring of holes) {
+            if (pointInPairRing(px, py, ring))
+                return false;
+        }
+    }
+    return true;
+}
+/** {@link pointInRing} over a `[[x, y], …]` ring (the hole convention), without splitting it. */
+export function pointInPairRing(px, py, ring) {
+    const n = ring.length;
+    if (n < 3)
+        return false;
+    let inside = false;
+    for (let i = 0, j = n - 1; i < n; j = i++) {
+        const xi = ring[i][0], yi = ring[i][1];
+        const xj = ring[j][0], yj = ring[j][1];
+        const intersect = yi > py !== yj > py && px < ((xj - xi) * (py - yi)) / (yj - yi) + xi;
+        if (intersect)
+            inside = !inside;
+    }
+    return inside;
+}
+/**
+ * Drop every vertex whose Euclidean distance to (cx, cy) is less than
+ * `radius`. Returns trimmed parallel arrays plus how many vertices were
+ * removed (so the caller can detect a no-op tick and skip the relayout).
+ */
+export function dropVerticesWithinRadius(xpoints, ypoints, cx, cy, radius) {
+    const r2 = radius * radius;
+    const xs = [];
+    const ys = [];
+    let removed = 0;
+    for (let i = 0; i < xpoints.length; i++) {
+        const dx = xpoints[i] - cx;
+        const dy = ypoints[i] - cy;
+        if (dx * dx + dy * dy < r2) {
+            removed++;
+            continue;
+        }
+        xs.push(xpoints[i]);
+        ys.push(ypoints[i]);
+    }
+    return { xpoints: xs, ypoints: ys, removed };
+}
+/**
+ * Erase the part of a ring that lies within `radius` of (cx, cy): every vertex
+ * inside the disc goes (as {@link dropVerticesWithinRadius}), and the outline is
+ * rejoined where it crosses the disc's edge rather than at the next vertex left
+ * standing. On a densely sampled outline the two agree; on a sparse one (a
+ * traced box is four corners) dropping a corner would otherwise cut straight to
+ * the neighbouring corners. A `closed` ring wraps; an open one (a polyline) has
+ * no edge from its last vertex back to its first. Returns the trimmed parallel
+ * arrays plus how many of the ring's own vertices were removed (0: unchanged).
+ */
+export function eraseRingWithinRadius(xpoints, ypoints, cx, cy, radius, closed = true) {
+    const n = xpoints.length;
+    // A vertex placed on the disc's edge by an earlier tick stays outside (no churn).
+    const r2 = radius * radius * (1 - 1e-9);
+    const inside = new Array(n);
+    let removed = 0;
+    for (let i = 0; i < n; i++) {
+        const dx = xpoints[i] - cx;
+        const dy = ypoints[i] - cy;
+        inside[i] = dx * dx + dy * dy < r2;
+        if (inside[i])
+            removed++;
+    }
+    if (removed === 0)
+        return { xpoints: xpoints.slice(), ypoints: ypoints.slice(), removed };
+    const xs = [];
+    const ys = [];
+    for (let i = 0; i < n; i++) {
+        if (!inside[i]) {
+            xs.push(xpoints[i]);
+            ys.push(ypoints[i]);
+        }
+        const j = i + 1 < n ? i + 1 : closed ? 0 : -1;
+        if (j < 0 || inside[i] === inside[j])
+            continue;
+        // The edge i→j crosses the disc's edge once: leaving it when i is inside.
+        const ax = xpoints[i] - cx, ay = ypoints[i] - cy;
+        const ex = xpoints[j] - xpoints[i], ey = ypoints[j] - ypoints[i];
+        const a = ex * ex + ey * ey;
+        const b = ax * ex + ay * ey;
+        const c = ax * ax + ay * ay - radius * radius;
+        const root = Math.sqrt(Math.max(0, b * b - a * c));
+        const t = (inside[i] ? -b + root : -b - root) / a;
+        // A crossing at the vertex left standing (one on the disc's edge) is that vertex.
+        if (inside[i] ? t >= 1 - 1e-12 : t <= 1e-12)
+            continue;
+        xs.push(xpoints[i] + t * ex);
+        ys.push(ypoints[i] + t * ey);
+    }
+    return { xpoints: xs, ypoints: ys, removed };
+}
+/** Perpendicular distance from (px,py) to the segment (ax,ay)–(bx,by). */
+function perpDistance(px, py, ax, ay, bx, by) {
+    const dx = bx - ax, dy = by - ay;
+    const len2 = dx * dx + dy * dy;
+    if (len2 === 0)
+        return Math.hypot(px - ax, py - ay);
+    let t = ((px - ax) * dx + (py - ay) * dy) / len2;
+    t = Math.max(0, Math.min(1, t));
+    return Math.hypot(px - (ax + t * dx), py - (ay + t * dy));
+}
+/**
+ * Simplify a closed ring with the Douglas–Peucker algorithm (jit-ui#85):
+ * drop every vertex that lies within `epsilon` pixels (the "altitude
+ * threshold") of the line between its kept neighbours. Anchored at vertex 0
+ * (the ring is closed by appending it). Rings of ≤ 3 vertices, or a
+ * non-positive epsilon, are returned unchanged. Output keeps the closed-ring
+ * convention (no repeated closing point).
+ */
+export function simplifyRing(xs, ys, epsilon) {
+    const n = xs.length;
+    if (n <= 3 || !(epsilon > 0))
+        return { xs: xs.slice(), ys: ys.slice() };
+    // Close the ring so both endpoints of the DP run are vertex 0.
+    const px = [...xs, xs[0]];
+    const py = [...ys, ys[0]];
+    const m = px.length;
+    const keep = new Array(m).fill(false);
+    keep[0] = true;
+    keep[m - 1] = true;
+    const stack = [[0, m - 1]];
+    while (stack.length) {
+        const [a, b] = stack.pop();
+        let maxD = -1, idx = -1;
+        for (let i = a + 1; i < b; i++) {
+            const d = perpDistance(px[i], py[i], px[a], py[a], px[b], py[b]);
+            if (d > maxD) {
+                maxD = d;
+                idx = i;
+            }
+        }
+        if (maxD > epsilon && idx > a) {
+            keep[idx] = true;
+            stack.push([a, idx], [idx, b]);
+        }
+    }
+    const outX = [], outY = [];
+    for (let i = 0; i < m - 1; i++) {
+        // drop the duplicated closing vertex
+        if (keep[i]) {
+            outX.push(px[i]);
+            outY.push(py[i]);
+        }
+    }
+    return { xs: outX, ys: outY };
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoicmluZy5qcyIsInNvdXJjZVJvb3QiOiIiLCJzb3VyY2VzIjpbIi4uLy4uLy4uLy4uL3NyYy9saWIvZ2VvbWV0cnkvcmluZy50cyJdLCJuYW1lcyI6W10sIm1hcHBpbmdzIjoiQUFBQTs7Ozs7O0dBTUc7QUFVSCxvREFBb0Q7QUFDcEQsTUFBTSxVQUFVLFVBQVUsQ0FBQyxFQUFxQixFQUFFLEVBQXFCO0lBQ3JFLElBQUksSUFBSSxHQUFHLFFBQVEsRUFDakIsSUFBSSxHQUFHLENBQUMsUUFBUSxFQUNoQixJQUFJLEdBQUcsUUFBUSxFQUNmLElBQUksR0FBRyxDQUFDLFFBQVEsQ0FBQztJQUNuQixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsRUFBRSxDQUFDLE1BQU0sRUFBRSxDQUFDLEVBQUUsRUFBRSxDQUFDO1FBQ25DLElBQUksRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUk7WUFBRSxJQUFJLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQy9CLElBQUksRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUk7WUFBRSxJQUFJLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQy9CLElBQUksRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUk7WUFBRSxJQUFJLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQy9CLElBQUksRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUk7WUFBRSxJQUFJLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ2pDLENBQUM7SUFDRCxPQUFPLEVBQUUsSUFBSSxFQUFFLElBQUksRUFBRSxJQUFJLEVBQUUsSUFBSSxFQUFFLENBQUM7QUFDcEMsQ0FBQztBQUVEOzs7O0dBSUc7QUFDSCxNQUFNLFVBQVUsV0FBVyxDQUFDLEVBQVUsRUFBRSxFQUFVLEVBQUUsRUFBcUIsRUFBRSxFQUFxQjtJQUM5RixNQUFNLENBQUMsR0FBRyxFQUFFLENBQUMsTUFBTSxDQUFDO0lBQ3BCLElBQUksQ0FBQyxHQUFHLENBQUM7UUFBRSxPQUFPLEtBQUssQ0FBQztJQUN4QixJQUFJLE1BQU0sR0FBRyxLQUFLLENBQUM7SUFDbkIsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxFQUFFLEVBQUUsQ0FBQztRQUMxQyxNQUFNLEVBQUUsR0FBRyxFQUFFLENBQUMsQ0FBQyxDQUFDLEVBQ2QsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNiLE1BQU0sRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDLENBQUMsRUFDZCxFQUFFLEdBQUcsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ2IsTUFBTSxTQUFTLEdBQUcsRUFBRSxHQUFHLEVBQUUsS0FBSyxFQUFFLEdBQUcsRUFBRSxJQUFJLEVBQUUsR0FBRyxDQUFDLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQyxHQUFHLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLEdBQUcsRUFBRSxDQUFDO1FBQ3ZGLElBQUksU0FBUztZQUFFLE1BQU0sR0FBRyxDQUFDLE1BQU0sQ0FBQztJQUNsQyxDQUFDO0lBQ0QsT0FBTyxNQUFNLENBQUM7QUFDaEIsQ0FBQztBQUVEOzs7O0dBSUc7QUFDSCxNQUFNLFVBQVUsdUJBQXVCLENBQ3JDLEVBQVUsRUFDVixFQUFVLEVBQ1YsRUFBcUIsRUFDckIsRUFBcUIsRUFDckIsS0FBb0I7SUFFcEIsSUFBSSxDQUFDLFdBQVcsQ0FBQyxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLENBQUM7UUFBRSxPQUFPLEtBQUssQ0FBQztJQUMvQyxJQUFJLEtBQUssRUFBRSxDQUFDO1FBQ1YsS0FBSyxNQUFNLElBQUksSUFBSSxLQUFLLEVBQUUsQ0FBQztZQUN6QixJQUFJLGVBQWUsQ0FBQyxFQUFFLEVBQUUsRUFBRSxFQUFFLElBQUksQ0FBQztnQkFBRSxPQUFPLEtBQUssQ0FBQztRQUNsRCxDQUFDO0lBQ0gsQ0FBQztJQUNELE9BQU8sSUFBSSxDQUFDO0FBQ2QsQ0FBQztBQUVELGlHQUFpRztBQUNqRyxNQUFNLFVBQVUsZUFBZSxDQUFDLEVBQVUsRUFBRSxFQUFVLEVBQUUsSUFBc0M7SUFDNUYsTUFBTSxDQUFDLEdBQUcsSUFBSSxDQUFDLE1BQU0sQ0FBQztJQUN0QixJQUFJLENBQUMsR0FBRyxDQUFDO1FBQUUsT0FBTyxLQUFLLENBQUM7SUFDeEIsSUFBSSxNQUFNLEdBQUcsS0FBSyxDQUFDO0lBQ25CLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsRUFBRSxFQUFFLENBQUM7UUFDMUMsTUFBTSxFQUFFLEdBQUcsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxFQUNuQixFQUFFLEdBQUcsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ2xCLE1BQU0sRUFBRSxHQUFHLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsRUFDbkIsRUFBRSxHQUFHLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNsQixNQUFNLFNBQVMsR0FBRyxFQUFFLEdBQUcsRUFBRSxLQUFLLEVBQUUsR0FBRyxFQUFFLElBQUksRUFBRSxHQUFHLENBQUMsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLEdBQUcsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLENBQUMsR0FBRyxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsR0FBRyxFQUFFLENBQUM7UUFDdkYsSUFBSSxTQUFTO1lBQUUsTUFBTSxHQUFHLENBQUMsTUFBTSxDQUFDO0lBQ2xDLENBQUM7SUFDRCxPQUFPLE1BQU0sQ0FBQztBQUNoQixDQUFDO0FBRUQ7Ozs7R0FJRztBQUNILE1BQU0sVUFBVSx3QkFBd0IsQ0FDdEMsT0FBaUIsRUFDakIsT0FBaUIsRUFDakIsRUFBVSxFQUNWLEVBQVUsRUFDVixNQUFjO0lBRWQsTUFBTSxFQUFFLEdBQUcsTUFBTSxHQUFHLE1BQU0sQ0FBQztJQUMzQixNQUFNLEVBQUUsR0FBYSxFQUFFLENBQUM7SUFDeEIsTUFBTSxFQUFFLEdBQWEsRUFBRSxDQUFDO0lBQ3hCLElBQUksT0FBTyxHQUFHLENBQUMsQ0FBQztJQUNoQixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsT0FBTyxDQUFDLE1BQU0sRUFBRSxDQUFDLEVBQUUsRUFBRSxDQUFDO1FBQ3hDLE1BQU0sRUFBRSxHQUFHLE9BQU8sQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUM7UUFDM0IsTUFBTSxFQUFFLEdBQUcsT0FBTyxDQUFDLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQztRQUMzQixJQUFJLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxFQUFFLEVBQUUsQ0FBQztZQUMzQixPQUFPLEVBQUUsQ0FBQztZQUNWLFNBQVM7UUFDWCxDQUFDO1FBQ0QsRUFBRSxDQUFDLElBQUksQ0FBQyxPQUFPLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNwQixFQUFFLENBQUMsSUFBSSxDQUFDLE9BQU8sQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ3RCLENBQUM7SUFDRCxPQUFPLEVBQUUsT0FBTyxFQUFFLEVBQUUsRUFBRSxPQUFPLEVBQUUsRUFBRSxFQUFFLE9BQU8sRUFBRSxDQUFDO0FBQy9DLENBQUM7QUFFRDs7Ozs7Ozs7O0dBU0c7QUFDSCxNQUFNLFVBQVUscUJBQXFCLENBQ25DLE9BQWlCLEVBQ2pCLE9BQWlCLEVBQ2pCLEVBQVUsRUFDVixFQUFVLEVBQ1YsTUFBYyxFQUNkLE1BQU0sR0FBRyxJQUFJO0lBRWIsTUFBTSxDQUFDLEdBQUcsT0FBTyxDQUFDLE1BQU0sQ0FBQztJQUN6QixrRkFBa0Y7SUFDbEYsTUFBTSxFQUFFLEdBQUcsTUFBTSxHQUFHLE1BQU0sR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsQ0FBQztJQUN4QyxNQUFNLE1BQU0sR0FBRyxJQUFJLEtBQUssQ0FBVSxDQUFDLENBQUMsQ0FBQztJQUNyQyxJQUFJLE9BQU8sR0FBRyxDQUFDLENBQUM7SUFDaEIsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEVBQUUsRUFBRSxDQUFDO1FBQzNCLE1BQU0sRUFBRSxHQUFHLE9BQU8sQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUM7UUFDM0IsTUFBTSxFQUFFLEdBQUcsT0FBTyxDQUFDLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQztRQUMzQixNQUFNLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQztRQUNuQyxJQUFJLE1BQU0sQ0FBQyxDQUFDLENBQUM7WUFBRSxPQUFPLEVBQUUsQ0FBQztJQUMzQixDQUFDO0lBQ0QsSUFBSSxPQUFPLEtBQUssQ0FBQztRQUFFLE9BQU8sRUFBRSxPQUFPLEVBQUUsT0FBTyxDQUFDLEtBQUssRUFBRSxFQUFFLE9BQU8sRUFBRSxPQUFPLENBQUMsS0FBSyxFQUFFLEVBQUUsT0FBTyxFQUFFLENBQUM7SUFDMUYsTUFBTSxFQUFFLEdBQWEsRUFBRSxDQUFDO0lBQ3hCLE1BQU0sRUFBRSxHQUFhLEVBQUUsQ0FBQztJQUN4QixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7UUFDM0IsSUFBSSxDQUFDLE1BQU0sQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDO1lBQ2YsRUFBRSxDQUFDLElBQUksQ0FBQyxPQUFPLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztZQUNwQixFQUFFLENBQUMsSUFBSSxDQUFDLE9BQU8sQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ3RCLENBQUM7UUFDRCxNQUFNLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQzlDLElBQUksQ0FBQyxHQUFHLENBQUMsSUFBSSxNQUFNLENBQUMsQ0FBQyxDQUFDLEtBQUssTUFBTSxDQUFDLENBQUMsQ0FBQztZQUFFLFNBQVM7UUFDL0MsMEVBQTBFO1FBQzFFLE1BQU0sRUFBRSxHQUFHLE9BQU8sQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFLEVBQ3hCLEVBQUUsR0FBRyxPQUFPLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxDQUFDO1FBQ3ZCLE1BQU0sRUFBRSxHQUFHLE9BQU8sQ0FBQyxDQUFDLENBQUMsR0FBRyxPQUFPLENBQUMsQ0FBQyxDQUFDLEVBQ2hDLEVBQUUsR0FBRyxPQUFPLENBQUMsQ0FBQyxDQUFDLEdBQUcsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQy9CLE1BQU0sQ0FBQyxHQUFHLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQztRQUM1QixNQUFNLENBQUMsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxFQUFFLENBQUM7UUFDNUIsTUFBTSxDQUFDLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLE1BQU0sR0FBRyxNQUFNLENBQUM7UUFDOUMsTUFBTSxJQUFJLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ25ELE1BQU0sQ0FBQyxHQUFHLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQztRQUNsRCxrRkFBa0Y7UUFDbEYsSUFBSSxNQUFNLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLEdBQUcsS0FBSyxDQUFDLENBQUMsQ0FBQyxDQUFDLElBQUksS0FBSztZQUFFLFNBQVM7UUFDdEQsRUFBRSxDQUFDLElBQUksQ0FBQyxPQUFPLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxDQUFDO1FBQzdCLEVBQUUsQ0FBQyxJQUFJLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxFQUFFLENBQUMsQ0FBQztJQUMvQixDQUFDO0lBQ0QsT0FBTyxFQUFFLE9BQU8sRUFBRSxFQUFFLEVBQUUsT0FBTyxFQUFFLEVBQUUsRUFBRSxPQUFPLEVBQUUsQ0FBQztBQUMvQyxDQUFDO0FBRUQsMEVBQTBFO0FBQzFFLFNBQVMsWUFBWSxDQUFDLEVBQVUsRUFBRSxFQUFVLEVBQUUsRUFBVSxFQUFFLEVBQVUsRUFBRSxFQUFVLEVBQUUsRUFBVTtJQUMxRixNQUFNLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxFQUNoQixFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQztJQUNmLE1BQU0sSUFBSSxHQUFHLEVBQUUsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQztJQUMvQixJQUFJLElBQUksS0FBSyxDQUFDO1FBQUUsT0FBTyxJQUFJLENBQUMsS0FBSyxDQUFDLEVBQUUsR0FBRyxFQUFFLEVBQUUsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDO0lBQ3BELElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLEdBQUcsRUFBRSxHQUFHLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxHQUFHLElBQUksQ0FBQztJQUNqRCxDQUFDLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUNoQyxPQUFPLElBQUksQ0FBQyxLQUFLLENBQUMsRUFBRSxHQUFHLENBQUMsRUFBRSxHQUFHLENBQUMsR0FBRyxFQUFFLENBQUMsRUFBRSxFQUFFLEdBQUcsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxDQUFDLENBQUM7QUFDNUQsQ0FBQztBQUVEOzs7Ozs7O0dBT0c7QUFDSCxNQUFNLFVBQVUsWUFBWSxDQUFDLEVBQVksRUFBRSxFQUFZLEVBQUUsT0FBZTtJQUN0RSxNQUFNLENBQUMsR0FBRyxFQUFFLENBQUMsTUFBTSxDQUFDO0lBQ3BCLElBQUksQ0FBQyxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUMsT0FBTyxHQUFHLENBQUMsQ0FBQztRQUFFLE9BQU8sRUFBRSxFQUFFLEVBQUUsRUFBRSxDQUFDLEtBQUssRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLENBQUMsS0FBSyxFQUFFLEVBQUUsQ0FBQztJQUV4RSwrREFBK0Q7SUFDL0QsTUFBTSxFQUFFLEdBQUcsQ0FBQyxHQUFHLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUMxQixNQUFNLEVBQUUsR0FBRyxDQUFDLEdBQUcsRUFBRSxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQzFCLE1BQU0sQ0FBQyxHQUFHLEVBQUUsQ0FBQyxNQUFNLENBQUM7SUFDcEIsTUFBTSxJQUFJLEdBQUcsSUFBSSxLQUFLLENBQVUsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLEtBQUssQ0FBQyxDQUFDO0lBQy9DLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUM7SUFDZixJQUFJLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQztJQUVuQixNQUFNLEtBQUssR0FBNEIsQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUNwRCxPQUFPLEtBQUssQ0FBQyxNQUFNLEVBQUUsQ0FBQztRQUNwQixNQUFNLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQyxHQUFHLEtBQUssQ0FBQyxHQUFHLEVBQXNCLENBQUM7UUFDL0MsSUFBSSxJQUFJLEdBQUcsQ0FBQyxDQUFDLEVBQ1gsR0FBRyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ1gsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQztZQUMvQixNQUFNLENBQUMsR0FBRyxZQUFZLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztZQUNqRSxJQUFJLENBQUMsR0FBRyxJQUFJLEVBQUUsQ0FBQztnQkFDYixJQUFJLEdBQUcsQ0FBQyxDQUFDO2dCQUNULEdBQUcsR0FBRyxDQUFDLENBQUM7WUFDVixDQUFDO1FBQ0gsQ0FBQztRQUNELElBQUksSUFBSSxHQUFHLE9BQU8sSUFBSSxHQUFHLEdBQUcsQ0FBQyxFQUFFLENBQUM7WUFDOUIsSUFBSSxDQUFDLEdBQUcsQ0FBQyxHQUFHLElBQUksQ0FBQztZQUNqQixLQUFLLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDakMsQ0FBQztJQUNILENBQUM7SUFFRCxNQUFNLElBQUksR0FBYSxFQUFFLEVBQ3ZCLElBQUksR0FBYSxFQUFFLENBQUM7SUFDdEIsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQztRQUMvQixxQ0FBcUM7UUFDckMsSUFBSSxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQztZQUNaLElBQUksQ0FBQyxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDakIsSUFBSSxDQUFDLElBQUksQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNuQixDQUFDO0lBQ0gsQ0FBQztJQUNELE9BQU8sRUFBRSxFQUFFLEVBQUUsSUFBSSxFQUFFLEVBQUUsRUFBRSxJQUFJLEVBQUUsQ0FBQztBQUNoQyxDQUFDIiwic291cmNlc0NvbnRlbnQiOlsiLyoqXG4gKiBQdXJlIHJpbmcgKGNsb3NlZCBwb2x5bGluZSkgaGVscGVyczogYm91bmRzLCBjb250YWlubWVudCwgdmVydGV4IGRyb3BwaW5nIGFuZFxuICogc2ltcGxpZmljYXRpb24uIE5vIERPTSwgbm8gQW5ndWxhciwgbm8gYHRoaXNgIOKAlCB1c2FibGUgZnJvbSB3b3JrZXJzLlxuICpcbiAqIEEgcmluZyBpcyBhIHBhaXIgb2YgcGFyYWxsZWwgYHhzYC9geXNgIGFycmF5cyB3aXRob3V0IGEgcmVwZWF0ZWQgY2xvc2luZ1xuICogdmVydGV4OyBob2xlIHJpbmdzIHVzZSB0aGUgYFtbeCwgeV0sIOKApl1gIGNvbnZlbnRpb24gb2Yge0BsaW5rIFBvbHlnb24uaG9sZXN9LlxuICovXG5cbi8qKiBBeGlzLWFsaWduZWQgYm91bmRzIG9mIGEgcmluZy4gSW5maW5pdGUgKG1pbiA+IG1heCkgZm9yIGFuIGVtcHR5IHJpbmcuICovXG5leHBvcnQgaW50ZXJmYWNlIFJpbmdCb3VuZHMge1xuICBtaW5YOiBudW1iZXI7XG4gIG1pblk6IG51bWJlcjtcbiAgbWF4WDogbnVtYmVyO1xuICBtYXhZOiBudW1iZXI7XG59XG5cbi8qKiBUaGUgYXhpcy1hbGlnbmVkIGJvdW5kcyBvZiBhIHJpbmcncyB2ZXJ0aWNlcy4gKi9cbmV4cG9ydCBmdW5jdGlvbiByaW5nQm91bmRzKHhzOiBBcnJheUxpa2U8bnVtYmVyPiwgeXM6IEFycmF5TGlrZTxudW1iZXI+KTogUmluZ0JvdW5kcyB7XG4gIGxldCBtaW5YID0gSW5maW5pdHksXG4gICAgbWF4WCA9IC1JbmZpbml0eSxcbiAgICBtaW5ZID0gSW5maW5pdHksXG4gICAgbWF4WSA9IC1JbmZpbml0eTtcbiAgZm9yIChsZXQgaSA9IDA7IGkgPCB4cy5sZW5ndGg7IGkrKykge1xuICAgIGlmICh4c1tpXSA8IG1pblgpIG1pblggPSB4c1tpXTtcbiAgICBpZiAoeHNbaV0gPiBtYXhYKSBtYXhYID0geHNbaV07XG4gICAgaWYgKHlzW2ldIDwgbWluWSkgbWluWSA9IHlzW2ldO1xuICAgIGlmICh5c1tpXSA+IG1heFkpIG1heFkgPSB5c1tpXTtcbiAgfVxuICByZXR1cm4geyBtaW5YLCBtaW5ZLCBtYXhYLCBtYXhZIH07XG59XG5cbi8qKlxuICogU3RhbmRhcmQgcmF5LWNhc3QgKGV2ZW4tb2RkKSBwb2ludC1pbi1yaW5nIHRlc3QuIChweCwgcHkpIGFuZCB0aGUgcmluZ1xuICogdmVydGljZXMgbXVzdCBiZSBpbiB0aGUgc2FtZSBjb29yZGluYXRlIHN5c3RlbS4gUmluZ3Mgd2l0aCBmZXdlciB0aGFuIHRocmVlXG4gKiB2ZXJ0aWNlcyBjb250YWluIG5vdGhpbmcuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBwb2ludEluUmluZyhweDogbnVtYmVyLCBweTogbnVtYmVyLCB4czogQXJyYXlMaWtlPG51bWJlcj4sIHlzOiBBcnJheUxpa2U8bnVtYmVyPik6IGJvb2xlYW4ge1xuICBjb25zdCBuID0geHMubGVuZ3RoO1xuICBpZiAobiA8IDMpIHJldHVybiBmYWxzZTtcbiAgbGV0IGluc2lkZSA9IGZhbHNlO1xuICBmb3IgKGxldCBpID0gMCwgaiA9IG4gLSAxOyBpIDwgbjsgaiA9IGkrKykge1xuICAgIGNvbnN0IHlpID0geXNbaV0sXG4gICAgICB5aiA9IHlzW2pdO1xuICAgIGNvbnN0IHhpID0geHNbaV0sXG4gICAgICB4aiA9IHhzW2pdO1xuICAgIGNvbnN0IGludGVyc2VjdCA9IHlpID4gcHkgIT09IHlqID4gcHkgJiYgcHggPCAoKHhqIC0geGkpICogKHB5IC0geWkpKSAvICh5aiAtIHlpKSArIHhpO1xuICAgIGlmIChpbnRlcnNlY3QpIGluc2lkZSA9ICFpbnNpZGU7XG4gIH1cbiAgcmV0dXJuIGluc2lkZTtcbn1cblxuLyoqXG4gKiBQb2ludC1pbi1yZWdpb24gdGVzdCB0aGF0IGhvbm91cnMgaW50ZXJpb3IgcmluZ3MgKGhvbGVzKTogaW5zaWRlIHRoZVxuICogZXh0ZXJpb3IgQU5EIG91dHNpZGUgZXZlcnkgaG9sZS4gQ29vcmRpbmF0ZXMgbXVzdCBzaGFyZSB0aGUgcG9seWdvbidzIGZyYW1lLlxuICogV2l0aCBubyBob2xlcyB0aGlzIGlzIGV4YWN0bHkge0BsaW5rIHBvaW50SW5SaW5nfS5cbiAqL1xuZXhwb3J0IGZ1bmN0aW9uIHBvaW50SW5Qb2x5Z29uV2l0aEhvbGVzKFxuICBweDogbnVtYmVyLFxuICBweTogbnVtYmVyLFxuICB4czogQXJyYXlMaWtlPG51bWJlcj4sXG4gIHlzOiBBcnJheUxpa2U8bnVtYmVyPixcbiAgaG9sZXM/OiBudW1iZXJbXVtdW10sXG4pOiBib29sZWFuIHtcbiAgaWYgKCFwb2ludEluUmluZyhweCwgcHksIHhzLCB5cykpIHJldHVybiBmYWxzZTtcbiAgaWYgKGhvbGVzKSB7XG4gICAgZm9yIChjb25zdCByaW5nIG9mIGhvbGVzKSB7XG4gICAgICBpZiAocG9pbnRJblBhaXJSaW5nKHB4LCBweSwgcmluZykpIHJldHVybiBmYWxzZTtcbiAgICB9XG4gIH1cbiAgcmV0dXJuIHRydWU7XG59XG5cbi8qKiB7QGxpbmsgcG9pbnRJblJpbmd9IG92ZXIgYSBgW1t4LCB5XSwg4oCmXWAgcmluZyAodGhlIGhvbGUgY29udmVudGlvbiksIHdpdGhvdXQgc3BsaXR0aW5nIGl0LiAqL1xuZXhwb3J0IGZ1bmN0aW9uIHBvaW50SW5QYWlyUmluZyhweDogbnVtYmVyLCBweTogbnVtYmVyLCByaW5nOiBSZWFkb25seUFycmF5PEFycmF5TGlrZTxudW1iZXI+Pik6IGJvb2xlYW4ge1xuICBjb25zdCBuID0gcmluZy5sZW5ndGg7XG4gIGlmIChuIDwgMykgcmV0dXJuIGZhbHNlO1xuICBsZXQgaW5zaWRlID0gZmFsc2U7XG4gIGZvciAobGV0IGkgPSAwLCBqID0gbiAtIDE7IGkgPCBuOyBqID0gaSsrKSB7XG4gICAgY29uc3QgeGkgPSByaW5nW2ldWzBdLFxuICAgICAgeWkgPSByaW5nW2ldWzFdO1xuICAgIGNvbnN0IHhqID0gcmluZ1tqXVswXSxcbiAgICAgIHlqID0gcmluZ1tqXVsxXTtcbiAgICBjb25zdCBpbnRlcnNlY3QgPSB5aSA+IHB5ICE9PSB5aiA+IHB5ICYmIHB4IDwgKCh4aiAtIHhpKSAqIChweSAtIHlpKSkgLyAoeWogLSB5aSkgKyB4aTtcbiAgICBpZiAoaW50ZXJzZWN0KSBpbnNpZGUgPSAhaW5zaWRlO1xuICB9XG4gIHJldHVybiBpbnNpZGU7XG59XG5cbi8qKlxuICogRHJvcCBldmVyeSB2ZXJ0ZXggd2hvc2UgRXVjbGlkZWFuIGRpc3RhbmNlIHRvIChjeCwgY3kpIGlzIGxlc3MgdGhhblxuICogYHJhZGl1c2AuIFJldHVybnMgdHJpbW1lZCBwYXJhbGxlbCBhcnJheXMgcGx1cyBob3cgbWFueSB2ZXJ0aWNlcyB3ZXJlXG4gKiByZW1vdmVkIChzbyB0aGUgY2FsbGVyIGNhbiBkZXRlY3QgYSBuby1vcCB0aWNrIGFuZCBza2lwIHRoZSByZWxheW91dCkuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBkcm9wVmVydGljZXNXaXRoaW5SYWRpdXMoXG4gIHhwb2ludHM6IG51bWJlcltdLFxuICB5cG9pbnRzOiBudW1iZXJbXSxcbiAgY3g6IG51bWJlcixcbiAgY3k6IG51bWJlcixcbiAgcmFkaXVzOiBudW1iZXIsXG4pOiB7IHhwb2ludHM6IG51bWJlcltdOyB5cG9pbnRzOiBudW1iZXJbXTsgcmVtb3ZlZDogbnVtYmVyIH0ge1xuICBjb25zdCByMiA9IHJhZGl1cyAqIHJhZGl1cztcbiAgY29uc3QgeHM6IG51bWJlcltdID0gW107XG4gIGNvbnN0IHlzOiBudW1iZXJbXSA9IFtdO1xuICBsZXQgcmVtb3ZlZCA9IDA7XG4gIGZvciAobGV0IGkgPSAwOyBpIDwgeHBvaW50cy5sZW5ndGg7IGkrKykge1xuICAgIGNvbnN0IGR4ID0geHBvaW50c1tpXSAtIGN4O1xuICAgIGNvbnN0IGR5ID0geXBvaW50c1tpXSAtIGN5O1xuICAgIGlmIChkeCAqIGR4ICsgZHkgKiBkeSA8IHIyKSB7XG4gICAgICByZW1vdmVkKys7XG4gICAgICBjb250aW51ZTtcbiAgICB9XG4gICAgeHMucHVzaCh4cG9pbnRzW2ldKTtcbiAgICB5cy5wdXNoKHlwb2ludHNbaV0pO1xuICB9XG4gIHJldHVybiB7IHhwb2ludHM6IHhzLCB5cG9pbnRzOiB5cywgcmVtb3ZlZCB9O1xufVxuXG4vKipcbiAqIEVyYXNlIHRoZSBwYXJ0IG9mIGEgcmluZyB0aGF0IGxpZXMgd2l0aGluIGByYWRpdXNgIG9mIChjeCwgY3kpOiBldmVyeSB2ZXJ0ZXhcbiAqIGluc2lkZSB0aGUgZGlzYyBnb2VzIChhcyB7QGxpbmsgZHJvcFZlcnRpY2VzV2l0aGluUmFkaXVzfSksIGFuZCB0aGUgb3V0bGluZSBpc1xuICogcmVqb2luZWQgd2hlcmUgaXQgY3Jvc3NlcyB0aGUgZGlzYydzIGVkZ2UgcmF0aGVyIHRoYW4gYXQgdGhlIG5leHQgdmVydGV4IGxlZnRcbiAqIHN0YW5kaW5nLiBPbiBhIGRlbnNlbHkgc2FtcGxlZCBvdXRsaW5lIHRoZSB0d28gYWdyZWU7IG9uIGEgc3BhcnNlIG9uZSAoYVxuICogdHJhY2VkIGJveCBpcyBmb3VyIGNvcm5lcnMpIGRyb3BwaW5nIGEgY29ybmVyIHdvdWxkIG90aGVyd2lzZSBjdXQgc3RyYWlnaHQgdG9cbiAqIHRoZSBuZWlnaGJvdXJpbmcgY29ybmVycy4gQSBgY2xvc2VkYCByaW5nIHdyYXBzOyBhbiBvcGVuIG9uZSAoYSBwb2x5bGluZSkgaGFzXG4gKiBubyBlZGdlIGZyb20gaXRzIGxhc3QgdmVydGV4IGJhY2sgdG8gaXRzIGZpcnN0LiBSZXR1cm5zIHRoZSB0cmltbWVkIHBhcmFsbGVsXG4gKiBhcnJheXMgcGx1cyBob3cgbWFueSBvZiB0aGUgcmluZydzIG93biB2ZXJ0aWNlcyB3ZXJlIHJlbW92ZWQgKDA6IHVuY2hhbmdlZCkuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBlcmFzZVJpbmdXaXRoaW5SYWRpdXMoXG4gIHhwb2ludHM6IG51bWJlcltdLFxuICB5cG9pbnRzOiBudW1iZXJbXSxcbiAgY3g6IG51bWJlcixcbiAgY3k6IG51bWJlcixcbiAgcmFkaXVzOiBudW1iZXIsXG4gIGNsb3NlZCA9IHRydWUsXG4pOiB7IHhwb2ludHM6IG51bWJlcltdOyB5cG9pbnRzOiBudW1iZXJbXTsgcmVtb3ZlZDogbnVtYmVyIH0ge1xuICBjb25zdCBuID0geHBvaW50cy5sZW5ndGg7XG4gIC8vIEEgdmVydGV4IHBsYWNlZCBvbiB0aGUgZGlzYydzIGVkZ2UgYnkgYW4gZWFybGllciB0aWNrIHN0YXlzIG91dHNpZGUgKG5vIGNodXJuKS5cbiAgY29uc3QgcjIgPSByYWRpdXMgKiByYWRpdXMgKiAoMSAtIDFlLTkpO1xuICBjb25zdCBpbnNpZGUgPSBuZXcgQXJyYXk8Ym9vbGVhbj4obik7XG4gIGxldCByZW1vdmVkID0gMDtcbiAgZm9yIChsZXQgaSA9IDA7IGkgPCBuOyBpKyspIHtcbiAgICBjb25zdCBkeCA9IHhwb2ludHNbaV0gLSBjeDtcbiAgICBjb25zdCBkeSA9IHlwb2ludHNbaV0gLSBjeTtcbiAgICBpbnNpZGVbaV0gPSBkeCAqIGR4ICsgZHkgKiBkeSA8IHIyO1xuICAgIGlmIChpbnNpZGVbaV0pIHJlbW92ZWQrKztcbiAgfVxuICBpZiAocmVtb3ZlZCA9PT0gMCkgcmV0dXJuIHsgeHBvaW50czogeHBvaW50cy5zbGljZSgpLCB5cG9pbnRzOiB5cG9pbnRzLnNsaWNlKCksIHJlbW92ZWQgfTtcbiAgY29uc3QgeHM6IG51bWJlcltdID0gW107XG4gIGNvbnN0IHlzOiBudW1iZXJbXSA9IFtdO1xuICBmb3IgKGxldCBpID0gMDsgaSA8IG47IGkrKykge1xuICAgIGlmICghaW5zaWRlW2ldKSB7XG4gICAgICB4cy5wdXNoKHhwb2ludHNbaV0pO1xuICAgICAgeXMucHVzaCh5cG9pbnRzW2ldKTtcbiAgICB9XG4gICAgY29uc3QgaiA9IGkgKyAxIDwgbiA/IGkgKyAxIDogY2xvc2VkID8gMCA6IC0xO1xuICAgIGlmIChqIDwgMCB8fCBpbnNpZGVbaV0gPT09IGluc2lkZVtqXSkgY29udGludWU7XG4gICAgLy8gVGhlIGVkZ2UgaeKGkmogY3Jvc3NlcyB0aGUgZGlzYydzIGVkZ2Ugb25jZTogbGVhdmluZyBpdCB3aGVuIGkgaXMgaW5zaWRlLlxuICAgIGNvbnN0IGF4ID0geHBvaW50c1tpXSAtIGN4LFxuICAgICAgYXkgPSB5cG9pbnRzW2ldIC0gY3k7XG4gICAgY29uc3QgZXggPSB4cG9pbnRzW2pdIC0geHBvaW50c1tpXSxcbiAgICAgIGV5ID0geXBvaW50c1tqXSAtIHlwb2ludHNbaV07XG4gICAgY29uc3QgYSA9IGV4ICogZXggKyBleSAqIGV5O1xuICAgIGNvbnN0IGIgPSBheCAqIGV4ICsgYXkgKiBleTtcbiAgICBjb25zdCBjID0gYXggKiBheCArIGF5ICogYXkgLSByYWRpdXMgKiByYWRpdXM7XG4gICAgY29uc3Qgcm9vdCA9IE1hdGguc3FydChNYXRoLm1heCgwLCBiICogYiAtIGEgKiBjKSk7XG4gICAgY29uc3QgdCA9IChpbnNpZGVbaV0gPyAtYiArIHJvb3QgOiAtYiAtIHJvb3QpIC8gYTtcbiAgICAvLyBBIGNyb3NzaW5nIGF0IHRoZSB2ZXJ0ZXggbGVmdCBzdGFuZGluZyAob25lIG9uIHRoZSBkaXNjJ3MgZWRnZSkgaXMgdGhhdCB2ZXJ0ZXguXG4gICAgaWYgKGluc2lkZVtpXSA/IHQgPj0gMSAtIDFlLTEyIDogdCA8PSAxZS0xMikgY29udGludWU7XG4gICAgeHMucHVzaCh4cG9pbnRzW2ldICsgdCAqIGV4KTtcbiAgICB5cy5wdXNoKHlwb2ludHNbaV0gKyB0ICogZXkpO1xuICB9XG4gIHJldHVybiB7IHhwb2ludHM6IHhzLCB5cG9pbnRzOiB5cywgcmVtb3ZlZCB9O1xufVxuXG4vKiogUGVycGVuZGljdWxhciBkaXN0YW5jZSBmcm9tIChweCxweSkgdG8gdGhlIHNlZ21lbnQgKGF4LGF5KeKAkyhieCxieSkuICovXG5mdW5jdGlvbiBwZXJwRGlzdGFuY2UocHg6IG51bWJlciwgcHk6IG51bWJlciwgYXg6IG51bWJlciwgYXk6IG51bWJlciwgYng6IG51bWJlciwgYnk6IG51bWJlcik6IG51bWJlciB7XG4gIGNvbnN0IGR4ID0gYnggLSBheCxcbiAgICBkeSA9IGJ5IC0gYXk7XG4gIGNvbnN0IGxlbjIgPSBkeCAqIGR4ICsgZHkgKiBkeTtcbiAgaWYgKGxlbjIgPT09IDApIHJldHVybiBNYXRoLmh5cG90KHB4IC0gYXgsIHB5IC0gYXkpO1xuICBsZXQgdCA9ICgocHggLSBheCkgKiBkeCArIChweSAtIGF5KSAqIGR5KSAvIGxlbjI7XG4gIHQgPSBNYXRoLm1heCgwLCBNYXRoLm1pbigxLCB0KSk7XG4gIHJldHVybiBNYXRoLmh5cG90KHB4IC0gKGF4ICsgdCAqIGR4KSwgcHkgLSAoYXkgKyB0ICogZHkpKTtcbn1cblxuLyoqXG4gKiBTaW1wbGlmeSBhIGNsb3NlZCByaW5nIHdpdGggdGhlIERvdWdsYXPigJNQZXVja2VyIGFsZ29yaXRobSAoaml0LXVpIzg1KTpcbiAqIGRyb3AgZXZlcnkgdmVydGV4IHRoYXQgbGllcyB3aXRoaW4gYGVwc2lsb25gIHBpeGVscyAodGhlIFwiYWx0aXR1ZGVcbiAqIHRocmVzaG9sZFwiKSBvZiB0aGUgbGluZSBiZXR3ZWVuIGl0cyBrZXB0IG5laWdoYm91cnMuIEFuY2hvcmVkIGF0IHZlcnRleCAwXG4gKiAodGhlIHJpbmcgaXMgY2xvc2VkIGJ5IGFwcGVuZGluZyBpdCkuIFJpbmdzIG9mIOKJpCAzIHZlcnRpY2VzLCBvciBhXG4gKiBub24tcG9zaXRpdmUgZXBzaWxvbiwgYXJlIHJldHVybmVkIHVuY2hhbmdlZC4gT3V0cHV0IGtlZXBzIHRoZSBjbG9zZWQtcmluZ1xuICogY29udmVudGlvbiAobm8gcmVwZWF0ZWQgY2xvc2luZyBwb2ludCkuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBzaW1wbGlmeVJpbmcoeHM6IG51bWJlcltdLCB5czogbnVtYmVyW10sIGVwc2lsb246IG51bWJlcik6IHsgeHM6IG51bWJlcltdOyB5czogbnVtYmVyW10gfSB7XG4gIGNvbnN0IG4gPSB4cy5sZW5ndGg7XG4gIGlmIChuIDw9IDMgfHwgIShlcHNpbG9uID4gMCkpIHJldHVybiB7IHhzOiB4cy5zbGljZSgpLCB5czogeXMuc2xpY2UoKSB9O1xuXG4gIC8vIENsb3NlIHRoZSByaW5nIHNvIGJvdGggZW5kcG9pbnRzIG9mIHRoZSBEUCBydW4gYXJlIHZlcnRleCAwLlxuICBjb25zdCBweCA9IFsuLi54cywgeHNbMF1dO1xuICBjb25zdCBweSA9IFsuLi55cywgeXNbMF1dO1xuICBjb25zdCBtID0gcHgubGVuZ3RoO1xuICBjb25zdCBrZWVwID0gbmV3IEFycmF5PGJvb2xlYW4+KG0pLmZpbGwoZmFsc2UpO1xuICBrZWVwWzBdID0gdHJ1ZTtcbiAga2VlcFttIC0gMV0gPSB0cnVlO1xuXG4gIGNvbnN0IHN0YWNrOiBBcnJheTxbbnVtYmVyLCBudW1iZXJdPiA9IFtbMCwgbSAtIDFdXTtcbiAgd2hpbGUgKHN0YWNrLmxlbmd0aCkge1xuICAgIGNvbnN0IFthLCBiXSA9IHN0YWNrLnBvcCgpIGFzIFtudW1iZXIsIG51bWJlcl07XG4gICAgbGV0IG1heEQgPSAtMSxcbiAgICAgIGlkeCA9IC0xO1xuICAgIGZvciAobGV0IGkgPSBhICsgMTsgaSA8IGI7IGkrKykge1xuICAgICAgY29uc3QgZCA9IHBlcnBEaXN0YW5jZShweFtpXSwgcHlbaV0sIHB4W2FdLCBweVthXSwgcHhbYl0sIHB5W2JdKTtcbiAgICAgIGlmIChkID4gbWF4RCkge1xuICAgICAgICBtYXhEID0gZDtcbiAgICAgICAgaWR4ID0gaTtcbiAgICAgIH1cbiAgICB9XG4gICAgaWYgKG1heEQgPiBlcHNpbG9uICYmIGlkeCA+IGEpIHtcbiAgICAgIGtlZXBbaWR4XSA9IHRydWU7XG4gICAgICBzdGFjay5wdXNoKFthLCBpZHhdLCBbaWR4LCBiXSk7XG4gICAgfVxuICB9XG5cbiAgY29uc3Qgb3V0WDogbnVtYmVyW10gPSBbXSxcbiAgICBvdXRZOiBudW1iZXJbXSA9IFtdO1xuICBmb3IgKGxldCBpID0gMDsgaSA8IG0gLSAxOyBpKyspIHtcbiAgICAvLyBkcm9wIHRoZSBkdXBsaWNhdGVkIGNsb3NpbmcgdmVydGV4XG4gICAgaWYgKGtlZXBbaV0pIHtcbiAgICAgIG91dFgucHVzaChweFtpXSk7XG4gICAgICBvdXRZLnB1c2gocHlbaV0pO1xuICAgIH1cbiAgfVxuICByZXR1cm4geyB4czogb3V0WCwgeXM6IG91dFkgfTtcbn1cbiJdfQ==