@performant-software/geospatial 3.1.28-beta.0 → 3.1.29

package/dist/index.cjs558.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs558.js","sources":["../../../node_modules/@allmaps/triangulate/node_modules/kdbush/index.js"],"sourcesContent":["\nconst ARRAY_TYPES = [\n    Int8Array, Uint8Array, Uint8ClampedArray, Int16Array, Uint16Array,\n    Int32Array, Uint32Array, Float32Array, Float64Array\n];\n\n/** @typedef {Int8ArrayConstructor | Uint8ArrayConstructor | Uint8ClampedArrayConstructor | Int16ArrayConstructor | Uint16ArrayConstructor | Int32ArrayConstructor | Uint32ArrayConstructor | Float32ArrayConstructor | Float64ArrayConstructor} TypedArrayConstructor */\n/** @typedef {Int8Array | Uint8Array | Uint8ClampedArray | Int16Array | Uint16Array | Int32Array | Uint32Array | Float32Array | Float64Array} TypedArray */\n\nconst VERSION = 1; // serialized format version\nconst HEADER_SIZE = 8;\n\n// Shared scratch stack for iterative DFS in range/within. Sized for the worst case:\n// 3 ints per frame * (treeHeight + 1), with treeHeight ≤ ceil(log2(2^32 / 3)) ≈ 31.\nconst STACK = new Uint32Array(96);\n\nexport default class KDBush {\n\n    /**\n     * Creates an index from raw `ArrayBuffer` data.\n     * @param {ArrayBufferLike} data\n     */\n    static from(data) {\n        // @ts-expect-error duck typing array buffers\n        if (!data || data.byteLength === undefined || data.buffer) {\n            throw new Error('Data must be an instance of ArrayBuffer or SharedArrayBuffer.');\n        }\n        const [magic, versionAndType] = new Uint8Array(data, 0, 2);\n        if (magic !== 0xdb) {\n            throw new Error('Data does not appear to be in a KDBush format.');\n        }\n        const version = versionAndType >> 4;\n        if (version !== VERSION) {\n            throw new Error(`Got v${version} data when expected v${VERSION}.`);\n        }\n        const ArrayType = ARRAY_TYPES[versionAndType & 0x0f];\n        if (!ArrayType) {\n            throw new Error('Unrecognized array type.');\n        }\n        const [nodeSize] = new Uint16Array(data, 2, 1);\n        const [numItems] = new Uint32Array(data, 4, 1);\n\n        return new KDBush(numItems, nodeSize, ArrayType, undefined, data);\n    }\n\n    /**\n     * Creates an index that will hold a given number of items.\n     * @param {number} numItems\n     * @param {number} [nodeSize=64] Size of the KD-tree node (64 by default).\n     * @param {TypedArrayConstructor} [ArrayType=Float64Array] The array type used for coordinates storage (`Float64Array` by default).\n     * @param {ArrayBufferConstructor | SharedArrayBufferConstructor} [ArrayBufferType=ArrayBuffer] The array buffer type used for storage (`ArrayBuffer` by default).\n     * @param {ArrayBufferLike} [data] (For internal use only)\n     */\n    constructor(numItems, nodeSize = 64, ArrayType = Float64Array, ArrayBufferType = ArrayBuffer, data) {\n        if (isNaN(numItems) || numItems < 0) throw new Error(`Unexpected numItems value: ${numItems}.`);\n\n        this.numItems = +numItems;\n        this.nodeSize = Math.min(Math.max(+nodeSize, 2), 65535);\n        this.ArrayType = ArrayType;\n        this.IndexArrayType = numItems < 65536 ? Uint16Array : Uint32Array;\n\n        const arrayTypeIndex = ARRAY_TYPES.indexOf(this.ArrayType);\n        const coordsByteSize = numItems * 2 * this.ArrayType.BYTES_PER_ELEMENT;\n        const idsByteSize = numItems * this.IndexArrayType.BYTES_PER_ELEMENT;\n        const padCoords = (8 - idsByteSize % 8) % 8;\n\n        if (arrayTypeIndex < 0) {\n            throw new Error(`Unexpected typed array class: ${ArrayType}.`);\n        }\n\n        if (data) { // reconstruct an index from a buffer\n            this.data = data;\n            // @ts-expect-error TS can't handle SharedArrayBuffer overloads\n            this.ids = new this.IndexArrayType(data, HEADER_SIZE, numItems);\n            // @ts-expect-error TS can't handle SharedArrayBuffer overloads\n            this.coords = new ArrayType(data, HEADER_SIZE + idsByteSize + padCoords, numItems * 2);\n            this._pos = numItems * 2;\n            this._finished = true;\n\n        } else { // initialize a new index\n            const data = this.data = new ArrayBufferType(HEADER_SIZE + coordsByteSize + idsByteSize + padCoords);\n            // @ts-expect-error TS can't handle SharedArrayBuffer overloads\n            this.ids = new this.IndexArrayType(data, HEADER_SIZE, numItems);\n            // @ts-expect-error TS can't handle SharedArrayBuffer overloads\n            this.coords = new ArrayType(data, HEADER_SIZE + idsByteSize + padCoords, numItems * 2);\n            this._pos = 0;\n            this._finished = false;\n\n            // set header\n            new Uint8Array(data, 0, 2).set([0xdb, (VERSION << 4) + arrayTypeIndex]);\n            new Uint16Array(data, 2, 1)[0] = nodeSize;\n            new Uint32Array(data, 4, 1)[0] = numItems;\n        }\n    }\n\n    /**\n     * Add a point to the index.\n     * @param {number} x\n     * @param {number} y\n     * @returns {number} An incremental index associated with the added item (starting from `0`).\n     */\n    add(x, y) {\n        const index = this._pos >> 1;\n        this.ids[index] = index;\n        this.coords[this._pos++] = x;\n        this.coords[this._pos++] = y;\n        return index;\n    }\n\n    /**\n     * Perform indexing of the added points.\n     */\n    finish() {\n        const numAdded = this._pos >> 1;\n        if (numAdded !== this.numItems) {\n            throw new Error(`Added ${numAdded} items when expected ${this.numItems}.`);\n        }\n        // kd-sort both arrays for efficient search\n        sort(this.ids, this.coords, this.nodeSize, 0, this.numItems - 1, 0);\n\n        this._finished = true;\n        return this;\n    }\n\n    /**\n     * Search the index for items within a given bounding box.\n     * @param {number} minX\n     * @param {number} minY\n     * @param {number} maxX\n     * @param {number} maxY\n     * @returns {number[]} An array of indices correponding to the found items.\n     */\n    range(minX, minY, maxX, maxY) {\n        if (!this._finished) throw new Error('Data not yet indexed - call index.finish().');\n\n        const {ids, coords, nodeSize} = this;\n        STACK[0] = 0;\n        STACK[1] = ids.length - 1;\n        STACK[2] = 0;\n        let sp = 3;\n        const result = [];\n\n        // recursively search for items in range in the kd-sorted arrays\n        while (sp > 0) {\n            const axis = STACK[--sp];\n            const right = STACK[--sp];\n            const left = STACK[--sp];\n\n            // if we reached \"tree node\", search linearly\n            if (right - left <= nodeSize) {\n                for (let i = left; i <= right; i++) {\n                    const x = coords[2 * i];\n                    const y = coords[2 * i + 1];\n                    if (x >= minX && x <= maxX && y >= minY && y <= maxY) result.push(ids[i]);\n                }\n                continue;\n            }\n\n            // otherwise find the middle index\n            const m = (left + right) >> 1;\n\n            // include the middle item if it's in range\n            const x = coords[2 * m];\n            const y = coords[2 * m + 1];\n            if (x >= minX && x <= maxX && y >= minY && y <= maxY) result.push(ids[m]);\n\n            // queue search in halves that intersect the query\n            if (axis === 0 ? minX <= x : minY <= y) {\n                STACK[sp++] = left;\n                STACK[sp++] = m - 1;\n                STACK[sp++] = 1 - axis;\n            }\n            if (axis === 0 ? maxX >= x : maxY >= y) {\n                STACK[sp++] = m + 1;\n                STACK[sp++] = right;\n                STACK[sp++] = 1 - axis;\n            }\n        }\n\n        return result;\n    }\n\n    /**\n     * Search the index for items within a given radius.\n     * @param {number} qx\n     * @param {number} qy\n     * @param {number} r Query radius.\n     * @returns {number[]} An array of indices correponding to the found items.\n     */\n    within(qx, qy, r) {\n        const result = /** @type {number[]} */ ([]);\n        this.withinInto(qx, qy, r, result);\n        return result;\n    }\n\n    /**\n     * Search the index for items within a given radius, writing matching ids into `out`\n     * via indexed assignment (`out[i] = id`). Accepts any indexed-writable container —\n     * a typed array sized to the expected upper bound (allocation-free, fast) or a plain\n     * `Array` (which will grow as needed). Returns the number of matches written.\n     * @param {number} qx\n     * @param {number} qy\n     * @param {number} r Query radius.\n     * @param {number[] | TypedArray} out Container to write matching ids into.\n     * @returns {number} The number of matches written to `out`.\n     */\n    withinInto(qx, qy, r, out) {\n        if (!this._finished) throw new Error('Data not yet indexed - call index.finish().');\n\n        const {ids, coords, nodeSize} = this;\n        STACK[0] = 0;\n        STACK[1] = ids.length - 1;\n        STACK[2] = 0;\n        let sp = 3;\n        let count = 0;\n        const r2 = r * r;\n\n        // recursively search for items within radius in the kd-sorted arrays\n        while (sp > 0) {\n            const axis = STACK[--sp];\n            const right = STACK[--sp];\n            const left = STACK[--sp];\n\n            // if we reached \"tree node\", search linearly\n            if (right - left <= nodeSize) {\n                for (let i = left; i <= right; i++) {\n                    if (sqDist(coords[2 * i], coords[2 * i + 1], qx, qy) <= r2) out[count++] = ids[i];\n                }\n                continue;\n            }\n\n            // otherwise find the middle index\n            const m = (left + right) >> 1;\n\n            // include the middle item if it's in range\n            const x = coords[2 * m];\n            const y = coords[2 * m + 1];\n            if (sqDist(x, y, qx, qy) <= r2) out[count++] = ids[m];\n\n            // queue search in halves that intersect the query\n            if (axis === 0 ? qx - r <= x : qy - r <= y) {\n                STACK[sp++] = left;\n                STACK[sp++] = m - 1;\n                STACK[sp++] = 1 - axis;\n            }\n            if (axis === 0 ? qx + r >= x : qy + r >= y) {\n                STACK[sp++] = m + 1;\n                STACK[sp++] = right;\n                STACK[sp++] = 1 - axis;\n            }\n        }\n\n        return count;\n    }\n}\n\n/**\n * @param {Uint16Array | Uint32Array} ids\n * @param {TypedArray} coords\n * @param {number} nodeSize\n * @param {number} left\n * @param {number} right\n * @param {number} axis\n */\nfunction sort(ids, coords, nodeSize, left, right, axis) {\n    if (right - left <= nodeSize) return;\n\n    const m = (left + right) >> 1; // middle index\n\n    // sort ids and coords around the middle index so that the halves lie\n    // either left/right or top/bottom correspondingly (taking turns)\n    select(ids, coords, m, left, right, axis);\n\n    // recursively kd-sort first half and second half on the opposite axis\n    sort(ids, coords, nodeSize, left, m - 1, 1 - axis);\n    sort(ids, coords, nodeSize, m + 1, right, 1 - axis);\n}\n\n/**\n * Custom Floyd-Rivest selection algorithm: sort ids and coords so that\n * [left..k-1] items are smaller than k-th item (on either x or y axis)\n * @param {Uint16Array | Uint32Array} ids\n * @param {TypedArray} coords\n * @param {number} k\n * @param {number} left\n * @param {number} right\n * @param {number} axis\n */\nfunction select(ids, coords, k, left, right, axis) {\n\n    while (right > left) {\n        if (right - left > 600) {\n            const n = right - left + 1;\n            const m = k - left + 1;\n            const z = Math.log(n);\n            const s = 0.5 * Math.exp(2 * z / 3);\n            const sd = 0.5 * Math.sqrt(z * s * (n - s) / n) * (m - n / 2 < 0 ? -1 : 1);\n            const newLeft = Math.max(left, Math.floor(k - m * s / n + sd));\n            const newRight = Math.min(right, Math.floor(k + (n - m) * s / n + sd));\n            select(ids, coords, k, newLeft, newRight, axis);\n        }\n\n        const t = coords[2 * k + axis];\n        let i = left;\n        let j = right;\n\n        swapItem(ids, coords, left, k);\n        if (coords[2 * right + axis] > t) swapItem(ids, coords, left, right);\n\n        while (i < j) {\n            swapItem(ids, coords, i, j);\n            i++;\n            j--;\n            while (coords[2 * i + axis] < t) i++;\n            while (coords[2 * j + axis] > t) j--;\n        }\n\n        if (coords[2 * left + axis] === t) swapItem(ids, coords, left, j);\n        else {\n            j++;\n            swapItem(ids, coords, j, right);\n        }\n\n        if (j <= k) left = j + 1;\n        if (k <= j) right = j - 1;\n    }\n}\n\n/**\n * @param {Uint16Array | Uint32Array} ids\n * @param {TypedArray} coords\n * @param {number} i\n * @param {number} j\n */\nfunction swapItem(ids, coords, i, j) {\n    swap(ids, i, j);\n    swap(coords, 2 * i, 2 * j);\n    swap(coords, 2 * i + 1, 2 * j + 1);\n}\n\n/**\n * @param {TypedArray} arr\n * @param {number} i\n * @param {number} j\n */\nfunction swap(arr, i, j) {\n    const tmp = arr[i];\n    arr[i] = arr[j];\n    arr[j] = tmp;\n}\n\n/**\n * @param {number} ax\n * @param {number} ay\n * @param {number} bx\n * @param {number} by\n */\nfunction sqDist(ax, ay, bx, by) {\n    const dx = ax - bx;\n    const dy = ay - by;\n    return dx * dx + dy * dy;\n}\n"],"names":["ARRAY_TYPES","VERSION","HEADER_SIZE","STACK","KDBush","data","magic","versionAndType","version","ArrayType","nodeSize","numItems","ArrayBufferType","arrayTypeIndex","coordsByteSize","idsByteSize","padCoords","x","y","index","numAdded","sort","minX","minY","maxX","maxY","ids","coords","sp","result","axis","right","left","i","m","qx","qy","r","out","count","r2","sqDist","select","k","n","z","s","sd","newLeft","newRight","t","j","swapItem","swap","arr","tmp","ax","ay","bx","by","dx","dy"],"mappings":"aACA,MAAMA,EAAc,CAChB,UAAW,WAAY,kBAAmB,WAAY,YACtD,WAAY,YAAa,aAAc,YAC3C,EAKMC,EAAU,EACVC,EAAc,EAIdC,EAAQ,IAAI,YAAY,EAAE,EAEhC,MAAqBC,CAAO,CAMxB,OAAO,KAAKC,EAAM,CAEd,GAAI,CAACA,GAAQA,EAAK,aAAe,QAAaA,EAAK,OAC/C,MAAM,IAAI,MAAM,+DAA+D,EAEnF,KAAM,CAACC,EAAOC,CAAc,EAAI,IAAI,WAAWF,EAAM,EAAG,CAAC,EACzD,GAAIC,IAAU,IACV,MAAM,IAAI,MAAM,gDAAgD,EAEpE,MAAME,EAAUD,GAAkB,EAClC,GAAIC,IAAYP,EACZ,MAAM,IAAI,MAAM,QAAQO,CAAO,wBAAwBP,CAAO,GAAG,EAErE,MAAMQ,EAAYT,EAAYO,EAAiB,EAAI,EACnD,GAAI,CAACE,EACD,MAAM,IAAI,MAAM,0BAA0B,EAE9C,KAAM,CAACC,CAAQ,EAAI,IAAI,YAAYL,EAAM,EAAG,CAAC,EACvC,CAACM,CAAQ,EAAI,IAAI,YAAYN,EAAM,EAAG,CAAC,EAE7C,OAAO,IAAID,EAAOO,EAAUD,EAAUD,EAAW,OAAWJ,CAAI,CACpE,CAUA,YAAYM,EAAUD,EAAW,GAAID,EAAY,aAAcG,EAAkB,YAAaP,EAAM,CAChG,GAAI,MAAMM,CAAQ,GAAKA,EAAW,QAAS,IAAI,MAAM,8BAA8BA,CAAQ,GAAG,EAE9F,KAAK,SAAW,CAACA,EACjB,KAAK,SAAW,KAAK,IAAI,KAAK,IAAI,CAACD,EAAU,CAAC,EAAG,KAAK,EACtD,KAAK,UAAYD,EACjB,KAAK,eAAiBE,EAAW,MAAQ,YAAc,YAEvD,MAAME,EAAiBb,EAAY,QAAQ,KAAK,SAAS,EACnDc,EAAiBH,EAAW,EAAI,KAAK,UAAU,kBAC/CI,EAAcJ,EAAW,KAAK,eAAe,kBAC7CK,GAAa,EAAID,EAAc,GAAK,EAE1C,GAAIF,EAAiB,EACjB,MAAM,IAAI,MAAM,iCAAiCJ,CAAS,GAAG,EAGjE,GAAIJ,EACA,KAAK,KAAOA,EAEZ,KAAK,IAAM,IAAI,KAAK,eAAeA,EAAMH,EAAaS,CAAQ,EAE9D,KAAK,OAAS,IAAIF,EAAUJ,EAAMH,EAAca,EAAcC,EAAWL,EAAW,CAAC,EACrF,KAAK,KAAOA,EAAW,EACvB,KAAK,UAAY,OAEd,CACH,MAAMN,EAAO,KAAK,KAAO,IAAIO,EAAgBV,EAAcY,EAAiBC,EAAcC,CAAS,EAEnG,KAAK,IAAM,IAAI,KAAK,eAAeX,EAAMH,EAAaS,CAAQ,EAE9D,KAAK,OAAS,IAAIF,EAAUJ,EAAMH,EAAca,EAAcC,EAAWL,EAAW,CAAC,EACrF,KAAK,KAAO,EACZ,KAAK,UAAY,GAGjB,IAAI,WAAWN,EAAM,EAAG,CAAC,EAAE,IAAI,CAAC,KAAOJ,GAAW,GAAKY,CAAc,CAAC,EACtE,IAAI,YAAYR,EAAM,EAAG,CAAC,EAAE,CAAC,EAAIK,EACjC,IAAI,YAAYL,EAAM,EAAG,CAAC,EAAE,CAAC,EAAIM,CACrC,CACJ,CAQA,IAAIM,EAAGC,EAAG,CACN,MAAMC,EAAQ,KAAK,MAAQ,EAC3B,YAAK,IAAIA,CAAK,EAAIA,EAClB,KAAK,OAAO,KAAK,MAAM,EAAIF,EAC3B,KAAK,OAAO,KAAK,MAAM,EAAIC,EACpBC,CACX,CAKA,QAAS,CACL,MAAMC,EAAW,KAAK,MAAQ,EAC9B,GAAIA,IAAa,KAAK,SAClB,MAAM,IAAI,MAAM,SAASA,CAAQ,wBAAwB,KAAK,QAAQ,GAAG,EAG7E,OAAAC,EAAK,KAAK,IAAK,KAAK,OAAQ,KAAK,SAAU,EAAG,KAAK,SAAW,EAAG,CAAC,EAElE,KAAK,UAAY,GACV,IACX,CAUA,MAAMC,EAAMC,EAAMC,EAAMC,EAAM,CAC1B,GAAI,CAAC,KAAK,UAAW,MAAM,IAAI,MAAM,6CAA6C,EAElF,KAAM,CAAC,IAAAC,EAAK,OAAAC,EAAQ,SAAAjB,CAAA,EAAY,KAChCP,EAAM,CAAC,EAAI,EACXA,EAAM,CAAC,EAAIuB,EAAI,OAAS,EACxBvB,EAAM,CAAC,EAAI,EACX,IAAIyB,EAAK,EACT,MAAMC,EAAS,CAAA,EAGf,KAAOD,EAAK,GAAG,CACX,MAAME,EAAO3B,EAAM,EAAEyB,CAAE,EACjBG,EAAQ5B,EAAM,EAAEyB,CAAE,EAClBI,EAAO7B,EAAM,EAAEyB,CAAE,EAGvB,GAAIG,EAAQC,GAAQtB,EAAU,CAC1B,QAASuB,EAAID,EAAMC,GAAKF,EAAOE,IAAK,CAChC,MAAMhB,EAAIU,EAAO,EAAIM,CAAC,EAChBf,EAAIS,EAAO,EAAIM,EAAI,CAAC,EACtBhB,GAAKK,GAAQL,GAAKO,GAAQN,GAAKK,GAAQL,GAAKO,GAAMI,EAAO,KAAKH,EAAIO,CAAC,CAAC,CAC5E,CACA,QACJ,CAGA,MAAMC,EAAKF,EAAOD,GAAU,EAGtBd,EAAIU,EAAO,EAAIO,CAAC,EAChBhB,EAAIS,EAAO,EAAIO,EAAI,CAAC,EACtBjB,GAAKK,GAAQL,GAAKO,GAAQN,GAAKK,GAAQL,GAAKO,GAAMI,EAAO,KAAKH,EAAIQ,CAAC,CAAC,GAGpEJ,IAAS,EAAIR,GAAQL,EAAIM,GAAQL,KACjCf,EAAMyB,GAAI,EAAII,EACd7B,EAAMyB,GAAI,EAAIM,EAAI,EAClB/B,EAAMyB,GAAI,EAAI,EAAIE,IAElBA,IAAS,EAAIN,GAAQP,EAAIQ,GAAQP,KACjCf,EAAMyB,GAAI,EAAIM,EAAI,EAClB/B,EAAMyB,GAAI,EAAIG,EACd5B,EAAMyB,GAAI,EAAI,EAAIE,EAE1B,CAEA,OAAOD,CACX,CASA,OAAOM,EAAIC,EAAIC,EAAG,CACd,MAAMR,EAAkC,CAAA,EACxC,YAAK,WAAWM,EAAIC,EAAIC,EAAGR,CAAM,EAC1BA,CACX,CAaA,WAAWM,EAAIC,EAAIC,EAAGC,EAAK,CACvB,GAAI,CAAC,KAAK,UAAW,MAAM,IAAI,MAAM,6CAA6C,EAElF,KAAM,CAAC,IAAAZ,EAAK,OAAAC,EAAQ,SAAAjB,CAAA,EAAY,KAChCP,EAAM,CAAC,EAAI,EACXA,EAAM,CAAC,EAAIuB,EAAI,OAAS,EACxBvB,EAAM,CAAC,EAAI,EACX,IAAIyB,EAAK,EACLW,EAAQ,EACZ,MAAMC,EAAKH,EAAIA,EAGf,KAAOT,EAAK,GAAG,CACX,MAAME,EAAO3B,EAAM,EAAEyB,CAAE,EACjBG,EAAQ5B,EAAM,EAAEyB,CAAE,EAClBI,EAAO7B,EAAM,EAAEyB,CAAE,EAGvB,GAAIG,EAAQC,GAAQtB,EAAU,CAC1B,QAASuB,EAAID,EAAMC,GAAKF,EAAOE,IACvBQ,EAAOd,EAAO,EAAIM,CAAC,EAAGN,EAAO,EAAIM,EAAI,CAAC,EAAGE,EAAIC,CAAE,GAAKI,MAAQD,GAAO,EAAIb,EAAIO,CAAC,GAEpF,QACJ,CAGA,MAAMC,EAAKF,EAAOD,GAAU,EAGtBd,EAAIU,EAAO,EAAIO,CAAC,EAChBhB,EAAIS,EAAO,EAAIO,EAAI,CAAC,EACtBO,EAAOxB,EAAGC,EAAGiB,EAAIC,CAAE,GAAKI,IAAIF,EAAIC,GAAO,EAAIb,EAAIQ,CAAC,IAGhDJ,IAAS,EAAIK,EAAKE,GAAKpB,EAAImB,EAAKC,GAAKnB,KACrCf,EAAMyB,GAAI,EAAII,EACd7B,EAAMyB,GAAI,EAAIM,EAAI,EAClB/B,EAAMyB,GAAI,EAAI,EAAIE,IAElBA,IAAS,EAAIK,EAAKE,GAAKpB,EAAImB,EAAKC,GAAKnB,KACrCf,EAAMyB,GAAI,EAAIM,EAAI,EAClB/B,EAAMyB,GAAI,EAAIG,EACd5B,EAAMyB,GAAI,EAAI,EAAIE,EAE1B,CAEA,OAAOS,CACX,CACJ,CAUA,SAASlB,EAAKK,EAAKC,EAAQjB,EAAUsB,EAAMD,EAAOD,EAAM,CACpD,GAAIC,EAAQC,GAAQtB,EAAU,OAE9B,MAAMwB,EAAKF,EAAOD,GAAU,EAI5BW,EAAOhB,EAAKC,EAAQO,EAAGF,EAAMD,EAAOD,CAAI,EAGxCT,EAAKK,EAAKC,EAAQjB,EAAUsB,EAAME,EAAI,EAAG,EAAIJ,CAAI,EACjDT,EAAKK,EAAKC,EAAQjB,EAAUwB,EAAI,EAAGH,EAAO,EAAID,CAAI,CACtD,CAYA,SAASY,EAAOhB,EAAKC,EAAQgB,EAAGX,EAAMD,EAAOD,EAAM,CAE/C,KAAOC,EAAQC,GAAM,CACjB,GAAID,EAAQC,EAAO,IAAK,CACpB,MAAMY,EAAIb,EAAQC,EAAO,EACnBE,EAAIS,EAAIX,EAAO,EACfa,EAAI,KAAK,IAAID,CAAC,EACdE,EAAI,GAAM,KAAK,IAAI,EAAID,EAAI,CAAC,EAC5BE,EAAK,GAAM,KAAK,KAAKF,EAAIC,GAAKF,EAAIE,GAAKF,CAAC,GAAKV,EAAIU,EAAI,EAAI,EAAI,GAAK,GAClEI,EAAU,KAAK,IAAIhB,EAAM,KAAK,MAAMW,EAAIT,EAAIY,EAAIF,EAAIG,CAAE,CAAC,EACvDE,EAAW,KAAK,IAAIlB,EAAO,KAAK,MAAMY,GAAKC,EAAIV,GAAKY,EAAIF,EAAIG,CAAE,CAAC,EACrEL,EAAOhB,EAAKC,EAAQgB,EAAGK,EAASC,EAAUnB,CAAI,CAClD,CAEA,MAAMoB,EAAIvB,EAAO,EAAIgB,EAAIb,CAAI,EAC7B,IAAIG,EAAID,EACJmB,EAAIpB,EAKR,IAHAqB,EAAS1B,EAAKC,EAAQK,EAAMW,CAAC,EACzBhB,EAAO,EAAII,EAAQD,CAAI,EAAIoB,GAAGE,EAAS1B,EAAKC,EAAQK,EAAMD,CAAK,EAE5DE,EAAIkB,GAAG,CAIV,IAHAC,EAAS1B,EAAKC,EAAQM,EAAGkB,CAAC,EAC1BlB,IACAkB,IACOxB,EAAO,EAAIM,EAAIH,CAAI,EAAIoB,GAAGjB,IACjC,KAAON,EAAO,EAAIwB,EAAIrB,CAAI,EAAIoB,GAAGC,GACrC,CAEIxB,EAAO,EAAIK,EAAOF,CAAI,IAAMoB,EAAGE,EAAS1B,EAAKC,EAAQK,EAAMmB,CAAC,GAE5DA,IACAC,EAAS1B,EAAKC,EAAQwB,EAAGpB,CAAK,GAG9BoB,GAAKR,IAAGX,EAAOmB,EAAI,GACnBR,GAAKQ,IAAGpB,EAAQoB,EAAI,EAC5B,CACJ,CAQA,SAASC,EAAS1B,EAAKC,EAAQM,EAAGkB,EAAG,CACjCE,EAAK3B,EAAKO,EAAGkB,CAAC,EACdE,EAAK1B,EAAQ,EAAIM,EAAG,EAAIkB,CAAC,EACzBE,EAAK1B,EAAQ,EAAIM,EAAI,EAAG,EAAIkB,EAAI,CAAC,CACrC,CAOA,SAASE,EAAKC,EAAKrB,EAAGkB,EAAG,CACrB,MAAMI,EAAMD,EAAIrB,CAAC,EACjBqB,EAAIrB,CAAC,EAAIqB,EAAIH,CAAC,EACdG,EAAIH,CAAC,EAAII,CACb,CAQA,SAASd,EAAOe,EAAIC,EAAIC,EAAIC,EAAI,CAC5B,MAAMC,EAAKJ,EAAKE,EACVG,EAAKJ,EAAKE,EAChB,OAAOC,EAAKA,EAAKC,EAAKA,CAC1B","x_google_ignoreList":[0]}

