@jax-data-science/sci-image-visualizer 0.8.3 → 0.10.0

package/esm2022/lib/geometry/reconstruct.mjs

@@ -0,0 +1,315 @@
+/**
+ * The 3D reconstruction of a region from its outlines on a stack of slices: each slice's outlines
+ * filled (each feature even-odd, the features as their union) into a voxel mask, the mask's
+ * boundary meshed by surface nets and smoothed; and its open lines lofted into a sheet.
+ *
+ * Pure: usable from workers.
+ */
+import { pointInRing, ringArea, smoothMesh, surfaceNets } from 'napari-js/geometry';
+import { loftLines } from './loft';
+/**
+ * The mesh of the region outlined by `slices`, in image pixels (x, y) and slice indices (z: slice
+ * k's plane at z = k): the filled outlines' surface, and the sheet through the lines (see
+ * {@link loftLines}). Empty when no slice has an area and no two slices have lines.
+ */
+export function reconstruct(slices, opts = {}) {
+    const filled = fill(slices, opts);
+    const lines = slices.filter((s) => s.lines?.length);
+    if (!lines.length)
+        return filled;
+    const sheet = loftLines(lines.map((s) => ({ z: s.z, lines: s.lines ?? [] })), opts.pixelStep ?? 1);
+    return mergeMeshes(filled, sheet);
+}
+/** The surface of the filled outlines. */
+function fill(slices, opts) {
+    let x0 = Infinity, y0 = Infinity, x1 = -Infinity, y1 = -Infinity, z0 = Infinity, z1 = -Infinity;
+    for (const s of slices) {
+        for (const r of s.rings) {
+            for (let i = 0; i + 1 < r.length; i += 2) {
+                if (r[i] < x0)
+                    x0 = r[i];
+                if (r[i] > x1)
+                    x1 = r[i];
+                if (r[i + 1] < y0)
+                    y0 = r[i + 1];
+                if (r[i + 1] > y1)
+                    y1 = r[i + 1];
+            }
+        }
+        if (s.rings.length) {
+            z0 = Math.min(z0, s.z);
+            z1 = Math.max(z1, s.z);
+        }
+    }
+    if (!(x1 > x0) || !(y1 > y0))
+        return { vertices: new Float32Array(0), faces: new Uint32Array(0) };
+    const nz = z1 - z0 + 1;
+    const maxVoxels = opts.maxVoxels ?? 32_000_000;
+    let step = Math.max(opts.pixelStep ?? 1, 1e-6);
+    while (Math.ceil((x1 - x0) / step) * Math.ceil((y1 - y0) / step) * nz > maxVoxels)
+        step *= 1.25;
+    const nx = Math.max(1, Math.ceil((x1 - x0) / step));
+    const ny = Math.max(1, Math.ceil((y1 - y0) / step));
+    const mask = new Uint8Array(nx * ny * nz);
+    for (const s of slices) {
+        for (const rings of featureRings(s.rings, s.ringGroups)) {
+            fillSlice(mask, nx, ny, (s.z - z0) * nx * ny, rings, x0, y0, step);
+        }
+    }
+    // Voxel i spanning [i, i + 1) (napari-js centres it on i): the scaling below maps cell edges to px.
+    const mesh = smoothMesh(surfaceNets(mask, [nx, ny, nz], { origin: [0.5, 0.5, 0.5] }), {
+        iterations: opts.smooth ?? 2,
+    });
+    const v = mesh.vertices;
+    for (let i = 0; i < v.length; i += 3) {
+        v[i] = x0 + v[i] * step;
+        v[i + 1] = y0 + v[i + 1] * step;
+        // Voxel k spans [k, k+1) with its centre (the slice) at k + 0.5.
+        v[i + 2] = z0 + v[i + 2] - 0.5;
+    }
+    return mesh;
+}
+/** `a` and `b` as one mesh. */
+function mergeMeshes(a, b) {
+    if (!b.faces.length)
+        return a;
+    if (!a.faces.length)
+        return b;
+    const vertices = new Float32Array(a.vertices.length + b.vertices.length);
+    vertices.set(a.vertices);
+    vertices.set(b.vertices, a.vertices.length);
+    const faces = new Uint32Array(a.faces.length + b.faces.length);
+    faces.set(a.faces);
+    const offset = a.vertices.length / 3;
+    for (let i = 0; i < b.faces.length; i++)
+        faces[a.faces.length + i] = b.faces[i] + offset;
+    return { vertices, faces };
+}
+/**
+ * Fill one feature into one slice of the mask: even-odd over its `rings`, sampled at voxel
+ * centres. Cells already set stay set, so the features of a slice fill as their union.
+ */
+function fillSlice(mask, nx, ny, offset, rings, x0, y0, step) {
+    const crossings = [];
+    // Only the rows the feature spans.
+    let top = Infinity, bottom = -Infinity;
+    for (const r of rings) {
+        for (let i = 1; i < r.length; i += 2) {
+            if (r[i] < top)
+                top = r[i];
+            if (r[i] > bottom)
+                bottom = r[i];
+        }
+    }
+    const j0 = Math.max(0, Math.floor((top - y0) / step - 0.5));
+    const j1 = Math.min(ny - 1, Math.ceil((bottom - y0) / step - 0.5));
+    for (let j = j0; j <= j1; j++) {
+        const y = y0 + (j + 0.5) * step;
+        crossings.length = 0;
+        for (const r of rings) {
+            const n = r.length >> 1;
+            for (let i = 0, k = n - 1; i < n; k = i++) {
+                const ya = r[2 * i + 1], yb = r[2 * k + 1];
+                if (ya > y !== yb > y) {
+                    const xa = r[2 * i], xb = r[2 * k];
+                    crossings.push(xa + ((y - ya) * (xb - xa)) / (yb - ya));
+                }
+            }
+        }
+        if (crossings.length < 2)
+            continue;
+        crossings.sort((a, b) => a - b);
+        const row = offset + j * nx;
+        for (let c = 0; c + 1 < crossings.length; c += 2) {
+            const from = Math.max(0, Math.ceil((crossings[c] - x0) / step - 0.5));
+            const to = Math.min(nx - 1, Math.floor((crossings[c + 1] - x0) / step - 0.5));
+            for (let i = from; i <= to; i++)
+                mask[row + i] = 1;
+        }
+    }
+}
+/**
+ * The area inside `rings` (interleaved `x, y`), even-odd: a ring inside an odd number of others is
+ * a hole. Exact for rings that do not cross (a polygon set's exteriors and holes).
+ */
+export function ringsArea(rings) {
+    const boxes = rings.map(ringBox);
+    let area = 0;
+    for (let i = 0; i < rings.length; i++) {
+        const r = rings[i];
+        if (r.length < 6)
+            continue;
+        const x = r[0], y = r[1];
+        let depth = 0;
+        for (let j = 0; j < rings.length; j++) {
+            if (j === i)
+                continue;
+            const b = boxes[j];
+            if (x < b[0] || x > b[2] || y < b[1] || y > b[3])
+                continue;
+            if (pointInRing(x, y, rings[j]))
+                depth++;
+        }
+        area += (depth % 2 ? -1 : 1) * Math.abs(ringArea(r));
+    }
+    return Math.max(0, area);
+}
+/**
+ * The rings of each feature (see {@link SliceOutlines.ringGroups}), in the order first met; one
+ * feature of every ring without `groups`.
+ */
+export function featureRings(rings, groups) {
+    if (!groups || !rings.length)
+        return rings.length ? [rings.slice()] : [];
+    const byGroup = new Map();
+    rings.forEach((r, i) => {
+        const g = groups[i] ?? 0;
+        const feature = byGroup.get(g);
+        if (feature)
+            feature.push(r);
+        else
+            byGroup.set(g, [r]);
+    });
+    return [...byGroup.values()];
+}
+/**
+ * The area of a slice's outlines (px²): each feature's (even-odd, see {@link ringsArea}), the
+ * features as their union, so a feature inside or over another of the slice adds only what it
+ * covers beyond it. Features whose extents overlap are measured together, by the length of their
+ * union along each row, exact between the rows through their vertices where no outline crosses
+ * another (and refined where one does).
+ */
+export function outlinesArea(rings, groups) {
+    const features = featureRings(rings, groups).map((f) => ({ rings: f, box: featureBox(f) }));
+    if (features.length < 2)
+        return features.length ? ringsArea(features[0].rings) : 0;
+    // Features whose extents overlap, joined (by x, so a feature meets only those it can overlap).
+    features.sort((a, b) => a.box[0] - b.box[0]);
+    const parent = features.map((_, i) => i);
+    const root = (i) => (parent[i] === i ? i : (parent[i] = root(parent[i])));
+    for (let i = 0; i < features.length; i++) {
+        const a = features[i].box;
+        for (let j = i + 1; j < features.length && features[j].box[0] < a[2]; j++) {
+            const b = features[j].box;
+            if (b[1] < a[3] && a[1] < b[3])
+                parent[root(j)] = root(i);
+        }
+    }
+    const clusters = new Map();
+    features.forEach((f, i) => {
+        const r = root(i);
+        const cluster = clusters.get(r);
+        if (cluster)
+            cluster.push(f.rings);
+        else
+            clusters.set(r, [f.rings]);
+    });
+    let area = 0;
+    for (const cluster of clusters.values())
+        area += cluster.length === 1 ? ringsArea(cluster[0]) : unionArea(cluster);
+    return area;
+}
+/**
+ * The volume of the region outlined by `slices`, in px² × slices (one slice thick per outline):
+ * the sum of each slice's area (see {@link outlinesArea}; Cavalieri's estimator). Times µm per px²
+ * and per slice, µm³.
+ */
+export function outlinedVolume(slices) {
+    let v = 0;
+    for (const s of slices)
+        v += outlinesArea(s.rings, s.ringGroups);
+    return v;
+}
+/** The area of the union of `features` (each even-odd over its rings), row band by row band. */
+function unionArea(features) {
+    const ys = [];
+    for (const f of features)
+        for (const r of f)
+            for (let i = 1; i < r.length; i += 2)
+                ys.push(r[i]);
+    ys.sort((a, b) => a - b);
+    let area = 0;
+    for (let k = 0; k + 1 < ys.length; k++)
+        if (ys[k + 1] > ys[k])
+            area += bandArea(features, ys[k], ys[k + 1], 0);
+    return area;
+}
+/**
+ * The union's area between rows `a` and `b`. Between two vertices' rows its width along a row is
+ * linear in y unless outlines cross there: then the band is halved until it is (or thin enough).
+ */
+function bandArea(features, a, b, depth) {
+    const h = b - a;
+    const at = (t) => unionWidth(features, a + t * h);
+    const mid = at(0.5);
+    if (depth >= 10)
+        return h * mid;
+    const e = 1e-6;
+    const tol = 1e-9 * (1 + Math.abs(mid));
+    const linear = Math.abs(at(0.25) + at(0.75) - 2 * mid) <= tol && Math.abs(at(e) + at(1 - e) - 2 * mid) <= tol;
+    if (linear)
+        return h * mid;
+    const m = a + h / 2;
+    return bandArea(features, a, m, depth + 1) + bandArea(features, m, b, depth + 1);
+}
+/** The length of the union of `features` along the row `y`. */
+function unionWidth(features, y) {
+    const spans = [];
+    const xs = [];
+    for (const f of features) {
+        xs.length = 0;
+        for (const r of f) {
+            const n = r.length >> 1;
+            for (let i = 0, k = n - 1; i < n; k = i++) {
+                const ya = r[2 * i + 1], yb = r[2 * k + 1];
+                if (ya > y !== yb > y)
+                    xs.push(r[2 * i] + ((y - ya) * (r[2 * k] - r[2 * i])) / (yb - ya));
+            }
+        }
+        xs.sort((p, q) => p - q);
+        for (let c = 0; c + 1 < xs.length; c += 2)
+            spans.push([xs[c], xs[c + 1]]);
+    }
+    spans.sort((p, q) => p[0] - q[0]);
+    let width = 0, from = -Infinity, to = -Infinity;
+    for (const [s0, s1] of spans) {
+        if (s0 > to) {
+            if (to > from)
+                width += to - from;
+            from = s0;
+            to = s1;
+        }
+        else if (s1 > to)
+            to = s1;
+    }
+    if (to > from)
+        width += to - from;
+    return width;
+}
+/** The extent of a feature's rings. */
+function featureBox(rings) {
+    const box = [Infinity, Infinity, -Infinity, -Infinity];
+    for (const r of rings) {
+        const [x0, y0, x1, y1] = ringBox(r);
+        box[0] = Math.min(box[0], x0);
+        box[1] = Math.min(box[1], y0);
+        box[2] = Math.max(box[2], x1);
+        box[3] = Math.max(box[3], y1);
+    }
+    return box;
+}
+function ringBox(r) {
+    let x0 = Infinity, y0 = Infinity, x1 = -Infinity, y1 = -Infinity;
+    for (let i = 0; i + 1 < r.length; i += 2) {
+        if (r[i] < x0)
+            x0 = r[i];
+        if (r[i] > x1)
+            x1 = r[i];
+        if (r[i + 1] < y0)
+            y0 = r[i + 1];
+        if (r[i + 1] > y1)
+            y1 = r[i + 1];
+    }
+    return [x0, y0, x1, y1];
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoicmVjb25zdHJ1Y3QuanMiLCJzb3VyY2VSb290IjoiIiwic291cmNlcyI6WyIuLi8uLi8uLi8uLi9zcmMvbGliL2dlb21ldHJ5L3JlY29uc3RydWN0LnRzIl0sIm5hbWVzIjpbXSwibWFwcGluZ3MiOiJBQUFBOzs7Ozs7R0FNRztBQUNILE9BQU8sRUFBRSxXQUFXLEVBQUUsUUFBUSxFQUFFLFVBQVUsRUFBRSxXQUFXLEVBQUUsTUFBTSxvQkFBb0IsQ0FBQztBQUdwRixPQUFPLEVBQUUsU0FBUyxFQUFFLE1BQU0sUUFBUSxDQUFDO0FBMkJuQzs7OztHQUlHO0FBQ0gsTUFBTSxVQUFVLFdBQVcsQ0FBQyxNQUFnQyxFQUFFLE9BQTJCLEVBQUU7SUFDekYsTUFBTSxNQUFNLEdBQUcsSUFBSSxDQUFDLE1BQU0sRUFBRSxJQUFJLENBQUMsQ0FBQztJQUNsQyxNQUFNLEtBQUssR0FBRyxNQUFNLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsS0FBSyxFQUFFLE1BQU0sQ0FBQyxDQUFDO0lBQ3BELElBQUksQ0FBQyxLQUFLLENBQUMsTUFBTTtRQUFFLE9BQU8sTUFBTSxDQUFDO0lBQ2pDLE1BQU0sS0FBSyxHQUFHLFNBQVMsQ0FDckIsS0FBSyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxFQUFFLEtBQUssRUFBRSxDQUFDLENBQUMsS0FBSyxJQUFJLEVBQUUsRUFBRSxDQUFDLENBQUMsRUFDcEQsSUFBSSxDQUFDLFNBQVMsSUFBSSxDQUFDLENBQ3BCLENBQUM7SUFDRixPQUFPLFdBQVcsQ0FBQyxNQUFNLEVBQUUsS0FBSyxDQUFDLENBQUM7QUFDcEMsQ0FBQztBQUVELDBDQUEwQztBQUMxQyxTQUFTLElBQUksQ0FBQyxNQUFnQyxFQUFFLElBQXdCO0lBQ3RFLElBQUksRUFBRSxHQUFHLFFBQVEsRUFDZixFQUFFLEdBQUcsUUFBUSxFQUNiLEVBQUUsR0FBRyxDQUFDLFFBQVEsRUFDZCxFQUFFLEdBQUcsQ0FBQyxRQUFRLEVBQ2QsRUFBRSxHQUFHLFFBQVEsRUFDYixFQUFFLEdBQUcsQ0FBQyxRQUFRLENBQUM7SUFDakIsS0FBSyxNQUFNLENBQUMsSUFBSSxNQUFNLEVBQUUsQ0FBQztRQUN2QixLQUFLLE1BQU0sQ0FBQyxJQUFJLENBQUMsQ0FBQyxLQUFLLEVBQUUsQ0FBQztZQUN4QixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxNQUFNLEVBQUUsQ0FBQyxJQUFJLENBQUMsRUFBRSxDQUFDO2dCQUN6QyxJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFO29CQUFFLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7Z0JBQ3pCLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLEVBQUU7b0JBQUUsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztnQkFDekIsSUFBSSxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLEVBQUU7b0JBQUUsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7Z0JBQ2pDLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxFQUFFO29CQUFFLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1lBQ25DLENBQUM7UUFDSCxDQUFDO1FBQ0QsSUFBSSxDQUFDLENBQUMsS0FBSyxDQUFDLE1BQU0sRUFBRSxDQUFDO1lBQ25CLEVBQUUsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDdkIsRUFBRSxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUN6QixDQUFDO0lBQ0gsQ0FBQztJQUNELElBQUksQ0FBQyxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsSUFBSSxDQUFDLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQztRQUFFLE9BQU8sRUFBRSxRQUFRLEVBQUUsSUFBSSxZQUFZLENBQUMsQ0FBQyxDQUFDLEVBQUUsS0FBSyxFQUFFLElBQUksV0FBVyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUM7SUFDbEcsTUFBTSxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxDQUFDLENBQUM7SUFDdkIsTUFBTSxTQUFTLEdBQUcsSUFBSSxDQUFDLFNBQVMsSUFBSSxVQUFVLENBQUM7SUFDL0MsSUFBSSxJQUFJLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxJQUFJLENBQUMsU0FBUyxJQUFJLENBQUMsRUFBRSxJQUFJLENBQUMsQ0FBQztJQUMvQyxPQUFPLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLEdBQUcsSUFBSSxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxFQUFFLEdBQUcsU0FBUztRQUFFLElBQUksSUFBSSxJQUFJLENBQUM7SUFDaEcsTUFBTSxFQUFFLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsR0FBRyxJQUFJLENBQUMsQ0FBQyxDQUFDO0lBQ3BELE1BQU0sRUFBRSxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxFQUFFLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxFQUFFLEdBQUcsRUFBRSxDQUFDLEdBQUcsSUFBSSxDQUFDLENBQUMsQ0FBQztJQUNwRCxNQUFNLElBQUksR0FBRyxJQUFJLFVBQVUsQ0FBQyxFQUFFLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDO0lBQzFDLEtBQUssTUFBTSxDQUFDLElBQUksTUFBTSxFQUFFLENBQUM7UUFDdkIsS0FBSyxNQUFNLEtBQUssSUFBSSxZQUFZLENBQUMsQ0FBQyxDQUFDLEtBQUssRUFBRSxDQUFDLENBQUMsVUFBVSxDQUFDLEVBQUUsQ0FBQztZQUN4RCxTQUFTLENBQUMsSUFBSSxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxHQUFHLEVBQUUsR0FBRyxFQUFFLEVBQUUsS0FBSyxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsSUFBSSxDQUFDLENBQUM7UUFDckUsQ0FBQztJQUNILENBQUM7SUFFRCxvR0FBb0c7SUFDcEcsTUFBTSxJQUFJLEdBQUcsVUFBVSxDQUFDLFdBQVcsQ0FBQyxJQUFJLEVBQUUsQ0FBQyxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsQ0FBQyxFQUFFLEVBQUUsTUFBTSxFQUFFLENBQUMsR0FBRyxFQUFFLEdBQUcsRUFBRSxHQUFHLENBQUMsRUFBRSxDQUFDLEVBQUU7UUFDcEYsVUFBVSxFQUFFLElBQUksQ0FBQyxNQUFNLElBQUksQ0FBQztLQUM3QixDQUFDLENBQUM7SUFDSCxNQUFNLENBQUMsR0FBRyxJQUFJLENBQUMsUUFBUSxDQUFDO0lBQ3hCLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLENBQUMsTUFBTSxFQUFFLENBQUMsSUFBSSxDQUFDLEVBQUUsQ0FBQztRQUNyQyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUM7UUFDeEIsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUM7UUFDaEMsaUVBQWlFO1FBQ2pFLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsR0FBRyxDQUFDO0lBQ2pDLENBQUM7SUFDRCxPQUFPLElBQUksQ0FBQztBQUNkLENBQUM7QUFFRCwrQkFBK0I7QUFDL0IsU0FBUyxXQUFXLENBQUMsQ0FBTyxFQUFFLENBQU87SUFDbkMsSUFBSSxDQUFDLENBQUMsQ0FBQyxLQUFLLENBQUMsTUFBTTtRQUFFLE9BQU8sQ0FBQyxDQUFDO0lBQzlCLElBQUksQ0FBQyxDQUFDLENBQUMsS0FBSyxDQUFDLE1BQU07UUFBRSxPQUFPLENBQUMsQ0FBQztJQUM5QixNQUFNLFFBQVEsR0FBRyxJQUFJLFlBQVksQ0FBQyxDQUFDLENBQUMsUUFBUSxDQUFDLE1BQU0sR0FBRyxDQUFDLENBQUMsUUFBUSxDQUFDLE1BQU0sQ0FBQyxDQUFDO0lBQ3pFLFFBQVEsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLFFBQVEsQ0FBQyxDQUFDO0lBQ3pCLFFBQVEsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLFFBQVEsRUFBRSxDQUFDLENBQUMsUUFBUSxDQUFDLE1BQU0sQ0FBQyxDQUFDO0lBQzVDLE1BQU0sS0FBSyxHQUFHLElBQUksV0FBVyxDQUFDLENBQUMsQ0FBQyxLQUFLLENBQUMsTUFBTSxHQUFHLENBQUMsQ0FBQyxLQUFLLENBQUMsTUFBTSxDQUFDLENBQUM7SUFDL0QsS0FBSyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsS0FBSyxDQUFDLENBQUM7SUFDbkIsTUFBTSxNQUFNLEdBQUcsQ0FBQyxDQUFDLFFBQVEsQ0FBQyxNQUFNLEdBQUcsQ0FBQyxDQUFDO0lBQ3JDLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLENBQUMsS0FBSyxDQUFDLE1BQU0sRUFBRSxDQUFDLEVBQUU7UUFBRSxLQUFLLENBQUMsQ0FBQyxDQUFDLEtBQUssQ0FBQyxNQUFNLEdBQUcsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsR0FBRyxNQUFNLENBQUM7SUFDekYsT0FBTyxFQUFFLFFBQVEsRUFBRSxLQUFLLEVBQUUsQ0FBQztBQUM3QixDQUFDO0FBRUQ7OztHQUdHO0FBQ0gsU0FBUyxTQUFTLENBQ2hCLElBQWdCLEVBQ2hCLEVBQVUsRUFDVixFQUFVLEVBQ1YsTUFBYyxFQUNkLEtBQW1DLEVBQ25DLEVBQVUsRUFDVixFQUFVLEVBQ1YsSUFBWTtJQUVaLE1BQU0sU0FBUyxHQUFhLEVBQUUsQ0FBQztJQUMvQixtQ0FBbUM7SUFDbkMsSUFBSSxHQUFHLEdBQUcsUUFBUSxFQUNoQixNQUFNLEdBQUcsQ0FBQyxRQUFRLENBQUM7SUFDckIsS0FBSyxNQUFNLENBQUMsSUFBSSxLQUFLLEVBQUUsQ0FBQztRQUN0QixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLE1BQU0sRUFBRSxDQUFDLElBQUksQ0FBQyxFQUFFLENBQUM7WUFDckMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsR0FBRztnQkFBRSxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1lBQzNCLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLE1BQU07Z0JBQUUsTUFBTSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNuQyxDQUFDO0lBQ0gsQ0FBQztJQUNELE1BQU0sRUFBRSxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxFQUFFLElBQUksQ0FBQyxLQUFLLENBQUMsQ0FBQyxHQUFHLEdBQUcsRUFBRSxDQUFDLEdBQUcsSUFBSSxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDNUQsTUFBTSxFQUFFLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxFQUFFLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxNQUFNLEdBQUcsRUFBRSxDQUFDLEdBQUcsSUFBSSxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDbkUsS0FBSyxJQUFJLENBQUMsR0FBRyxFQUFFLEVBQUUsQ0FBQyxJQUFJLEVBQUUsRUFBRSxDQUFDLEVBQUUsRUFBRSxDQUFDO1FBQzlCLE1BQU0sQ0FBQyxHQUFHLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxHQUFHLENBQUMsR0FBRyxJQUFJLENBQUM7UUFDaEMsU0FBUyxDQUFDLE1BQU0sR0FBRyxDQUFDLENBQUM7UUFDckIsS0FBSyxNQUFNLENBQUMsSUFBSSxLQUFLLEVBQUUsQ0FBQztZQUN0QixNQUFNLENBQUMsR0FBRyxDQUFDLENBQUMsTUFBTSxJQUFJLENBQUMsQ0FBQztZQUN4QixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsRUFBRSxDQUFDO2dCQUMxQyxNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsRUFDckIsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO2dCQUNwQixJQUFJLEVBQUUsR0FBRyxDQUFDLEtBQUssRUFBRSxHQUFHLENBQUMsRUFBRSxDQUFDO29CQUN0QixNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxFQUNqQixFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztvQkFDaEIsU0FBUyxDQUFDLElBQUksQ0FBQyxFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsQ0FBQyxHQUFHLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDLENBQUM7Z0JBQzFELENBQUM7WUFDSCxDQUFDO1FBQ0gsQ0FBQztRQUNELElBQUksU0FBUyxDQUFDLE1BQU0sR0FBRyxDQUFDO1lBQUUsU0FBUztRQUNuQyxTQUFTLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ2hDLE1BQU0sR0FBRyxHQUFHLE1BQU0sR0FBRyxDQUFDLEdBQUcsRUFBRSxDQUFDO1FBQzVCLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEdBQUcsU0FBUyxDQUFDLE1BQU0sRUFBRSxDQUFDLElBQUksQ0FBQyxFQUFFLENBQUM7WUFDakQsTUFBTSxJQUFJLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLFNBQVMsQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUMsR0FBRyxJQUFJLEdBQUcsR0FBRyxDQUFDLENBQUMsQ0FBQztZQUN0RSxNQUFNLEVBQUUsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsR0FBRyxDQUFDLEVBQUUsSUFBSSxDQUFDLEtBQUssQ0FBQyxDQUFDLFNBQVMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsRUFBRSxDQUFDLEdBQUcsSUFBSSxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUM7WUFDOUUsS0FBSyxJQUFJLENBQUMsR0FBRyxJQUFJLEVBQUUsQ0FBQyxJQUFJLEVBQUUsRUFBRSxDQUFDLEVBQUU7Z0JBQUUsSUFBSSxDQUFDLEdBQUcsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUM7UUFDckQsQ0FBQztJQUNILENBQUM7QUFDSCxDQUFDO0FBRUQ7OztHQUdHO0FBQ0gsTUFBTSxVQUFVLFNBQVMsQ0FBQyxLQUFtQztJQUMzRCxNQUFNLEtBQUssR0FBRyxLQUFLLENBQUMsR0FBRyxDQUFDLE9BQU8sQ0FBQyxDQUFDO0lBQ2pDLElBQUksSUFBSSxHQUFHLENBQUMsQ0FBQztJQUNiLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxLQUFLLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7UUFDdEMsTUFBTSxDQUFDLEdBQUcsS0FBSyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ25CLElBQUksQ0FBQyxDQUFDLE1BQU0sR0FBRyxDQUFDO1lBQUUsU0FBUztRQUMzQixNQUFNLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQ1osQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNYLElBQUksS0FBSyxHQUFHLENBQUMsQ0FBQztRQUNkLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxLQUFLLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7WUFDdEMsSUFBSSxDQUFDLEtBQUssQ0FBQztnQkFBRSxTQUFTO1lBQ3RCLE1BQU0sQ0FBQyxHQUFHLEtBQUssQ0FBQyxDQUFDLENBQUMsQ0FBQztZQUNuQixJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDO2dCQUFFLFNBQVM7WUFDM0QsSUFBSSxXQUFXLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxLQUFLLENBQUMsQ0FBQyxDQUFDLENBQUM7Z0JBQUUsS0FBSyxFQUFFLENBQUM7UUFDM0MsQ0FBQztRQUNELElBQUksSUFBSSxDQUFDLEtBQUssR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLFFBQVEsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ3ZELENBQUM7SUFDRCxPQUFPLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxFQUFFLElBQUksQ0FBQyxDQUFDO0FBQzNCLENBQUM7QUFFRDs7O0dBR0c7QUFDSCxNQUFNLFVBQVUsWUFBWSxDQUMxQixLQUFtQyxFQUNuQyxNQUEwQjtJQUUxQixJQUFJLENBQUMsTUFBTSxJQUFJLENBQUMsS0FBSyxDQUFDLE1BQU07UUFBRSxPQUFPLEtBQUssQ0FBQyxNQUFNLENBQUMsQ0FBQyxDQUFDLENBQUMsS0FBSyxDQUFDLEtBQUssRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQztJQUN6RSxNQUFNLE9BQU8sR0FBRyxJQUFJLEdBQUcsRUFBK0IsQ0FBQztJQUN2RCxLQUFLLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFO1FBQ3JCLE1BQU0sQ0FBQyxHQUFHLE1BQU0sQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUM7UUFDekIsTUFBTSxPQUFPLEdBQUcsT0FBTyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUMvQixJQUFJLE9BQU87WUFBRSxPQUFPLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDOztZQUN4QixPQUFPLENBQUMsR0FBRyxDQUFDLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDM0IsQ0FBQyxDQUFDLENBQUM7SUFDSCxPQUFPLENBQUMsR0FBRyxPQUFPLENBQUMsTUFBTSxFQUFFLENBQUMsQ0FBQztBQUMvQixDQUFDO0FBRUQ7Ozs7OztHQU1HO0FBQ0gsTUFBTSxVQUFVLFlBQVksQ0FBQyxLQUFtQyxFQUFFLE1BQTBCO0lBQzFGLE1BQU0sUUFBUSxHQUFHLFlBQVksQ0FBQyxLQUFLLEVBQUUsTUFBTSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLEVBQUUsS0FBSyxFQUFFLENBQUMsRUFBRSxHQUFHLEVBQUUsVUFBVSxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDO0lBQzVGLElBQUksUUFBUSxDQUFDLE1BQU0sR0FBRyxDQUFDO1FBQUUsT0FBTyxRQUFRLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxTQUFTLENBQUMsUUFBUSxDQUFDLENBQUMsQ0FBQyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDbkYsK0ZBQStGO0lBQy9GLFFBQVEsQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUM3QyxNQUFNLE1BQU0sR0FBRyxRQUFRLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDekMsTUFBTSxJQUFJLEdBQUcsQ0FBQyxDQUFTLEVBQVUsRUFBRSxDQUFDLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxLQUFLLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLE1BQU0sQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsTUFBTSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQzFGLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxRQUFRLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7UUFDekMsTUFBTSxDQUFDLEdBQUcsUUFBUSxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQztRQUMxQixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLFFBQVEsQ0FBQyxNQUFNLElBQUksUUFBUSxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQztZQUMxRSxNQUFNLENBQUMsR0FBRyxRQUFRLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDO1lBQzFCLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztnQkFBRSxNQUFNLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQzVELENBQUM7SUFDSCxDQUFDO0lBQ0QsTUFBTSxRQUFRLEdBQUcsSUFBSSxHQUFHLEVBQWlDLENBQUM7SUFDMUQsUUFBUSxDQUFDLE9BQU8sQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDLEVBQUUsRUFBRTtRQUN4QixNQUFNLENBQUMsR0FBRyxJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDbEIsTUFBTSxPQUFPLEdBQUcsUUFBUSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNoQyxJQUFJLE9BQU87WUFBRSxPQUFPLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxLQUFLLENBQUMsQ0FBQzs7WUFDOUIsUUFBUSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUMsS0FBSyxDQUFDLENBQUMsQ0FBQztJQUNsQyxDQUFDLENBQUMsQ0FBQztJQUNILElBQUksSUFBSSxHQUFHLENBQUMsQ0FBQztJQUNiLEtBQUssTUFBTSxPQUFPLElBQUksUUFBUSxDQUFDLE1BQU0sRUFBRTtRQUNyQyxJQUFJLElBQUksT0FBTyxDQUFDLE1BQU0sS0FBSyxDQUFDLENBQUMsQ0FBQyxDQUFDLFNBQVMsQ0FBQyxPQUFPLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsU0FBUyxDQUFDLE9BQU8sQ0FBQyxDQUFDO0lBQzVFLE9BQU8sSUFBSSxDQUFDO0FBQ2QsQ0FBQztBQUVEOzs7O0dBSUc7QUFDSCxNQUFNLFVBQVUsY0FBYyxDQUFDLE1BQWdDO0lBQzdELElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQztJQUNWLEtBQUssTUFBTSxDQUFDLElBQUksTUFBTTtRQUFFLENBQUMsSUFBSSxZQUFZLENBQUMsQ0FBQyxDQUFDLEtBQUssRUFBRSxDQUFDLENBQUMsVUFBVSxDQUFDLENBQUM7SUFDakUsT0FBTyxDQUFDLENBQUM7QUFDWCxDQUFDO0FBRUQsZ0dBQWdHO0FBQ2hHLFNBQVMsU0FBUyxDQUFDLFFBQXdDO0lBQ3pELE1BQU0sRUFBRSxHQUFhLEVBQUUsQ0FBQztJQUN4QixLQUFLLE1BQU0sQ0FBQyxJQUFJLFFBQVE7UUFBRSxLQUFLLE1BQU0sQ0FBQyxJQUFJLENBQUM7WUFBRSxLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLE1BQU0sRUFBRSxDQUFDLElBQUksQ0FBQztnQkFBRSxFQUFFLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ2pHLEVBQUUsQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDekIsSUFBSSxJQUFJLEdBQUcsQ0FBQyxDQUFDO0lBQ2IsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsR0FBRyxFQUFFLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRTtRQUFFLElBQUksRUFBRSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxFQUFFLENBQUMsQ0FBQyxDQUFDO1lBQUUsSUFBSSxJQUFJLFFBQVEsQ0FBQyxRQUFRLEVBQUUsRUFBRSxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUM7SUFDL0csT0FBTyxJQUFJLENBQUM7QUFDZCxDQUFDO0FBRUQ7OztHQUdHO0FBQ0gsU0FBUyxRQUFRLENBQUMsUUFBd0MsRUFBRSxDQUFTLEVBQUUsQ0FBUyxFQUFFLEtBQWE7SUFDN0YsTUFBTSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQztJQUNoQixNQUFNLEVBQUUsR0FBRyxDQUFDLENBQVMsRUFBRSxFQUFFLENBQUMsVUFBVSxDQUFDLFFBQVEsRUFBRSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO0lBQzFELE1BQU0sR0FBRyxHQUFHLEVBQUUsQ0FBQyxHQUFHLENBQUMsQ0FBQztJQUNwQixJQUFJLEtBQUssSUFBSSxFQUFFO1FBQUUsT0FBTyxDQUFDLEdBQUcsR0FBRyxDQUFDO0lBQ2hDLE1BQU0sQ0FBQyxHQUFHLElBQUksQ0FBQztJQUNmLE1BQU0sR0FBRyxHQUFHLElBQUksR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDdkMsTUFBTSxNQUFNLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsSUFBSSxDQUFDLEdBQUcsRUFBRSxDQUFDLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxHQUFHLENBQUMsSUFBSSxHQUFHLElBQUksSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLEdBQUcsR0FBRyxDQUFDLElBQUksR0FBRyxDQUFDO0lBQzlHLElBQUksTUFBTTtRQUFFLE9BQU8sQ0FBQyxHQUFHLEdBQUcsQ0FBQztJQUMzQixNQUFNLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQztJQUNwQixPQUFPLFFBQVEsQ0FBQyxRQUFRLEVBQUUsQ0FBQyxFQUFFLENBQUMsRUFBRSxLQUFLLEdBQUcsQ0FBQyxDQUFDLEdBQUcsUUFBUSxDQUFDLFFBQVEsRUFBRSxDQUFDLEVBQUUsQ0FBQyxFQUFFLEtBQUssR0FBRyxDQUFDLENBQUMsQ0FBQztBQUNuRixDQUFDO0FBRUQsK0RBQStEO0FBQy9ELFNBQVMsVUFBVSxDQUFDLFFBQXdDLEVBQUUsQ0FBUztJQUNyRSxNQUFNLEtBQUssR0FBNEIsRUFBRSxDQUFDO0lBQzFDLE1BQU0sRUFBRSxHQUFhLEVBQUUsQ0FBQztJQUN4QixLQUFLLE1BQU0sQ0FBQyxJQUFJLFFBQVEsRUFBRSxDQUFDO1FBQ3pCLEVBQUUsQ0FBQyxNQUFNLEdBQUcsQ0FBQyxDQUFDO1FBQ2QsS0FBSyxNQUFNLENBQUMsSUFBSSxDQUFDLEVBQUUsQ0FBQztZQUNsQixNQUFNLENBQUMsR0FBRyxDQUFDLENBQUMsTUFBTSxJQUFJLENBQUMsQ0FBQztZQUN4QixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsRUFBRSxDQUFDO2dCQUMxQyxNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsRUFDckIsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO2dCQUNwQixJQUFJLEVBQUUsR0FBRyxDQUFDLEtBQUssRUFBRSxHQUFHLENBQUM7b0JBQUUsRUFBRSxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsRUFBRSxHQUFHLEVBQUUsQ0FBQyxDQUFDLENBQUM7WUFDNUYsQ0FBQztRQUNILENBQUM7UUFDRCxFQUFFLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ3pCLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEdBQUcsRUFBRSxDQUFDLE1BQU0sRUFBRSxDQUFDLElBQUksQ0FBQztZQUFFLEtBQUssQ0FBQyxJQUFJLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDNUUsQ0FBQztJQUNELEtBQUssQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDbEMsSUFBSSxLQUFLLEdBQUcsQ0FBQyxFQUNYLElBQUksR0FBRyxDQUFDLFFBQVEsRUFDaEIsRUFBRSxHQUFHLENBQUMsUUFBUSxDQUFDO0lBQ2pCLEtBQUssTUFBTSxDQUFDLEVBQUUsRUFBRSxFQUFFLENBQUMsSUFBSSxLQUFLLEVBQUUsQ0FBQztRQUM3QixJQUFJLEVBQUUsR0FBRyxFQUFFLEVBQUUsQ0FBQztZQUNaLElBQUksRUFBRSxHQUFHLElBQUk7Z0JBQUUsS0FBSyxJQUFJLEVBQUUsR0FBRyxJQUFJLENBQUM7WUFDbEMsSUFBSSxHQUFHLEVBQUUsQ0FBQztZQUNWLEVBQUUsR0FBRyxFQUFFLENBQUM7UUFDVixDQUFDO2FBQU0sSUFBSSxFQUFFLEdBQUcsRUFBRTtZQUFFLEVBQUUsR0FBRyxFQUFFLENBQUM7SUFDOUIsQ0FBQztJQUNELElBQUksRUFBRSxHQUFHLElBQUk7UUFBRSxLQUFLLElBQUksRUFBRSxHQUFHLElBQUksQ0FBQztJQUNsQyxPQUFPLEtBQUssQ0FBQztBQUNmLENBQUM7QUFFRCx1Q0FBdUM7QUFDdkMsU0FBUyxVQUFVLENBQUMsS0FBbUM7SUFDckQsTUFBTSxHQUFHLEdBQXFDLENBQUMsUUFBUSxFQUFFLFFBQVEsRUFBRSxDQUFDLFFBQVEsRUFBRSxDQUFDLFFBQVEsQ0FBQyxDQUFDO0lBQ3pGLEtBQUssTUFBTSxDQUFDLElBQUksS0FBSyxFQUFFLENBQUM7UUFDdEIsTUFBTSxDQUFDLEVBQUUsRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsQ0FBQyxHQUFHLE9BQU8sQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNwQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUM7UUFDOUIsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDO1FBQzlCLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQztRQUM5QixHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUM7SUFDaEMsQ0FBQztJQUNELE9BQU8sR0FBRyxDQUFDO0FBQ2IsQ0FBQztBQUVELFNBQVMsT0FBTyxDQUFDLENBQW9CO0lBQ25DLElBQUksRUFBRSxHQUFHLFFBQVEsRUFDZixFQUFFLEdBQUcsUUFBUSxFQUNiLEVBQUUsR0FBRyxDQUFDLFFBQVEsRUFDZCxFQUFFLEdBQUcsQ0FBQyxRQUFRLENBQUM7SUFDakIsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsTUFBTSxFQUFFLENBQUMsSUFBSSxDQUFDLEVBQUUsQ0FBQztRQUN6QyxJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFO1lBQUUsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUN6QixJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxFQUFFO1lBQUUsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUN6QixJQUFJLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsRUFBRTtZQUFFLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ2pDLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxFQUFFO1lBQUUsRUFBRSxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDbkMsQ0FBQztJQUNELE9BQU8sQ0FBQyxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQztBQUMxQixDQUFDIiwic291cmNlc0NvbnRlbnQiOlsiLyoqXG4gKiBUaGUgM0QgcmVjb25zdHJ1Y3Rpb24gb2YgYSByZWdpb24gZnJvbSBpdHMgb3V0bGluZXMgb24gYSBzdGFjayBvZiBzbGljZXM6IGVhY2ggc2xpY2UncyBvdXRsaW5lc1xuICogZmlsbGVkIChlYWNoIGZlYXR1cmUgZXZlbi1vZGQsIHRoZSBmZWF0dXJlcyBhcyB0aGVpciB1bmlvbikgaW50byBhIHZveGVsIG1hc2ssIHRoZSBtYXNrJ3NcbiAqIGJvdW5kYXJ5IG1lc2hlZCBieSBzdXJmYWNlIG5ldHMgYW5kIHNtb290aGVkOyBhbmQgaXRzIG9wZW4gbGluZXMgbG9mdGVkIGludG8gYSBzaGVldC5cbiAqXG4gKiBQdXJlOiB1c2FibGUgZnJvbSB3b3JrZXJzLlxuICovXG5pbXBvcnQgeyBwb2ludEluUmluZywgcmluZ0FyZWEsIHNtb290aE1lc2gsIHN1cmZhY2VOZXRzIH0gZnJvbSAnbmFwYXJpLWpzL2dlb21ldHJ5JztcbmltcG9ydCB0eXBlIHsgTWVzaCB9IGZyb20gJ25hcGFyaS1qcy9nZW9tZXRyeSc7XG5cbmltcG9ydCB7IGxvZnRMaW5lcyB9IGZyb20gJy4vbG9mdCc7XG5cbi8qKlxuICogQSByZWdpb24ncyBvdXRsaW5lcyBvbiBvbmUgc2xpY2U6IHJpbmdzIGFzIGludGVybGVhdmVkIGB4LCB5YCBpbWFnZSBwaXhlbHMgKGhvbGVzIGluY2x1ZGVkKSwgYW5kXG4gKiBvcGVuIHBvbHlsaW5lcyAoYSBsYXllciBib3VuZGFyeSksIGxpa2V3aXNlLlxuICovXG5leHBvcnQgaW50ZXJmYWNlIFNsaWNlT3V0bGluZXMge1xuICB6OiBudW1iZXI7XG4gIHJpbmdzOiBBcnJheUxpa2U8bnVtYmVyPltdO1xuICAvKipcbiAgICogVGhlIGZlYXR1cmUgZWFjaCByaW5nIGlzIHBhcnQgb2YgKGVxdWFsIGlkcywgb25lIGZlYXR1cmUpOiBhIGZlYXR1cmUncyByaW5ncyBmaWxsIGV2ZW4tb2RkIChpdHNcbiAgICogaG9sZXMpLCB0aGUgZmVhdHVyZXMgYXMgdGhlaXIgdW5pb24sIHNvIG9uZSBpbnNpZGUgb3Igb3ZlciBhbm90aGVyIGN1dHMgbm8gaG9sZSBpbiBpdC4gQWJzZW50OlxuICAgKiB0aGUgcmluZ3MgYXJlIG9uZSBmZWF0dXJlLlxuICAgKi9cbiAgcmluZ0dyb3Vwcz86IEFycmF5TGlrZTxudW1iZXI+O1xuICBsaW5lcz86IEFycmF5TGlrZTxudW1iZXI+W107XG59XG5cbmV4cG9ydCBpbnRlcmZhY2UgUmVjb25zdHJ1Y3RPcHRpb25zIHtcbiAgLyoqIEluLXBsYW5lIHZveGVsIGVkZ2UsIGltYWdlIHB4IChkZWZhdWx0IDEpLiBSYWlzZWQgdG8ga2VlcCB0aGUgbWFzayB3aXRoaW4gYG1heFZveGVsc2AuICovXG4gIHBpeGVsU3RlcD86IG51bWJlcjtcbiAgLyoqIEF0IG1vc3QgdGhpcyBtYW55IHZveGVscyAoZGVmYXVsdCAzMiBNKS4gKi9cbiAgbWF4Vm94ZWxzPzogbnVtYmVyO1xuICAvKiogTGFwbGFjaWFuIHNtb290aGluZyBwYXNzZXMgKGRlZmF1bHQgMikuICovXG4gIHNtb290aD86IG51bWJlcjtcbn1cblxuLyoqXG4gKiBUaGUgbWVzaCBvZiB0aGUgcmVnaW9uIG91dGxpbmVkIGJ5IGBzbGljZXNgLCBpbiBpbWFnZSBwaXhlbHMgKHgsIHkpIGFuZCBzbGljZSBpbmRpY2VzICh6OiBzbGljZVxuICogaydzIHBsYW5lIGF0IHogPSBrKTogdGhlIGZpbGxlZCBvdXRsaW5lcycgc3VyZmFjZSwgYW5kIHRoZSBzaGVldCB0aHJvdWdoIHRoZSBsaW5lcyAoc2VlXG4gKiB7QGxpbmsgbG9mdExpbmVzfSkuIEVtcHR5IHdoZW4gbm8gc2xpY2UgaGFzIGFuIGFyZWEgYW5kIG5vIHR3byBzbGljZXMgaGF2ZSBsaW5lcy5cbiAqL1xuZXhwb3J0IGZ1bmN0aW9uIHJlY29uc3RydWN0KHNsaWNlczogcmVhZG9ubHkgU2xpY2VPdXRsaW5lc1tdLCBvcHRzOiBSZWNvbnN0cnVjdE9wdGlvbnMgPSB7fSk6IE1lc2gge1xuICBjb25zdCBmaWxsZWQgPSBmaWxsKHNsaWNlcywgb3B0cyk7XG4gIGNvbnN0IGxpbmVzID0gc2xpY2VzLmZpbHRlcigocykgPT4gcy5saW5lcz8ubGVuZ3RoKTtcbiAgaWYgKCFsaW5lcy5sZW5ndGgpIHJldHVybiBmaWxsZWQ7XG4gIGNvbnN0IHNoZWV0ID0gbG9mdExpbmVzKFxuICAgIGxpbmVzLm1hcCgocykgPT4gKHsgejogcy56LCBsaW5lczogcy5saW5lcyA/PyBbXSB9KSksXG4gICAgb3B0cy5waXhlbFN0ZXAgPz8gMSxcbiAgKTtcbiAgcmV0dXJuIG1lcmdlTWVzaGVzKGZpbGxlZCwgc2hlZXQpO1xufVxuXG4vKiogVGhlIHN1cmZhY2Ugb2YgdGhlIGZpbGxlZCBvdXRsaW5lcy4gKi9cbmZ1bmN0aW9uIGZpbGwoc2xpY2VzOiByZWFkb25seSBTbGljZU91dGxpbmVzW10sIG9wdHM6IFJlY29uc3RydWN0T3B0aW9ucyk6IE1lc2gge1xuICBsZXQgeDAgPSBJbmZpbml0eSxcbiAgICB5MCA9IEluZmluaXR5LFxuICAgIHgxID0gLUluZmluaXR5LFxuICAgIHkxID0gLUluZmluaXR5LFxuICAgIHowID0gSW5maW5pdHksXG4gICAgejEgPSAtSW5maW5pdHk7XG4gIGZvciAoY29uc3QgcyBvZiBzbGljZXMpIHtcbiAgICBmb3IgKGNvbnN0IHIgb2Ygcy5yaW5ncykge1xuICAgICAgZm9yIChsZXQgaSA9IDA7IGkgKyAxIDwgci5sZW5ndGg7IGkgKz0gMikge1xuICAgICAgICBpZiAocltpXSA8IHgwKSB4MCA9IHJbaV07XG4gICAgICAgIGlmIChyW2ldID4geDEpIHgxID0gcltpXTtcbiAgICAgICAgaWYgKHJbaSArIDFdIDwgeTApIHkwID0gcltpICsgMV07XG4gICAgICAgIGlmIChyW2kgKyAxXSA+IHkxKSB5MSA9IHJbaSArIDFdO1xuICAgICAgfVxuICAgIH1cbiAgICBpZiAocy5yaW5ncy5sZW5ndGgpIHtcbiAgICAgIHowID0gTWF0aC5taW4oejAsIHMueik7XG4gICAgICB6MSA9IE1hdGgubWF4KHoxLCBzLnopO1xuICAgIH1cbiAgfVxuICBpZiAoISh4MSA+IHgwKSB8fCAhKHkxID4geTApKSByZXR1cm4geyB2ZXJ0aWNlczogbmV3IEZsb2F0MzJBcnJheSgwKSwgZmFjZXM6IG5ldyBVaW50MzJBcnJheSgwKSB9O1xuICBjb25zdCBueiA9IHoxIC0gejAgKyAxO1xuICBjb25zdCBtYXhWb3hlbHMgPSBvcHRzLm1heFZveGVscyA/PyAzMl8wMDBfMDAwO1xuICBsZXQgc3RlcCA9IE1hdGgubWF4KG9wdHMucGl4ZWxTdGVwID8/IDEsIDFlLTYpO1xuICB3aGlsZSAoTWF0aC5jZWlsKCh4MSAtIHgwKSAvIHN0ZXApICogTWF0aC5jZWlsKCh5MSAtIHkwKSAvIHN0ZXApICogbnogPiBtYXhWb3hlbHMpIHN0ZXAgKj0gMS4yNTtcbiAgY29uc3QgbnggPSBNYXRoLm1heCgxLCBNYXRoLmNlaWwoKHgxIC0geDApIC8gc3RlcCkpO1xuICBjb25zdCBueSA9IE1hdGgubWF4KDEsIE1hdGguY2VpbCgoeTEgLSB5MCkgLyBzdGVwKSk7XG4gIGNvbnN0IG1hc2sgPSBuZXcgVWludDhBcnJheShueCAqIG55ICogbnopO1xuICBmb3IgKGNvbnN0IHMgb2Ygc2xpY2VzKSB7XG4gICAgZm9yIChjb25zdCByaW5ncyBvZiBmZWF0dXJlUmluZ3Mocy5yaW5ncywgcy5yaW5nR3JvdXBzKSkge1xuICAgICAgZmlsbFNsaWNlKG1hc2ssIG54LCBueSwgKHMueiAtIHowKSAqIG54ICogbnksIHJpbmdzLCB4MCwgeTAsIHN0ZXApO1xuICAgIH1cbiAgfVxuXG4gIC8vIFZveGVsIGkgc3Bhbm5pbmcgW2ksIGkgKyAxKSAobmFwYXJpLWpzIGNlbnRyZXMgaXQgb24gaSk6IHRoZSBzY2FsaW5nIGJlbG93IG1hcHMgY2VsbCBlZGdlcyB0byBweC5cbiAgY29uc3QgbWVzaCA9IHNtb290aE1lc2goc3VyZmFjZU5ldHMobWFzaywgW254LCBueSwgbnpdLCB7IG9yaWdpbjogWzAuNSwgMC41LCAwLjVdIH0pLCB7XG4gICAgaXRlcmF0aW9uczogb3B0cy5zbW9vdGggPz8gMixcbiAgfSk7XG4gIGNvbnN0IHYgPSBtZXNoLnZlcnRpY2VzO1xuICBmb3IgKGxldCBpID0gMDsgaSA8IHYubGVuZ3RoOyBpICs9IDMpIHtcbiAgICB2W2ldID0geDAgKyB2W2ldICogc3RlcDtcbiAgICB2W2kgKyAxXSA9IHkwICsgdltpICsgMV0gKiBzdGVwO1xuICAgIC8vIFZveGVsIGsgc3BhbnMgW2ssIGsrMSkgd2l0aCBpdHMgY2VudHJlICh0aGUgc2xpY2UpIGF0IGsgKyAwLjUuXG4gICAgdltpICsgMl0gPSB6MCArIHZbaSArIDJdIC0gMC41O1xuICB9XG4gIHJldHVybiBtZXNoO1xufVxuXG4vKiogYGFgIGFuZCBgYmAgYXMgb25lIG1lc2guICovXG5mdW5jdGlvbiBtZXJnZU1lc2hlcyhhOiBNZXNoLCBiOiBNZXNoKTogTWVzaCB7XG4gIGlmICghYi5mYWNlcy5sZW5ndGgpIHJldHVybiBhO1xuICBpZiAoIWEuZmFjZXMubGVuZ3RoKSByZXR1cm4gYjtcbiAgY29uc3QgdmVydGljZXMgPSBuZXcgRmxvYXQzMkFycmF5KGEudmVydGljZXMubGVuZ3RoICsgYi52ZXJ0aWNlcy5sZW5ndGgpO1xuICB2ZXJ0aWNlcy5zZXQoYS52ZXJ0aWNlcyk7XG4gIHZlcnRpY2VzLnNldChiLnZlcnRpY2VzLCBhLnZlcnRpY2VzLmxlbmd0aCk7XG4gIGNvbnN0IGZhY2VzID0gbmV3IFVpbnQzMkFycmF5KGEuZmFjZXMubGVuZ3RoICsgYi5mYWNlcy5sZW5ndGgpO1xuICBmYWNlcy5zZXQoYS5mYWNlcyk7XG4gIGNvbnN0IG9mZnNldCA9IGEudmVydGljZXMubGVuZ3RoIC8gMztcbiAgZm9yIChsZXQgaSA9IDA7IGkgPCBiLmZhY2VzLmxlbmd0aDsgaSsrKSBmYWNlc1thLmZhY2VzLmxlbmd0aCArIGldID0gYi5mYWNlc1tpXSArIG9mZnNldDtcbiAgcmV0dXJuIHsgdmVydGljZXMsIGZhY2VzIH07XG59XG5cbi8qKlxuICogRmlsbCBvbmUgZmVhdHVyZSBpbnRvIG9uZSBzbGljZSBvZiB0aGUgbWFzazogZXZlbi1vZGQgb3ZlciBpdHMgYHJpbmdzYCwgc2FtcGxlZCBhdCB2b3hlbFxuICogY2VudHJlcy4gQ2VsbHMgYWxyZWFkeSBzZXQgc3RheSBzZXQsIHNvIHRoZSBmZWF0dXJlcyBvZiBhIHNsaWNlIGZpbGwgYXMgdGhlaXIgdW5pb24uXG4gKi9cbmZ1bmN0aW9uIGZpbGxTbGljZShcbiAgbWFzazogVWludDhBcnJheSxcbiAgbng6IG51bWJlcixcbiAgbnk6IG51bWJlcixcbiAgb2Zmc2V0OiBudW1iZXIsXG4gIHJpbmdzOiByZWFkb25seSBBcnJheUxpa2U8bnVtYmVyPltdLFxuICB4MDogbnVtYmVyLFxuICB5MDogbnVtYmVyLFxuICBzdGVwOiBudW1iZXIsXG4pOiB2b2lkIHtcbiAgY29uc3QgY3Jvc3NpbmdzOiBudW1iZXJbXSA9IFtdO1xuICAvLyBPbmx5IHRoZSByb3dzIHRoZSBmZWF0dXJlIHNwYW5zLlxuICBsZXQgdG9wID0gSW5maW5pdHksXG4gICAgYm90dG9tID0gLUluZmluaXR5O1xuICBmb3IgKGNvbnN0IHIgb2YgcmluZ3MpIHtcbiAgICBmb3IgKGxldCBpID0gMTsgaSA8IHIubGVuZ3RoOyBpICs9IDIpIHtcbiAgICAgIGlmIChyW2ldIDwgdG9wKSB0b3AgPSByW2ldO1xuICAgICAgaWYgKHJbaV0gPiBib3R0b20pIGJvdHRvbSA9IHJbaV07XG4gICAgfVxuICB9XG4gIGNvbnN0IGowID0gTWF0aC5tYXgoMCwgTWF0aC5mbG9vcigodG9wIC0geTApIC8gc3RlcCAtIDAuNSkpO1xuICBjb25zdCBqMSA9IE1hdGgubWluKG55IC0gMSwgTWF0aC5jZWlsKChib3R0b20gLSB5MCkgLyBzdGVwIC0gMC41KSk7XG4gIGZvciAobGV0IGogPSBqMDsgaiA8PSBqMTsgaisrKSB7XG4gICAgY29uc3QgeSA9IHkwICsgKGogKyAwLjUpICogc3RlcDtcbiAgICBjcm9zc2luZ3MubGVuZ3RoID0gMDtcbiAgICBmb3IgKGNvbnN0IHIgb2YgcmluZ3MpIHtcbiAgICAgIGNvbnN0IG4gPSByLmxlbmd0aCA+PiAxO1xuICAgICAgZm9yIChsZXQgaSA9IDAsIGsgPSBuIC0gMTsgaSA8IG47IGsgPSBpKyspIHtcbiAgICAgICAgY29uc3QgeWEgPSByWzIgKiBpICsgMV0sXG4gICAgICAgICAgeWIgPSByWzIgKiBrICsgMV07XG4gICAgICAgIGlmICh5YSA+IHkgIT09IHliID4geSkge1xuICAgICAgICAgIGNvbnN0IHhhID0gclsyICogaV0sXG4gICAgICAgICAgICB4YiA9IHJbMiAqIGtdO1xuICAgICAgICAgIGNyb3NzaW5ncy5wdXNoKHhhICsgKCh5IC0geWEpICogKHhiIC0geGEpKSAvICh5YiAtIHlhKSk7XG4gICAgICAgIH1cbiAgICAgIH1cbiAgICB9XG4gICAgaWYgKGNyb3NzaW5ncy5sZW5ndGggPCAyKSBjb250aW51ZTtcbiAgICBjcm9zc2luZ3Muc29ydCgoYSwgYikgPT4gYSAtIGIpO1xuICAgIGNvbnN0IHJvdyA9IG9mZnNldCArIGogKiBueDtcbiAgICBmb3IgKGxldCBjID0gMDsgYyArIDEgPCBjcm9zc2luZ3MubGVuZ3RoOyBjICs9IDIpIHtcbiAgICAgIGNvbnN0IGZyb20gPSBNYXRoLm1heCgwLCBNYXRoLmNlaWwoKGNyb3NzaW5nc1tjXSAtIHgwKSAvIHN0ZXAgLSAwLjUpKTtcbiAgICAgIGNvbnN0IHRvID0gTWF0aC5taW4obnggLSAxLCBNYXRoLmZsb29yKChjcm9zc2luZ3NbYyArIDFdIC0geDApIC8gc3RlcCAtIDAuNSkpO1xuICAgICAgZm9yIChsZXQgaSA9IGZyb207IGkgPD0gdG87IGkrKykgbWFza1tyb3cgKyBpXSA9IDE7XG4gICAgfVxuICB9XG59XG5cbi8qKlxuICogVGhlIGFyZWEgaW5zaWRlIGByaW5nc2AgKGludGVybGVhdmVkIGB4LCB5YCksIGV2ZW4tb2RkOiBhIHJpbmcgaW5zaWRlIGFuIG9kZCBudW1iZXIgb2Ygb3RoZXJzIGlzXG4gKiBhIGhvbGUuIEV4YWN0IGZvciByaW5ncyB0aGF0IGRvIG5vdCBjcm9zcyAoYSBwb2x5Z29uIHNldCdzIGV4dGVyaW9ycyBhbmQgaG9sZXMpLlxuICovXG5leHBvcnQgZnVuY3Rpb24gcmluZ3NBcmVhKHJpbmdzOiByZWFkb25seSBBcnJheUxpa2U8bnVtYmVyPltdKTogbnVtYmVyIHtcbiAgY29uc3QgYm94ZXMgPSByaW5ncy5tYXAocmluZ0JveCk7XG4gIGxldCBhcmVhID0gMDtcbiAgZm9yIChsZXQgaSA9IDA7IGkgPCByaW5ncy5sZW5ndGg7IGkrKykge1xuICAgIGNvbnN0IHIgPSByaW5nc1tpXTtcbiAgICBpZiAoci5sZW5ndGggPCA2KSBjb250aW51ZTtcbiAgICBjb25zdCB4ID0gclswXSxcbiAgICAgIHkgPSByWzFdO1xuICAgIGxldCBkZXB0aCA9IDA7XG4gICAgZm9yIChsZXQgaiA9IDA7IGogPCByaW5ncy5sZW5ndGg7IGorKykge1xuICAgICAgaWYgKGogPT09IGkpIGNvbnRpbnVlO1xuICAgICAgY29uc3QgYiA9IGJveGVzW2pdO1xuICAgICAgaWYgKHggPCBiWzBdIHx8IHggPiBiWzJdIHx8IHkgPCBiWzFdIHx8IHkgPiBiWzNdKSBjb250aW51ZTtcbiAgICAgIGlmIChwb2ludEluUmluZyh4LCB5LCByaW5nc1tqXSkpIGRlcHRoKys7XG4gICAgfVxuICAgIGFyZWEgKz0gKGRlcHRoICUgMiA/IC0xIDogMSkgKiBNYXRoLmFicyhyaW5nQXJlYShyKSk7XG4gIH1cbiAgcmV0dXJuIE1hdGgubWF4KDAsIGFyZWEpO1xufVxuXG4vKipcbiAqIFRoZSByaW5ncyBvZiBlYWNoIGZlYXR1cmUgKHNlZSB7QGxpbmsgU2xpY2VPdXRsaW5lcy5yaW5nR3JvdXBzfSksIGluIHRoZSBvcmRlciBmaXJzdCBtZXQ7IG9uZVxuICogZmVhdHVyZSBvZiBldmVyeSByaW5nIHdpdGhvdXQgYGdyb3Vwc2AuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBmZWF0dXJlUmluZ3MoXG4gIHJpbmdzOiByZWFkb25seSBBcnJheUxpa2U8bnVtYmVyPltdLFxuICBncm91cHM/OiBBcnJheUxpa2U8bnVtYmVyPixcbik6IEFycmF5TGlrZTxudW1iZXI+W11bXSB7XG4gIGlmICghZ3JvdXBzIHx8ICFyaW5ncy5sZW5ndGgpIHJldHVybiByaW5ncy5sZW5ndGggPyBbcmluZ3Muc2xpY2UoKV0gOiBbXTtcbiAgY29uc3QgYnlHcm91cCA9IG5ldyBNYXA8bnVtYmVyLCBBcnJheUxpa2U8bnVtYmVyPltdPigpO1xuICByaW5ncy5mb3JFYWNoKChyLCBpKSA9PiB7XG4gICAgY29uc3QgZyA9IGdyb3Vwc1tpXSA/PyAwO1xuICAgIGNvbnN0IGZlYXR1cmUgPSBieUdyb3VwLmdldChnKTtcbiAgICBpZiAoZmVhdHVyZSkgZmVhdHVyZS5wdXNoKHIpO1xuICAgIGVsc2UgYnlHcm91cC5zZXQoZywgW3JdKTtcbiAgfSk7XG4gIHJldHVybiBbLi4uYnlHcm91cC52YWx1ZXMoKV07XG59XG5cbi8qKlxuICogVGhlIGFyZWEgb2YgYSBzbGljZSdzIG91dGxpbmVzIChweMKyKTogZWFjaCBmZWF0dXJlJ3MgKGV2ZW4tb2RkLCBzZWUge0BsaW5rIHJpbmdzQXJlYX0pLCB0aGVcbiAqIGZlYXR1cmVzIGFzIHRoZWlyIHVuaW9uLCBzbyBhIGZlYXR1cmUgaW5zaWRlIG9yIG92ZXIgYW5vdGhlciBvZiB0aGUgc2xpY2UgYWRkcyBvbmx5IHdoYXQgaXRcbiAqIGNvdmVycyBiZXlvbmQgaXQuIEZlYXR1cmVzIHdob3NlIGV4dGVudHMgb3ZlcmxhcCBhcmUgbWVhc3VyZWQgdG9nZXRoZXIsIGJ5IHRoZSBsZW5ndGggb2YgdGhlaXJcbiAqIHVuaW9uIGFsb25nIGVhY2ggcm93LCBleGFjdCBiZXR3ZWVuIHRoZSByb3dzIHRocm91Z2ggdGhlaXIgdmVydGljZXMgd2hlcmUgbm8gb3V0bGluZSBjcm9zc2VzXG4gKiBhbm90aGVyIChhbmQgcmVmaW5lZCB3aGVyZSBvbmUgZG9lcykuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBvdXRsaW5lc0FyZWEocmluZ3M6IHJlYWRvbmx5IEFycmF5TGlrZTxudW1iZXI+W10sIGdyb3Vwcz86IEFycmF5TGlrZTxudW1iZXI+KTogbnVtYmVyIHtcbiAgY29uc3QgZmVhdHVyZXMgPSBmZWF0dXJlUmluZ3MocmluZ3MsIGdyb3VwcykubWFwKChmKSA9PiAoeyByaW5nczogZiwgYm94OiBmZWF0dXJlQm94KGYpIH0pKTtcbiAgaWYgKGZlYXR1cmVzLmxlbmd0aCA8IDIpIHJldHVybiBmZWF0dXJlcy5sZW5ndGggPyByaW5nc0FyZWEoZmVhdHVyZXNbMF0ucmluZ3MpIDogMDtcbiAgLy8gRmVhdHVyZXMgd2hvc2UgZXh0ZW50cyBvdmVybGFwLCBqb2luZWQgKGJ5IHgsIHNvIGEgZmVhdHVyZSBtZWV0cyBvbmx5IHRob3NlIGl0IGNhbiBvdmVybGFwKS5cbiAgZmVhdHVyZXMuc29ydCgoYSwgYikgPT4gYS5ib3hbMF0gLSBiLmJveFswXSk7XG4gIGNvbnN0IHBhcmVudCA9IGZlYXR1cmVzLm1hcCgoXywgaSkgPT4gaSk7XG4gIGNvbnN0IHJvb3QgPSAoaTogbnVtYmVyKTogbnVtYmVyID0+IChwYXJlbnRbaV0gPT09IGkgPyBpIDogKHBhcmVudFtpXSA9IHJvb3QocGFyZW50W2ldKSkpO1xuICBmb3IgKGxldCBpID0gMDsgaSA8IGZlYXR1cmVzLmxlbmd0aDsgaSsrKSB7XG4gICAgY29uc3QgYSA9IGZlYXR1cmVzW2ldLmJveDtcbiAgICBmb3IgKGxldCBqID0gaSArIDE7IGogPCBmZWF0dXJlcy5sZW5ndGggJiYgZmVhdHVyZXNbal0uYm94WzBdIDwgYVsyXTsgaisrKSB7XG4gICAgICBjb25zdCBiID0gZmVhdHVyZXNbal0uYm94O1xuICAgICAgaWYgKGJbMV0gPCBhWzNdICYmIGFbMV0gPCBiWzNdKSBwYXJlbnRbcm9vdChqKV0gPSByb290KGkpO1xuICAgIH1cbiAgfVxuICBjb25zdCBjbHVzdGVycyA9IG5ldyBNYXA8bnVtYmVyLCBBcnJheUxpa2U8bnVtYmVyPltdW10+KCk7XG4gIGZlYXR1cmVzLmZvckVhY2goKGYsIGkpID0+IHtcbiAgICBjb25zdCByID0gcm9vdChpKTtcbiAgICBjb25zdCBjbHVzdGVyID0gY2x1c3RlcnMuZ2V0KHIpO1xuICAgIGlmIChjbHVzdGVyKSBjbHVzdGVyLnB1c2goZi5yaW5ncyk7XG4gICAgZWxzZSBjbHVzdGVycy5zZXQociwgW2YucmluZ3NdKTtcbiAgfSk7XG4gIGxldCBhcmVhID0gMDtcbiAgZm9yIChjb25zdCBjbHVzdGVyIG9mIGNsdXN0ZXJzLnZhbHVlcygpKVxuICAgIGFyZWEgKz0gY2x1c3Rlci5sZW5ndGggPT09IDEgPyByaW5nc0FyZWEoY2x1c3RlclswXSkgOiB1bmlvbkFyZWEoY2x1c3Rlcik7XG4gIHJldHVybiBhcmVhO1xufVxuXG4vKipcbiAqIFRoZSB2b2x1bWUgb2YgdGhlIHJlZ2lvbiBvdXRsaW5lZCBieSBgc2xpY2VzYCwgaW4gcHjCsiDDlyBzbGljZXMgKG9uZSBzbGljZSB0aGljayBwZXIgb3V0bGluZSk6XG4gKiB0aGUgc3VtIG9mIGVhY2ggc2xpY2UncyBhcmVhIChzZWUge0BsaW5rIG91dGxpbmVzQXJlYX07IENhdmFsaWVyaSdzIGVzdGltYXRvcikuIFRpbWVzIMK1bSBwZXIgcHjCslxuICogYW5kIHBlciBzbGljZSwgwrVtwrMuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiBvdXRsaW5lZFZvbHVtZShzbGljZXM6IHJlYWRvbmx5IFNsaWNlT3V0bGluZXNbXSk6IG51bWJlciB7XG4gIGxldCB2ID0gMDtcbiAgZm9yIChjb25zdCBzIG9mIHNsaWNlcykgdiArPSBvdXRsaW5lc0FyZWEocy5yaW5ncywgcy5yaW5nR3JvdXBzKTtcbiAgcmV0dXJuIHY7XG59XG5cbi8qKiBUaGUgYXJlYSBvZiB0aGUgdW5pb24gb2YgYGZlYXR1cmVzYCAoZWFjaCBldmVuLW9kZCBvdmVyIGl0cyByaW5ncyksIHJvdyBiYW5kIGJ5IHJvdyBiYW5kLiAqL1xuZnVuY3Rpb24gdW5pb25BcmVhKGZlYXR1cmVzOiByZWFkb25seSBBcnJheUxpa2U8bnVtYmVyPltdW10pOiBudW1iZXIge1xuICBjb25zdCB5czogbnVtYmVyW10gPSBbXTtcbiAgZm9yIChjb25zdCBmIG9mIGZlYXR1cmVzKSBmb3IgKGNvbnN0IHIgb2YgZikgZm9yIChsZXQgaSA9IDE7IGkgPCByLmxlbmd0aDsgaSArPSAyKSB5cy5wdXNoKHJbaV0pO1xuICB5cy5zb3J0KChhLCBiKSA9PiBhIC0gYik7XG4gIGxldCBhcmVhID0gMDtcbiAgZm9yIChsZXQgayA9IDA7IGsgKyAxIDwgeXMubGVuZ3RoOyBrKyspIGlmICh5c1trICsgMV0gPiB5c1trXSkgYXJlYSArPSBiYW5kQXJlYShmZWF0dXJlcywgeXNba10sIHlzW2sgKyAxXSwgMCk7XG4gIHJldHVybiBhcmVhO1xufVxuXG4vKipcbiAqIFRoZSB1bmlvbidzIGFyZWEgYmV0d2VlbiByb3dzIGBhYCBhbmQgYGJgLiBCZXR3ZWVuIHR3byB2ZXJ0aWNlcycgcm93cyBpdHMgd2lkdGggYWxvbmcgYSByb3cgaXNcbiAqIGxpbmVhciBpbiB5IHVubGVzcyBvdXRsaW5lcyBjcm9zcyB0aGVyZTogdGhlbiB0aGUgYmFuZCBpcyBoYWx2ZWQgdW50aWwgaXQgaXMgKG9yIHRoaW4gZW5vdWdoKS5cbiAqL1xuZnVuY3Rpb24gYmFuZEFyZWEoZmVhdHVyZXM6IHJlYWRvbmx5IEFycmF5TGlrZTxudW1iZXI+W11bXSwgYTogbnVtYmVyLCBiOiBudW1iZXIsIGRlcHRoOiBudW1iZXIpOiBudW1iZXIge1xuICBjb25zdCBoID0gYiAtIGE7XG4gIGNvbnN0IGF0ID0gKHQ6IG51bWJlcikgPT4gdW5pb25XaWR0aChmZWF0dXJlcywgYSArIHQgKiBoKTtcbiAgY29uc3QgbWlkID0gYXQoMC41KTtcbiAgaWYgKGRlcHRoID49IDEwKSByZXR1cm4gaCAqIG1pZDtcbiAgY29uc3QgZSA9IDFlLTY7XG4gIGNvbnN0IHRvbCA9IDFlLTkgKiAoMSArIE1hdGguYWJzKG1pZCkpO1xuICBjb25zdCBsaW5lYXIgPSBNYXRoLmFicyhhdCgwLjI1KSArIGF0KDAuNzUpIC0gMiAqIG1pZCkgPD0gdG9sICYmIE1hdGguYWJzKGF0KGUpICsgYXQoMSAtIGUpIC0gMiAqIG1pZCkgPD0gdG9sO1xuICBpZiAobGluZWFyKSByZXR1cm4gaCAqIG1pZDtcbiAgY29uc3QgbSA9IGEgKyBoIC8gMjtcbiAgcmV0dXJuIGJhbmRBcmVhKGZlYXR1cmVzLCBhLCBtLCBkZXB0aCArIDEpICsgYmFuZEFyZWEoZmVhdHVyZXMsIG0sIGIsIGRlcHRoICsgMSk7XG59XG5cbi8qKiBUaGUgbGVuZ3RoIG9mIHRoZSB1bmlvbiBvZiBgZmVhdHVyZXNgIGFsb25nIHRoZSByb3cgYHlgLiAqL1xuZnVuY3Rpb24gdW5pb25XaWR0aChmZWF0dXJlczogcmVhZG9ubHkgQXJyYXlMaWtlPG51bWJlcj5bXVtdLCB5OiBudW1iZXIpOiBudW1iZXIge1xuICBjb25zdCBzcGFuczogQXJyYXk8W251bWJlciwgbnVtYmVyXT4gPSBbXTtcbiAgY29uc3QgeHM6IG51bWJlcltdID0gW107XG4gIGZvciAoY29uc3QgZiBvZiBmZWF0dXJlcykge1xuICAgIHhzLmxlbmd0aCA9IDA7XG4gICAgZm9yIChjb25zdCByIG9mIGYpIHtcbiAgICAgIGNvbnN0IG4gPSByLmxlbmd0aCA+PiAxO1xuICAgICAgZm9yIChsZXQgaSA9IDAsIGsgPSBuIC0gMTsgaSA8IG47IGsgPSBpKyspIHtcbiAgICAgICAgY29uc3QgeWEgPSByWzIgKiBpICsgMV0sXG4gICAgICAgICAgeWIgPSByWzIgKiBrICsgMV07XG4gICAgICAgIGlmICh5YSA+IHkgIT09IHliID4geSkgeHMucHVzaChyWzIgKiBpXSArICgoeSAtIHlhKSAqIChyWzIgKiBrXSAtIHJbMiAqIGldKSkgLyAoeWIgLSB5YSkpO1xuICAgICAgfVxuICAgIH1cbiAgICB4cy5zb3J0KChwLCBxKSA9PiBwIC0gcSk7XG4gICAgZm9yIChsZXQgYyA9IDA7IGMgKyAxIDwgeHMubGVuZ3RoOyBjICs9IDIpIHNwYW5zLnB1c2goW3hzW2NdLCB4c1tjICsgMV1dKTtcbiAgfVxuICBzcGFucy5zb3J0KChwLCBxKSA9PiBwWzBdIC0gcVswXSk7XG4gIGxldCB3aWR0aCA9IDAsXG4gICAgZnJvbSA9IC1JbmZpbml0eSxcbiAgICB0byA9IC1JbmZpbml0eTtcbiAgZm9yIChjb25zdCBbczAsIHMxXSBvZiBzcGFucykge1xuICAgIGlmIChzMCA+IHRvKSB7XG4gICAgICBpZiAodG8gPiBmcm9tKSB3aWR0aCArPSB0byAtIGZyb207XG4gICAgICBmcm9tID0gczA7XG4gICAgICB0byA9IHMxO1xuICAgIH0gZWxzZSBpZiAoczEgPiB0bykgdG8gPSBzMTtcbiAgfVxuICBpZiAodG8gPiBmcm9tKSB3aWR0aCArPSB0byAtIGZyb207XG4gIHJldHVybiB3aWR0aDtcbn1cblxuLyoqIFRoZSBleHRlbnQgb2YgYSBmZWF0dXJlJ3MgcmluZ3MuICovXG5mdW5jdGlvbiBmZWF0dXJlQm94KHJpbmdzOiByZWFkb25seSBBcnJheUxpa2U8bnVtYmVyPltdKTogW251bWJlciwgbnVtYmVyLCBudW1iZXIsIG51bWJlcl0ge1xuICBjb25zdCBib3g6IFtudW1iZXIsIG51bWJlciwgbnVtYmVyLCBudW1iZXJdID0gW0luZmluaXR5LCBJbmZpbml0eSwgLUluZmluaXR5LCAtSW5maW5pdHldO1xuICBmb3IgKGNvbnN0IHIgb2YgcmluZ3MpIHtcbiAgICBjb25zdCBbeDAsIHkwLCB4MSwgeTFdID0gcmluZ0JveChyKTtcbiAgICBib3hbMF0gPSBNYXRoLm1pbihib3hbMF0sIHgwKTtcbiAgICBib3hbMV0gPSBNYXRoLm1pbihib3hbMV0sIHkwKTtcbiAgICBib3hbMl0gPSBNYXRoLm1heChib3hbMl0sIHgxKTtcbiAgICBib3hbM10gPSBNYXRoLm1heChib3hbM10sIHkxKTtcbiAgfVxuICByZXR1cm4gYm94O1xufVxuXG5mdW5jdGlvbiByaW5nQm94KHI6IEFycmF5TGlrZTxudW1iZXI+KTogW251bWJlciwgbnVtYmVyLCBudW1iZXIsIG51bWJlcl0ge1xuICBsZXQgeDAgPSBJbmZpbml0eSxcbiAgICB5MCA9IEluZmluaXR5LFxuICAgIHgxID0gLUluZmluaXR5LFxuICAgIHkxID0gLUluZmluaXR5O1xuICBmb3IgKGxldCBpID0gMDsgaSArIDEgPCByLmxlbmd0aDsgaSArPSAyKSB7XG4gICAgaWYgKHJbaV0gPCB4MCkgeDAgPSByW2ldO1xuICAgIGlmIChyW2ldID4geDEpIHgxID0gcltpXTtcbiAgICBpZiAocltpICsgMV0gPCB5MCkgeTAgPSByW2kgKyAxXTtcbiAgICBpZiAocltpICsgMV0gPiB5MSkgeTEgPSByW2kgKyAxXTtcbiAgfVxuICByZXR1cm4gW3gwLCB5MCwgeDEsIHkxXTtcbn1cbiJdfQ==

package/esm2022/lib/geometry/reslice-image.mjs

@@ -0,0 +1,57 @@
+/**
+ * An image of a plane through a volume (the X or Y section beside the 3D view), cut from the
+ * volume's voxels and displayed the way napari-js displays the volume: each channel windowed,
+ * inverted and gamma-corrected (napari-js's `windowGamma`) and coloured by its colormap's LUT
+ * (`buildLut`), the channels added (an additive composite, clamped).
+ *
+ * Pure: no DOM, no GPU.
+ */
+import { buildLut, windowGamma } from 'napari-js';
+/**
+ * The image of the plane `axis = index` (a voxel index): through x, slices across and y down;
+ * through y, x across and slices down — the section coordinates of `geometry/section`.
+ */
+export function resliceImage(vol, axis, index) {
+    const { width: w, height: h, depth: d } = vol;
+    const at = Math.min(Math.max(0, Math.round(index)), (axis === 'x' ? w : h) - 1);
+    const outW = axis === 'x' ? d : w;
+    const outH = axis === 'x' ? h : d;
+    const rgb = new Float32Array(outW * outH * 3);
+    for (const ch of vol.channels) {
+        if (!ch.visible)
+            continue;
+        const lut = channelLut(ch, vol.invert);
+        for (let row = 0; row < outH; row++) {
+            for (let col = 0; col < outW; col++) {
+                const voxel = axis === 'x' ? at + w * (row + h * col) : col + w * (at + h * row);
+                const k = (row * outW + col) * 3;
+                const v = ch.data[voxel] * 3;
+                rgb[k] += lut[v];
+                rgb[k + 1] += lut[v + 1];
+                rgb[k + 2] += lut[v + 2];
+            }
+        }
+    }
+    const data = new Uint8ClampedArray(outW * outH * 4);
+    for (let i = 0, j = 0; i < rgb.length; i += 3, j += 4) {
+        data[j] = rgb[i];
+        data[j + 1] = rgb[i + 1];
+        data[j + 2] = rgb[i + 2];
+        data[j + 3] = 255;
+    }
+    return { width: outW, height: outH, data };
+}
+/** Voxel value (0..255) → displayed RGB (0..255 each), as 256 RGB triples. */
+function channelLut(ch, invert) {
+    const colors = buildLut(ch.colormap); // RGBA8 × 256
+    const last = colors.length / 4 - 1;
+    const lut = new Float32Array(256 * 3);
+    for (let v = 0; v < 256; v++) {
+        const i = 4 * Math.round(windowGamma(v, ch.min, ch.max, ch.gamma, invert) * last);
+        lut[3 * v] = colors[i];
+        lut[3 * v + 1] = colors[i + 1];
+        lut[3 * v + 2] = colors[i + 2];
+    }
+    return lut;
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoicmVzbGljZS1pbWFnZS5qcyIsInNvdXJjZVJvb3QiOiIiLCJzb3VyY2VzIjpbIi4uLy4uLy4uLy4uL3NyYy9saWIvZ2VvbWV0cnkvcmVzbGljZS1pbWFnZS50cyJdLCJuYW1lcyI6W10sIm1hcHBpbmdzIjoiQUFBQTs7Ozs7OztHQU9HO0FBQ0gsT0FBTyxFQUFFLFFBQVEsRUFBRSxXQUFXLEVBQUUsTUFBTSxXQUFXLENBQUM7QUFrQ2xEOzs7R0FHRztBQUNILE1BQU0sVUFBVSxZQUFZLENBQUMsR0FBaUIsRUFBRSxJQUFpQixFQUFFLEtBQWE7SUFDOUUsTUFBTSxFQUFFLEtBQUssRUFBRSxDQUFDLEVBQUUsTUFBTSxFQUFFLENBQUMsRUFBRSxLQUFLLEVBQUUsQ0FBQyxFQUFFLEdBQUcsR0FBRyxDQUFDO0lBQzlDLE1BQU0sRUFBRSxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEVBQUUsSUFBSSxDQUFDLEtBQUssQ0FBQyxLQUFLLENBQUMsQ0FBQyxFQUFFLENBQUMsSUFBSSxLQUFLLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztJQUNoRixNQUFNLElBQUksR0FBRyxJQUFJLEtBQUssR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUNsQyxNQUFNLElBQUksR0FBRyxJQUFJLEtBQUssR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUNsQyxNQUFNLEdBQUcsR0FBRyxJQUFJLFlBQVksQ0FBQyxJQUFJLEdBQUcsSUFBSSxHQUFHLENBQUMsQ0FBQyxDQUFDO0lBQzlDLEtBQUssTUFBTSxFQUFFLElBQUksR0FBRyxDQUFDLFFBQVEsRUFBRSxDQUFDO1FBQzlCLElBQUksQ0FBQyxFQUFFLENBQUMsT0FBTztZQUFFLFNBQVM7UUFDMUIsTUFBTSxHQUFHLEdBQUcsVUFBVSxDQUFDLEVBQUUsRUFBRSxHQUFHLENBQUMsTUFBTSxDQUFDLENBQUM7UUFDdkMsS0FBSyxJQUFJLEdBQUcsR0FBRyxDQUFDLEVBQUUsR0FBRyxHQUFHLElBQUksRUFBRSxHQUFHLEVBQUUsRUFBRSxDQUFDO1lBQ3BDLEtBQUssSUFBSSxHQUFHLEdBQUcsQ0FBQyxFQUFFLEdBQUcsR0FBRyxJQUFJLEVBQUUsR0FBRyxFQUFFLEVBQUUsQ0FBQztnQkFDcEMsTUFBTSxLQUFLLEdBQUcsSUFBSSxLQUFLLEdBQUcsQ0FBQyxDQUFDLENBQUMsRUFBRSxHQUFHLENBQUMsR0FBRyxDQUFDLEdBQUcsR0FBRyxDQUFDLEdBQUcsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsR0FBRyxDQUFDLEdBQUcsQ0FBQyxFQUFFLEdBQUcsQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDO2dCQUNqRixNQUFNLENBQUMsR0FBRyxDQUFDLEdBQUcsR0FBRyxJQUFJLEdBQUcsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDO2dCQUNqQyxNQUFNLENBQUMsR0FBRyxFQUFFLENBQUMsSUFBSSxDQUFDLEtBQUssQ0FBQyxHQUFHLENBQUMsQ0FBQztnQkFDN0IsR0FBRyxDQUFDLENBQUMsQ0FBQyxJQUFJLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztnQkFDakIsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsSUFBSSxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO2dCQUN6QixHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxJQUFJLEdBQUcsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7WUFDM0IsQ0FBQztRQUNILENBQUM7SUFDSCxDQUFDO0lBQ0QsTUFBTSxJQUFJLEdBQUcsSUFBSSxpQkFBaUIsQ0FBQyxJQUFJLEdBQUcsSUFBSSxHQUFHLENBQUMsQ0FBQyxDQUFDO0lBQ3BELEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLEdBQUcsQ0FBQyxNQUFNLEVBQUUsQ0FBQyxJQUFJLENBQUMsRUFBRSxDQUFDLElBQUksQ0FBQyxFQUFFLENBQUM7UUFDdEQsSUFBSSxDQUFDLENBQUMsQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNqQixJQUFJLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUM7UUFDekIsSUFBSSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ3pCLElBQUksQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsR0FBRyxDQUFDO0lBQ3BCLENBQUM7SUFDRCxPQUFPLEVBQUUsS0FBSyxFQUFFLElBQUksRUFBRSxNQUFNLEVBQUUsSUFBSSxFQUFFLElBQUksRUFBRSxDQUFDO0FBQzdDLENBQUM7QUFFRCw4RUFBOEU7QUFDOUUsU0FBUyxVQUFVLENBQUMsRUFBdUIsRUFBRSxNQUFlO0lBQzFELE1BQU0sTUFBTSxHQUFHLFFBQVEsQ0FBQyxFQUFFLENBQUMsUUFBUSxDQUFDLENBQUMsQ0FBQyxjQUFjO0lBQ3BELE1BQU0sSUFBSSxHQUFHLE1BQU0sQ0FBQyxNQUFNLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQztJQUNuQyxNQUFNLEdBQUcsR0FBRyxJQUFJLFlBQVksQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUM7SUFDdEMsS0FBSyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLEdBQUcsRUFBRSxDQUFDLEVBQUUsRUFBRSxDQUFDO1FBQzdCLE1BQU0sQ0FBQyxHQUFHLENBQUMsR0FBRyxJQUFJLENBQUMsS0FBSyxDQUFDLFdBQVcsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLEdBQUcsRUFBRSxFQUFFLENBQUMsR0FBRyxFQUFFLEVBQUUsQ0FBQyxLQUFLLEVBQUUsTUFBTSxDQUFDLEdBQUcsSUFBSSxDQUFDLENBQUM7UUFDbEYsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxNQUFNLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDdkIsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsTUFBTSxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztRQUMvQixHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxNQUFNLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO0lBQ2pDLENBQUM7SUFDRCxPQUFPLEdBQUcsQ0FBQztBQUNiLENBQUMiLCJzb3VyY2VzQ29udGVudCI6WyIvKipcbiAqIEFuIGltYWdlIG9mIGEgcGxhbmUgdGhyb3VnaCBhIHZvbHVtZSAodGhlIFggb3IgWSBzZWN0aW9uIGJlc2lkZSB0aGUgM0QgdmlldyksIGN1dCBmcm9tIHRoZVxuICogdm9sdW1lJ3Mgdm94ZWxzIGFuZCBkaXNwbGF5ZWQgdGhlIHdheSBuYXBhcmktanMgZGlzcGxheXMgdGhlIHZvbHVtZTogZWFjaCBjaGFubmVsIHdpbmRvd2VkLFxuICogaW52ZXJ0ZWQgYW5kIGdhbW1hLWNvcnJlY3RlZCAobmFwYXJpLWpzJ3MgYHdpbmRvd0dhbW1hYCkgYW5kIGNvbG91cmVkIGJ5IGl0cyBjb2xvcm1hcCdzIExVVFxuICogKGBidWlsZEx1dGApLCB0aGUgY2hhbm5lbHMgYWRkZWQgKGFuIGFkZGl0aXZlIGNvbXBvc2l0ZSwgY2xhbXBlZCkuXG4gKlxuICogUHVyZTogbm8gRE9NLCBubyBHUFUuXG4gKi9cbmltcG9ydCB7IGJ1aWxkTHV0LCB3aW5kb3dHYW1tYSB9IGZyb20gJ25hcGFyaS1qcyc7XG5pbXBvcnQgdHlwZSB7IENvbG9ybWFwIH0gZnJvbSAnbmFwYXJpLWpzJztcblxuaW1wb3J0IHR5cGUgeyBTZWN0aW9uQXhpcyB9IGZyb20gJy4vc2VjdGlvbic7XG5cbi8qKiBPbmUgY2hhbm5lbCBvZiBhIHZvbHVtZSBhcyBpdCBpcyBkaXNwbGF5ZWQuICovXG5leHBvcnQgaW50ZXJmYWNlIFZvbHVtZUNoYW5uZWxWb3hlbHMge1xuICAvKiogdWludDggdm94ZWxzLCB4IGZhc3Rlc3QsIHRoZW4geSwgdGhlbiB6LiAqL1xuICBkYXRhOiBVaW50OEFycmF5O1xuICAvKiogVGhlIGNvbG9ybWFwIHRoZSAzRCB2aWV3IGRyYXdzIHRoZSBjaGFubmVsIHdpdGggKGEgdGludCwgb3IgdGhlIGRpc3BsYXkgY29sb3JtYXApLiAqL1xuICBjb2xvcm1hcDogQ29sb3JtYXA7XG4gIC8qKiBEaXNwbGF5IHdpbmRvdyBvdmVyIDAuLjI1NS4gKi9cbiAgbWluOiBudW1iZXI7XG4gIG1heDogbnVtYmVyO1xuICAvKiogbmFwYXJpLWpzIGdhbW1hIChzZWUgYHRvTmFwYXJpR2FtbWFgKTogdGhlIHdpbmRvd2VkIHZhbHVlIHQgc2hvd3MgYXMgdF7Osy4gKi9cbiAgZ2FtbWE6IG51bWJlcjtcbiAgdmlzaWJsZTogYm9vbGVhbjtcbn1cblxuZXhwb3J0IGludGVyZmFjZSBWb2x1bWVWb3hlbHMge1xuICB3aWR0aDogbnVtYmVyO1xuICBoZWlnaHQ6IG51bWJlcjtcbiAgZGVwdGg6IG51bWJlcjtcbiAgY2hhbm5lbHM6IHJlYWRvbmx5IFZvbHVtZUNoYW5uZWxWb3hlbHNbXTtcbiAgaW52ZXJ0OiBib29sZWFuO1xufVxuXG4vKiogQW4gUkdCQSBpbWFnZSwgcm93IGJ5IHJvdy4gKi9cbmV4cG9ydCBpbnRlcmZhY2UgUmdiYUltYWdlIHtcbiAgd2lkdGg6IG51bWJlcjtcbiAgaGVpZ2h0OiBudW1iZXI7XG4gIGRhdGE6IFVpbnQ4Q2xhbXBlZEFycmF5O1xufVxuXG4vKipcbiAqIFRoZSBpbWFnZSBvZiB0aGUgcGxhbmUgYGF4aXMgPSBpbmRleGAgKGEgdm94ZWwgaW5kZXgpOiB0aHJvdWdoIHgsIHNsaWNlcyBhY3Jvc3MgYW5kIHkgZG93bjtcbiAqIHRocm91Z2ggeSwgeCBhY3Jvc3MgYW5kIHNsaWNlcyBkb3duIOKAlCB0aGUgc2VjdGlvbiBjb29yZGluYXRlcyBvZiBgZ2VvbWV0cnkvc2VjdGlvbmAuXG4gKi9cbmV4cG9ydCBmdW5jdGlvbiByZXNsaWNlSW1hZ2Uodm9sOiBWb2x1bWVWb3hlbHMsIGF4aXM6IFNlY3Rpb25BeGlzLCBpbmRleDogbnVtYmVyKTogUmdiYUltYWdlIHtcbiAgY29uc3QgeyB3aWR0aDogdywgaGVpZ2h0OiBoLCBkZXB0aDogZCB9ID0gdm9sO1xuICBjb25zdCBhdCA9IE1hdGgubWluKE1hdGgubWF4KDAsIE1hdGgucm91bmQoaW5kZXgpKSwgKGF4aXMgPT09ICd4JyA/IHcgOiBoKSAtIDEpO1xuICBjb25zdCBvdXRXID0gYXhpcyA9PT0gJ3gnID8gZCA6IHc7XG4gIGNvbnN0IG91dEggPSBheGlzID09PSAneCcgPyBoIDogZDtcbiAgY29uc3QgcmdiID0gbmV3IEZsb2F0MzJBcnJheShvdXRXICogb3V0SCAqIDMpO1xuICBmb3IgKGNvbnN0IGNoIG9mIHZvbC5jaGFubmVscykge1xuICAgIGlmICghY2gudmlzaWJsZSkgY29udGludWU7XG4gICAgY29uc3QgbHV0ID0gY2hhbm5lbEx1dChjaCwgdm9sLmludmVydCk7XG4gICAgZm9yIChsZXQgcm93ID0gMDsgcm93IDwgb3V0SDsgcm93KyspIHtcbiAgICAgIGZvciAobGV0IGNvbCA9IDA7IGNvbCA8IG91dFc7IGNvbCsrKSB7XG4gICAgICAgIGNvbnN0IHZveGVsID0gYXhpcyA9PT0gJ3gnID8gYXQgKyB3ICogKHJvdyArIGggKiBjb2wpIDogY29sICsgdyAqIChhdCArIGggKiByb3cpO1xuICAgICAgICBjb25zdCBrID0gKHJvdyAqIG91dFcgKyBjb2wpICogMztcbiAgICAgICAgY29uc3QgdiA9IGNoLmRhdGFbdm94ZWxdICogMztcbiAgICAgICAgcmdiW2tdICs9IGx1dFt2XTtcbiAgICAgICAgcmdiW2sgKyAxXSArPSBsdXRbdiArIDFdO1xuICAgICAgICByZ2JbayArIDJdICs9IGx1dFt2ICsgMl07XG4gICAgICB9XG4gICAgfVxuICB9XG4gIGNvbnN0IGRhdGEgPSBuZXcgVWludDhDbGFtcGVkQXJyYXkob3V0VyAqIG91dEggKiA0KTtcbiAgZm9yIChsZXQgaSA9IDAsIGogPSAwOyBpIDwgcmdiLmxlbmd0aDsgaSArPSAzLCBqICs9IDQpIHtcbiAgICBkYXRhW2pdID0gcmdiW2ldO1xuICAgIGRhdGFbaiArIDFdID0gcmdiW2kgKyAxXTtcbiAgICBkYXRhW2ogKyAyXSA9IHJnYltpICsgMl07XG4gICAgZGF0YVtqICsgM10gPSAyNTU7XG4gIH1cbiAgcmV0dXJuIHsgd2lkdGg6IG91dFcsIGhlaWdodDogb3V0SCwgZGF0YSB9O1xufVxuXG4vKiogVm94ZWwgdmFsdWUgKDAuLjI1NSkg4oaSIGRpc3BsYXllZCBSR0IgKDAuLjI1NSBlYWNoKSwgYXMgMjU2IFJHQiB0cmlwbGVzLiAqL1xuZnVuY3Rpb24gY2hhbm5lbEx1dChjaDogVm9sdW1lQ2hhbm5lbFZveGVscywgaW52ZXJ0OiBib29sZWFuKTogRmxvYXQzMkFycmF5IHtcbiAgY29uc3QgY29sb3JzID0gYnVpbGRMdXQoY2guY29sb3JtYXApOyAvLyBSR0JBOCDDlyAyNTZcbiAgY29uc3QgbGFzdCA9IGNvbG9ycy5sZW5ndGggLyA0IC0gMTtcbiAgY29uc3QgbHV0ID0gbmV3IEZsb2F0MzJBcnJheSgyNTYgKiAzKTtcbiAgZm9yIChsZXQgdiA9IDA7IHYgPCAyNTY7IHYrKykge1xuICAgIGNvbnN0IGkgPSA0ICogTWF0aC5yb3VuZCh3aW5kb3dHYW1tYSh2LCBjaC5taW4sIGNoLm1heCwgY2guZ2FtbWEsIGludmVydCkgKiBsYXN0KTtcbiAgICBsdXRbMyAqIHZdID0gY29sb3JzW2ldO1xuICAgIGx1dFszICogdiArIDFdID0gY29sb3JzW2kgKyAxXTtcbiAgICBsdXRbMyAqIHYgKyAyXSA9IGNvbG9yc1tpICsgMl07XG4gIH1cbiAgcmV0dXJuIGx1dDtcbn1cbiJdfQ==