package/dist/index.cjs559.js

@@ -0,0 +1,2 @@
+"use strict";Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const n=require("./index.cjs611.js"),K=(3+16*n.epsilon)*n.epsilon,L=(2+12*n.epsilon)*n.epsilon,N=(9+64*n.epsilon)*n.epsilon*n.epsilon,T=n.vec(4),E=n.vec(8),F=n.vec(12),G=n.vec(16),f=n.vec(4);function Q(q,A,O,P,D,S,k){let v,h,m,M,t,e,r,o,i,l,s,c,u,p,a,d,C,b;const g=q-D,w=O-D,_=A-S,j=P-S;p=g*j,e=n.splitter*g,r=e-(e-g),o=g-r,e=n.splitter*j,i=e-(e-j),l=j-i,a=o*l-(p-r*i-o*i-r*l),d=_*w,e=n.splitter*_,r=e-(e-_),o=_-r,e=n.splitter*w,i=e-(e-w),l=w-i,C=o*l-(d-r*i-o*i-r*l),s=a-C,t=a-s,T[0]=a-(s+t)+(t-C),c=p+s,t=c-p,u=p-(c-t)+(s-t),s=u-d,t=u-s,T[1]=u-(s+t)+(t-d),b=c+s,t=b-c,T[2]=c-(b-t)+(s-t),T[3]=b;let B=n.estimate(4,T),z=L*k;if(B>=z||-B>=z||(t=q-g,v=q-(g+t)+(t-D),t=O-w,m=O-(w+t)+(t-D),t=A-_,h=A-(_+t)+(t-S),t=P-j,M=P-(j+t)+(t-S),v===0&&h===0&&m===0&&M===0)||(z=N*k+n.resulterrbound*Math.abs(B),B+=g*M+j*v-(_*m+w*h),B>=z||-B>=z))return B;p=v*j,e=n.splitter*v,r=e-(e-v),o=v-r,e=n.splitter*j,i=e-(e-j),l=j-i,a=o*l-(p-r*i-o*i-r*l),d=h*w,e=n.splitter*h,r=e-(e-h),o=h-r,e=n.splitter*w,i=e-(e-w),l=w-i,C=o*l-(d-r*i-o*i-r*l),s=a-C,t=a-s,f[0]=a-(s+t)+(t-C),c=p+s,t=c-p,u=p-(c-t)+(s-t),s=u-d,t=u-s,f[1]=u-(s+t)+(t-d),b=c+s,t=b-c,f[2]=c-(b-t)+(s-t),f[3]=b;const H=n.sum(4,T,4,f,E);p=g*M,e=n.splitter*g,r=e-(e-g),o=g-r,e=n.splitter*M,i=e-(e-M),l=M-i,a=o*l-(p-r*i-o*i-r*l),d=_*m,e=n.splitter*_,r=e-(e-_),o=_-r,e=n.splitter*m,i=e-(e-m),l=m-i,C=o*l-(d-r*i-o*i-r*l),s=a-C,t=a-s,f[0]=a-(s+t)+(t-C),c=p+s,t=c-p,u=p-(c-t)+(s-t),s=u-d,t=u-s,f[1]=u-(s+t)+(t-d),b=c+s,t=b-c,f[2]=c-(b-t)+(s-t),f[3]=b;const I=n.sum(H,E,4,f,F);p=v*M,e=n.splitter*v,r=e-(e-v),o=v-r,e=n.splitter*M,i=e-(e-M),l=M-i,a=o*l-(p-r*i-o*i-r*l),d=h*m,e=n.splitter*h,r=e-(e-h),o=h-r,e=n.splitter*m,i=e-(e-m),l=m-i,C=o*l-(d-r*i-o*i-r*l),s=a-C,t=a-s,f[0]=a-(s+t)+(t-C),c=p+s,t=c-p,u=p-(c-t)+(s-t),s=u-d,t=u-s,f[1]=u-(s+t)+(t-d),b=c+s,t=b-c,f[2]=c-(b-t)+(s-t),f[3]=b;const J=n.sum(I,F,4,f,G);return G[J-1]}function R(q,A,O,P,D,S){const k=(A-S)*(O-D),v=(q-D)*(P-S),h=k-v,m=Math.abs(k+v);return Math.abs(h)>=K*m?h:-Q(q,A,O,P,D,S,m)}exports.orient2d=R;
+//# sourceMappingURL=index.cjs559.js.map

package/dist/index.cjs559.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs559.js","sources":["../../../node_modules/@kninnug/constrainautor/node_modules/robust-predicates/esm/orient2d.js"],"sourcesContent":["import {epsilon, splitter, resulterrbound, estimate, vec, sum} from './util.js';\n\nconst ccwerrboundA = (3 + 16 * epsilon) * epsilon;\nconst ccwerrboundB = (2 + 12 * epsilon) * epsilon;\nconst ccwerrboundC = (9 + 64 * epsilon) * epsilon * epsilon;\n\nconst B = vec(4);\nconst C1 = vec(8);\nconst C2 = vec(12);\nconst D = vec(16);\nconst u = vec(4);\n\nfunction orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {\n    let acxtail, acytail, bcxtail, bcytail;\n    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;\n\n    const acx = ax - cx;\n    const bcx = bx - cx;\n    const acy = ay - cy;\n    const bcy = by - cy;\n\n    s1 = acx * bcy;\n    c = splitter * acx;\n    ahi = c - (c - acx);\n    alo = acx - ahi;\n    c = splitter * bcy;\n    bhi = c - (c - bcy);\n    blo = bcy - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = acy * bcx;\n    c = splitter * acy;\n    ahi = c - (c - acy);\n    alo = acy - ahi;\n    c = splitter * bcx;\n    bhi = c - (c - bcx);\n    blo = bcx - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    B[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    B[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    B[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    B[3] = u3;\n\n    let det = estimate(4, B);\n    let errbound = ccwerrboundB * detsum;\n    if (det >= errbound || -det >= errbound) {\n        return det;\n    }\n\n    bvirt = ax - acx;\n    acxtail = ax - (acx + bvirt) + (bvirt - cx);\n    bvirt = bx - bcx;\n    bcxtail = bx - (bcx + bvirt) + (bvirt - cx);\n    bvirt = ay - acy;\n    acytail = ay - (acy + bvirt) + (bvirt - cy);\n    bvirt = by - bcy;\n    bcytail = by - (bcy + bvirt) + (bvirt - cy);\n\n    if (acxtail === 0 && acytail === 0 && bcxtail === 0 && bcytail === 0) {\n        return det;\n    }\n\n    errbound = ccwerrboundC * detsum + resulterrbound * Math.abs(det);\n    det += (acx * bcytail + bcy * acxtail) - (acy * bcxtail + bcx * acytail);\n    if (det >= errbound || -det >= errbound) return det;\n\n    s1 = acxtail * bcy;\n    c = splitter * acxtail;\n    ahi = c - (c - acxtail);\n    alo = acxtail - ahi;\n    c = splitter * bcy;\n    bhi = c - (c - bcy);\n    blo = bcy - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = acytail * bcx;\n    c = splitter * acytail;\n    ahi = c - (c - acytail);\n    alo = acytail - ahi;\n    c = splitter * bcx;\n    bhi = c - (c - bcx);\n    blo = bcx - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    u[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    u[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    u[3] = u3;\n    const C1len = sum(4, B, 4, u, C1);\n\n    s1 = acx * bcytail;\n    c = splitter * acx;\n    ahi = c - (c - acx);\n    alo = acx - ahi;\n    c = splitter * bcytail;\n    bhi = c - (c - bcytail);\n    blo = bcytail - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = acy * bcxtail;\n    c = splitter * acy;\n    ahi = c - (c - acy);\n    alo = acy - ahi;\n    c = splitter * bcxtail;\n    bhi = c - (c - bcxtail);\n    blo = bcxtail - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    u[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    u[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    u[3] = u3;\n    const C2len = sum(C1len, C1, 4, u, C2);\n\n    s1 = acxtail * bcytail;\n    c = splitter * acxtail;\n    ahi = c - (c - acxtail);\n    alo = acxtail - ahi;\n    c = splitter * bcytail;\n    bhi = c - (c - bcytail);\n    blo = bcytail - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = acytail * bcxtail;\n    c = splitter * acytail;\n    ahi = c - (c - acytail);\n    alo = acytail - ahi;\n    c = splitter * bcxtail;\n    bhi = c - (c - bcxtail);\n    blo = bcxtail - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    u[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    u[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    u[3] = u3;\n    const Dlen = sum(C2len, C2, 4, u, D);\n\n    return D[Dlen - 1];\n}\n\nexport function orient2d(ax, ay, bx, by, cx, cy) {\n    const detleft = (ay - cy) * (bx - cx);\n    const detright = (ax - cx) * (by - cy);\n    const det = detleft - detright;\n\n    const detsum = Math.abs(detleft + detright);\n    if (Math.abs(det) >= ccwerrboundA * detsum) return det;\n\n    return -orient2dadapt(ax, ay, bx, by, cx, cy, detsum);\n}\n\nexport function orient2dfast(ax, ay, bx, by, cx, cy) {\n    return (ay - cy) * (bx - cx) - (ax - cx) * (by - cy);\n}\n"],"names":["ccwerrboundA","epsilon","ccwerrboundB","ccwerrboundC","B","vec","C1","C2","D","u","orient2dadapt","ax","ay","bx","by","cx","cy","detsum","acxtail","acytail","bcxtail","bcytail","bvirt","c","ahi","alo","bhi","blo","_i","_j","_0","s1","s0","t1","t0","u3","acx","bcx","acy","bcy","splitter","det","estimate","errbound","resulterrbound","C1len","sum","C2len","Dlen","orient2d","detleft","detright"],"mappings":"qHAEMA,GAAgB,EAAI,GAAKC,EAAAA,SAAWA,EAAAA,QACpCC,GAAgB,EAAI,GAAKD,EAAAA,SAAWA,EAAAA,QACpCE,GAAgB,EAAI,GAAKF,EAAAA,SAAWA,EAAAA,QAAUA,EAAAA,QAE9CG,EAAIC,EAAAA,IAAI,CAAC,EACTC,EAAKD,EAAAA,IAAI,CAAC,EACVE,EAAKF,EAAAA,IAAI,EAAE,EACXG,EAAIH,EAAAA,IAAI,EAAE,EACVI,EAAIJ,EAAAA,IAAI,CAAC,EAEf,SAASK,EAAcC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAQ,CACnD,IAAIC,EAASC,EAASC,EAASC,EAC3BC,EAAOC,EAAGC,EAAKC,EAAKC,EAAKC,EAAKC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAE9D,MAAMC,EAAMzB,EAAKI,EACXsB,EAAMxB,EAAKE,EACXuB,EAAM1B,EAAKI,EACXuB,EAAMzB,EAAKE,EAEjBe,EAAKK,EAAMG,EACXhB,EAAIiB,EAAAA,SAAWJ,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAIiB,EAAAA,SAAWD,EACfb,EAAMH,GAAKA,EAAIgB,GACfZ,EAAMY,EAAMb,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKK,EAAMD,EACXd,EAAIiB,EAAAA,SAAWF,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAIiB,EAAAA,SAAWH,EACfX,EAAMH,GAAKA,EAAIc,GACfV,EAAMU,EAAMX,EACZQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbxB,EAAE,CAAC,EAAI4B,GAAMJ,EAAKN,IAAUA,EAAQY,GACpCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbxB,EAAE,CAAC,EAAI0B,GAAMF,EAAKN,IAAUA,EAAQW,GACpCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbzB,EAAE,CAAC,EAAIyB,GAAMM,EAAKb,IAAUM,EAAKN,GACjClB,EAAE,CAAC,EAAI+B,EAEP,IAAIM,EAAMC,EAAAA,SAAS,EAAGtC,CAAC,EACnBuC,EAAWzC,EAAee,EAoB9B,GAnBIwB,GAAOE,GAAY,CAACF,GAAOE,IAI/BrB,EAAQX,EAAKyB,EACblB,EAAUP,GAAMyB,EAAMd,IAAUA,EAAQP,GACxCO,EAAQT,EAAKwB,EACbjB,EAAUP,GAAMwB,EAAMf,IAAUA,EAAQP,GACxCO,EAAQV,EAAK0B,EACbnB,EAAUP,GAAM0B,EAAMhB,IAAUA,EAAQN,GACxCM,EAAQR,EAAKyB,EACblB,EAAUP,GAAMyB,EAAMjB,IAAUA,EAAQN,GAEpCE,IAAY,GAAKC,IAAY,GAAKC,IAAY,GAAKC,IAAY,KAInEsB,EAAWxC,EAAec,EAAS2B,EAAAA,eAAiB,KAAK,IAAIH,CAAG,EAChEA,GAAQL,EAAMf,EAAUkB,EAAMrB,GAAYoB,EAAMlB,EAAUiB,EAAMlB,GAC5DsB,GAAOE,GAAY,CAACF,GAAOE,GAAU,OAAOF,EAEhDV,EAAKb,EAAUqB,EACfhB,EAAIiB,EAAAA,SAAWtB,EACfM,EAAMD,GAAKA,EAAIL,GACfO,EAAMP,EAAUM,EAChBD,EAAIiB,EAAAA,SAAWD,EACfb,EAAMH,GAAKA,EAAIgB,GACfZ,EAAMY,EAAMb,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKd,EAAUkB,EACfd,EAAIiB,EAAAA,SAAWrB,EACfK,EAAMD,GAAKA,EAAIJ,GACfM,EAAMN,EAAUK,EAChBD,EAAIiB,EAAAA,SAAWH,EACfX,EAAMH,GAAKA,EAAIc,GACfV,EAAMU,EAAMX,EACZQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbnB,EAAE,CAAC,EAAIuB,GAAMJ,EAAKN,IAAUA,EAAQY,GACpCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbnB,EAAE,CAAC,EAAIqB,GAAMF,EAAKN,IAAUA,EAAQW,GACpCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbpB,EAAE,CAAC,EAAIoB,GAAMM,EAAKb,IAAUM,EAAKN,GACjCb,EAAE,CAAC,EAAI0B,EACP,MAAMU,EAAQC,EAAAA,IAAI,EAAG1C,EAAG,EAAGK,EAAGH,CAAE,EAEhCyB,EAAKK,EAAMf,EACXE,EAAIiB,EAAAA,SAAWJ,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAIiB,EAAAA,SAAWnB,EACfK,EAAMH,GAAKA,EAAIF,GACfM,EAAMN,EAAUK,EAChBM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKK,EAAMlB,EACXG,EAAIiB,EAAAA,SAAWF,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAIiB,EAAAA,SAAWpB,EACfM,EAAMH,GAAKA,EAAIH,GACfO,EAAMP,EAAUM,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbnB,EAAE,CAAC,EAAIuB,GAAMJ,EAAKN,IAAUA,EAAQY,GACpCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbnB,EAAE,CAAC,EAAIqB,GAAMF,EAAKN,IAAUA,EAAQW,GACpCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbpB,EAAE,CAAC,EAAIoB,GAAMM,EAAKb,IAAUM,EAAKN,GACjCb,EAAE,CAAC,EAAI0B,EACP,MAAMY,EAAQD,EAAAA,IAAID,EAAOvC,EAAI,EAAGG,EAAGF,CAAE,EAErCwB,EAAKb,EAAUG,EACfE,EAAIiB,EAAAA,SAAWtB,EACfM,EAAMD,GAAKA,EAAIL,GACfO,EAAMP,EAAUM,EAChBD,EAAIiB,EAAAA,SAAWnB,EACfK,EAAMH,GAAKA,EAAIF,GACfM,EAAMN,EAAUK,EAChBM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKd,EAAUC,EACfG,EAAIiB,EAAAA,SAAWrB,EACfK,EAAMD,GAAKA,EAAIJ,GACfM,EAAMN,EAAUK,EAChBD,EAAIiB,EAAAA,SAAWpB,EACfM,EAAMH,GAAKA,EAAIH,GACfO,EAAMP,EAAUM,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbnB,EAAE,CAAC,EAAIuB,GAAMJ,EAAKN,IAAUA,EAAQY,GACpCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbnB,EAAE,CAAC,EAAIqB,GAAMF,EAAKN,IAAUA,EAAQW,GACpCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbpB,EAAE,CAAC,EAAIoB,GAAMM,EAAKb,IAAUM,EAAKN,GACjCb,EAAE,CAAC,EAAI0B,EACP,MAAMa,EAAOF,EAAAA,IAAIC,EAAOxC,EAAI,EAAGE,EAAGD,CAAC,EAEnC,OAAOA,EAAEwC,EAAO,CAAC,CACrB,CAEO,SAASC,EAAStC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAI,CAC7C,MAAMkC,GAAWtC,EAAKI,IAAOH,EAAKE,GAC5BoC,GAAYxC,EAAKI,IAAOD,EAAKE,GAC7ByB,EAAMS,EAAUC,EAEhBlC,EAAS,KAAK,IAAIiC,EAAUC,CAAQ,EAC1C,OAAI,KAAK,IAAIV,CAAG,GAAKzC,EAAeiB,EAAewB,EAE5C,CAAC/B,EAAcC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,CAAM,CACxD","x_google_ignoreList":[0]}

package/dist/index.cjs56.js

@@ -1,2 +1,2 @@
-"use strict";const t=require("./index.cjs94.js"),c=require("./index.cjs38.js"),l=require("./index.cjs161.js"),f=require("./index.cjs229.js"),q=require("./index.cjs228.js"),v=t(function(s,e,o){var i,u;return c(e)?u=e:(e=q(e),i=e.slice(0,-1),e=e[e.length-1]),l(s,function(r){var n=u;if(!n){if(i&&i.length&&(r=f(r,i)),r==null)return;n=r[e]}return n==null?n:n.apply(r,o)})});module.exports=v;
+"use strict";const t=require("./index.cjs94.js"),c=require("./index.cjs38.js"),l=require("./index.cjs161.js"),f=require("./index.cjs246.js"),q=require("./index.cjs245.js"),v=t(function(s,e,o){var i,u;return c(e)?u=e:(e=q(e),i=e.slice(0,-1),e=e[e.length-1]),l(s,function(r){var n=u;if(!n){if(i&&i.length&&(r=f(r,i)),r==null)return;n=r[e]}return n==null?n:n.apply(r,o)})});module.exports=v;
 //# sourceMappingURL=index.cjs56.js.map

package/dist/index.cjs560.js

@@ -0,0 +1,2 @@
+"use strict";Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const e=require("./index.cjs611.js"),ke=(10+96*e.epsilon)*e.epsilon,we=(4+48*e.epsilon)*e.epsilon,ze=(44+576*e.epsilon)*e.epsilon*e.epsilon,$=e.vec(4),ee=e.vec(4),se=e.vec(4),Q=e.vec(4),R=e.vec(4),U=e.vec(4),F=e.vec(4),G=e.vec(4),qe=e.vec(8),Ae=e.vec(8),Be=e.vec(8),Ce=e.vec(8),Oe=e.vec(8),Pe=e.vec(8),ge=e.vec(8),je=e.vec(8),Se=e.vec(8),ie=e.vec(4),re=e.vec(4),oe=e.vec(4),p=e.vec(8),d=e.vec(16),T=e.vec(16),k=e.vec(16),P=e.vec(32),te=e.vec(32),w=e.vec(48),H=e.vec(64);let he=e.vec(1152),Te=e.vec(1152);function z(K,W,X){K=e.sum(K,he,W,X,Te);const Y=he;return he=Te,Te=Y,K}function De(K,W,X,Y,pe,ue,I,J,ce){let m,h,M,y,g,x,j,be,fe,ve,me,de,_e,le,Z,ae,ne,ye,xe,L,N,s,t,l,a,n,o,c,i,u,r,b,f,_,v;const S=K-I,q=X-I,A=pe-I,C=W-J,B=Y-J,O=ue-J;r=q*O,t=e.splitter*q,l=t-(t-q),a=q-l,t=e.splitter*O,n=t-(t-O),o=O-n,b=a*o-(r-l*n-a*n-l*o),f=A*B,t=e.splitter*A,l=t-(t-A),a=A-l,t=e.splitter*B,n=t-(t-B),o=B-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,$[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,$[1]=u-(c+s)+(s-f),v=i+c,s=v-i,$[2]=i-(v-s)+(c-s),$[3]=v,r=A*C,t=e.splitter*A,l=t-(t-A),a=A-l,t=e.splitter*C,n=t-(t-C),o=C-n,b=a*o-(r-l*n-a*n-l*o),f=S*O,t=e.splitter*S,l=t-(t-S),a=S-l,t=e.splitter*O,n=t-(t-O),o=O-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,ee[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,ee[1]=u-(c+s)+(s-f),v=i+c,s=v-i,ee[2]=i-(v-s)+(c-s),ee[3]=v,r=S*B,t=e.splitter*S,l=t-(t-S),a=S-l,t=e.splitter*B,n=t-(t-B),o=B-n,b=a*o-(r-l*n-a*n-l*o),f=q*C,t=e.splitter*q,l=t-(t-q),a=q-l,t=e.splitter*C,n=t-(t-C),o=C-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,se[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,se[1]=u-(c+s)+(s-f),v=i+c,s=v-i,se[2]=i-(v-s)+(c-s),se[3]=v,m=e.sum(e.sum(e.sum(e.scale(e.scale(4,$,S,p),p,S,d),d,e.scale(e.scale(4,$,C,p),p,C,T),T,P),P,e.sum(e.scale(e.scale(4,ee,q,p),p,q,d),d,e.scale(e.scale(4,ee,B,p),p,B,T),T,te),te,H),H,e.sum(e.scale(e.scale(4,se,A,p),p,A,d),d,e.scale(e.scale(4,se,O,p),p,O,T),T,P),P,he);let V=e.estimate(m,he),Me=we*ce;if(V>=Me||-V>=Me||(s=K-S,h=K-(S+s)+(s-I),s=W-C,g=W-(C+s)+(s-J),s=X-q,M=X-(q+s)+(s-I),s=Y-B,x=Y-(B+s)+(s-J),s=pe-A,y=pe-(A+s)+(s-I),s=ue-O,j=ue-(O+s)+(s-J),h===0&&M===0&&y===0&&g===0&&x===0&&j===0)||(Me=ze*ce+e.resulterrbound*Math.abs(V),V+=(S*S+C*C)*(q*j+O*M-(B*y+A*x))+2*(S*h+C*g)*(q*O-B*A)+((q*q+B*B)*(A*g+C*y-(O*h+S*j))+2*(q*M+B*x)*(A*C-O*S))+((A*A+O*O)*(S*x+B*h-(C*M+q*g))+2*(A*y+O*j)*(S*B-C*q)),V>=Me||-V>=Me))return V;if((M!==0||x!==0||y!==0||j!==0)&&(r=S*S,t=e.splitter*S,l=t-(t-S),a=S-l,b=a*a-(r-l*l-(l+l)*a),f=C*C,t=e.splitter*C,l=t-(t-C),a=C-l,_=a*a-(f-l*l-(l+l)*a),c=b+_,s=c-b,Q[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,Q[1]=u-(c-s)+(f-s),v=i+c,s=v-i,Q[2]=i-(v-s)+(c-s),Q[3]=v),(y!==0||j!==0||h!==0||g!==0)&&(r=q*q,t=e.splitter*q,l=t-(t-q),a=q-l,b=a*a-(r-l*l-(l+l)*a),f=B*B,t=e.splitter*B,l=t-(t-B),a=B-l,_=a*a-(f-l*l-(l+l)*a),c=b+_,s=c-b,R[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,R[1]=u-(c-s)+(f-s),v=i+c,s=v-i,R[2]=i-(v-s)+(c-s),R[3]=v),(h!==0||g!==0||M!==0||x!==0)&&(r=A*A,t=e.splitter*A,l=t-(t-A),a=A-l,b=a*a-(r-l*l-(l+l)*a),f=O*O,t=e.splitter*O,l=t-(t-O),a=O-l,_=a*a-(f-l*l-(l+l)*a),c=b+_,s=c-b,U[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,U[1]=u-(c-s)+(f-s),v=i+c,s=v-i,U[2]=i-(v-s)+(c-s),U[3]=v),h!==0&&(be=e.scale(4,$,h,qe),m=z(m,e.sum_three(e.scale(be,qe,2*S,d),d,e.scale(e.scale(4,U,h,p),p,B,T),T,e.scale(e.scale(4,R,h,p),p,-O,k),k,P,w),w)),g!==0&&(fe=e.scale(4,$,g,Ae),m=z(m,e.sum_three(e.scale(fe,Ae,2*C,d),d,e.scale(e.scale(4,R,g,p),p,A,T),T,e.scale(e.scale(4,U,g,p),p,-q,k),k,P,w),w)),M!==0&&(ve=e.scale(4,ee,M,Be),m=z(m,e.sum_three(e.scale(ve,Be,2*q,d),d,e.scale(e.scale(4,Q,M,p),p,O,T),T,e.scale(e.scale(4,U,M,p),p,-C,k),k,P,w),w)),x!==0&&(me=e.scale(4,ee,x,Ce),m=z(m,e.sum_three(e.scale(me,Ce,2*B,d),d,e.scale(e.scale(4,U,x,p),p,S,T),T,e.scale(e.scale(4,Q,x,p),p,-A,k),k,P,w),w)),y!==0&&(de=e.scale(4,se,y,Oe),m=z(m,e.sum_three(e.scale(de,Oe,2*A,d),d,e.scale(e.scale(4,R,y,p),p,C,T),T,e.scale(e.scale(4,Q,y,p),p,-B,k),k,P,w),w)),j!==0&&(_e=e.scale(4,se,j,Pe),m=z(m,e.sum_three(e.scale(_e,Pe,2*O,d),d,e.scale(e.scale(4,Q,j,p),p,q,T),T,e.scale(e.scale(4,R,j,p),p,-S,k),k,P,w),w)),h!==0||g!==0){if(M!==0||x!==0||y!==0||j!==0?(r=M*O,t=e.splitter*M,l=t-(t-M),a=M-l,t=e.splitter*O,n=t-(t-O),o=O-n,b=a*o-(r-l*n-a*n-l*o),f=q*j,t=e.splitter*q,l=t-(t-q),a=q-l,t=e.splitter*j,n=t-(t-j),o=j-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,F[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,F[1]=u-(c-s)+(f-s),v=i+c,s=v-i,F[2]=i-(v-s)+(c-s),F[3]=v,r=y*-B,t=e.splitter*y,l=t-(t-y),a=y-l,t=e.splitter*-B,n=t-(t- -B),o=-B-n,b=a*o-(r-l*n-a*n-l*o),f=A*-x,t=e.splitter*A,l=t-(t-A),a=A-l,t=e.splitter*-x,n=t-(t- -x),o=-x-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,G[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,G[1]=u-(c-s)+(f-s),v=i+c,s=v-i,G[2]=i-(v-s)+(c-s),G[3]=v,Z=e.sum(4,F,4,G,je),r=M*j,t=e.splitter*M,l=t-(t-M),a=M-l,t=e.splitter*j,n=t-(t-j),o=j-n,b=a*o-(r-l*n-a*n-l*o),f=y*x,t=e.splitter*y,l=t-(t-y),a=y-l,t=e.splitter*x,n=t-(t-x),o=x-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,re[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,re[1]=u-(c+s)+(s-f),v=i+c,s=v-i,re[2]=i-(v-s)+(c-s),re[3]=v,ye=4):(je[0]=0,Z=1,re[0]=0,ye=1),h!==0){const D=e.scale(Z,je,h,k);m=z(m,e.sum(e.scale(be,qe,h,d),d,e.scale(D,k,2*S,P),P,w),w);const E=e.scale(ye,re,h,p);m=z(m,e.sum_three(e.scale(E,p,2*S,d),d,e.scale(E,p,h,T),T,e.scale(D,k,h,P),P,te,H),H),x!==0&&(m=z(m,e.scale(e.scale(4,U,h,p),p,x,d),d)),j!==0&&(m=z(m,e.scale(e.scale(4,R,-h,p),p,j,d),d))}if(g!==0){const D=e.scale(Z,je,g,k);m=z(m,e.sum(e.scale(fe,Ae,g,d),d,e.scale(D,k,2*C,P),P,w),w);const E=e.scale(ye,re,g,p);m=z(m,e.sum_three(e.scale(E,p,2*C,d),d,e.scale(E,p,g,T),T,e.scale(D,k,g,P),P,te,H),H)}}if(M!==0||x!==0){if(y!==0||j!==0||h!==0||g!==0?(r=y*C,t=e.splitter*y,l=t-(t-y),a=y-l,t=e.splitter*C,n=t-(t-C),o=C-n,b=a*o-(r-l*n-a*n-l*o),f=A*g,t=e.splitter*A,l=t-(t-A),a=A-l,t=e.splitter*g,n=t-(t-g),o=g-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,F[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,F[1]=u-(c-s)+(f-s),v=i+c,s=v-i,F[2]=i-(v-s)+(c-s),F[3]=v,L=-O,N=-j,r=h*L,t=e.splitter*h,l=t-(t-h),a=h-l,t=e.splitter*L,n=t-(t-L),o=L-n,b=a*o-(r-l*n-a*n-l*o),f=S*N,t=e.splitter*S,l=t-(t-S),a=S-l,t=e.splitter*N,n=t-(t-N),o=N-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,G[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,G[1]=u-(c-s)+(f-s),v=i+c,s=v-i,G[2]=i-(v-s)+(c-s),G[3]=v,ae=e.sum(4,F,4,G,Se),r=y*g,t=e.splitter*y,l=t-(t-y),a=y-l,t=e.splitter*g,n=t-(t-g),o=g-n,b=a*o-(r-l*n-a*n-l*o),f=h*j,t=e.splitter*h,l=t-(t-h),a=h-l,t=e.splitter*j,n=t-(t-j),o=j-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,oe[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,oe[1]=u-(c+s)+(s-f),v=i+c,s=v-i,oe[2]=i-(v-s)+(c-s),oe[3]=v,xe=4):(Se[0]=0,ae=1,oe[0]=0,xe=1),M!==0){const D=e.scale(ae,Se,M,k);m=z(m,e.sum(e.scale(ve,Be,M,d),d,e.scale(D,k,2*q,P),P,w),w);const E=e.scale(xe,oe,M,p);m=z(m,e.sum_three(e.scale(E,p,2*q,d),d,e.scale(E,p,M,T),T,e.scale(D,k,M,P),P,te,H),H),j!==0&&(m=z(m,e.scale(e.scale(4,Q,M,p),p,j,d),d)),g!==0&&(m=z(m,e.scale(e.scale(4,U,-M,p),p,g,d),d))}if(x!==0){const D=e.scale(ae,Se,x,k);m=z(m,e.sum(e.scale(me,Ce,x,d),d,e.scale(D,k,2*B,P),P,w),w);const E=e.scale(xe,oe,x,p);m=z(m,e.sum_three(e.scale(E,p,2*B,d),d,e.scale(E,p,x,T),T,e.scale(D,k,x,P),P,te,H),H)}}if(y!==0||j!==0){if(h!==0||g!==0||M!==0||x!==0?(r=h*B,t=e.splitter*h,l=t-(t-h),a=h-l,t=e.splitter*B,n=t-(t-B),o=B-n,b=a*o-(r-l*n-a*n-l*o),f=S*x,t=e.splitter*S,l=t-(t-S),a=S-l,t=e.splitter*x,n=t-(t-x),o=x-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,F[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,F[1]=u-(c-s)+(f-s),v=i+c,s=v-i,F[2]=i-(v-s)+(c-s),F[3]=v,L=-C,N=-g,r=M*L,t=e.splitter*M,l=t-(t-M),a=M-l,t=e.splitter*L,n=t-(t-L),o=L-n,b=a*o-(r-l*n-a*n-l*o),f=q*N,t=e.splitter*q,l=t-(t-q),a=q-l,t=e.splitter*N,n=t-(t-N),o=N-n,_=a*o-(f-l*n-a*n-l*o),c=b+_,s=c-b,G[0]=b-(c-s)+(_-s),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u+f,s=c-u,G[1]=u-(c-s)+(f-s),v=i+c,s=v-i,G[2]=i-(v-s)+(c-s),G[3]=v,le=e.sum(4,F,4,G,ge),r=h*x,t=e.splitter*h,l=t-(t-h),a=h-l,t=e.splitter*x,n=t-(t-x),o=x-n,b=a*o-(r-l*n-a*n-l*o),f=M*g,t=e.splitter*M,l=t-(t-M),a=M-l,t=e.splitter*g,n=t-(t-g),o=g-n,_=a*o-(f-l*n-a*n-l*o),c=b-_,s=b-c,ie[0]=b-(c+s)+(s-_),i=r+c,s=i-r,u=r-(i-s)+(c-s),c=u-f,s=u-c,ie[1]=u-(c+s)+(s-f),v=i+c,s=v-i,ie[2]=i-(v-s)+(c-s),ie[3]=v,ne=4):(ge[0]=0,le=1,ie[0]=0,ne=1),y!==0){const D=e.scale(le,ge,y,k);m=z(m,e.sum(e.scale(de,Oe,y,d),d,e.scale(D,k,2*A,P),P,w),w);const E=e.scale(ne,ie,y,p);m=z(m,e.sum_three(e.scale(E,p,2*A,d),d,e.scale(E,p,y,T),T,e.scale(D,k,y,P),P,te,H),H),g!==0&&(m=z(m,e.scale(e.scale(4,R,y,p),p,g,d),d)),x!==0&&(m=z(m,e.scale(e.scale(4,Q,-y,p),p,x,d),d))}if(j!==0){const D=e.scale(le,ge,j,k);m=z(m,e.sum(e.scale(_e,Pe,j,d),d,e.scale(D,k,2*O,P),P,w),w);const E=e.scale(ne,ie,j,p);m=z(m,e.sum_three(e.scale(E,p,2*O,d),d,e.scale(E,p,j,T),T,e.scale(D,k,j,P),P,te,H),H)}}return he[m-1]}function Ee(K,W,X,Y,pe,ue,I,J){const ce=K-I,m=X-I,h=pe-I,M=W-J,y=Y-J,g=ue-J,x=m*g,j=h*y,be=ce*ce+M*M,fe=h*M,ve=ce*g,me=m*m+y*y,de=ce*y,_e=m*M,le=h*h+g*g,Z=be*(x-j)+me*(fe-ve)+le*(de-_e),ae=(Math.abs(x)+Math.abs(j))*be+(Math.abs(fe)+Math.abs(ve))*me+(Math.abs(de)+Math.abs(_e))*le,ne=ke*ae;return Z>ne||-Z>ne?Z:De(K,W,X,Y,pe,ue,I,J,ae)}exports.incircle=Ee;
+//# sourceMappingURL=index.cjs560.js.map

package/dist/index.cjs560.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.cjs560.js","sources":["../../../node_modules/@kninnug/constrainautor/node_modules/robust-predicates/esm/incircle.js"],"sourcesContent":["import {epsilon, splitter, resulterrbound, estimate, vec, sum, sum_three, scale} from './util.js';\n\nconst iccerrboundA = (10 + 96 * epsilon) * epsilon;\nconst iccerrboundB = (4 + 48 * epsilon) * epsilon;\nconst iccerrboundC = (44 + 576 * epsilon) * epsilon * epsilon;\n\nconst bc = vec(4);\nconst ca = vec(4);\nconst ab = vec(4);\nconst aa = vec(4);\nconst bb = vec(4);\nconst cc = vec(4);\nconst u = vec(4);\nconst v = vec(4);\nconst axtbc = vec(8);\nconst aytbc = vec(8);\nconst bxtca = vec(8);\nconst bytca = vec(8);\nconst cxtab = vec(8);\nconst cytab = vec(8);\nconst abt = vec(8);\nconst bct = vec(8);\nconst cat = vec(8);\nconst abtt = vec(4);\nconst bctt = vec(4);\nconst catt = vec(4);\n\nconst _8 = vec(8);\nconst _16 = vec(16);\nconst _16b = vec(16);\nconst _16c = vec(16);\nconst _32 = vec(32);\nconst _32b = vec(32);\nconst _48 = vec(48);\nconst _64 = vec(64);\n\nlet fin = vec(1152);\nlet fin2 = vec(1152);\n\nfunction finadd(finlen, a, alen) {\n    finlen = sum(finlen, fin, a, alen, fin2);\n    const tmp = fin; fin = fin2; fin2 = tmp;\n    return finlen;\n}\n\nfunction incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent) {\n    let finlen;\n    let adxtail, bdxtail, cdxtail, adytail, bdytail, cdytail;\n    let axtbclen, aytbclen, bxtcalen, bytcalen, cxtablen, cytablen;\n    let abtlen, bctlen, catlen;\n    let abttlen, bcttlen, cattlen;\n    let n1, n0;\n\n    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;\n\n    const adx = ax - dx;\n    const bdx = bx - dx;\n    const cdx = cx - dx;\n    const ady = ay - dy;\n    const bdy = by - dy;\n    const cdy = cy - dy;\n\n    s1 = bdx * cdy;\n    c = splitter * bdx;\n    ahi = c - (c - bdx);\n    alo = bdx - ahi;\n    c = splitter * cdy;\n    bhi = c - (c - cdy);\n    blo = cdy - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = cdx * bdy;\n    c = splitter * cdx;\n    ahi = c - (c - cdx);\n    alo = cdx - ahi;\n    c = splitter * bdy;\n    bhi = c - (c - bdy);\n    blo = bdy - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    bc[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    bc[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    bc[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    bc[3] = u3;\n    s1 = cdx * ady;\n    c = splitter * cdx;\n    ahi = c - (c - cdx);\n    alo = cdx - ahi;\n    c = splitter * ady;\n    bhi = c - (c - ady);\n    blo = ady - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = adx * cdy;\n    c = splitter * adx;\n    ahi = c - (c - adx);\n    alo = adx - ahi;\n    c = splitter * cdy;\n    bhi = c - (c - cdy);\n    blo = cdy - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    ca[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    ca[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    ca[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    ca[3] = u3;\n    s1 = adx * bdy;\n    c = splitter * adx;\n    ahi = c - (c - adx);\n    alo = adx - ahi;\n    c = splitter * bdy;\n    bhi = c - (c - bdy);\n    blo = bdy - bhi;\n    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n    t1 = bdx * ady;\n    c = splitter * bdx;\n    ahi = c - (c - bdx);\n    alo = bdx - ahi;\n    c = splitter * ady;\n    bhi = c - (c - ady);\n    blo = ady - bhi;\n    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n    _i = s0 - t0;\n    bvirt = s0 - _i;\n    ab[0] = s0 - (_i + bvirt) + (bvirt - t0);\n    _j = s1 + _i;\n    bvirt = _j - s1;\n    _0 = s1 - (_j - bvirt) + (_i - bvirt);\n    _i = _0 - t1;\n    bvirt = _0 - _i;\n    ab[1] = _0 - (_i + bvirt) + (bvirt - t1);\n    u3 = _j + _i;\n    bvirt = u3 - _j;\n    ab[2] = _j - (u3 - bvirt) + (_i - bvirt);\n    ab[3] = u3;\n\n    finlen = sum(\n        sum(\n            sum(\n                scale(scale(4, bc, adx, _8), _8, adx, _16), _16,\n                scale(scale(4, bc, ady, _8), _8, ady, _16b), _16b, _32), _32,\n            sum(\n                scale(scale(4, ca, bdx, _8), _8, bdx, _16), _16,\n                scale(scale(4, ca, bdy, _8), _8, bdy, _16b), _16b, _32b), _32b, _64), _64,\n        sum(\n            scale(scale(4, ab, cdx, _8), _8, cdx, _16), _16,\n            scale(scale(4, ab, cdy, _8), _8, cdy, _16b), _16b, _32), _32, fin);\n\n    let det = estimate(finlen, fin);\n    let errbound = iccerrboundB * permanent;\n    if (det >= errbound || -det >= errbound) {\n        return det;\n    }\n\n    bvirt = ax - adx;\n    adxtail = ax - (adx + bvirt) + (bvirt - dx);\n    bvirt = ay - ady;\n    adytail = ay - (ady + bvirt) + (bvirt - dy);\n    bvirt = bx - bdx;\n    bdxtail = bx - (bdx + bvirt) + (bvirt - dx);\n    bvirt = by - bdy;\n    bdytail = by - (bdy + bvirt) + (bvirt - dy);\n    bvirt = cx - cdx;\n    cdxtail = cx - (cdx + bvirt) + (bvirt - dx);\n    bvirt = cy - cdy;\n    cdytail = cy - (cdy + bvirt) + (bvirt - dy);\n    if (adxtail === 0 && bdxtail === 0 && cdxtail === 0 && adytail === 0 && bdytail === 0 && cdytail === 0) {\n        return det;\n    }\n\n    errbound = iccerrboundC * permanent + resulterrbound * Math.abs(det);\n    det += ((adx * adx + ady * ady) * ((bdx * cdytail + cdy * bdxtail) - (bdy * cdxtail + cdx * bdytail)) +\n        2 * (adx * adxtail + ady * adytail) * (bdx * cdy - bdy * cdx)) +\n        ((bdx * bdx + bdy * bdy) * ((cdx * adytail + ady * cdxtail) - (cdy * adxtail + adx * cdytail)) +\n        2 * (bdx * bdxtail + bdy * bdytail) * (cdx * ady - cdy * adx)) +\n        ((cdx * cdx + cdy * cdy) * ((adx * bdytail + bdy * adxtail) - (ady * bdxtail + bdx * adytail)) +\n        2 * (cdx * cdxtail + cdy * cdytail) * (adx * bdy - ady * bdx));\n\n    if (det >= errbound || -det >= errbound) {\n        return det;\n    }\n\n    if (bdxtail !== 0 || bdytail !== 0 || cdxtail !== 0 || cdytail !== 0) {\n        s1 = adx * adx;\n        c = splitter * adx;\n        ahi = c - (c - adx);\n        alo = adx - ahi;\n        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);\n        t1 = ady * ady;\n        c = splitter * ady;\n        ahi = c - (c - ady);\n        alo = ady - ahi;\n        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);\n        _i = s0 + t0;\n        bvirt = _i - s0;\n        aa[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n        _j = s1 + _i;\n        bvirt = _j - s1;\n        _0 = s1 - (_j - bvirt) + (_i - bvirt);\n        _i = _0 + t1;\n        bvirt = _i - _0;\n        aa[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n        u3 = _j + _i;\n        bvirt = u3 - _j;\n        aa[2] = _j - (u3 - bvirt) + (_i - bvirt);\n        aa[3] = u3;\n    }\n    if (cdxtail !== 0 || cdytail !== 0 || adxtail !== 0 || adytail !== 0) {\n        s1 = bdx * bdx;\n        c = splitter * bdx;\n        ahi = c - (c - bdx);\n        alo = bdx - ahi;\n        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);\n        t1 = bdy * bdy;\n        c = splitter * bdy;\n        ahi = c - (c - bdy);\n        alo = bdy - ahi;\n        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);\n        _i = s0 + t0;\n        bvirt = _i - s0;\n        bb[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n        _j = s1 + _i;\n        bvirt = _j - s1;\n        _0 = s1 - (_j - bvirt) + (_i - bvirt);\n        _i = _0 + t1;\n        bvirt = _i - _0;\n        bb[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n        u3 = _j + _i;\n        bvirt = u3 - _j;\n        bb[2] = _j - (u3 - bvirt) + (_i - bvirt);\n        bb[3] = u3;\n    }\n    if (adxtail !== 0 || adytail !== 0 || bdxtail !== 0 || bdytail !== 0) {\n        s1 = cdx * cdx;\n        c = splitter * cdx;\n        ahi = c - (c - cdx);\n        alo = cdx - ahi;\n        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);\n        t1 = cdy * cdy;\n        c = splitter * cdy;\n        ahi = c - (c - cdy);\n        alo = cdy - ahi;\n        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);\n        _i = s0 + t0;\n        bvirt = _i - s0;\n        cc[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n        _j = s1 + _i;\n        bvirt = _j - s1;\n        _0 = s1 - (_j - bvirt) + (_i - bvirt);\n        _i = _0 + t1;\n        bvirt = _i - _0;\n        cc[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n        u3 = _j + _i;\n        bvirt = u3 - _j;\n        cc[2] = _j - (u3 - bvirt) + (_i - bvirt);\n        cc[3] = u3;\n    }\n\n    if (adxtail !== 0) {\n        axtbclen = scale(4, bc, adxtail, axtbc);\n        finlen = finadd(finlen, sum_three(\n            scale(axtbclen, axtbc, 2 * adx, _16), _16,\n            scale(scale(4, cc, adxtail, _8), _8, bdy, _16b), _16b,\n            scale(scale(4, bb, adxtail, _8), _8, -cdy, _16c), _16c, _32, _48), _48);\n    }\n    if (adytail !== 0) {\n        aytbclen = scale(4, bc, adytail, aytbc);\n        finlen = finadd(finlen, sum_three(\n            scale(aytbclen, aytbc, 2 * ady, _16), _16,\n            scale(scale(4, bb, adytail, _8), _8, cdx, _16b), _16b,\n            scale(scale(4, cc, adytail, _8), _8, -bdx, _16c), _16c, _32, _48), _48);\n    }\n    if (bdxtail !== 0) {\n        bxtcalen = scale(4, ca, bdxtail, bxtca);\n        finlen = finadd(finlen, sum_three(\n            scale(bxtcalen, bxtca, 2 * bdx, _16), _16,\n            scale(scale(4, aa, bdxtail, _8), _8, cdy, _16b), _16b,\n            scale(scale(4, cc, bdxtail, _8), _8, -ady, _16c), _16c, _32, _48), _48);\n    }\n    if (bdytail !== 0) {\n        bytcalen = scale(4, ca, bdytail, bytca);\n        finlen = finadd(finlen, sum_three(\n            scale(bytcalen, bytca, 2 * bdy, _16), _16,\n            scale(scale(4, cc, bdytail, _8), _8, adx, _16b), _16b,\n            scale(scale(4, aa, bdytail, _8), _8, -cdx, _16c), _16c, _32, _48), _48);\n    }\n    if (cdxtail !== 0) {\n        cxtablen = scale(4, ab, cdxtail, cxtab);\n        finlen = finadd(finlen, sum_three(\n            scale(cxtablen, cxtab, 2 * cdx, _16), _16,\n            scale(scale(4, bb, cdxtail, _8), _8, ady, _16b), _16b,\n            scale(scale(4, aa, cdxtail, _8), _8, -bdy, _16c), _16c, _32, _48), _48);\n    }\n    if (cdytail !== 0) {\n        cytablen = scale(4, ab, cdytail, cytab);\n        finlen = finadd(finlen, sum_three(\n            scale(cytablen, cytab, 2 * cdy, _16), _16,\n            scale(scale(4, aa, cdytail, _8), _8, bdx, _16b), _16b,\n            scale(scale(4, bb, cdytail, _8), _8, -adx, _16c), _16c, _32, _48), _48);\n    }\n\n    if (adxtail !== 0 || adytail !== 0) {\n        if (bdxtail !== 0 || bdytail !== 0 || cdxtail !== 0 || cdytail !== 0) {\n            s1 = bdxtail * cdy;\n            c = splitter * bdxtail;\n            ahi = c - (c - bdxtail);\n            alo = bdxtail - ahi;\n            c = splitter * cdy;\n            bhi = c - (c - cdy);\n            blo = cdy - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = bdx * cdytail;\n            c = splitter * bdx;\n            ahi = c - (c - bdx);\n            alo = bdx - ahi;\n            c = splitter * cdytail;\n            bhi = c - (c - cdytail);\n            blo = cdytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            u[3] = u3;\n            s1 = cdxtail * -bdy;\n            c = splitter * cdxtail;\n            ahi = c - (c - cdxtail);\n            alo = cdxtail - ahi;\n            c = splitter * -bdy;\n            bhi = c - (c - -bdy);\n            blo = -bdy - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = cdx * -bdytail;\n            c = splitter * cdx;\n            ahi = c - (c - cdx);\n            alo = cdx - ahi;\n            c = splitter * -bdytail;\n            bhi = c - (c - -bdytail);\n            blo = -bdytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            v[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            v[3] = u3;\n            bctlen = sum(4, u, 4, v, bct);\n            s1 = bdxtail * cdytail;\n            c = splitter * bdxtail;\n            ahi = c - (c - bdxtail);\n            alo = bdxtail - ahi;\n            c = splitter * cdytail;\n            bhi = c - (c - cdytail);\n            blo = cdytail - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = cdxtail * bdytail;\n            c = splitter * cdxtail;\n            ahi = c - (c - cdxtail);\n            alo = cdxtail - ahi;\n            c = splitter * bdytail;\n            bhi = c - (c - bdytail);\n            blo = bdytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 - t0;\n            bvirt = s0 - _i;\n            bctt[0] = s0 - (_i + bvirt) + (bvirt - t0);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 - t1;\n            bvirt = _0 - _i;\n            bctt[1] = _0 - (_i + bvirt) + (bvirt - t1);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            bctt[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            bctt[3] = u3;\n            bcttlen = 4;\n        } else {\n            bct[0] = 0;\n            bctlen = 1;\n            bctt[0] = 0;\n            bcttlen = 1;\n        }\n        if (adxtail !== 0) {\n            const len = scale(bctlen, bct, adxtail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(axtbclen, axtbc, adxtail, _16), _16,\n                scale(len, _16c, 2 * adx, _32), _32, _48), _48);\n\n            const len2 = scale(bcttlen, bctt, adxtail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * adx, _16), _16,\n                scale(len2, _8, adxtail, _16b), _16b,\n                scale(len, _16c, adxtail, _32), _32, _32b, _64), _64);\n\n            if (bdytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, cc, adxtail, _8), _8, bdytail, _16), _16);\n            }\n            if (cdytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, bb, -adxtail, _8), _8, cdytail, _16), _16);\n            }\n        }\n        if (adytail !== 0) {\n            const len = scale(bctlen, bct, adytail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(aytbclen, aytbc, adytail, _16), _16,\n                scale(len, _16c, 2 * ady, _32), _32, _48), _48);\n\n            const len2 = scale(bcttlen, bctt, adytail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * ady, _16), _16,\n                scale(len2, _8, adytail, _16b), _16b,\n                scale(len, _16c, adytail, _32), _32, _32b, _64), _64);\n        }\n    }\n    if (bdxtail !== 0 || bdytail !== 0) {\n        if (cdxtail !== 0 || cdytail !== 0 || adxtail !== 0 || adytail !== 0) {\n            s1 = cdxtail * ady;\n            c = splitter * cdxtail;\n            ahi = c - (c - cdxtail);\n            alo = cdxtail - ahi;\n            c = splitter * ady;\n            bhi = c - (c - ady);\n            blo = ady - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = cdx * adytail;\n            c = splitter * cdx;\n            ahi = c - (c - cdx);\n            alo = cdx - ahi;\n            c = splitter * adytail;\n            bhi = c - (c - adytail);\n            blo = adytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            u[3] = u3;\n            n1 = -cdy;\n            n0 = -cdytail;\n            s1 = adxtail * n1;\n            c = splitter * adxtail;\n            ahi = c - (c - adxtail);\n            alo = adxtail - ahi;\n            c = splitter * n1;\n            bhi = c - (c - n1);\n            blo = n1 - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = adx * n0;\n            c = splitter * adx;\n            ahi = c - (c - adx);\n            alo = adx - ahi;\n            c = splitter * n0;\n            bhi = c - (c - n0);\n            blo = n0 - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            v[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            v[3] = u3;\n            catlen = sum(4, u, 4, v, cat);\n            s1 = cdxtail * adytail;\n            c = splitter * cdxtail;\n            ahi = c - (c - cdxtail);\n            alo = cdxtail - ahi;\n            c = splitter * adytail;\n            bhi = c - (c - adytail);\n            blo = adytail - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = adxtail * cdytail;\n            c = splitter * adxtail;\n            ahi = c - (c - adxtail);\n            alo = adxtail - ahi;\n            c = splitter * cdytail;\n            bhi = c - (c - cdytail);\n            blo = cdytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 - t0;\n            bvirt = s0 - _i;\n            catt[0] = s0 - (_i + bvirt) + (bvirt - t0);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 - t1;\n            bvirt = _0 - _i;\n            catt[1] = _0 - (_i + bvirt) + (bvirt - t1);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            catt[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            catt[3] = u3;\n            cattlen = 4;\n        } else {\n            cat[0] = 0;\n            catlen = 1;\n            catt[0] = 0;\n            cattlen = 1;\n        }\n        if (bdxtail !== 0) {\n            const len = scale(catlen, cat, bdxtail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(bxtcalen, bxtca, bdxtail, _16), _16,\n                scale(len, _16c, 2 * bdx, _32), _32, _48), _48);\n\n            const len2 = scale(cattlen, catt, bdxtail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * bdx, _16), _16,\n                scale(len2, _8, bdxtail, _16b), _16b,\n                scale(len, _16c, bdxtail, _32), _32, _32b, _64), _64);\n\n            if (cdytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, aa, bdxtail, _8), _8, cdytail, _16), _16);\n            }\n            if (adytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, cc, -bdxtail, _8), _8, adytail, _16), _16);\n            }\n        }\n        if (bdytail !== 0) {\n            const len = scale(catlen, cat, bdytail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(bytcalen, bytca, bdytail, _16), _16,\n                scale(len, _16c, 2 * bdy, _32), _32, _48), _48);\n\n            const len2 = scale(cattlen, catt, bdytail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * bdy, _16), _16,\n                scale(len2, _8, bdytail, _16b), _16b,\n                scale(len, _16c, bdytail, _32), _32,  _32b, _64), _64);\n        }\n    }\n    if (cdxtail !== 0 || cdytail !== 0) {\n        if (adxtail !== 0 || adytail !== 0 || bdxtail !== 0 || bdytail !== 0) {\n            s1 = adxtail * bdy;\n            c = splitter * adxtail;\n            ahi = c - (c - adxtail);\n            alo = adxtail - ahi;\n            c = splitter * bdy;\n            bhi = c - (c - bdy);\n            blo = bdy - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = adx * bdytail;\n            c = splitter * adx;\n            ahi = c - (c - adx);\n            alo = adx - ahi;\n            c = splitter * bdytail;\n            bhi = c - (c - bdytail);\n            blo = bdytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            u[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            u[3] = u3;\n            n1 = -ady;\n            n0 = -adytail;\n            s1 = bdxtail * n1;\n            c = splitter * bdxtail;\n            ahi = c - (c - bdxtail);\n            alo = bdxtail - ahi;\n            c = splitter * n1;\n            bhi = c - (c - n1);\n            blo = n1 - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = bdx * n0;\n            c = splitter * bdx;\n            ahi = c - (c - bdx);\n            alo = bdx - ahi;\n            c = splitter * n0;\n            bhi = c - (c - n0);\n            blo = n0 - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 + t0;\n            bvirt = _i - s0;\n            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 + t1;\n            bvirt = _i - _0;\n            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            v[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            v[3] = u3;\n            abtlen = sum(4, u, 4, v, abt);\n            s1 = adxtail * bdytail;\n            c = splitter * adxtail;\n            ahi = c - (c - adxtail);\n            alo = adxtail - ahi;\n            c = splitter * bdytail;\n            bhi = c - (c - bdytail);\n            blo = bdytail - bhi;\n            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);\n            t1 = bdxtail * adytail;\n            c = splitter * bdxtail;\n            ahi = c - (c - bdxtail);\n            alo = bdxtail - ahi;\n            c = splitter * adytail;\n            bhi = c - (c - adytail);\n            blo = adytail - bhi;\n            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);\n            _i = s0 - t0;\n            bvirt = s0 - _i;\n            abtt[0] = s0 - (_i + bvirt) + (bvirt - t0);\n            _j = s1 + _i;\n            bvirt = _j - s1;\n            _0 = s1 - (_j - bvirt) + (_i - bvirt);\n            _i = _0 - t1;\n            bvirt = _0 - _i;\n            abtt[1] = _0 - (_i + bvirt) + (bvirt - t1);\n            u3 = _j + _i;\n            bvirt = u3 - _j;\n            abtt[2] = _j - (u3 - bvirt) + (_i - bvirt);\n            abtt[3] = u3;\n            abttlen = 4;\n        } else {\n            abt[0] = 0;\n            abtlen = 1;\n            abtt[0] = 0;\n            abttlen = 1;\n        }\n        if (cdxtail !== 0) {\n            const len = scale(abtlen, abt, cdxtail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(cxtablen, cxtab, cdxtail, _16), _16,\n                scale(len, _16c, 2 * cdx, _32), _32, _48), _48);\n\n            const len2 = scale(abttlen, abtt, cdxtail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * cdx, _16), _16,\n                scale(len2, _8, cdxtail, _16b), _16b,\n                scale(len, _16c, cdxtail, _32), _32, _32b, _64), _64);\n\n            if (adytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, bb, cdxtail, _8), _8, adytail, _16), _16);\n            }\n            if (bdytail !== 0) {\n                finlen = finadd(finlen, scale(scale(4, aa, -cdxtail, _8), _8, bdytail, _16), _16);\n            }\n        }\n        if (cdytail !== 0) {\n            const len = scale(abtlen, abt, cdytail, _16c);\n            finlen = finadd(finlen, sum(\n                scale(cytablen, cytab, cdytail, _16), _16,\n                scale(len, _16c, 2 * cdy, _32), _32, _48), _48);\n\n            const len2 = scale(abttlen, abtt, cdytail, _8);\n            finlen = finadd(finlen, sum_three(\n                scale(len2, _8, 2 * cdy, _16), _16,\n                scale(len2, _8, cdytail, _16b), _16b,\n                scale(len, _16c, cdytail, _32), _32, _32b, _64), _64);\n        }\n    }\n\n    return fin[finlen - 1];\n}\n\nexport function incircle(ax, ay, bx, by, cx, cy, dx, dy) {\n    const adx = ax - dx;\n    const bdx = bx - dx;\n    const cdx = cx - dx;\n    const ady = ay - dy;\n    const bdy = by - dy;\n    const cdy = cy - dy;\n\n    const bdxcdy = bdx * cdy;\n    const cdxbdy = cdx * bdy;\n    const alift = adx * adx + ady * ady;\n\n    const cdxady = cdx * ady;\n    const adxcdy = adx * cdy;\n    const blift = bdx * bdx + bdy * bdy;\n\n    const adxbdy = adx * bdy;\n    const bdxady = bdx * ady;\n    const clift = cdx * cdx + cdy * cdy;\n\n    const det =\n        alift * (bdxcdy - cdxbdy) +\n        blift * (cdxady - adxcdy) +\n        clift * (adxbdy - bdxady);\n\n    const permanent =\n        (Math.abs(bdxcdy) + Math.abs(cdxbdy)) * alift +\n        (Math.abs(cdxady) + Math.abs(adxcdy)) * blift +\n        (Math.abs(adxbdy) + Math.abs(bdxady)) * clift;\n\n    const errbound = iccerrboundA * permanent;\n\n    if (det > errbound || -det > errbound) {\n        return det;\n    }\n    return incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent);\n}\n\nexport function incirclefast(ax, ay, bx, by, cx, cy, dx, dy) {\n    const adx = ax - dx;\n    const ady = ay - dy;\n    const bdx = bx - dx;\n    const bdy = by - dy;\n    const cdx = cx - dx;\n    const cdy = cy - dy;\n\n    const abdet = adx * bdy - bdx * ady;\n    const bcdet = bdx * cdy - cdx * bdy;\n    const cadet = cdx * ady - adx * cdy;\n    const alift = adx * adx + ady * ady;\n    const blift = bdx * bdx + bdy * bdy;\n    const clift = cdx * cdx + cdy * cdy;\n\n    return alift * bcdet + blift * cadet + clift * abdet;\n}\n"],"names":["iccerrboundA","epsilon","iccerrboundB","iccerrboundC","bc","vec","ca","ab","aa","bb","cc","u","v","axtbc","aytbc","bxtca","bytca","cxtab","cytab","abt","bct","cat","abtt","bctt","catt","_8","_16","_16b","_16c","_32","_32b","_48","_64","fin","fin2","finadd","finlen","a","alen","sum","tmp","incircleadapt","ax","ay","bx","by","cx","cy","dx","dy","permanent","adxtail","bdxtail","cdxtail","adytail","bdytail","cdytail","axtbclen","aytbclen","bxtcalen","bytcalen","cxtablen","cytablen","abtlen","bctlen","catlen","abttlen","bcttlen","cattlen","n1","n0","bvirt","c","ahi","alo","bhi","blo","_i","_j","_0","s1","s0","t1","t0","u3","adx","bdx","cdx","ady","bdy","cdy","splitter","scale","det","estimate","errbound","resulterrbound","sum_three","len","len2","incircle","bdxcdy","cdxbdy","alift","cdxady","adxcdy","blift","adxbdy","bdxady","clift"],"mappings":"qHAEMA,IAAgB,GAAK,GAAKC,EAAAA,SAAWA,EAAAA,QACrCC,IAAgB,EAAI,GAAKD,EAAAA,SAAWA,EAAAA,QACpCE,IAAgB,GAAK,IAAMF,EAAAA,SAAWA,EAAAA,QAAUA,EAAAA,QAEhDG,EAAKC,EAAAA,IAAI,CAAC,EACVC,GAAKD,EAAAA,IAAI,CAAC,EACVE,GAAKF,EAAAA,IAAI,CAAC,EACVG,EAAKH,EAAAA,IAAI,CAAC,EACVI,EAAKJ,EAAAA,IAAI,CAAC,EACVK,EAAKL,EAAAA,IAAI,CAAC,EACVM,EAAIN,EAAAA,IAAI,CAAC,EACTO,EAAIP,EAAAA,IAAI,CAAC,EACTQ,GAAQR,EAAAA,IAAI,CAAC,EACbS,GAAQT,EAAAA,IAAI,CAAC,EACbU,GAAQV,EAAAA,IAAI,CAAC,EACbW,GAAQX,EAAAA,IAAI,CAAC,EACbY,GAAQZ,EAAAA,IAAI,CAAC,EACba,GAAQb,EAAAA,IAAI,CAAC,EACbc,GAAMd,EAAAA,IAAI,CAAC,EACXe,GAAMf,EAAAA,IAAI,CAAC,EACXgB,GAAMhB,EAAAA,IAAI,CAAC,EACXiB,GAAOjB,EAAAA,IAAI,CAAC,EACZkB,GAAOlB,EAAAA,IAAI,CAAC,EACZmB,GAAOnB,EAAAA,IAAI,CAAC,EAEZoB,EAAKpB,EAAAA,IAAI,CAAC,EACVqB,EAAMrB,EAAAA,IAAI,EAAE,EACZsB,EAAOtB,EAAAA,IAAI,EAAE,EACbuB,EAAOvB,EAAAA,IAAI,EAAE,EACbwB,EAAMxB,EAAAA,IAAI,EAAE,EACZyB,GAAOzB,EAAAA,IAAI,EAAE,EACb0B,EAAM1B,EAAAA,IAAI,EAAE,EACZ2B,EAAM3B,EAAAA,IAAI,EAAE,EAElB,IAAI4B,GAAM5B,EAAAA,IAAI,IAAI,EACd6B,GAAO7B,EAAAA,IAAI,IAAI,EAEnB,SAAS8B,EAAOC,EAAQC,EAAGC,EAAM,CAC7BF,EAASG,EAAAA,IAAIH,EAAQH,GAAKI,EAAGC,EAAMJ,EAAI,EACvC,MAAMM,EAAMP,GAAK,OAAAA,GAAMC,GAAMA,GAAOM,EAC7BJ,CACX,CAEA,SAASK,GAAcC,EAAIC,EAAIC,EAAIC,EAAIC,GAAIC,GAAIC,EAAIC,EAAIC,GAAW,CAC9D,IAAId,EACAe,EAASC,EAASC,EAASC,EAASC,EAASC,EAC7CC,GAAUC,GAAUC,GAAUC,GAAUC,GAAUC,GAClDC,GAAQC,EAAQC,GAChBC,GAASC,GAASC,GAClBC,EAAIC,EAEJC,EAAOC,EAAGC,EAAKC,EAAKC,EAAKC,EAAKC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAAIC,EAE9D,MAAMC,EAAM3C,EAAKM,EACXsC,EAAM1C,EAAKI,EACXuC,EAAMzC,GAAKE,EACXwC,EAAM7C,EAAKM,EACXwC,EAAM5C,EAAKI,EACXyC,EAAM3C,GAAKE,EAEjB+B,EAAKM,EAAMI,EACXlB,EAAImB,EAAAA,SAAWL,EACfb,EAAMD,GAAKA,EAAIc,GACfZ,EAAMY,EAAMb,EACZD,EAAImB,EAAAA,SAAWD,EACff,EAAMH,GAAKA,EAAIkB,GACfd,EAAMc,EAAMf,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKK,EAAME,EACXjB,EAAImB,EAAAA,SAAWJ,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAImB,EAAAA,SAAWF,EACfd,EAAMH,GAAKA,EAAIiB,GACfb,EAAMa,EAAMd,EACZQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbzE,EAAG,CAAC,EAAI6E,GAAMJ,EAAKN,IAAUA,EAAQY,GACrCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbzE,EAAG,CAAC,EAAI2E,GAAMF,EAAKN,IAAUA,EAAQW,GACrCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACb1E,EAAG,CAAC,EAAI0E,GAAMM,EAAKb,IAAUM,EAAKN,GAClCnE,EAAG,CAAC,EAAIgF,EACRJ,EAAKO,EAAMC,EACXhB,EAAImB,EAAAA,SAAWJ,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAImB,EAAAA,SAAWH,EACfb,EAAMH,GAAKA,EAAIgB,GACfZ,EAAMY,EAAMb,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKG,EAAMK,EACXlB,EAAImB,EAAAA,SAAWN,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAImB,EAAAA,SAAWD,EACff,EAAMH,GAAKA,EAAIkB,GACfd,EAAMc,EAAMf,EACZQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbvE,GAAG,CAAC,EAAI2E,GAAMJ,EAAKN,IAAUA,EAAQY,GACrCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbvE,GAAG,CAAC,EAAIyE,GAAMF,EAAKN,IAAUA,EAAQW,GACrCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbxE,GAAG,CAAC,EAAIwE,GAAMM,EAAKb,IAAUM,EAAKN,GAClCjE,GAAG,CAAC,EAAI8E,EACRJ,EAAKK,EAAMI,EACXjB,EAAImB,EAAAA,SAAWN,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAImB,EAAAA,SAAWF,EACfd,EAAMH,GAAKA,EAAIiB,GACfb,EAAMa,EAAMd,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKI,EAAME,EACXhB,EAAImB,EAAAA,SAAWL,EACfb,EAAMD,GAAKA,EAAIc,GACfZ,EAAMY,EAAMb,EACZD,EAAImB,EAAAA,SAAWH,EACfb,EAAMH,GAAKA,EAAIgB,GACfZ,EAAMY,EAAMb,EACZQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbtE,GAAG,CAAC,EAAI0E,GAAMJ,EAAKN,IAAUA,EAAQY,GACrCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbtE,GAAG,CAAC,EAAIwE,GAAMF,EAAKN,IAAUA,EAAQW,GACrCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbvE,GAAG,CAAC,EAAIuE,GAAMM,EAAKb,IAAUM,EAAKN,GAClChE,GAAG,CAAC,EAAI6E,EAERhD,EAASG,EAAAA,IACLA,EAAAA,IACIA,EAAAA,IACIqD,QAAMA,EAAAA,MAAM,EAAGxF,EAAIiF,EAAK5D,CAAE,EAAGA,EAAI4D,EAAK3D,CAAG,EAAGA,EAC5CkE,QAAMA,EAAAA,MAAM,EAAGxF,EAAIoF,EAAK/D,CAAE,EAAGA,EAAI+D,EAAK7D,CAAI,EAAGA,EAAME,CAAA,EAAMA,EAC7DU,EAAAA,IACIqD,QAAMA,EAAAA,MAAM,EAAGtF,GAAIgF,EAAK7D,CAAE,EAAGA,EAAI6D,EAAK5D,CAAG,EAAGA,EAC5CkE,QAAMA,EAAAA,MAAM,EAAGtF,GAAImF,EAAKhE,CAAE,EAAGA,EAAIgE,EAAK9D,CAAI,EAAGA,EAAMG,EAAA,EAAOA,GAAME,CAAA,EAAMA,EAC9EO,EAAAA,IACIqD,QAAMA,EAAAA,MAAM,EAAGrF,GAAIgF,EAAK9D,CAAE,EAAGA,EAAI8D,EAAK7D,CAAG,EAAGA,EAC5CkE,QAAMA,EAAAA,MAAM,EAAGrF,GAAImF,EAAKjE,CAAE,EAAGA,EAAIiE,EAAK/D,CAAI,EAAGA,EAAME,CAAA,EAAMA,EAAKI,EAAA,EAEtE,IAAI4D,EAAMC,EAAAA,SAAS1D,EAAQH,EAAG,EAC1B8D,GAAW7F,GAAegD,GA6B9B,GA5BI2C,GAAOE,IAAY,CAACF,GAAOE,KAI/BxB,EAAQ7B,EAAK2C,EACblC,EAAUT,GAAM2C,EAAMd,IAAUA,EAAQvB,GACxCuB,EAAQ5B,EAAK6C,EACblC,EAAUX,GAAM6C,EAAMjB,IAAUA,EAAQtB,GACxCsB,EAAQ3B,EAAK0C,EACblC,EAAUR,GAAM0C,EAAMf,IAAUA,EAAQvB,GACxCuB,EAAQ1B,EAAK4C,EACblC,EAAUV,GAAM4C,EAAMlB,IAAUA,EAAQtB,GACxCsB,EAAQzB,GAAKyC,EACblC,EAAUP,IAAMyC,EAAMhB,IAAUA,EAAQvB,GACxCuB,EAAQxB,GAAK2C,EACblC,EAAUT,IAAM2C,EAAMnB,IAAUA,EAAQtB,GACpCE,IAAY,GAAKC,IAAY,GAAKC,IAAY,GAAKC,IAAY,GAAKC,IAAY,GAAKC,IAAY,KAIrGuC,GAAW5F,GAAe+C,GAAY8C,EAAAA,eAAiB,KAAK,IAAIH,CAAG,EACnEA,IAASR,EAAMA,EAAMG,EAAMA,IAASF,EAAM9B,EAAUkC,EAAMtC,GAAYqC,EAAMpC,EAAUkC,EAAMhC,IACxF,GAAK8B,EAAMlC,EAAUqC,EAAMlC,IAAYgC,EAAMI,EAAMD,EAAMF,KACvDD,EAAMA,EAAMG,EAAMA,IAASF,EAAMjC,EAAUkC,EAAMnC,GAAYqC,EAAMvC,EAAUkC,EAAM7B,IACrF,GAAK8B,EAAMlC,EAAUqC,EAAMlC,IAAYgC,EAAMC,EAAME,EAAML,MACvDE,EAAMA,EAAMG,EAAMA,IAASL,EAAM9B,EAAUkC,EAAMtC,GAAYqC,EAAMpC,EAAUkC,EAAMhC,IACrF,GAAKiC,EAAMlC,EAAUqC,EAAMlC,IAAY6B,EAAMI,EAAMD,EAAMF,IAEzDO,GAAOE,IAAY,CAACF,GAAOE,IAC3B,OAAOF,EA0HX,IAvHIzC,IAAY,GAAKG,IAAY,GAAKF,IAAY,GAAKG,IAAY,KAC/DwB,EAAKK,EAAMA,EACXb,EAAImB,EAAAA,SAAWN,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZQ,EAAKP,EAAMA,GAAOM,EAAKP,EAAMA,GAAOA,EAAMA,GAAOC,GACjDQ,EAAKM,EAAMA,EACXhB,EAAImB,EAAAA,SAAWH,EACff,EAAMD,GAAKA,EAAIgB,GACfd,EAAMc,EAAMf,EACZU,EAAKT,EAAMA,GAAOQ,EAAKT,EAAMA,GAAOA,EAAMA,GAAOC,GACjDG,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbzE,EAAG,CAAC,EAAIyE,GAAMJ,EAAKN,IAAUY,EAAKZ,GAClCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbvE,EAAG,CAAC,EAAIuE,GAAMF,EAAKN,IAAUW,EAAKX,GAClCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbtE,EAAG,CAAC,EAAIsE,GAAMM,EAAKb,IAAUM,EAAKN,GAClC/D,EAAG,CAAC,EAAI4E,IAER/B,IAAY,GAAKG,IAAY,GAAKL,IAAY,GAAKG,IAAY,KAC/D0B,EAAKM,EAAMA,EACXd,EAAImB,EAAAA,SAAWL,EACfb,EAAMD,GAAKA,EAAIc,GACfZ,EAAMY,EAAMb,EACZQ,EAAKP,EAAMA,GAAOM,EAAKP,EAAMA,GAAOA,EAAMA,GAAOC,GACjDQ,EAAKO,EAAMA,EACXjB,EAAImB,EAAAA,SAAWF,EACfhB,EAAMD,GAAKA,EAAIiB,GACff,EAAMe,EAAMhB,EACZU,EAAKT,EAAMA,GAAOQ,EAAKT,EAAMA,GAAOA,EAAMA,GAAOC,GACjDG,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbxE,EAAG,CAAC,EAAIwE,GAAMJ,EAAKN,IAAUY,EAAKZ,GAClCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbtE,EAAG,CAAC,EAAIsE,GAAMF,EAAKN,IAAUW,EAAKX,GAClCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbrE,EAAG,CAAC,EAAIqE,GAAMM,EAAKb,IAAUM,EAAKN,GAClC9D,EAAG,CAAC,EAAI2E,IAERjC,IAAY,GAAKG,IAAY,GAAKF,IAAY,GAAKG,IAAY,KAC/DyB,EAAKO,EAAMA,EACXf,EAAImB,EAAAA,SAAWJ,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZQ,EAAKP,EAAMA,GAAOM,EAAKP,EAAMA,GAAOA,EAAMA,GAAOC,GACjDQ,EAAKQ,EAAMA,EACXlB,EAAImB,EAAAA,SAAWD,EACfjB,EAAMD,GAAKA,EAAIkB,GACfhB,EAAMgB,EAAMjB,EACZU,EAAKT,EAAMA,GAAOQ,EAAKT,EAAMA,GAAOA,EAAMA,GAAOC,GACjDG,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbvE,EAAG,CAAC,EAAIuE,GAAMJ,EAAKN,IAAUY,EAAKZ,GAClCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbrE,EAAG,CAAC,EAAIqE,GAAMF,EAAKN,IAAUW,EAAKX,GAClCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbpE,EAAG,CAAC,EAAIoE,GAAMM,EAAKb,IAAUM,EAAKN,GAClC7D,EAAG,CAAC,EAAI0E,GAGRjC,IAAY,IACZM,GAAWmC,EAAAA,MAAM,EAAGxF,EAAI+C,EAAStC,EAAK,EACtCuB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMnC,GAAU5C,GAAO,EAAIwE,EAAK3D,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGlF,EAAIyC,EAAS1B,CAAE,EAAGA,EAAIgE,EAAK9D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGnF,EAAI0C,EAAS1B,CAAE,EAAGA,EAAI,CAACiE,EAAK9D,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAE1EuB,IAAY,IACZI,GAAWkC,EAAAA,MAAM,EAAGxF,EAAIkD,EAASxC,EAAK,EACtCsB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMlC,GAAU5C,GAAO,EAAI0E,EAAK9D,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGnF,EAAI6C,EAAS7B,CAAE,EAAGA,EAAI8D,EAAK5D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGlF,EAAI4C,EAAS7B,CAAE,EAAGA,EAAI,CAAC6D,EAAK1D,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAE1EqB,IAAY,IACZO,GAAWiC,EAAAA,MAAM,EAAGtF,GAAI8C,EAASrC,EAAK,EACtCqB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMjC,GAAU5C,GAAO,EAAIuE,EAAK5D,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGpF,EAAI4C,EAAS3B,CAAE,EAAGA,EAAIiE,EAAK/D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGlF,EAAI0C,EAAS3B,CAAE,EAAGA,EAAI,CAAC+D,EAAK5D,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAE1EwB,IAAY,IACZK,GAAWgC,EAAAA,MAAM,EAAGtF,GAAIiD,EAASvC,EAAK,EACtCoB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMhC,GAAU5C,GAAO,EAAIyE,EAAK/D,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGlF,EAAI6C,EAAS9B,CAAE,EAAGA,EAAI4D,EAAK1D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGpF,EAAI+C,EAAS9B,CAAE,EAAGA,EAAI,CAAC8D,EAAK3D,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAE1EsB,IAAY,IACZQ,GAAW+B,EAAAA,MAAM,EAAGrF,GAAI8C,EAASpC,EAAK,EACtCmB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAM/B,GAAU5C,GAAO,EAAIsE,EAAK7D,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGnF,EAAI4C,EAAS5B,CAAE,EAAGA,EAAI+D,EAAK7D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGpF,EAAI6C,EAAS5B,CAAE,EAAGA,EAAI,CAACgE,EAAK7D,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAE1EyB,IAAY,IACZM,GAAW8B,EAAAA,MAAM,EAAGrF,GAAIiD,EAAStC,EAAK,EACtCkB,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAM9B,GAAU5C,GAAO,EAAIwE,EAAKhE,CAAG,EAAGA,EACtCkE,QAAMA,EAAAA,MAAM,EAAGpF,EAAIgD,EAAS/B,CAAE,EAAGA,EAAI6D,EAAK3D,CAAI,EAAGA,EACjDiE,QAAMA,EAAAA,MAAM,EAAGnF,EAAI+C,EAAS/B,CAAE,EAAGA,EAAI,CAAC4D,EAAKzD,CAAI,EAAGA,EAAMC,EAAKE,CAAA,EAAMA,CAAG,GAG1EoB,IAAY,GAAKG,IAAY,EAAG,CAiGhC,GAhGIF,IAAY,GAAKG,IAAY,GAAKF,IAAY,GAAKG,IAAY,GAC/DwB,EAAK5B,EAAUsC,EACflB,EAAImB,EAAAA,SAAWvC,EACfqB,EAAMD,GAAKA,EAAIpB,GACfsB,EAAMtB,EAAUqB,EAChBD,EAAImB,EAAAA,SAAWD,EACff,EAAMH,GAAKA,EAAIkB,GACfd,EAAMc,EAAMf,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKI,EAAM9B,EACXgB,EAAImB,EAAAA,SAAWL,EACfb,EAAMD,GAAKA,EAAIc,GACfZ,EAAMY,EAAMb,EACZD,EAAImB,EAAAA,SAAWnC,EACfmB,EAAMH,GAAKA,EAAIhB,GACfoB,EAAMpB,EAAUmB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbtE,EAAE,CAAC,EAAIsE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbpE,EAAE,CAAC,EAAIoE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbnE,EAAE,CAAC,EAAImE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC5D,EAAE,CAAC,EAAIyE,EACPJ,EAAK3B,EAAU,CAACoC,EAChBjB,EAAImB,EAAAA,SAAWtC,EACfoB,EAAMD,GAAKA,EAAInB,GACfqB,EAAMrB,EAAUoB,EAChBD,EAAImB,EAAAA,SAAW,CAACF,EAChBd,EAAMH,GAAKA,GAAI,CAACiB,GAChBb,EAAM,CAACa,EAAMd,EACbM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKK,EAAM,CAAChC,EACZiB,EAAImB,EAAAA,SAAWJ,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAImB,EAAAA,SAAW,CAACpC,EAChBoB,EAAMH,GAAKA,GAAI,CAACjB,GAChBqB,EAAM,CAACrB,EAAUoB,EACjBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbrE,EAAE,CAAC,EAAIqE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbnE,EAAE,CAAC,EAAImE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACblE,EAAE,CAAC,EAAIkE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC3D,EAAE,CAAC,EAAIwE,EACPpB,EAASzB,EAAAA,IAAI,EAAG5B,EAAG,EAAGC,EAAGQ,EAAG,EAC5B4D,EAAK5B,EAAUI,EACfgB,EAAImB,EAAAA,SAAWvC,EACfqB,EAAMD,GAAKA,EAAIpB,GACfsB,EAAMtB,EAAUqB,EAChBD,EAAImB,EAAAA,SAAWnC,EACfmB,EAAMH,GAAKA,EAAIhB,GACfoB,EAAMpB,EAAUmB,EAChBM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAK7B,EAAUE,EACfiB,EAAImB,EAAAA,SAAWtC,EACfoB,EAAMD,GAAKA,EAAInB,GACfqB,EAAMrB,EAAUoB,EAChBD,EAAImB,EAAAA,SAAWpC,EACfoB,EAAMH,GAAKA,EAAIjB,GACfqB,EAAMrB,EAAUoB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbtD,GAAK,CAAC,EAAI0D,GAAMJ,EAAKN,IAAUA,EAAQY,GACvCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbtD,GAAK,CAAC,EAAIwD,GAAMF,EAAKN,IAAUA,EAAQW,GACvCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbvD,GAAK,CAAC,EAAIuD,GAAMM,EAAKb,IAAUM,EAAKN,GACpChD,GAAK,CAAC,EAAI6D,EACVjB,GAAU,IAEV/C,GAAI,CAAC,EAAI,EACT4C,EAAS,EACTzC,GAAK,CAAC,EAAI,EACV4C,GAAU,GAEVhB,IAAY,EAAG,CACf,MAAM+C,EAAMN,EAAAA,MAAM5B,EAAQ5C,GAAK+B,EAASvB,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAMnC,GAAU5C,GAAOsC,EAASzB,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAIyD,EAAKxD,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAMzB,GAAS5C,GAAM4B,EAAS1B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAI4D,EAAK3D,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI0B,EAASxB,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAMuB,EAAStB,CAAG,EAAGA,EAAKC,GAAME,CAAA,EAAMA,CAAG,EAEpDuB,IAAY,IACZnB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGlF,EAAIyC,EAAS1B,CAAE,EAAGA,EAAI8B,EAAS7B,CAAG,EAAGA,CAAG,GAE/E8B,IAAY,IACZpB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGnF,EAAI,CAAC0C,EAAS1B,CAAE,EAAGA,EAAI+B,EAAS9B,CAAG,EAAGA,CAAG,EAExF,CACA,GAAI4B,IAAY,EAAG,CACf,MAAM4C,EAAMN,EAAAA,MAAM5B,EAAQ5C,GAAKkC,EAAS1B,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAMlC,GAAU5C,GAAOwC,EAAS5B,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAI4D,EAAK3D,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAMzB,GAAS5C,GAAM+B,EAAS7B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAI+D,EAAK9D,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI6B,EAAS3B,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAM0B,EAASzB,CAAG,EAAGA,EAAKC,GAAME,CAAA,EAAMA,CAAG,CAC5D,CACJ,CACA,GAAIoB,IAAY,GAAKG,IAAY,EAAG,CAmGhC,GAlGIF,IAAY,GAAKG,IAAY,GAAKL,IAAY,GAAKG,IAAY,GAC/D0B,EAAK3B,EAAUmC,EACfhB,EAAImB,EAAAA,SAAWtC,EACfoB,EAAMD,GAAKA,EAAInB,GACfqB,EAAMrB,EAAUoB,EAChBD,EAAImB,EAAAA,SAAWH,EACfb,EAAMH,GAAKA,EAAIgB,GACfZ,EAAMY,EAAMb,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKK,EAAMjC,EACXkB,EAAImB,EAAAA,SAAWJ,EACfd,EAAMD,GAAKA,EAAIe,GACfb,EAAMa,EAAMd,EACZD,EAAImB,EAAAA,SAAWrC,EACfqB,EAAMH,GAAKA,EAAIlB,GACfsB,EAAMtB,EAAUqB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbtE,EAAE,CAAC,EAAIsE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbpE,EAAE,CAAC,EAAIoE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbnE,EAAE,CAAC,EAAImE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC5D,EAAE,CAAC,EAAIyE,EACPf,EAAK,CAACqB,EACNpB,EAAK,CAACd,EACNwB,EAAK7B,EAAUkB,EACfG,EAAImB,EAAAA,SAAWxC,EACfsB,EAAMD,GAAKA,EAAIrB,GACfuB,EAAMvB,EAAUsB,EAChBD,EAAImB,EAAAA,SAAWtB,EACfM,EAAMH,GAAKA,EAAIH,GACfO,EAAMP,EAAKM,EACXM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKG,EAAMf,EACXE,EAAImB,EAAAA,SAAWN,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAImB,EAAAA,SAAWrB,EACfK,EAAMH,GAAKA,EAAIF,GACfM,EAAMN,EAAKK,EACXQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbrE,EAAE,CAAC,EAAIqE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbnE,EAAE,CAAC,EAAImE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACblE,EAAE,CAAC,EAAIkE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC3D,EAAE,CAAC,EAAIwE,EACPnB,GAAS1B,EAAAA,IAAI,EAAG5B,EAAG,EAAGC,EAAGS,EAAG,EAC5B2D,EAAK3B,EAAUC,EACfkB,EAAImB,EAAAA,SAAWtC,EACfoB,EAAMD,GAAKA,EAAInB,GACfqB,EAAMrB,EAAUoB,EAChBD,EAAImB,EAAAA,SAAWrC,EACfqB,EAAMH,GAAKA,EAAIlB,GACfsB,EAAMtB,EAAUqB,EAChBM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAK/B,EAAUK,EACfgB,EAAImB,EAAAA,SAAWxC,EACfsB,EAAMD,GAAKA,EAAIrB,GACfuB,EAAMvB,EAAUsB,EAChBD,EAAImB,EAAAA,SAAWnC,EACfmB,EAAMH,GAAKA,EAAIhB,GACfoB,EAAMpB,EAAUmB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbrD,GAAK,CAAC,EAAIyD,GAAMJ,EAAKN,IAAUA,EAAQY,GACvCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbrD,GAAK,CAAC,EAAIuD,GAAMF,EAAKN,IAAUA,EAAQW,GACvCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbtD,GAAK,CAAC,EAAIsD,GAAMM,EAAKb,IAAUM,EAAKN,GACpC/C,GAAK,CAAC,EAAI4D,EACVhB,GAAU,IAEV/C,GAAI,CAAC,EAAI,EACT4C,GAAS,EACTzC,GAAK,CAAC,EAAI,EACV4C,GAAU,GAEVhB,IAAY,EAAG,CACf,MAAM8C,EAAMN,EAAAA,MAAM3B,GAAQ5C,GAAK+B,EAASxB,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAMjC,GAAU5C,GAAOqC,EAAS1B,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAI0D,EAAKzD,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAMxB,GAAS5C,GAAM4B,EAAS3B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAI6D,EAAK5D,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI2B,EAASzB,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAMwB,EAASvB,CAAG,EAAGA,EAAKC,GAAME,CAAA,EAAMA,CAAG,EAEpDwB,IAAY,IACZpB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGpF,EAAI4C,EAAS3B,CAAE,EAAGA,EAAI+B,EAAS9B,CAAG,EAAGA,CAAG,GAE/E4B,IAAY,IACZlB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGlF,EAAI,CAAC0C,EAAS3B,CAAE,EAAGA,EAAI6B,EAAS5B,CAAG,EAAGA,CAAG,EAExF,CACA,GAAI6B,IAAY,EAAG,CACf,MAAM2C,EAAMN,EAAAA,MAAM3B,GAAQ5C,GAAKkC,EAAS3B,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAMhC,GAAU5C,GAAOuC,EAAS7B,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAI6D,EAAK5D,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAMxB,GAAS5C,GAAM+B,EAAS9B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAIgE,EAAK/D,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI8B,EAAS5B,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAM2B,EAAS1B,CAAG,EAAGA,EAAMC,GAAME,CAAA,EAAMA,CAAG,CAC7D,CACJ,CACA,GAAIqB,IAAY,GAAKG,IAAY,EAAG,CAmGhC,GAlGIL,IAAY,GAAKG,IAAY,GAAKF,IAAY,GAAKG,IAAY,GAC/DyB,EAAK7B,EAAUsC,EACfjB,EAAImB,EAAAA,SAAWxC,EACfsB,EAAMD,GAAKA,EAAIrB,GACfuB,EAAMvB,EAAUsB,EAChBD,EAAImB,EAAAA,SAAWF,EACfd,EAAMH,GAAKA,EAAIiB,GACfb,EAAMa,EAAMd,EACZM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKG,EAAM9B,EACXiB,EAAImB,EAAAA,SAAWN,EACfZ,EAAMD,GAAKA,EAAIa,GACfX,EAAMW,EAAMZ,EACZD,EAAImB,EAAAA,SAAWpC,EACfoB,EAAMH,GAAKA,EAAIjB,GACfqB,EAAMrB,EAAUoB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbtE,EAAE,CAAC,EAAIsE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbpE,EAAE,CAAC,EAAIoE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbnE,EAAE,CAAC,EAAImE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC5D,EAAE,CAAC,EAAIyE,EACPf,EAAK,CAACmB,EACNlB,EAAK,CAAChB,EACN0B,EAAK5B,EAAUiB,EACfG,EAAImB,EAAAA,SAAWvC,EACfqB,EAAMD,GAAKA,EAAIpB,GACfsB,EAAMtB,EAAUqB,EAChBD,EAAImB,EAAAA,SAAWtB,EACfM,EAAMH,GAAKA,EAAIH,GACfO,EAAMP,EAAKM,EACXM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAKI,EAAMhB,EACXE,EAAImB,EAAAA,SAAWL,EACfb,EAAMD,GAAKA,EAAIc,GACfZ,EAAMY,EAAMb,EACZD,EAAImB,EAAAA,SAAWrB,EACfK,EAAMH,GAAKA,EAAIF,GACfM,EAAMN,EAAKK,EACXQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQM,EAAKI,EACbrE,EAAE,CAAC,EAAIqE,GAAMJ,EAAKN,IAAUY,EAAKZ,GACjCO,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQM,EAAKE,EACbnE,EAAE,CAAC,EAAImE,GAAMF,EAAKN,IAAUW,EAAKX,GACjCa,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACblE,EAAE,CAAC,EAAIkE,GAAMM,EAAKb,IAAUM,EAAKN,GACjC3D,EAAE,CAAC,EAAIwE,EACPrB,GAASxB,EAAAA,IAAI,EAAG5B,EAAG,EAAGC,EAAGO,EAAG,EAC5B6D,EAAK7B,EAAUI,EACfiB,EAAImB,EAAAA,SAAWxC,EACfsB,EAAMD,GAAKA,EAAIrB,GACfuB,EAAMvB,EAAUsB,EAChBD,EAAImB,EAAAA,SAAWpC,EACfoB,EAAMH,GAAKA,EAAIjB,GACfqB,EAAMrB,EAAUoB,EAChBM,EAAKP,EAAME,GAAOI,EAAKP,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDM,EAAK9B,EAAUE,EACfkB,EAAImB,EAAAA,SAAWvC,EACfqB,EAAMD,GAAKA,EAAIpB,GACfsB,EAAMtB,EAAUqB,EAChBD,EAAImB,EAAAA,SAAWrC,EACfqB,EAAMH,GAAKA,EAAIlB,GACfsB,EAAMtB,EAAUqB,EAChBQ,EAAKT,EAAME,GAAOM,EAAKT,EAAME,EAAMD,EAAMC,EAAMF,EAAMG,GACrDC,EAAKI,EAAKE,EACVZ,EAAQU,EAAKJ,EACbvD,GAAK,CAAC,EAAI2D,GAAMJ,EAAKN,IAAUA,EAAQY,GACvCL,EAAKE,EAAKH,EACVN,EAAQO,EAAKE,EACbD,EAAKC,GAAMF,EAAKP,IAAUM,EAAKN,GAC/BM,EAAKE,EAAKG,EACVX,EAAQQ,EAAKF,EACbvD,GAAK,CAAC,EAAIyD,GAAMF,EAAKN,IAAUA,EAAQW,GACvCE,EAAKN,EAAKD,EACVN,EAAQa,EAAKN,EACbxD,GAAK,CAAC,EAAIwD,GAAMM,EAAKb,IAAUM,EAAKN,GACpCjD,GAAK,CAAC,EAAI8D,EACVlB,GAAU,IAEV/C,GAAI,CAAC,EAAI,EACT4C,GAAS,EACTzC,GAAK,CAAC,EAAI,EACV4C,GAAU,GAEVb,IAAY,EAAG,CACf,MAAM6C,EAAMN,EAAAA,MAAM7B,GAAQ5C,GAAKkC,EAASzB,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAM/B,GAAU5C,GAAOoC,EAAS3B,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAI2D,EAAK1D,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAM1B,GAAS5C,GAAM+B,EAAS5B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAI8D,EAAK7D,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI4B,EAAS1B,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAMyB,EAASxB,CAAG,EAAGA,EAAKC,GAAME,CAAA,EAAMA,CAAG,EAEpDsB,IAAY,IACZlB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGnF,EAAI4C,EAAS5B,CAAE,EAAGA,EAAI6B,EAAS5B,CAAG,EAAGA,CAAG,GAE/E6B,IAAY,IACZnB,EAASD,EAAOC,EAAQwD,EAAAA,MAAMA,EAAAA,MAAM,EAAGpF,EAAI,CAAC6C,EAAS5B,CAAE,EAAGA,EAAI8B,EAAS7B,CAAG,EAAGA,CAAG,EAExF,CACA,GAAI8B,IAAY,EAAG,CACf,MAAM0C,EAAMN,EAAAA,MAAM7B,GAAQ5C,GAAKqC,EAAS5B,CAAI,EAC5CQ,EAASD,EAAOC,EAAQG,EAAAA,IACpBqD,EAAAA,MAAM9B,GAAU5C,GAAOsC,EAAS9B,CAAG,EAAGA,EACtCkE,EAAAA,MAAMM,EAAKtE,EAAM,EAAI8D,EAAK7D,CAAG,EAAGA,EAAKE,CAAA,EAAMA,CAAG,EAElD,MAAMoE,EAAOP,EAAAA,MAAM1B,GAAS5C,GAAMkC,EAAS/B,CAAE,EAC7CW,EAASD,EAAOC,EAAQ6D,EAAAA,UACpBL,EAAAA,MAAMO,EAAM1E,EAAI,EAAIiE,EAAKhE,CAAG,EAAGA,EAC/BkE,EAAAA,MAAMO,EAAM1E,EAAI+B,EAAS7B,CAAI,EAAGA,EAChCiE,EAAAA,MAAMM,EAAKtE,EAAM4B,EAAS3B,CAAG,EAAGA,EAAKC,GAAME,CAAA,EAAMA,CAAG,CAC5D,CACJ,CAEA,OAAOC,GAAIG,EAAS,CAAC,CACzB,CAEO,SAASgE,GAAS1D,EAAIC,EAAIC,EAAIC,EAAIC,GAAIC,GAAIC,EAAIC,EAAI,CACrD,MAAMoC,GAAM3C,EAAKM,EACXsC,EAAM1C,EAAKI,EACXuC,EAAMzC,GAAKE,EACXwC,EAAM7C,EAAKM,EACXwC,EAAM5C,EAAKI,EACXyC,EAAM3C,GAAKE,EAEXoD,EAASf,EAAMI,EACfY,EAASf,EAAME,EACfc,GAAQlB,GAAMA,GAAMG,EAAMA,EAE1BgB,GAASjB,EAAMC,EACfiB,GAASpB,GAAMK,EACfgB,GAAQpB,EAAMA,EAAMG,EAAMA,EAE1BkB,GAAStB,GAAMI,EACfmB,GAAStB,EAAME,EACfqB,GAAQtB,EAAMA,EAAMG,EAAMA,EAE1BG,EACFU,IAASF,EAASC,GAClBI,IAASF,GAASC,IAClBI,IAASF,GAASC,IAEhB1D,IACD,KAAK,IAAImD,CAAM,EAAI,KAAK,IAAIC,CAAM,GAAKC,IACvC,KAAK,IAAIC,EAAM,EAAI,KAAK,IAAIC,EAAM,GAAKC,IACvC,KAAK,IAAIC,EAAM,EAAI,KAAK,IAAIC,EAAM,GAAKC,GAEtCd,GAAW/F,GAAekD,GAEhC,OAAI2C,EAAME,IAAY,CAACF,EAAME,GAClBF,EAEJpD,GAAcC,EAAIC,EAAIC,EAAIC,EAAIC,GAAIC,GAAIC,EAAIC,EAAIC,EAAS,CAClE","x_google_ignoreList":[0]}

package/dist/index.cjs561.js

@@ -0,0 +1,9 @@
+"use strict";Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const q=require("./index.cjs562.js");var It;function Te(){if(It)return q.__exports;It=1,Object.defineProperty(q.__exports,"__esModule",{value:!0});const Nt=Object.prototype.toString;function B(i){const e=Nt.call(i);return e.endsWith("Array]")&&!e.includes("Big")}function Tt(i,e={}){if(!B(i))throw new TypeError("input must be an array");if(i.length===0)throw new TypeError("input must not be empty");const{fromIndex:s=0,toIndex:t=i.length}=e;if(s<0||s>=i.length||!Number.isInteger(s))throw new Error("fromIndex must be a positive integer smaller than length");if(t<=s||t>i.length||!Number.isInteger(t))throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");let o=i[s];for(let r=s+1;r<t;r++)i[r]>o&&(o=i[r]);return o}function qt(i,e={}){if(!B(i))throw new TypeError("input must be an array");if(i.length===0)throw new TypeError("input must not be empty");const{fromIndex:s=0,toIndex:t=i.length}=e;if(s<0||s>=i.length||!Number.isInteger(s))throw new Error("fromIndex must be a positive integer smaller than length");if(t<=s||t>i.length||!Number.isInteger(t))throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");let o=i[s];for(let r=s+1;r<t;r++)i[r]<o&&(o=i[r]);return o}function dt(i,e={}){if(B(i)){if(i.length===0)throw new TypeError("input must not be empty")}else throw new TypeError("input must be an array");let s;if(e.output!==void 0){if(!B(e.output))throw new TypeError("output option must be an array if specified");s=e.output}else s=new Array(i.length);const t=qt(i),o=Tt(i);if(t===o)throw new RangeError("minimum and maximum input values are equal. Cannot rescale a constant array");const{min:r=e.autoMinMax?t:0,max:l=e.autoMinMax?o:1}=e;if(r>=l)throw new RangeError("min option must be smaller than max option");const n=(l-r)/(o-t);for(let h=0;h<i.length;h++)s[h]=(i[h]-t)*n+r;return s}const rt=" ".repeat(2),yt=" ".repeat(4);function _t(){return Mt(this)}function Mt(i,e={}){const{maxRows:s=15,maxColumns:t=10,maxNumSize:o=8,padMinus:r="auto"}=e;return`${i.constructor.name} {
+${rt}[
+${yt}${Dt(i,s,t,o,r)}
+${rt}]
+${rt}rows: ${i.rows}
+${rt}columns: ${i.columns}
+}`}function Dt(i,e,s,t,o){const{rows:r,columns:l}=i,n=Math.min(r,e),h=Math.min(l,s),u=[];if(o==="auto"){o=!1;t:for(let w=0;w<n;w++)for(let f=0;f<h;f++)if(i.get(w,f)<0){o=!0;break t}}for(let w=0;w<n;w++){let f=[];for(let a=0;a<h;a++)f.push(Ft(i.get(w,a),t,o));u.push(`${f.join(" ")}`)}return h!==l&&(u[u.length-1]+=` ... ${l-s} more columns`),n!==r&&u.push(`... ${r-e} more rows`),u.join(`
+${yt}`)}function Ft(i,e,s){return(i>=0&&s?` ${jt(i,e-1)}`:jt(i,e)).padEnd(e)}function jt(i,e){let s=i.toString();if(s.length<=e)return s;let t=i.toFixed(e);if(t.length>e&&(t=i.toFixed(Math.max(0,e-(t.length-e)))),t.length<=e&&!t.startsWith("0.000")&&!t.startsWith("-0.000"))return t;let o=i.toExponential(e);return o.length>e&&(o=i.toExponential(Math.max(0,e-(o.length-e)))),o.slice(0)}function zt(i,e){i.prototype.add=function(t){return typeof t=="number"?this.addS(t):this.addM(t)},i.prototype.addS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)+t);return this},i.prototype.addM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)+t.get(o,r));return this},i.add=function(t,o){return new e(t).add(o)},i.prototype.sub=function(t){return typeof t=="number"?this.subS(t):this.subM(t)},i.prototype.subS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)-t);return this},i.prototype.subM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)-t.get(o,r));return this},i.sub=function(t,o){return new e(t).sub(o)},i.prototype.subtract=i.prototype.sub,i.prototype.subtractS=i.prototype.subS,i.prototype.subtractM=i.prototype.subM,i.subtract=i.sub,i.prototype.mul=function(t){return typeof t=="number"?this.mulS(t):this.mulM(t)},i.prototype.mulS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)*t);return this},i.prototype.mulM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)*t.get(o,r));return this},i.mul=function(t,o){return new e(t).mul(o)},i.prototype.multiply=i.prototype.mul,i.prototype.multiplyS=i.prototype.mulS,i.prototype.multiplyM=i.prototype.mulM,i.multiply=i.mul,i.prototype.div=function(t){return typeof t=="number"?this.divS(t):this.divM(t)},i.prototype.divS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)/t);return this},i.prototype.divM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)/t.get(o,r));return this},i.div=function(t,o){return new e(t).div(o)},i.prototype.divide=i.prototype.div,i.prototype.divideS=i.prototype.divS,i.prototype.divideM=i.prototype.divM,i.divide=i.div,i.prototype.mod=function(t){return typeof t=="number"?this.modS(t):this.modM(t)},i.prototype.modS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)%t);return this},i.prototype.modM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)%t.get(o,r));return this},i.mod=function(t,o){return new e(t).mod(o)},i.prototype.modulus=i.prototype.mod,i.prototype.modulusS=i.prototype.modS,i.prototype.modulusM=i.prototype.modM,i.modulus=i.mod,i.prototype.and=function(t){return typeof t=="number"?this.andS(t):this.andM(t)},i.prototype.andS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)&t);return this},i.prototype.andM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)&t.get(o,r));return this},i.and=function(t,o){return new e(t).and(o)},i.prototype.or=function(t){return typeof t=="number"?this.orS(t):this.orM(t)},i.prototype.orS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)|t);return this},i.prototype.orM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)|t.get(o,r));return this},i.or=function(t,o){return new e(t).or(o)},i.prototype.xor=function(t){return typeof t=="number"?this.xorS(t):this.xorM(t)},i.prototype.xorS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)^t);return this},i.prototype.xorM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)^t.get(o,r));return this},i.xor=function(t,o){return new e(t).xor(o)},i.prototype.leftShift=function(t){return typeof t=="number"?this.leftShiftS(t):this.leftShiftM(t)},i.prototype.leftShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)<<t);return this},i.prototype.leftShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)<<t.get(o,r));return this},i.leftShift=function(t,o){return new e(t).leftShift(o)},i.prototype.signPropagatingRightShift=function(t){return typeof t=="number"?this.signPropagatingRightShiftS(t):this.signPropagatingRightShiftM(t)},i.prototype.signPropagatingRightShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>t);return this},i.prototype.signPropagatingRightShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>t.get(o,r));return this},i.signPropagatingRightShift=function(t,o){return new e(t).signPropagatingRightShift(o)},i.prototype.rightShift=function(t){return typeof t=="number"?this.rightShiftS(t):this.rightShiftM(t)},i.prototype.rightShiftS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>>t);return this},i.prototype.rightShiftM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)>>>t.get(o,r));return this},i.rightShift=function(t,o){return new e(t).rightShift(o)},i.prototype.zeroFillRightShift=i.prototype.rightShift,i.prototype.zeroFillRightShiftS=i.prototype.rightShiftS,i.prototype.zeroFillRightShiftM=i.prototype.rightShiftM,i.zeroFillRightShift=i.rightShift,i.prototype.not=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,~this.get(t,o));return this},i.not=function(t){return new e(t).not()},i.prototype.abs=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.abs(this.get(t,o)));return this},i.abs=function(t){return new e(t).abs()},i.prototype.acos=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.acos(this.get(t,o)));return this},i.acos=function(t){return new e(t).acos()},i.prototype.acosh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.acosh(this.get(t,o)));return this},i.acosh=function(t){return new e(t).acosh()},i.prototype.asin=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.asin(this.get(t,o)));return this},i.asin=function(t){return new e(t).asin()},i.prototype.asinh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.asinh(this.get(t,o)));return this},i.asinh=function(t){return new e(t).asinh()},i.prototype.atan=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.atan(this.get(t,o)));return this},i.atan=function(t){return new e(t).atan()},i.prototype.atanh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.atanh(this.get(t,o)));return this},i.atanh=function(t){return new e(t).atanh()},i.prototype.cbrt=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cbrt(this.get(t,o)));return this},i.cbrt=function(t){return new e(t).cbrt()},i.prototype.ceil=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.ceil(this.get(t,o)));return this},i.ceil=function(t){return new e(t).ceil()},i.prototype.clz32=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.clz32(this.get(t,o)));return this},i.clz32=function(t){return new e(t).clz32()},i.prototype.cos=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cos(this.get(t,o)));return this},i.cos=function(t){return new e(t).cos()},i.prototype.cosh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.cosh(this.get(t,o)));return this},i.cosh=function(t){return new e(t).cosh()},i.prototype.exp=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.exp(this.get(t,o)));return this},i.exp=function(t){return new e(t).exp()},i.prototype.expm1=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.expm1(this.get(t,o)));return this},i.expm1=function(t){return new e(t).expm1()},i.prototype.floor=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.floor(this.get(t,o)));return this},i.floor=function(t){return new e(t).floor()},i.prototype.fround=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.fround(this.get(t,o)));return this},i.fround=function(t){return new e(t).fround()},i.prototype.log=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log(this.get(t,o)));return this},i.log=function(t){return new e(t).log()},i.prototype.log1p=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log1p(this.get(t,o)));return this},i.log1p=function(t){return new e(t).log1p()},i.prototype.log10=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log10(this.get(t,o)));return this},i.log10=function(t){return new e(t).log10()},i.prototype.log2=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.log2(this.get(t,o)));return this},i.log2=function(t){return new e(t).log2()},i.prototype.round=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.round(this.get(t,o)));return this},i.round=function(t){return new e(t).round()},i.prototype.sign=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sign(this.get(t,o)));return this},i.sign=function(t){return new e(t).sign()},i.prototype.sin=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sin(this.get(t,o)));return this},i.sin=function(t){return new e(t).sin()},i.prototype.sinh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sinh(this.get(t,o)));return this},i.sinh=function(t){return new e(t).sinh()},i.prototype.sqrt=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.sqrt(this.get(t,o)));return this},i.sqrt=function(t){return new e(t).sqrt()},i.prototype.tan=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.tan(this.get(t,o)));return this},i.tan=function(t){return new e(t).tan()},i.prototype.tanh=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.tanh(this.get(t,o)));return this},i.tanh=function(t){return new e(t).tanh()},i.prototype.trunc=function(){for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.set(t,o,Math.trunc(this.get(t,o)));return this},i.trunc=function(t){return new e(t).trunc()},i.pow=function(t,o){return new e(t).pow(o)},i.prototype.pow=function(t){return typeof t=="number"?this.powS(t):this.powM(t)},i.prototype.powS=function(t){for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)**t);return this},i.prototype.powM=function(t){if(t=e.checkMatrix(t),this.rows!==t.rows||this.columns!==t.columns)throw new RangeError("Matrices dimensions must be equal");for(let o=0;o<this.rows;o++)for(let r=0;r<this.columns;r++)this.set(o,r,this.get(o,r)**t.get(o,r));return this}}function J(i,e,s){let t=s?i.rows:i.rows-1;if(e<0||e>t)throw new RangeError("Row index out of range")}function G(i,e,s){let t=s?i.columns:i.columns-1;if(e<0||e>t)throw new RangeError("Column index out of range")}function tt(i,e){if(e.to1DArray&&(e=e.to1DArray()),e.length!==i.columns)throw new RangeError("vector size must be the same as the number of columns");return e}function et(i,e){if(e.to1DArray&&(e=e.to1DArray()),e.length!==i.rows)throw new RangeError("vector size must be the same as the number of rows");return e}function gt(i,e){if(!B(e))throw new TypeError("row indices must be an array");for(let s=0;s<e.length;s++)if(e[s]<0||e[s]>=i.rows)throw new RangeError("row indices are out of range")}function wt(i,e){if(!B(e))throw new TypeError("column indices must be an array");for(let s=0;s<e.length;s++)if(e[s]<0||e[s]>=i.columns)throw new RangeError("column indices are out of range")}function mt(i,e,s,t,o){if(arguments.length!==5)throw new RangeError("expected 4 arguments");if(nt("startRow",e),nt("endRow",s),nt("startColumn",t),nt("endColumn",o),e>s||t>o||e<0||e>=i.rows||s<0||s>=i.rows||t<0||t>=i.columns||o<0||o>=i.columns)throw new RangeError("Submatrix indices are out of range")}function it(i,e=0){let s=[];for(let t=0;t<i;t++)s.push(e);return s}function nt(i,e){if(typeof e!="number")throw new TypeError(`${i} must be a number`)}function st(i){if(i.isEmpty())throw new Error("Empty matrix has no elements to index")}function Vt(i){let e=it(i.rows);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[s]+=i.get(s,t);return e}function Ct(i){let e=it(i.columns);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[t]+=i.get(s,t);return e}function Pt(i){let e=0;for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)e+=i.get(s,t);return e}function $t(i){let e=it(i.rows,1);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[s]*=i.get(s,t);return e}function Lt(i){let e=it(i.columns,1);for(let s=0;s<i.rows;++s)for(let t=0;t<i.columns;++t)e[t]*=i.get(s,t);return e}function Bt(i){let e=1;for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)e*=i.get(s,t);return e}function Ut(i,e,s){const t=i.rows,o=i.columns,r=[];for(let l=0;l<t;l++){let n=0,h=0,u=0;for(let w=0;w<o;w++)u=i.get(l,w)-s[l],n+=u,h+=u*u;e?r.push((h-n*n/o)/(o-1)):r.push((h-n*n/o)/o)}return r}function Ot(i,e,s){const t=i.rows,o=i.columns,r=[];for(let l=0;l<o;l++){let n=0,h=0,u=0;for(let w=0;w<t;w++)u=i.get(w,l)-s[l],n+=u,h+=u*u;e?r.push((h-n*n/t)/(t-1)):r.push((h-n*n/t)/t)}return r}function Qt(i,e,s){const t=i.rows,o=i.columns,r=t*o;let l=0,n=0,h=0;for(let u=0;u<t;u++)for(let w=0;w<o;w++)h=i.get(u,w)-s,l+=h,n+=h*h;return e?(n-l*l/r)/(r-1):(n-l*l/r)/r}function Wt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e[s])}function Yt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e[t])}function Jt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)-e)}function Gt(i){const e=[];for(let s=0;s<i.rows;s++){let t=0;for(let o=0;o<i.columns;o++)t+=i.get(s,o)**2/(i.columns-1);e.push(Math.sqrt(t))}return e}function Kt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e[s])}function Xt(i){const e=[];for(let s=0;s<i.columns;s++){let t=0;for(let o=0;o<i.rows;o++)t+=i.get(o,s)**2/(i.rows-1);e.push(Math.sqrt(t))}return e}function Zt(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e[t])}function xt(i){const e=i.size-1;let s=0;for(let t=0;t<i.columns;t++)for(let o=0;o<i.rows;o++)s+=i.get(o,t)**2/e;return Math.sqrt(s)}function At(i,e){for(let s=0;s<i.rows;s++)for(let t=0;t<i.columns;t++)i.set(s,t,i.get(s,t)/e)}class T{static from1DArray(e,s,t){if(e*s!==t.length)throw new RangeError("data length does not match given dimensions");let r=new S(e,s);for(let l=0;l<e;l++)for(let n=0;n<s;n++)r.set(l,n,t[l*s+n]);return r}static rowVector(e){let s=new S(1,e.length);for(let t=0;t<e.length;t++)s.set(0,t,e[t]);return s}static columnVector(e){let s=new S(e.length,1);for(let t=0;t<e.length;t++)s.set(t,0,e[t]);return s}static zeros(e,s){return new S(e,s)}static ones(e,s){return new S(e,s).fill(1)}static rand(e,s,t={}){if(typeof t!="object")throw new TypeError("options must be an object");const{random:o=Math.random}=t;let r=new S(e,s);for(let l=0;l<e;l++)for(let n=0;n<s;n++)r.set(l,n,o());return r}static randInt(e,s,t={}){if(typeof t!="object")throw new TypeError("options must be an object");const{min:o=0,max:r=1e3,random:l=Math.random}=t;if(!Number.isInteger(o))throw new TypeError("min must be an integer");if(!Number.isInteger(r))throw new TypeError("max must be an integer");if(o>=r)throw new RangeError("min must be smaller than max");let n=r-o,h=new S(e,s);for(let u=0;u<e;u++)for(let w=0;w<s;w++){let f=o+Math.round(l()*n);h.set(u,w,f)}return h}static eye(e,s,t){s===void 0&&(s=e),t===void 0&&(t=1);let o=Math.min(e,s),r=this.zeros(e,s);for(let l=0;l<o;l++)r.set(l,l,t);return r}static diag(e,s,t){let o=e.length;s===void 0&&(s=o),t===void 0&&(t=s);let r=Math.min(o,s,t),l=this.zeros(s,t);for(let n=0;n<r;n++)l.set(n,n,e[n]);return l}static min(e,s){e=this.checkMatrix(e),s=this.checkMatrix(s);let t=e.rows,o=e.columns,r=new S(t,o);for(let l=0;l<t;l++)for(let n=0;n<o;n++)r.set(l,n,Math.min(e.get(l,n),s.get(l,n)));return r}static max(e,s){e=this.checkMatrix(e),s=this.checkMatrix(s);let t=e.rows,o=e.columns,r=new this(t,o);for(let l=0;l<t;l++)for(let n=0;n<o;n++)r.set(l,n,Math.max(e.get(l,n),s.get(l,n)));return r}static checkMatrix(e){return T.isMatrix(e)?e:new S(e)}static isMatrix(e){return e!=null&&e.klass==="Matrix"}get size(){return this.rows*this.columns}apply(e){if(typeof e!="function")throw new TypeError("callback must be a function");for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.call(this,s,t);return this}to1DArray(){let e=[];for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.push(this.get(s,t));return e}to2DArray(){let e=[];for(let s=0;s<this.rows;s++){e.push([]);for(let t=0;t<this.columns;t++)e[s].push(this.get(s,t))}return e}toJSON(){return this.to2DArray()}isRowVector(){return this.rows===1}isColumnVector(){return this.columns===1}isVector(){return this.rows===1||this.columns===1}isSquare(){return this.rows===this.columns}isEmpty(){return this.rows===0||this.columns===0}isSymmetric(){if(this.isSquare()){for(let e=0;e<this.rows;e++)for(let s=0;s<=e;s++)if(this.get(e,s)!==this.get(s,e))return!1;return!0}return!1}isDistance(){if(!this.isSymmetric())return!1;for(let e=0;e<this.rows;e++)if(this.get(e,e)!==0)return!1;return!0}isEchelonForm(){let e=0,s=0,t=-1,o=!0,r=!1;for(;e<this.rows&&o;){for(s=0,r=!1;s<this.columns&&r===!1;)this.get(e,s)===0?s++:this.get(e,s)===1&&s>t?(r=!0,t=s):(o=!1,r=!0);e++}return o}isReducedEchelonForm(){let e=0,s=0,t=-1,o=!0,r=!1;for(;e<this.rows&&o;){for(s=0,r=!1;s<this.columns&&r===!1;)this.get(e,s)===0?s++:this.get(e,s)===1&&s>t?(r=!0,t=s):(o=!1,r=!0);for(let l=s+1;l<this.rows;l++)this.get(e,l)!==0&&(o=!1);e++}return o}echelonForm(){let e=this.clone(),s=0,t=0;for(;s<e.rows&&t<e.columns;){let o=s;for(let r=s;r<e.rows;r++)e.get(r,t)>e.get(o,t)&&(o=r);if(e.get(o,t)===0)t++;else{e.swapRows(s,o);let r=e.get(s,t);for(let l=t;l<e.columns;l++)e.set(s,l,e.get(s,l)/r);for(let l=s+1;l<e.rows;l++){let n=e.get(l,t)/e.get(s,t);e.set(l,t,0);for(let h=t+1;h<e.columns;h++)e.set(l,h,e.get(l,h)-e.get(s,h)*n)}s++,t++}}return e}reducedEchelonForm(){let e=this.echelonForm(),s=e.columns,t=e.rows,o=t-1;for(;o>=0;)if(e.maxRow(o)===0)o--;else{let r=0,l=!1;for(;r<t&&l===!1;)e.get(o,r)===1?l=!0:r++;for(let n=0;n<o;n++){let h=e.get(n,r);for(let u=r;u<s;u++){let w=e.get(n,u)-h*e.get(o,u);e.set(n,u,w)}}o--}return e}set(){throw new Error("set method is unimplemented")}get(){throw new Error("get method is unimplemented")}repeat(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{rows:s=1,columns:t=1}=e;if(!Number.isInteger(s)||s<=0)throw new TypeError("rows must be a positive integer");if(!Number.isInteger(t)||t<=0)throw new TypeError("columns must be a positive integer");let o=new S(this.rows*s,this.columns*t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)o.setSubMatrix(this,this.rows*r,this.columns*l);return o}fill(e){for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,e);return this}neg(){return this.mulS(-1)}getRow(e){J(this,e);let s=[];for(let t=0;t<this.columns;t++)s.push(this.get(e,t));return s}getRowVector(e){return S.rowVector(this.getRow(e))}setRow(e,s){J(this,e),s=tt(this,s);for(let t=0;t<this.columns;t++)this.set(e,t,s[t]);return this}swapRows(e,s){J(this,e),J(this,s);for(let t=0;t<this.columns;t++){let o=this.get(e,t);this.set(e,t,this.get(s,t)),this.set(s,t,o)}return this}getColumn(e){G(this,e);let s=[];for(let t=0;t<this.rows;t++)s.push(this.get(t,e));return s}getColumnVector(e){return S.columnVector(this.getColumn(e))}setColumn(e,s){G(this,e),s=et(this,s);for(let t=0;t<this.rows;t++)this.set(t,e,s[t]);return this}swapColumns(e,s){G(this,e),G(this,s);for(let t=0;t<this.rows;t++){let o=this.get(t,e);this.set(t,e,this.get(t,s)),this.set(t,s,o)}return this}addRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)+e[t]);return this}subRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)-e[t]);return this}mulRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)*e[t]);return this}divRowVector(e){e=tt(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)/e[t]);return this}addColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)+e[s]);return this}subColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)-e[s]);return this}mulColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)*e[s]);return this}divColumnVector(e){e=et(this,e);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)this.set(s,t,this.get(s,t)/e[s]);return this}mulRow(e,s){J(this,e);for(let t=0;t<this.columns;t++)this.set(e,t,this.get(e,t)*s);return this}mulColumn(e,s){G(this,e);for(let t=0;t<this.rows;t++)this.set(t,e,this.get(t,e)*s);return this}max(e){if(this.isEmpty())return NaN;switch(e){case"row":{const s=new Array(this.rows).fill(Number.NEGATIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s[t]&&(s[t]=this.get(t,o));return s}case"column":{const s=new Array(this.columns).fill(Number.NEGATIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s[o]&&(s[o]=this.get(t,o));return s}case void 0:{let s=this.get(0,0);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>s&&(s=this.get(t,o));return s}default:throw new Error(`invalid option: ${e}`)}}maxIndex(){st(this);let e=this.get(0,0),s=[0,0];for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)>e&&(e=this.get(t,o),s[0]=t,s[1]=o);return s}min(e){if(this.isEmpty())return NaN;switch(e){case"row":{const s=new Array(this.rows).fill(Number.POSITIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s[t]&&(s[t]=this.get(t,o));return s}case"column":{const s=new Array(this.columns).fill(Number.POSITIVE_INFINITY);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s[o]&&(s[o]=this.get(t,o));return s}case void 0:{let s=this.get(0,0);for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<s&&(s=this.get(t,o));return s}default:throw new Error(`invalid option: ${e}`)}}minIndex(){st(this);let e=this.get(0,0),s=[0,0];for(let t=0;t<this.rows;t++)for(let o=0;o<this.columns;o++)this.get(t,o)<e&&(e=this.get(t,o),s[0]=t,s[1]=o);return s}maxRow(e){if(J(this,e),this.isEmpty())return NaN;let s=this.get(e,0);for(let t=1;t<this.columns;t++)this.get(e,t)>s&&(s=this.get(e,t));return s}maxRowIndex(e){J(this,e),st(this);let s=this.get(e,0),t=[e,0];for(let o=1;o<this.columns;o++)this.get(e,o)>s&&(s=this.get(e,o),t[1]=o);return t}minRow(e){if(J(this,e),this.isEmpty())return NaN;let s=this.get(e,0);for(let t=1;t<this.columns;t++)this.get(e,t)<s&&(s=this.get(e,t));return s}minRowIndex(e){J(this,e),st(this);let s=this.get(e,0),t=[e,0];for(let o=1;o<this.columns;o++)this.get(e,o)<s&&(s=this.get(e,o),t[1]=o);return t}maxColumn(e){if(G(this,e),this.isEmpty())return NaN;let s=this.get(0,e);for(let t=1;t<this.rows;t++)this.get(t,e)>s&&(s=this.get(t,e));return s}maxColumnIndex(e){G(this,e),st(this);let s=this.get(0,e),t=[0,e];for(let o=1;o<this.rows;o++)this.get(o,e)>s&&(s=this.get(o,e),t[0]=o);return t}minColumn(e){if(G(this,e),this.isEmpty())return NaN;let s=this.get(0,e);for(let t=1;t<this.rows;t++)this.get(t,e)<s&&(s=this.get(t,e));return s}minColumnIndex(e){G(this,e),st(this);let s=this.get(0,e),t=[0,e];for(let o=1;o<this.rows;o++)this.get(o,e)<s&&(s=this.get(o,e),t[0]=o);return t}diag(){let e=Math.min(this.rows,this.columns),s=[];for(let t=0;t<e;t++)s.push(this.get(t,t));return s}norm(e="frobenius"){switch(e){case"max":return this.max();case"frobenius":return Math.sqrt(this.dot(this));default:throw new RangeError(`unknown norm type: ${e}`)}}cumulativeSum(){let e=0;for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e+=this.get(s,t),this.set(s,t,e);return this}dot(e){T.isMatrix(e)&&(e=e.to1DArray());let s=this.to1DArray();if(s.length!==e.length)throw new RangeError("vectors do not have the same size");let t=0;for(let o=0;o<s.length;o++)t+=s[o]*e[o];return t}mmul(e){e=S.checkMatrix(e);let s=this.rows,t=this.columns,o=e.columns,r=new S(s,o),l=new Float64Array(t);for(let n=0;n<o;n++){for(let h=0;h<t;h++)l[h]=e.get(h,n);for(let h=0;h<s;h++){let u=0;for(let w=0;w<t;w++)u+=this.get(h,w)*l[w];r.set(h,n,u)}}return r}gram(){const e=this.rows,s=this.columns,t=new Float64Array(s*s);for(let r=0;r<e;r++)for(let l=0;l<s;l++){const n=this.get(r,l);if(n===0)continue;const h=l*s;for(let u=l;u<s;u++)t[h+u]+=n*this.get(r,u)}const o=new S(s,s);for(let r=0;r<s;r++){const l=r*s;for(let n=r;n<s;n++){const h=t[l+n];o.set(r,n,h),o.set(n,r,h)}}return o}mmulByTranspose(e){let s=this.rows,t=this.columns;if(e!==void 0&&e.length!==t)throw new RangeError("scale must have one value per column");let o=new S(s,s),r=new Float64Array(t);for(let l=0;l<s;l++){if(e===void 0)for(let n=0;n<t;n++)r[n]=this.get(l,n);else for(let n=0;n<t;n++)r[n]=e[n]*this.get(l,n);for(let n=l;n<s;n++){let h=0;for(let u=0;u<t;u++)h+=this.get(n,u)*r[u];o.set(n,l,h),o.set(l,n,h)}}return o}mpow(e){if(!this.isSquare())throw new RangeError("Matrix must be square");if(!Number.isInteger(e)||e<0)throw new RangeError("Exponent must be a non-negative integer");let s=S.eye(this.rows),t=this;for(let o=e;o>=1;o/=2)(o&1)!==0&&(s=s.mmul(t)),t=t.mmul(t);return s}strassen2x2(e){e=S.checkMatrix(e);let s=new S(2,2);const t=this.get(0,0),o=e.get(0,0),r=this.get(0,1),l=e.get(0,1),n=this.get(1,0),h=e.get(1,0),u=this.get(1,1),w=e.get(1,1),f=(t+u)*(o+w),a=(n+u)*o,j=t*(l-w),m=u*(h-o),p=(t+r)*w,E=(n-t)*(o+l),c=(r-u)*(h+w),M=f+m-p+c,k=j+p,I=a+m,N=f-a+j+E;return s.set(0,0,M),s.set(0,1,k),s.set(1,0,I),s.set(1,1,N),s}strassen3x3(e){e=S.checkMatrix(e);let s=new S(3,3);const t=this.get(0,0),o=this.get(0,1),r=this.get(0,2),l=this.get(1,0),n=this.get(1,1),h=this.get(1,2),u=this.get(2,0),w=this.get(2,1),f=this.get(2,2),a=e.get(0,0),j=e.get(0,1),m=e.get(0,2),p=e.get(1,0),E=e.get(1,1),c=e.get(1,2),M=e.get(2,0),k=e.get(2,1),I=e.get(2,2),N=(t+o+r-l-n-w-f)*E,V=(t-l)*(-j+E),v=n*(-a+j+p-E-c-M+I),R=(-t+l+n)*(a-j+E),F=(l+n)*(-a+j),g=t*a,d=(-t+u+w)*(a-m+c),b=(-t+u)*(m-c),y=(u+w)*(-a+m),C=(t+o+r-n-h-u-w)*c,D=w*(-a+m+p-E-c-M+k),z=(-r+w+f)*(E+M-k),P=(r-f)*(E-k),U=r*M,K=(w+f)*(-M+k),$=(-r+n+h)*(c+M-I),X=(r-h)*(c-I),A=(n+h)*(-M+I),_=o*p,O=h*k,W=l*m,Y=u*j,L=f*I,je=g+U+_,Se=N+R+F+g+z+U+K,Ee=g+d+y+C+U+$+A,be=V+v+R+g+U+$+X,ke=V+R+F+g+O,ve=U+$+X+A+W,Re=g+d+b+D+z+P+U,Ie=z+P+U+K+Y,Ne=g+d+b+y+L;return s.set(0,0,je),s.set(0,1,Se),s.set(0,2,Ee),s.set(1,0,be),s.set(1,1,ke),s.set(1,2,ve),s.set(2,0,Re),s.set(2,1,Ie),s.set(2,2,Ne),s}mmulStrassen(e){e=S.checkMatrix(e);let s=this.clone(),t=s.rows,o=s.columns,r=e.rows,l=e.columns;o!==r&&console.warn(`Multiplying ${t} x ${o} and ${r} x ${l} matrix: dimensions do not match.`);function n(f,a,j){let m=f.rows,p=f.columns;if(m===a&&p===j)return f;{let E=T.zeros(a,j);return E=E.setSubMatrix(f,0,0),E}}let h=Math.max(t,r),u=Math.max(o,l);s=n(s,h,u),e=n(e,h,u);function w(f,a,j,m){if(j<=512||m<=512)return f.mmul(a);j%2===1&&m%2===1?(f=n(f,j+1,m+1),a=n(a,j+1,m+1)):j%2===1?(f=n(f,j+1,m),a=n(a,j+1,m)):m%2===1&&(f=n(f,j,m+1),a=n(a,j,m+1));let p=parseInt(f.rows/2,10),E=parseInt(f.columns/2,10),c=f.subMatrix(0,p-1,0,E-1),M=a.subMatrix(0,p-1,0,E-1),k=f.subMatrix(0,p-1,E,f.columns-1),I=a.subMatrix(0,p-1,E,a.columns-1),N=f.subMatrix(p,f.rows-1,0,E-1),V=a.subMatrix(p,a.rows-1,0,E-1),v=f.subMatrix(p,f.rows-1,E,f.columns-1),R=a.subMatrix(p,a.rows-1,E,a.columns-1),F=w(T.add(c,v),T.add(M,R),p,E),g=w(T.add(N,v),M,p,E),d=w(c,T.sub(I,R),p,E),b=w(v,T.sub(V,M),p,E),y=w(T.add(c,k),R,p,E),C=w(T.sub(N,c),T.add(M,I),p,E),D=w(T.sub(k,v),T.add(V,R),p,E),z=T.add(F,b);z.sub(y),z.add(D);let P=T.add(d,y),U=T.add(g,b),K=T.sub(F,g);K.add(d),K.add(C);let $=T.zeros(2*z.rows,2*z.columns);return $=$.setSubMatrix(z,0,0),$=$.setSubMatrix(P,z.rows,0),$=$.setSubMatrix(U,0,z.columns),$=$.setSubMatrix(K,z.rows,z.columns),$.subMatrix(0,j-1,0,m-1)}return w(s,e,h,u)}scaleRows(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{min:s=0,max:t=1}=e;if(!Number.isFinite(s))throw new TypeError("min must be a number");if(!Number.isFinite(t))throw new TypeError("max must be a number");if(s>=t)throw new RangeError("min must be smaller than max");let o=new S(this.rows,this.columns);for(let r=0;r<this.rows;r++){const l=this.getRow(r);l.length>0&&dt(l,{min:s,max:t,output:l}),o.setRow(r,l)}return o}scaleColumns(e={}){if(typeof e!="object")throw new TypeError("options must be an object");const{min:s=0,max:t=1}=e;if(!Number.isFinite(s))throw new TypeError("min must be a number");if(!Number.isFinite(t))throw new TypeError("max must be a number");if(s>=t)throw new RangeError("min must be smaller than max");let o=new S(this.rows,this.columns);for(let r=0;r<this.columns;r++){const l=this.getColumn(r);l.length&&dt(l,{min:s,max:t,output:l}),o.setColumn(r,l)}return o}flipRows(){const e=Math.ceil(this.columns/2);for(let s=0;s<this.rows;s++)for(let t=0;t<e;t++){let o=this.get(s,t),r=this.get(s,this.columns-1-t);this.set(s,t,r),this.set(s,this.columns-1-t,o)}return this}flipColumns(){const e=Math.ceil(this.rows/2);for(let s=0;s<this.columns;s++)for(let t=0;t<e;t++){let o=this.get(t,s),r=this.get(this.rows-1-t,s);this.set(t,s,r),this.set(this.rows-1-t,s,o)}return this}kroneckerProduct(e){e=S.checkMatrix(e);let s=this.rows,t=this.columns,o=e.rows,r=e.columns,l=new S(s*o,t*r);for(let n=0;n<s;n++)for(let h=0;h<t;h++)for(let u=0;u<o;u++)for(let w=0;w<r;w++)l.set(o*n+u,r*h+w,this.get(n,h)*e.get(u,w));return l}kroneckerSum(e){if(e=S.checkMatrix(e),!this.isSquare()||!e.isSquare())throw new Error("Kronecker Sum needs two Square Matrices");let s=this.rows,t=e.rows,o=this.kroneckerProduct(S.eye(t,t)),r=S.eye(s,s).kroneckerProduct(e);return o.add(r)}transpose(){let e=new S(this.columns,this.rows);for(let s=0;s<this.rows;s++)for(let t=0;t<this.columns;t++)e.set(t,s,this.get(s,t));return e}sortRows(e=St){for(let s=0;s<this.rows;s++)this.setRow(s,this.getRow(s).sort(e));return this}sortColumns(e=St){for(let s=0;s<this.columns;s++)this.setColumn(s,this.getColumn(s).sort(e));return this}subMatrix(e,s,t,o){mt(this,e,s,t,o);let r=new S(s-e+1,o-t+1);for(let l=e;l<=s;l++)for(let n=t;n<=o;n++)r.set(l-e,n-t,this.get(l,n));return r}subMatrixRow(e,s,t){if(s===void 0&&(s=0),t===void 0&&(t=this.columns-1),s>t||s<0||s>=this.columns||t<0||t>=this.columns)throw new RangeError("Argument out of range");let o=new S(e.length,t-s+1);for(let r=0;r<e.length;r++)for(let l=s;l<=t;l++){if(e[r]<0||e[r]>=this.rows)throw new RangeError(`Row index out of range: ${e[r]}`);o.set(r,l-s,this.get(e[r],l))}return o}subMatrixColumn(e,s,t){if(s===void 0&&(s=0),t===void 0&&(t=this.rows-1),s>t||s<0||s>=this.rows||t<0||t>=this.rows)throw new RangeError("Argument out of range");let o=new S(t-s+1,e.length);for(let r=0;r<e.length;r++)for(let l=s;l<=t;l++){if(e[r]<0||e[r]>=this.columns)throw new RangeError(`Column index out of range: ${e[r]}`);o.set(l-s,r,this.get(l,e[r]))}return o}setSubMatrix(e,s,t){if(e=S.checkMatrix(e),e.isEmpty())return this;let o=s+e.rows-1,r=t+e.columns-1;mt(this,s,o,t,r);for(let l=0;l<e.rows;l++)for(let n=0;n<e.columns;n++)this.set(s+l,t+n,e.get(l,n));return this}selection(e,s){gt(this,e),wt(this,s);let t=new S(e.length,s.length);for(let o=0;o<e.length;o++){let r=e[o];for(let l=0;l<s.length;l++){let n=s[l];t.set(o,l,this.get(r,n))}}return t}trace(){let e=Math.min(this.rows,this.columns),s=0;for(let t=0;t<e;t++)s+=this.get(t,t);return s}clone(){return this.constructor.copy(this,new S(this.rows,this.columns))}static copy(e,s){for(const[t,o,r]of e.entries())s.set(t,o,r);return s}sum(e){switch(e){case"row":return Vt(this);case"column":return Ct(this);case void 0:return Pt(this);default:throw new Error(`invalid option: ${e}`)}}product(e){switch(e){case"row":return $t(this);case"column":return Lt(this);case void 0:return Bt(this);default:throw new Error(`invalid option: ${e}`)}}mean(e){const s=this.sum(e);switch(e){case"row":{for(let t=0;t<this.rows;t++)s[t]/=this.columns;return s}case"column":{for(let t=0;t<this.columns;t++)s[t]/=this.rows;return s}case void 0:return s/this.size;default:throw new Error(`invalid option: ${e}`)}}variance(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");const{unbiased:t=!0,mean:o=this.mean(e)}=s;if(typeof t!="boolean")throw new TypeError("unbiased must be a boolean");switch(e){case"row":{if(!B(o))throw new TypeError("mean must be an array");return Ut(this,t,o)}case"column":{if(!B(o))throw new TypeError("mean must be an array");return Ot(this,t,o)}case void 0:{if(typeof o!="number")throw new TypeError("mean must be a number");return Qt(this,t,o)}default:throw new Error(`invalid option: ${e}`)}}standardDeviation(e,s){typeof e=="object"&&(s=e,e=void 0);const t=this.variance(e,s);if(e===void 0)return Math.sqrt(t);for(let o=0;o<t.length;o++)t[o]=Math.sqrt(t[o]);return t}center(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");const{center:t=this.mean(e)}=s;switch(e){case"row":{if(!B(t))throw new TypeError("center must be an array");return Wt(this,t),this}case"column":{if(!B(t))throw new TypeError("center must be an array");return Yt(this,t),this}case void 0:{if(typeof t!="number")throw new TypeError("center must be a number");return Jt(this,t),this}default:throw new Error(`invalid option: ${e}`)}}scale(e,s={}){if(typeof e=="object"&&(s=e,e=void 0),typeof s!="object")throw new TypeError("options must be an object");let t=s.scale;switch(e){case"row":{if(t===void 0)t=Gt(this);else if(!B(t))throw new TypeError("scale must be an array");return Kt(this,t),this}case"column":{if(t===void 0)t=Xt(this);else if(!B(t))throw new TypeError("scale must be an array");return Zt(this,t),this}case void 0:{if(t===void 0)t=xt(this);else if(typeof t!="number")throw new TypeError("scale must be a number");return At(this,t),this}default:throw new Error(`invalid option: ${e}`)}}toString(e){return Mt(this,e)}[Symbol.iterator](){return this.entries()}*entries(){for(let e=0;e<this.rows;e++)for(let s=0;s<this.columns;s++)yield[e,s,this.get(e,s)]}*values(){for(let e=0;e<this.rows;e++)for(let s=0;s<this.columns;s++)yield this.get(e,s)}}T.prototype.klass="Matrix",typeof Symbol<"u"&&(T.prototype[Symbol.for("nodejs.util.inspect.custom")]=_t);function St(i,e){return i-e}function Ht(i){return i.every(e=>typeof e=="number")}T.random=T.rand,T.randomInt=T.randInt,T.diagonal=T.diag,T.prototype.diagonal=T.prototype.diag,T.identity=T.eye,T.prototype.negate=T.prototype.neg,T.prototype.tensorProduct=T.prototype.kroneckerProduct;class S extends T{data;#t(e,s){if(this.data=[],Number.isInteger(s)&&s>=0)for(let t=0;t<e;t++)this.data.push(new Float64Array(s));else throw new TypeError("nColumns must be a positive integer");this.rows=e,this.columns=s}constructor(e,s){if(super(),S.isMatrix(e))this.#t(e.rows,e.columns),S.copy(e,this);else if(Number.isInteger(e)&&e>=0)this.#t(e,s);else if(B(e)){const t=e;if(e=t.length,s=e?t[0].length:0,typeof s!="number")throw new TypeError("Data must be a 2D array with at least one element");this.data=[];for(let o=0;o<e;o++){if(t[o].length!==s)throw new RangeError("Inconsistent array dimensions");if(!Ht(t[o]))throw new TypeError("Input data contains non-numeric values");this.data.push(Float64Array.from(t[o]))}this.rows=e,this.columns=s}else throw new TypeError("First argument must be a positive number or an array")}set(e,s,t){return this.data[e][s]=t,this}get(e,s){return this.data[e][s]}removeRow(e){return J(this,e),this.data.splice(e,1),this.rows-=1,this}addRow(e,s){return s===void 0&&(s=e,e=this.rows),J(this,e,!0),s=Float64Array.from(tt(this,s)),this.data.splice(e,0,s),this.rows+=1,this}removeColumn(e){G(this,e);for(let s=0;s<this.rows;s++){const t=new Float64Array(this.columns-1);for(let o=0;o<e;o++)t[o]=this.data[s][o];for(let o=e+1;o<this.columns;o++)t[o-1]=this.data[s][o];this.data[s]=t}return this.columns-=1,this}addColumn(e,s){typeof s>"u"&&(s=e,e=this.columns),G(this,e,!0),s=et(this,s);for(let t=0;t<this.rows;t++){const o=new Float64Array(this.columns+1);let r=0;for(;r<e;r++)o[r]=this.data[t][r];for(o[r++]=s[t];r<this.columns+1;r++)o[r]=this.data[t][r-1];this.data[t]=o}return this.columns+=1,this}}zt(T,S);class H extends T{#t;get size(){return this.#t.size}get rows(){return this.#t.rows}get columns(){return this.#t.columns}get diagonalSize(){return this.rows}static isSymmetricMatrix(e){return S.isMatrix(e)&&e.klassType==="SymmetricMatrix"}static zeros(e){return new this(e)}static ones(e){return new this(e).fill(1)}constructor(e){if(super(),S.isMatrix(e)){if(!e.isSymmetric())throw new TypeError("not symmetric data");this.#t=S.copy(e,new S(e.rows,e.rows))}else if(Number.isInteger(e)&&e>=0)this.#t=new S(e,e);else if(this.#t=new S(e),!this.isSymmetric())throw new TypeError("not symmetric data")}clone(){const e=new H(this.diagonalSize);for(const[s,t,o]of this.upperRightEntries())e.set(s,t,o);return e}toMatrix(){return new S(this)}get(e,s){return this.#t.get(e,s)}set(e,s,t){return this.#t.set(e,s,t),this.#t.set(s,e,t),this}removeCross(e){return this.#t.removeRow(e),this.#t.removeColumn(e),this}addCross(e,s){s===void 0&&(s=e,e=this.diagonalSize);const t=s.slice();return t.splice(e,1),this.#t.addRow(e,t),this.#t.addColumn(e,s),this}applyMask(e){if(e.length!==this.diagonalSize)throw new RangeError("Mask size do not match with matrix size");const s=[];for(const[t,o]of e.entries())o||s.push(t);s.reverse();for(const t of s)this.removeCross(t);return this}toCompact(){const{diagonalSize:e}=this,s=new Array(e*(e+1)/2);for(let t=0,o=0,r=0;r<s.length;r++)s[r]=this.get(o,t),++t>=e&&(t=++o);return s}static fromCompact(e){const s=e.length,t=(Math.sqrt(8*s+1)-1)/2;if(!Number.isInteger(t))throw new TypeError(`This array is not a compact representation of a Symmetric Matrix, ${JSON.stringify(e)}`);const o=new H(t);for(let r=0,l=0,n=0;n<s;n++)o.set(r,l,e[n]),++r>=t&&(r=++l);return o}*upperRightEntries(){for(let e=0,s=0;e<this.diagonalSize;void 0){const t=this.get(e,s);yield[e,s,t],++s>=this.diagonalSize&&(s=++e)}}*upperRightValues(){for(let e=0,s=0;e<this.diagonalSize;void 0)yield this.get(e,s),++s>=this.diagonalSize&&(s=++e)}}H.prototype.klassType="SymmetricMatrix";class lt extends H{static isDistanceMatrix(e){return H.isSymmetricMatrix(e)&&e.klassSubType==="DistanceMatrix"}constructor(e){if(super(e),!this.isDistance())throw new TypeError("Provided arguments do no produce a distance matrix")}set(e,s,t){return e===s&&(t=0),super.set(e,s,t)}addCross(e,s){return s===void 0&&(s=e,e=this.diagonalSize),s=s.slice(),s[e]=0,super.addCross(e,s)}toSymmetricMatrix(){return new H(this)}clone(){const e=new lt(this.diagonalSize);for(const[s,t,o]of this.upperRightEntries())s!==t&&e.set(s,t,o);return e}toCompact(){const{diagonalSize:e}=this,s=(e-1)*e/2,t=new Array(s);for(let o=1,r=0,l=0;l<t.length;l++)t[l]=this.get(r,o),++o>=e&&(o=++r+1);return t}static fromCompact(e){const s=e.length;if(s===0)return new this(0);const t=(Math.sqrt(8*s+1)+1)/2;if(!Number.isInteger(t))throw new TypeError(`This array is not a compact representation of a DistanceMatrix, ${JSON.stringify(e)}`);const o=new this(t);for(let r=1,l=0,n=0;n<s;n++)o.set(r,l,e[n]),++r>=t&&(r=++l+1);return o}}lt.prototype.klassSubType="DistanceMatrix";class Z extends T{constructor(e,s,t){super(),this.matrix=e,this.rows=s,this.columns=t}}class te extends Z{constructor(e,s){G(e,s),super(e,e.rows,1),this.column=s}set(e,s,t){return this.matrix.set(e,this.column,t),this}get(e){return this.matrix.get(e,this.column)}}class ee extends Z{constructor(e,s){wt(e,s),super(e,e.rows,s.length),this.columnIndices=s}set(e,s,t){return this.matrix.set(e,this.columnIndices[s],t),this}get(e,s){return this.matrix.get(e,this.columnIndices[s])}}class se extends Z{constructor(e){super(e,e.rows,e.columns)}set(e,s,t){return this.matrix.set(e,this.columns-s-1,t),this}get(e,s){return this.matrix.get(e,this.columns-s-1)}}class oe extends Z{constructor(e){super(e,e.rows,e.columns)}set(e,s,t){return this.matrix.set(this.rows-e-1,s,t),this}get(e,s){return this.matrix.get(this.rows-e-1,s)}}class re extends Z{constructor(e,s){J(e,s),super(e,1,e.columns),this.row=s}set(e,s,t){return this.matrix.set(this.row,s,t),this}get(e,s){return this.matrix.get(this.row,s)}}class ie extends Z{constructor(e,s){gt(e,s),super(e,s.length,e.columns),this.rowIndices=s}set(e,s,t){return this.matrix.set(this.rowIndices[e],s,t),this}get(e,s){return this.matrix.get(this.rowIndices[e],s)}}class ht extends Z{constructor(e,s,t){gt(e,s),wt(e,t),super(e,s.length,t.length),this.rowIndices=s,this.columnIndices=t}set(e,s,t){return this.matrix.set(this.rowIndices[e],this.columnIndices[s],t),this}get(e,s){return this.matrix.get(this.rowIndices[e],this.columnIndices[s])}}class ne extends Z{constructor(e,s,t,o,r){mt(e,s,t,o,r),super(e,t-s+1,r-o+1),this.startRow=s,this.startColumn=o}set(e,s,t){return this.matrix.set(this.startRow+e,this.startColumn+s,t),this}get(e,s){return this.matrix.get(this.startRow+e,this.startColumn+s)}}class le extends Z{constructor(e){super(e,e.columns,e.rows)}set(e,s,t){return this.matrix.set(s,e,t),this}get(e,s){return this.matrix.get(s,e)}}class Et extends T{constructor(e,s={}){const{rows:t=1}=s;if(e.length%t!==0)throw new Error("the data length is not divisible by the number of rows");super(),this.rows=t,this.columns=e.length/t,this.data=e}set(e,s,t){let o=this._calculateIndex(e,s);return this.data[o]=t,this}get(e,s){let t=this._calculateIndex(e,s);return this.data[t]}_calculateIndex(e,s){return e*this.columns+s}}class Q extends T{constructor(e){super(),this.data=e,this.rows=e.length,this.columns=e[0].length}set(e,s,t){return this.data[e][s]=t,this}get(e,s){return this.data[e][s]}}function he(i,e){if(B(i))return i[0]&&B(i[0])?new Q(i):new Et(i,e);throw new Error("the argument is not an array")}class ut{constructor(e){e=Q.checkMatrix(e);let s=e.clone(),t=s.rows,o=s.columns,r=new Float64Array(t),l=1,n,h,u,w,f,a,j,m,p;for(n=0;n<t;n++)r[n]=n;for(m=new Float64Array(t),h=0;h<o;h++){for(n=0;n<t;n++)m[n]=s.get(n,h);for(n=0;n<t;n++){for(p=Math.min(n,h),f=0,u=0;u<p;u++)f+=s.get(n,u)*m[u];m[n]-=f,s.set(n,h,m[n])}for(w=h,n=h+1;n<t;n++)Math.abs(m[n])>Math.abs(m[w])&&(w=n);if(w!==h){for(u=0;u<o;u++)a=s.get(w,u),s.set(w,u,s.get(h,u)),s.set(h,u,a);j=r[w],r[w]=r[h],r[h]=j,l=-l}if(h<t&&s.get(h,h)!==0)for(n=h+1;n<t;n++)s.set(n,h,s.get(n,h)/s.get(h,h))}this.LU=s,this.pivotVector=r,this.pivotSign=l}isSingular(){let e=this.LU,s=e.columns;for(let t=0;t<s;t++)if(e.get(t,t)===0)return!0;return!1}solve(e){e=S.checkMatrix(e);let s=this.LU;if(s.rows!==e.rows)throw new Error("Invalid matrix dimensions");if(this.isSingular())throw new Error("LU matrix is singular");let o=e.columns,r=e.subMatrixRow(this.pivotVector,0,o-1),l=s.columns,n,h,u;for(u=0;u<l;u++)for(n=u+1;n<l;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u));for(u=l-1;u>=0;u--){for(h=0;h<o;h++)r.set(u,h,r.get(u,h)/s.get(u,u));for(n=0;n<u;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u))}return r}get determinant(){let e=this.LU;if(!e.isSquare())throw new Error("Matrix must be square");let s=this.pivotSign,t=e.columns;for(let o=0;o<t;o++)s*=e.get(o,o);return s}get lowerTriangularMatrix(){let e=this.LU,s=e.rows,t=e.columns,o=new S(s,t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)r>l?o.set(r,l,e.get(r,l)):r===l?o.set(r,l,1):o.set(r,l,0);return o}get upperTriangularMatrix(){let e=this.LU,s=e.rows,t=e.columns,o=new S(s,t);for(let r=0;r<s;r++)for(let l=0;l<t;l++)r<=l?o.set(r,l,e.get(r,l)):o.set(r,l,0);return o}get pivotPermutationVector(){return Array.from(this.pivotVector)}}function at(i){const e=i.data,s=i.rows;for(let t=0;t<s;t++){const o=e[t];for(let r=t+1;r<s;r++){const l=o[r];o[r]=e[r][t],e[r][t]=l}}return i}function x(i,e){let s=0;return Math.abs(i)>Math.abs(e)?(s=e/i,Math.abs(i)*Math.sqrt(1+s*s)):e!==0?(s=i/e,Math.abs(e)*Math.sqrt(1+s*s)):0}class pt{constructor(e){e=Q.checkMatrix(e);let s=e.clone(),t=e.rows,o=e.columns,r=new Float64Array(o),l,n,h,u;for(h=0;h<o;h++){let w=0;for(l=h;l<t;l++)w=x(w,s.get(l,h));if(w!==0){for(s.get(h,h)<0&&(w=-w),l=h;l<t;l++)s.set(l,h,s.get(l,h)/w);for(s.set(h,h,s.get(h,h)+1),n=h+1;n<o;n++){for(u=0,l=h;l<t;l++)u+=s.get(l,h)*s.get(l,n);for(u=-u/s.get(h,h),l=h;l<t;l++)s.set(l,n,s.get(l,n)+u*s.get(l,h))}}r[h]=-w}this.QR=s,this.Rdiag=r}solve(e){e=S.checkMatrix(e);let s=this.QR,t=s.rows;if(e.rows!==t)throw new Error("Matrix row dimensions must agree");if(!this.isFullRank())throw new Error("Matrix is rank deficient");let o=e.columns,r=e.clone(),l=s.columns,n,h,u,w;for(u=0;u<l;u++)for(h=0;h<o;h++){for(w=0,n=u;n<t;n++)w+=s.get(n,u)*r.get(n,h);for(w=-w/s.get(u,u),n=u;n<t;n++)r.set(n,h,r.get(n,h)+w*s.get(n,u))}for(u=l-1;u>=0;u--){for(h=0;h<o;h++)r.set(u,h,r.get(u,h)/this.Rdiag[u]);for(n=0;n<u;n++)for(h=0;h<o;h++)r.set(n,h,r.get(n,h)-r.get(u,h)*s.get(n,u))}return r.subMatrix(0,l-1,0,o-1)}isFullRank(){let e=this.QR.columns;for(let s=0;s<e;s++)if(this.Rdiag[s]===0)return!1;return!0}get upperTriangularMatrix(){let e=this.QR,s=e.columns,t=new S(s,s),o,r;for(o=0;o<s;o++)for(r=0;r<s;r++)o<r?t.set(o,r,e.get(o,r)):o===r?t.set(o,r,this.Rdiag[o]):t.set(o,r,0);return t}get orthogonalMatrix(){let e=this.QR,s=e.rows,t=e.columns,o=new S(s,t),r,l,n,h;for(n=t-1;n>=0;n--){for(r=0;r<s;r++)o.set(r,n,0);for(o.set(n,n,1),l=n;l<t;l++)if(e.get(n,n)!==0){for(h=0,r=n;r<s;r++)h+=e.get(r,n)*o.get(r,l);for(h=-h/e.get(n,n),r=n;r<s;r++)o.set(r,l,o.get(r,l)+h*e.get(r,n))}}return o}}class ot{constructor(e,s={}){if(e=Q.checkMatrix(e),e.isEmpty())throw new Error("Matrix must be non-empty");let t=e.rows,o=e.columns;const{computeLeftSingularVectors:r=!0,computeRightSingularVectors:l=!0,autoTranspose:n=!1}=s;let h=!!r,u=!!l,w=!1,f;if(t<o)if(!n)console.warn("Computing SVD on a matrix with more columns than rows. Consider enabling autoTranspose"),f=e.transpose();else{f=e.clone(),t=e.columns,o=e.rows,w=!0;let g=h;h=u,u=g}else f=e.transpose();let a=Math.min(t,o),j=Math.min(t+1,o),m=new Float64Array(j),p=new S(a,t),E=new S(o,o),c=new Float64Array(o),M=new Float64Array(t),k=new Float64Array(j);for(let g=0;g<j;g++)k[g]=g;let I=Math.min(t-1,o),N=Math.max(0,Math.min(o-2,t)),V=Math.max(I,N);for(let g=0;g<V;g++){if(g<I){m[g]=0;for(let d=g;d<t;d++)m[g]=x(m[g],f.get(g,d));if(m[g]!==0){f.get(g,g)<0&&(m[g]=-m[g]);for(let d=g;d<t;d++)f.set(g,d,f.get(g,d)/m[g]);f.set(g,g,f.get(g,g)+1)}m[g]=-m[g]}for(let d=g+1;d<o;d++){if(g<I&&m[g]!==0){let b=0;for(let y=g;y<t;y++)b+=f.get(g,y)*f.get(d,y);b=-b/f.get(g,g);for(let y=g;y<t;y++)f.set(d,y,f.get(d,y)+b*f.get(g,y))}c[d]=f.get(d,g)}if(h&&g<I)for(let d=g;d<t;d++)p.set(g,d,f.get(g,d));if(g<N){c[g]=0;for(let d=g+1;d<o;d++)c[g]=x(c[g],c[d]);if(c[g]!==0){c[g+1]<0&&(c[g]=0-c[g]);for(let d=g+1;d<o;d++)c[d]/=c[g];c[g+1]+=1}if(c[g]=-c[g],g+1<t&&c[g]!==0){for(let d=g+1;d<t;d++)M[d]=0;for(let d=g+1;d<t;d++)for(let b=g+1;b<o;b++)M[d]+=c[b]*f.get(b,d);for(let d=g+1;d<o;d++){let b=-c[d]/c[g+1];for(let y=g+1;y<t;y++)f.set(d,y,f.get(d,y)+b*M[y])}}if(u)for(let d=g+1;d<o;d++)E.set(g,d,c[d])}}let v=Math.min(o,t+1);if(I<o&&(m[I]=f.get(I,I)),t<v&&(m[v-1]=0),N+1<v&&(c[N]=f.get(v-1,N)),c[v-1]=0,h){for(let g=I;g<a;g++){for(let d=0;d<t;d++)p.set(g,d,0);p.set(g,g,1)}for(let g=I-1;g>=0;g--)if(m[g]!==0){for(let d=g+1;d<a;d++){let b=0;for(let y=g;y<t;y++)b+=p.get(g,y)*p.get(d,y);b=-b/p.get(g,g);for(let y=g;y<t;y++)p.set(d,y,p.get(d,y)+b*p.get(g,y))}for(let d=g;d<t;d++)p.set(g,d,-p.get(g,d));p.set(g,g,1+p.get(g,g));for(let d=0;d<g-1;d++)p.set(g,d,0)}else{for(let d=0;d<t;d++)p.set(g,d,0);p.set(g,g,1)}}if(u)for(let g=o-1;g>=0;g--){if(g<N&&c[g]!==0)for(let d=g+1;d<o;d++){let b=0;for(let y=g+1;y<o;y++)b+=E.get(g,y)*E.get(d,y);b=-b/E.get(g,g+1);for(let y=g+1;y<o;y++)E.set(d,y,E.get(d,y)+b*E.get(g,y))}for(let d=0;d<o;d++)E.set(g,d,0);E.set(g,g,1)}let R=v-1,F=Number.EPSILON;for(;v>0;){let g,d;for(g=v-2;g>=-1&&g!==-1;g--){const b=Number.MIN_VALUE+F*Math.abs(m[g]+Math.abs(m[g+1]));if(Math.abs(c[g])<=b||Number.isNaN(c[g])){c[g]=0;break}}if(g===v-2)d=4;else{let b;for(b=v-1;b>=g&&b!==g;b--){let y=(b!==v?Math.abs(c[b]):0)+(b!==g+1?Math.abs(c[b-1]):0);if(Math.abs(m[b])<=F*y){m[b]=0;break}}b===g?d=3:b===v-1?d=1:(d=2,g=b)}switch(g++,d){case 1:{let b=c[v-2];c[v-2]=0;for(let y=v-2;y>=g;y--){let C=x(m[y],b),D=m[y]/C,z=b/C;if(m[y]=C,y!==g&&(b=-z*c[y-1],c[y-1]=D*c[y-1]),u)for(let P=0;P<o;P++)C=D*E.get(y,P)+z*E.get(v-1,P),E.set(v-1,P,-z*E.get(y,P)+D*E.get(v-1,P)),E.set(y,P,C)}break}case 2:{let b=c[g-1];c[g-1]=0;for(let y=g;y<v;y++){let C=x(m[y],b),D=m[y]/C,z=b/C;if(m[y]=C,b=-z*c[y],c[y]=D*c[y],h)for(let P=0;P<t;P++)C=D*p.get(y,P)+z*p.get(g-1,P),p.set(g-1,P,-z*p.get(y,P)+D*p.get(g-1,P)),p.set(y,P,C)}break}case 3:{const b=Math.max(Math.abs(m[v-1]),Math.abs(m[v-2]),Math.abs(c[v-2]),Math.abs(m[g]),Math.abs(c[g])),y=m[v-1]/b,C=m[v-2]/b,D=c[v-2]/b,z=m[g]/b,P=c[g]/b,U=((C+y)*(C-y)+D*D)/2,K=y*D*(y*D);let $=0;(U!==0||K!==0)&&(U<0?$=0-Math.sqrt(U*U+K):$=Math.sqrt(U*U+K),$=K/(U+$));let X=(z+y)*(z-y)+$,A=z*P;for(let _=g;_<v-1;_++){let O=x(X,A);O===0&&(O=Number.MIN_VALUE);let W=X/O,Y=A/O;if(_!==g&&(c[_-1]=O),X=W*m[_]+Y*c[_],c[_]=W*c[_]-Y*m[_],A=Y*m[_+1],m[_+1]=W*m[_+1],u)for(let L=0;L<o;L++)O=W*E.get(_,L)+Y*E.get(_+1,L),E.set(_+1,L,-Y*E.get(_,L)+W*E.get(_+1,L)),E.set(_,L,O);if(O=x(X,A),O===0&&(O=Number.MIN_VALUE),W=X/O,Y=A/O,m[_]=O,X=W*c[_]+Y*m[_+1],m[_+1]=-Y*c[_]+W*m[_+1],A=Y*c[_+1],c[_+1]=W*c[_+1],h&&_<t-1)for(let L=0;L<t;L++)O=W*p.get(_,L)+Y*p.get(_+1,L),p.set(_+1,L,-Y*p.get(_,L)+W*p.get(_+1,L)),p.set(_,L,O)}c[v-2]=X;break}case 4:{if(m[g]<=0&&(m[g]=m[g]<0?-m[g]:0,u))for(let b=0;b<=R;b++)E.set(g,b,-E.get(g,b));for(;g<R&&!(m[g]>=m[g+1]);){let b=m[g];if(m[g]=m[g+1],m[g+1]=b,u&&g<o-1)for(let y=0;y<o;y++)b=E.get(g+1,y),E.set(g+1,y,E.get(g,y)),E.set(g,y,b);if(h&&g<t-1)for(let y=0;y<t;y++)b=p.get(g+1,y),p.set(g+1,y,p.get(g,y)),p.set(g,y,b);g++}v--;break}}}if(p=p.isSquare()?at(p):p.transpose(),E=at(E),w){let g=E;E=p,p=g}this.m=t,this.n=o,this.s=m,this.U=p,this.V=E}solve(e){let s=e,t=this.threshold,o=this.s.length,r=S.zeros(o,o);for(let a=0;a<o;a++)Math.abs(this.s[a])<=t?r.set(a,a,0):r.set(a,a,1/this.s[a]);let l=this.U,n=this.rightSingularVectors,h=n.mmul(r),u=n.rows,w=l.rows,f=S.zeros(u,w);for(let a=0;a<u;a++)for(let j=0;j<w;j++){let m=0;for(let p=0;p<o;p++)m+=h.get(a,p)*l.get(j,p);f.set(a,j,m)}return f.mmul(s)}solveForDiagonal(e){return this.solve(S.diag(e))}inverse(){let e=this.V,s=this.threshold,t=e.rows,o=e.columns,r=new S(t,this.s.length);for(let w=0;w<t;w++)for(let f=0;f<o;f++)Math.abs(this.s[f])>s&&r.set(w,f,e.get(w,f)/this.s[f]);let l=this.U,n=l.rows,h=l.columns,u=new S(t,n);for(let w=0;w<t;w++)for(let f=0;f<n;f++){let a=0;for(let j=0;j<h;j++)a+=r.get(w,j)*l.get(f,j);u.set(w,f,a)}return u}get condition(){return this.s[0]/this.s[Math.min(this.m,this.n)-1]}get norm2(){return this.s[0]}get rank(){let e=Math.max(this.m,this.n)*this.s[0]*Number.EPSILON,s=0,t=this.s;for(let o=0,r=t.length;o<r;o++)t[o]>e&&s++;return s}get diagonal(){return Array.from(this.s)}get threshold(){return Number.EPSILON/2*Math.max(this.m,this.n)*this.s[0]}get leftSingularVectors(){return this.U}get rightSingularVectors(){return this.V}get diagonalMatrix(){return S.diag(this.s)}}function ue(i,e=!1){return i=Q.checkMatrix(i),e?new ot(i).inverse():bt(i,S.eye(i.rows))}function bt(i,e,s=!1){return i=Q.checkMatrix(i),e=Q.checkMatrix(e),s?new ot(i).solve(e):i.isSquare()?new ut(i).solve(e):new pt(i).solve(e)}function ft(i){if(i=S.checkMatrix(i),i.isSquare()){if(i.columns===0)return 1;let e,s,t,o;if(i.columns===2)return e=i.get(0,0),s=i.get(0,1),t=i.get(1,0),o=i.get(1,1),e*o-s*t;if(i.columns===3){let r,l,n;return r=new ht(i,[1,2],[1,2]),l=new ht(i,[1,2],[0,2]),n=new ht(i,[1,2],[0,1]),e=i.get(0,0),s=i.get(0,1),t=i.get(0,2),e*ft(r)-s*ft(l)+t*ft(n)}else return new ut(i).determinant}else throw Error("determinant can only be calculated for a square matrix")}function fe(i,e){let s=[];for(let t=0;t<i;t++)t!==e&&s.push(t);return s}function ce(i,e,s,t=1e-9,o=1e-9){if(i>o)return new Array(e.rows+1).fill(0);{let r=e.addRow(s,[0]);for(let l=0;l<r.rows;l++)Math.abs(r.get(l,0))<t&&r.set(l,0,0);return r.to1DArray()}}function ge(i,e={}){const{thresholdValue:s=1e-9,thresholdError:t=1e-9}=e;i=S.checkMatrix(i);let o=i.rows,r=new S(o,o);for(let l=0;l<o;l++){let n=S.columnVector(i.getRow(l)),h=i.subMatrixRow(fe(o,l)).transpose(),w=new ot(h).solve(n),f=S.sub(n,h.mmul(w)).abs().max();r.setRow(l,ce(f,w,l,s,t))}return r}function we(i,e=Number.EPSILON){if(i=S.checkMatrix(i),i.isEmpty())return i.transpose();let s=new ot(i,{autoTranspose:!0}),t=s.leftSingularVectors,o=s.rightSingularVectors,r=s.diagonal;for(let l=0;l<r.length;l++)Math.abs(r[l])>e?r[l]=1/r[l]:r[l]=0;return o.mmul(S.diag(r).mmul(t.transpose()))}function me(i,e=i,s={}){i=new S(i);let t=!1;if(typeof e=="object"&&!S.isMatrix(e)&&!B(e)?(s=e,e=i,t=!0):e=new S(e),i.rows!==e.rows)throw new TypeError("Both matrices must have the same number of rows");const{center:o=!0}=s;o&&(i=i.center("column"),t||(e=e.center("column")));const r=i.transpose().mmul(e);for(let l=0;l<r.rows;l++)for(let n=0;n<r.columns;n++)r.set(l,n,r.get(l,n)*(1/(i.rows-1)));return r}function ae(i,e=i,s={}){i=new S(i);let t=!1;if(typeof e=="object"&&!S.isMatrix(e)&&!B(e)?(s=e,e=i,t=!0):e=new S(e),i.rows!==e.rows)throw new TypeError("Both matrices must have the same number of rows");const{center:o=!0,scale:r=!0}=s;o&&(i.center("column"),t||e.center("column")),r&&(i.scale("column"),t||e.scale("column"));const l=i.standardDeviation("column",{unbiased:!0}),n=t?l:e.standardDeviation("column",{unbiased:!0}),h=i.transpose().mmul(e);for(let u=0;u<h.rows;u++)for(let w=0;w<h.columns;w++)h.set(u,w,h.get(u,w)*(1/(l[u]*n[w]))*(1/(i.rows-1)));return h}class kt{constructor(e,s={}){const{assumeSymmetric:t=!1}=s;if(e=Q.checkMatrix(e),!e.isSquare())throw new Error("Matrix is not a square matrix");if(e.isEmpty())throw new Error("Matrix must be non-empty");let o=e.columns,r=new S(o,o),l=new Float64Array(o),n=new Float64Array(o),h=e,u,w,f=!1;if(t?f=!0:f=e.isSymmetric(),f){for(u=0;u<o;u++)for(w=0;w<o;w++)r.set(w,u,h.get(u,w));pe(o,n,l,r),de(o,n,l,r),at(r)}else{let a=new S(o,o),j=new Float64Array(o);for(w=0;w<o;w++)for(u=0;u<o;u++)a.set(u,w,h.get(u,w));ye(o,a,j,r),Me(o,n,l,r,a)}this.n=o,this.e=n,this.d=l,this.V=r}get realEigenvalues(){return Array.from(this.d)}get imaginaryEigenvalues(){return Array.from(this.e)}get eigenvectorMatrix(){return this.V}get diagonalMatrix(){let e=this.n,s=this.e,t=this.d,o=new S(e,e),r,l;for(r=0;r<e;r++){for(l=0;l<e;l++)o.set(r,l,0);o.set(r,r,t[r]),s[r]>0?o.set(r,r+1,s[r]):s[r]<0&&o.set(r,r-1,s[r])}return o}}function pe(i,e,s,t){let o,r,l,n,h,u,w,f;for(h=0;h<i;h++)s[h]=t.get(h,i-1);for(n=i-1;n>0;n--){for(f=0,l=0,u=0;u<n;u++)f=f+Math.abs(s[u]);if(f===0)for(e[n]=s[n-1],h=0;h<n;h++)s[h]=t.get(h,n-1),t.set(h,n,0),t.set(n,h,0);else{for(u=0;u<n;u++)s[u]/=f,l+=s[u]*s[u];for(o=s[n-1],r=Math.sqrt(l),o>0&&(r=-r),e[n]=f*r,l=l-o*r,s[n-1]=o-r,h=0;h<n;h++)e[h]=0;for(h=0;h<n;h++){for(o=s[h],t.set(n,h,o),r=e[h]+t.get(h,h)*o,u=h+1;u<=n-1;u++)r+=t.get(h,u)*s[u],e[u]+=t.get(h,u)*o;e[h]=r}for(o=0,h=0;h<n;h++)e[h]/=l,o+=e[h]*s[h];for(w=o/(l+l),h=0;h<n;h++)e[h]-=w*s[h];for(h=0;h<n;h++){for(o=s[h],r=e[h],u=h;u<=n-1;u++)t.set(h,u,t.get(h,u)-(o*e[u]+r*s[u]));s[h]=t.get(h,n-1),t.set(h,n,0)}}s[n]=l}for(n=0;n<i-1;n++){if(t.set(n,i-1,t.get(n,n)),t.set(n,n,1),l=s[n+1],l!==0){for(u=0;u<=n;u++)s[u]=t.get(n+1,u)/l;for(h=0;h<=n;h++){for(r=0,u=0;u<=n;u++)r+=t.get(n+1,u)*t.get(h,u);for(u=0;u<=n;u++)t.set(h,u,t.get(h,u)-r*s[u])}}for(u=0;u<=n;u++)t.set(n+1,u,0)}for(h=0;h<i;h++)s[h]=t.get(h,i-1),t.set(h,i-1,0);t.set(i-1,i-1,1),e[0]=0}function de(i,e,s,t){let o,r,l,n,h,u,w,f,a,j,m,p,E,c,M,k;for(l=1;l<i;l++)e[l-1]=e[l];e[i-1]=0;let I=0,N=0,V=Number.EPSILON;for(u=0;u<i;u++){for(N=Math.max(N,Math.abs(s[u])+Math.abs(e[u])),w=u;w<i&&!(Math.abs(e[w])<=V*N);)w++;if(w>u)do{for(o=s[u],f=(s[u+1]-o)/(2*e[u]),a=x(f,1),f<0&&(a=-a),s[u]=e[u]/(f+a),s[u+1]=e[u]*(f+a),j=s[u+1],r=o-s[u],l=u+2;l<i;l++)s[l]-=r;for(I=I+r,f=s[w],m=1,p=m,E=m,c=e[u+1],M=0,k=0,l=w-1;l>=u;l--)for(E=p,p=m,k=M,o=m*e[l],r=m*f,a=x(f,e[l]),e[l+1]=M*a,M=e[l]/a,m=f/a,f=m*s[l]-M*o,s[l+1]=r+M*(m*o+M*s[l]),h=0;h<i;h++)r=t.get(l+1,h),t.set(l+1,h,M*t.get(l,h)+m*r),t.set(l,h,m*t.get(l,h)-M*r);f=-M*k*E*c*e[u]/j,e[u]=M*f,s[u]=m*f}while(Math.abs(e[u])>V*N);s[u]=s[u]+I,e[u]=0}for(l=0;l<i-1;l++){for(h=l,f=s[l],n=l+1;n<i;n++)s[n]<f&&(h=n,f=s[n]);if(h!==l)for(s[h]=s[l],s[l]=f,n=0;n<i;n++)f=t.get(l,n),t.set(l,n,t.get(h,n)),t.set(h,n,f)}}function ye(i,e,s,t){let o=0,r=i-1,l,n,h,u,w,f,a;for(f=o+1;f<=r-1;f++){for(a=0,u=f;u<=r;u++)a=a+Math.abs(e.get(u,f-1));if(a!==0){for(h=0,u=r;u>=f;u--)s[u]=e.get(u,f-1)/a,h+=s[u]*s[u];for(n=Math.sqrt(h),s[f]>0&&(n=-n),h=h-s[f]*n,s[f]=s[f]-n,w=f;w<i;w++){for(l=0,u=r;u>=f;u--)l+=s[u]*e.get(u,w);for(l=l/h,u=f;u<=r;u++)e.set(u,w,e.get(u,w)-l*s[u])}for(u=0;u<=r;u++){for(l=0,w=r;w>=f;w--)l+=s[w]*e.get(u,w);for(l=l/h,w=f;w<=r;w++)e.set(u,w,e.get(u,w)-l*s[w])}s[f]=a*s[f],e.set(f,f-1,a*n)}}for(u=0;u<i;u++)for(w=0;w<i;w++)t.set(u,w,u===w?1:0);for(f=r-1;f>=o+1;f--)if(e.get(f,f-1)!==0){for(u=f+1;u<=r;u++)s[u]=e.get(u,f-1);for(w=f;w<=r;w++){for(n=0,u=f;u<=r;u++)n+=s[u]*t.get(u,w);for(n=n/s[f]/e.get(f,f-1),u=f;u<=r;u++)t.set(u,w,t.get(u,w)+n*s[u])}}}function Me(i,e,s,t,o){let r=i-1,l=0,n=i-1,h=Number.EPSILON,u=0,w=0,f=0,a=0,j=0,m=0,p=0,E=0,c,M,k,I,N,V,v,R,F,g,d,b,y,C,D;for(c=0;c<i;c++)for((c<l||c>n)&&(s[c]=o.get(c,c),e[c]=0),M=Math.max(c-1,0);M<i;M++)w=w+Math.abs(o.get(c,M));for(;r>=l;){for(I=r;I>l&&(m=Math.abs(o.get(I-1,I-1))+Math.abs(o.get(I,I)),m===0&&(m=w),!(Math.abs(o.get(I,I-1))<h*m));)I--;if(I===r)o.set(r,r,o.get(r,r)+u),s[r]=o.get(r,r),e[r]=0,r--,E=0;else if(I===r-1){if(v=o.get(r,r-1)*o.get(r-1,r),f=(o.get(r-1,r-1)-o.get(r,r))/2,a=f*f+v,p=Math.sqrt(Math.abs(a)),o.set(r,r,o.get(r,r)+u),o.set(r-1,r-1,o.get(r-1,r-1)+u),R=o.get(r,r),a>=0){for(p=f>=0?f+p:f-p,s[r-1]=R+p,s[r]=s[r-1],p!==0&&(s[r]=R-v/p),e[r-1]=0,e[r]=0,R=o.get(r,r-1),m=Math.abs(R)+Math.abs(p),f=R/m,a=p/m,j=Math.sqrt(f*f+a*a),f=f/j,a=a/j,M=r-1;M<i;M++)p=o.get(r-1,M),o.set(r-1,M,a*p+f*o.get(r,M)),o.set(r,M,a*o.get(r,M)-f*p);for(c=0;c<=r;c++)p=o.get(c,r-1),o.set(c,r-1,a*p+f*o.get(c,r)),o.set(c,r,a*o.get(c,r)-f*p);for(c=l;c<=n;c++)p=t.get(c,r-1),t.set(c,r-1,a*p+f*t.get(c,r)),t.set(c,r,a*t.get(c,r)-f*p)}else s[r-1]=R+f,s[r]=R+f,e[r-1]=p,e[r]=-p;r=r-2,E=0}else{if(R=o.get(r,r),F=0,v=0,I<r&&(F=o.get(r-1,r-1),v=o.get(r,r-1)*o.get(r-1,r)),E===10){for(u+=R,c=l;c<=r;c++)o.set(c,c,o.get(c,c)-R);m=Math.abs(o.get(r,r-1))+Math.abs(o.get(r-1,r-2)),R=F=.75*m,v=-.4375*m*m}if(E===30&&(m=(F-R)/2,m=m*m+v,m>0)){for(m=Math.sqrt(m),F<R&&(m=-m),m=R-v/((F-R)/2+m),c=l;c<=r;c++)o.set(c,c,o.get(c,c)-m);u+=m,R=F=v=.964}for(E=E+1,N=r-2;N>=I&&(p=o.get(N,N),j=R-p,m=F-p,f=(j*m-v)/o.get(N+1,N)+o.get(N,N+1),a=o.get(N+1,N+1)-p-j-m,j=o.get(N+2,N+1),m=Math.abs(f)+Math.abs(a)+Math.abs(j),f=f/m,a=a/m,j=j/m,!(N===I||Math.abs(o.get(N,N-1))*(Math.abs(a)+Math.abs(j))<h*(Math.abs(f)*(Math.abs(o.get(N-1,N-1))+Math.abs(p)+Math.abs(o.get(N+1,N+1))))));)N--;for(c=N+2;c<=r;c++)o.set(c,c-2,0),c>N+2&&o.set(c,c-3,0);for(k=N;k<=r-1&&(C=k!==r-1,k!==N&&(f=o.get(k,k-1),a=o.get(k+1,k-1),j=C?o.get(k+2,k-1):0,R=Math.abs(f)+Math.abs(a)+Math.abs(j),R!==0&&(f=f/R,a=a/R,j=j/R)),R!==0);k++)if(m=Math.sqrt(f*f+a*a+j*j),f<0&&(m=-m),m!==0){for(k!==N?o.set(k,k-1,-m*R):I!==N&&o.set(k,k-1,-o.get(k,k-1)),f=f+m,R=f/m,F=a/m,p=j/m,a=a/f,j=j/f,M=k;M<i;M++)f=o.get(k,M)+a*o.get(k+1,M),C&&(f=f+j*o.get(k+2,M),o.set(k+2,M,o.get(k+2,M)-f*p)),o.set(k,M,o.get(k,M)-f*R),o.set(k+1,M,o.get(k+1,M)-f*F);for(c=0;c<=Math.min(r,k+3);c++)f=R*o.get(c,k)+F*o.get(c,k+1),C&&(f=f+p*o.get(c,k+2),o.set(c,k+2,o.get(c,k+2)-f*j)),o.set(c,k,o.get(c,k)-f),o.set(c,k+1,o.get(c,k+1)-f*a);for(c=l;c<=n;c++)f=R*t.get(c,k)+F*t.get(c,k+1),C&&(f=f+p*t.get(c,k+2),t.set(c,k+2,t.get(c,k+2)-f*j)),t.set(c,k,t.get(c,k)-f),t.set(c,k+1,t.get(c,k+1)-f*a)}}}if(w!==0){for(r=i-1;r>=0;r--)if(f=s[r],a=e[r],a===0)for(I=r,o.set(r,r,1),c=r-1;c>=0;c--){for(v=o.get(c,c)-f,j=0,M=I;M<=r;M++)j=j+o.get(c,M)*o.get(M,r);if(e[c]<0)p=v,m=j;else if(I=c,e[c]===0?o.set(c,r,v!==0?-j/v:-j/(h*w)):(R=o.get(c,c+1),F=o.get(c+1,c),a=(s[c]-f)*(s[c]-f)+e[c]*e[c],V=(R*m-p*j)/a,o.set(c,r,V),o.set(c+1,r,Math.abs(R)>Math.abs(p)?(-j-v*V)/R:(-m-F*V)/p)),V=Math.abs(o.get(c,r)),h*V*V>1)for(M=c;M<=r;M++)o.set(M,r,o.get(M,r)/V)}else if(a<0)for(I=r-1,Math.abs(o.get(r,r-1))>Math.abs(o.get(r-1,r))?(o.set(r-1,r-1,a/o.get(r,r-1)),o.set(r-1,r,-(o.get(r,r)-f)/o.get(r,r-1))):(D=ct(0,-o.get(r-1,r),o.get(r-1,r-1)-f,a),o.set(r-1,r-1,D[0]),o.set(r-1,r,D[1])),o.set(r,r-1,0),o.set(r,r,1),c=r-2;c>=0;c--){for(g=0,d=0,M=I;M<=r;M++)g=g+o.get(c,M)*o.get(M,r-1),d=d+o.get(c,M)*o.get(M,r);if(v=o.get(c,c)-f,e[c]<0)p=v,j=g,m=d;else if(I=c,e[c]===0?(D=ct(-g,-d,v,a),o.set(c,r-1,D[0]),o.set(c,r,D[1])):(R=o.get(c,c+1),F=o.get(c+1,c),b=(s[c]-f)*(s[c]-f)+e[c]*e[c]-a*a,y=(s[c]-f)*2*a,b===0&&y===0&&(b=h*w*(Math.abs(v)+Math.abs(a)+Math.abs(R)+Math.abs(F)+Math.abs(p))),D=ct(R*j-p*g+a*d,R*m-p*d-a*g,b,y),o.set(c,r-1,D[0]),o.set(c,r,D[1]),Math.abs(R)>Math.abs(p)+Math.abs(a)?(o.set(c+1,r-1,(-g-v*o.get(c,r-1)+a*o.get(c,r))/R),o.set(c+1,r,(-d-v*o.get(c,r)-a*o.get(c,r-1))/R)):(D=ct(-j-F*o.get(c,r-1),-m-F*o.get(c,r),p,a),o.set(c+1,r-1,D[0]),o.set(c+1,r,D[1]))),V=Math.max(Math.abs(o.get(c,r-1)),Math.abs(o.get(c,r))),h*V*V>1)for(M=c;M<=r;M++)o.set(M,r-1,o.get(M,r-1)/V),o.set(M,r,o.get(M,r)/V)}for(c=0;c<i;c++)if(c<l||c>n)for(M=c;M<i;M++)t.set(c,M,o.get(c,M));for(M=i-1;M>=l;M--)for(c=l;c<=n;c++){for(p=0,k=l;k<=Math.min(M,n);k++)p=p+t.get(c,k)*o.get(k,M);t.set(c,M,p)}}}function ct(i,e,s,t){let o,r;return Math.abs(s)>Math.abs(t)?(o=t/s,r=s+o*t,[(i+o*e)/r,(e-o*i)/r]):(o=s/t,r=t+o*s,[(o*i+e)/r,(o*e-i)/r])}class vt{constructor(e){if(e=Q.checkMatrix(e),!e.isSymmetric())throw new Error("Matrix is not symmetric");let s=e,t=s.rows,o=new S(t,t),r=!0,l,n,h;for(n=0;n<t;n++){let u=0;for(h=0;h<n;h++){let w=0;for(l=0;l<h;l++)w+=o.get(h,l)*o.get(n,l);w=(s.get(n,h)-w)/o.get(h,h),o.set(n,h,w),u=u+w*w}for(u=s.get(n,n)-u,r&&=u>0,o.set(n,n,Math.sqrt(Math.max(u,0))),h=n+1;h<t;h++)o.set(n,h,0)}this.L=o,this.positiveDefinite=r}isPositiveDefinite(){return this.positiveDefinite}solve(e){e=Q.checkMatrix(e);let s=this.L,t=s.rows;if(e.rows!==t)throw new Error("Matrix dimensions do not match");if(this.isPositiveDefinite()===!1)throw new Error("Matrix is not positive definite");let o=e.columns,r=e.clone(),l,n,h;for(h=0;h<t;h++)for(n=0;n<o;n++){for(l=0;l<h;l++)r.set(h,n,r.get(h,n)-r.get(l,n)*s.get(h,l));r.set(h,n,r.get(h,n)/s.get(h,h))}for(h=t-1;h>=0;h--)for(n=0;n<o;n++){for(l=h+1;l<t;l++)r.set(h,n,r.get(h,n)-r.get(l,n)*s.get(l,h));r.set(h,n,r.get(h,n)/s.get(h,h))}return r}get lowerTriangularMatrix(){return this.L}}class Rt{constructor(e,s={}){e=Q.checkMatrix(e);let{Y:t}=s;const{scaleScores:o=!1,maxIterations:r=1e3,terminationCriteria:l=1e-10}=s;let n;if(t){if(B(t)&&typeof t[0]=="number"?t=S.columnVector(t):t=Q.checkMatrix(t),t.rows!==e.rows)throw new Error("Y should have the same number of rows as X");n=t.getColumnVector(0)}else n=e.getColumnVector(0);let h=1,u,w,f,a;for(let j=0;j<r&&h>l;j++)f=e.transpose().mmul(n).div(n.transpose().mmul(n).get(0,0)),f=f.div(f.norm()),u=e.mmul(f).div(f.transpose().mmul(f).get(0,0)),j>0&&(h=u.clone().sub(a).pow(2).sum()),a=u.clone(),t?(w=t.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0)),w=w.div(w.norm()),n=t.mmul(w).div(w.transpose().mmul(w).get(0,0))):n=u;if(t){let j=e.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0));j=j.div(j.norm());let m=e.clone().sub(u.clone().mmul(j.transpose())),p=n.transpose().mmul(u).div(u.transpose().mmul(u).get(0,0)),E=t.clone().sub(u.clone().mulS(p.get(0,0)).mmul(w.transpose()));this.t=u,this.p=j.transpose(),this.w=f.transpose(),this.q=w,this.u=n,this.s=u.transpose().mmul(u),this.xResidual=m,this.yResidual=E,this.betas=p}else this.w=f.transpose(),this.s=u.transpose().mmul(u).sqrt(),o?this.t=u.clone().div(this.s.get(0,0)):this.t=u,this.xResidual=e.sub(u.mmul(f.transpose()))}}return q.__exports.AbstractMatrix=T,q.__exports.CHO=vt,q.__exports.CholeskyDecomposition=vt,q.__exports.DistanceMatrix=lt,q.__exports.EVD=kt,q.__exports.EigenvalueDecomposition=kt,q.__exports.LU=ut,q.__exports.LuDecomposition=ut,q.__exports.Matrix=S,q.__exports.MatrixColumnSelectionView=ee,q.__exports.MatrixColumnView=te,q.__exports.MatrixFlipColumnView=se,q.__exports.MatrixFlipRowView=oe,q.__exports.MatrixRowSelectionView=ie,q.__exports.MatrixRowView=re,q.__exports.MatrixSelectionView=ht,q.__exports.MatrixSubView=ne,q.__exports.MatrixTransposeView=le,q.__exports.NIPALS=Rt,q.__exports.Nipals=Rt,q.__exports.QR=pt,q.__exports.QrDecomposition=pt,q.__exports.SVD=ot,q.__exports.SingularValueDecomposition=ot,q.__exports.SymmetricMatrix=H,q.__exports.WrapperMatrix1D=Et,q.__exports.WrapperMatrix2D=Q,q.__exports.correlation=ae,q.__exports.covariance=me,q.__exports.default=S,q.__exports.determinant=ft,q.__exports.inverse=ue,q.__exports.linearDependencies=ge,q.__exports.pseudoInverse=we,q.__exports.solve=bt,q.__exports.wrap=he,q.__exports}exports.__require=Te;
+//# sourceMappingURL=index.cjs561.js.map