{"version":3,"file":"spectrum-generator.min.js","sources":["../node_modules/ml-peak-shape-generator/lib-esm/util/constants.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/1d/gaussian/Gaussian.js","../node_modules/ml-peak-shape-generator/lib-esm/util/erfinv.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/1d/lorentzian/Lorentzian.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/1d/pseudoVoigt/PseudoVoigt.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/1d/generalizedLorentzian/GeneralizedLorentzian.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/2d/gaussian2D/Gaussian2D.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/1d/getShape1D.js","../node_modules/ml-peak-shape-generator/lib-esm/shapes/2d/getShape2D.js","../node_modules/is-any-array/lib-esm/index.js","../node_modules/ml-spectra-processing/lib-esm/x/xFindClosestIndex.js","../node_modules/ml-spectra-processing/lib-esm/utils/createArray.js","../node_modules/ml-xsadd/lib-esm/XSadd.js","../node_modules/ml-spectra-processing/lib-esm/utils/createRandomArray.js","../node_modules/ml-spectra-processing/lib-esm/x/xMaxValue.js","../node_modules/ml-spectra-processing/lib-esm/x/xCheck.js","../node_modules/ml-spectra-processing/lib-esm/x/xGetFromToIndex.js","../node_modules/ml-spectra-processing/lib-esm/matrix/matrixCreateEmpty.js","../node_modules/ml-spectra-processing/lib-esm/matrix/matrixMinMaxZ.js","../node_modules/ml-spectra-processing/lib-esm/matrix/matrixCheck.js","../lib-esm/util/addNoise.js","../node_modules/ml-spectra-processing/lib-esm/x/xAdd.js","../lib-esm/SpectrumGenerator.js","../lib-esm/util/addBaseline.js","../lib-esm/Spectrum2DGenerator.js"],"sourcesContent":["export const GAUSSIAN_EXP_FACTOR = -4 * Math.LN2;\nexport const ROOT_PI_OVER_LN2 = Math.sqrt(Math.PI / Math.LN2);\nexport const ROOT_THREE = Math.sqrt(3);\nexport const ROOT_2LN2 = Math.sqrt(2 * Math.LN2);\nexport const ROOT_2LN2_MINUS_ONE = Math.sqrt(2 * Math.LN2) - 1;\n//# sourceMappingURL=constants.js.map","import { ROOT_2LN2, GAUSSIAN_EXP_FACTOR, ROOT_PI_OVER_LN2, } from '../../../util/constants';\nimport erfinv from '../../../util/erfinv';\nexport class Gaussian {\n    constructor(options = {}) {\n        const { fwhm = 500, sd } = options;\n        this.fwhm = sd ? gaussianWidthToFWHM(2 * sd) : fwhm;\n    }\n    fwhmToWidth(fwhm = this.fwhm) {\n        return gaussianFwhmToWidth(fwhm);\n    }\n    widthToFWHM(width) {\n        return gaussianWidthToFWHM(width);\n    }\n    fct(x) {\n        return gaussianFct(x, this.fwhm);\n    }\n    getArea(height = calculateGaussianHeight({ fwhm: this.fwhm })) {\n        return getGaussianArea({ fwhm: this.fwhm, height });\n    }\n    getFactor(area) {\n        return getGaussianFactor(area);\n    }\n    getData(options = {}) {\n        return getGaussianData(this, options);\n    }\n    calculateHeight(area = 1) {\n        return calculateGaussianHeight({ fwhm: this.fwhm, area });\n    }\n    getParameters() {\n        return ['fwhm'];\n    }\n}\nexport function calculateGaussianHeight(options) {\n    let { fwhm = 500, area = 1, sd } = options;\n    if (sd)\n        fwhm = gaussianWidthToFWHM(2 * sd);\n    return (2 * area) / ROOT_PI_OVER_LN2 / fwhm;\n}\n/**\n * Calculate the height of the gaussian function of a specific width (fwhm) at a speicifc\n * x position (the gaussian is centered on x=0)\n * @param x\n * @param fwhm\n * @returns y\n */\nexport function gaussianFct(x, fwhm) {\n    return Math.exp(GAUSSIAN_EXP_FACTOR * Math.pow(x / fwhm, 2));\n}\nexport function gaussianWidthToFWHM(width) {\n    return width * ROOT_2LN2;\n}\nexport function gaussianFwhmToWidth(fwhm) {\n    return fwhm / ROOT_2LN2;\n}\nexport function getGaussianArea(options) {\n    let { fwhm = 500, sd, height = 1 } = options;\n    if (sd)\n        fwhm = gaussianWidthToFWHM(2 * sd);\n    return (height * ROOT_PI_OVER_LN2 * fwhm) / 2;\n}\nexport function getGaussianFactor(area = 0.9999) {\n    return Math.sqrt(2) * erfinv(area);\n}\nexport function getGaussianData(shape = {}, options = {}) {\n    let { fwhm = 500, sd } = shape;\n    if (sd)\n        fwhm = gaussianWidthToFWHM(2 * sd);\n    let { length, factor = getGaussianFactor(), height = calculateGaussianHeight({ fwhm }), } = options;\n    if (!length) {\n        length = Math.min(Math.ceil(fwhm * factor), Math.pow(2, 25) - 1);\n        if (length % 2 === 0)\n            length++;\n    }\n    const center = (length - 1) / 2;\n    const data = new Float64Array(length);\n    for (let i = 0; i <= center; i++) {\n        data[i] = gaussianFct(i - center, fwhm) * height;\n        data[length - 1 - i] = data[i];\n    }\n    return data;\n}\n//# sourceMappingURL=Gaussian.js.map","// https://en.wikipedia.org/wiki/Error_function#Inverse_functions\n// This code yields to a good approximation\n// If needed a better implementation using polynomial can be found on https://en.wikipedia.org/wiki/Error_function#Inverse_functions\nexport default function erfinv(x) {\n    let a = 0.147;\n    if (x === 0)\n        return 0;\n    let ln1MinusXSqrd = Math.log(1 - x * x);\n    let lnEtcBy2Plus2 = ln1MinusXSqrd / 2 + 2 / (Math.PI * a);\n    let firstSqrt = Math.sqrt(lnEtcBy2Plus2 ** 2 - ln1MinusXSqrd / a);\n    let secondSqrt = Math.sqrt(firstSqrt - lnEtcBy2Plus2);\n    return secondSqrt * (x > 0 ? 1 : -1);\n}\n//# sourceMappingURL=erfinv.js.map","import { ROOT_THREE } from '../../../util/constants';\nexport class Lorentzian {\n    constructor(options = {}) {\n        const { fwhm = 500 } = options;\n        this.fwhm = fwhm;\n    }\n    fwhmToWidth(fwhm = this.fwhm) {\n        return lorentzianFwhmToWidth(fwhm);\n    }\n    widthToFWHM(width) {\n        return lorentzianWidthToFWHM(width);\n    }\n    fct(x) {\n        return lorentzianFct(x, this.fwhm);\n    }\n    getArea(height = 1) {\n        return getLorentzianArea({ fwhm: this.fwhm, height });\n    }\n    getFactor(area) {\n        return getLorentzianFactor(area);\n    }\n    getData(options = {}) {\n        return getLorentzianData(this, options);\n    }\n    calculateHeight(area = 1) {\n        return calculateLorentzianHeight({ fwhm: this.fwhm, area });\n    }\n    getParameters() {\n        return ['fwhm'];\n    }\n}\nexport const calculateLorentzianHeight = ({ fwhm = 1, area = 1 }) => {\n    return (2 * area) / Math.PI / fwhm;\n};\nexport const getLorentzianArea = (options) => {\n    const { fwhm = 500, height = 1 } = options;\n    return (height * Math.PI * fwhm) / 2;\n};\nexport const lorentzianFct = (x, fwhm) => {\n    return fwhm ** 2 / (4 * x ** 2 + fwhm ** 2);\n};\nexport const lorentzianWidthToFWHM = (width) => {\n    return width * ROOT_THREE;\n};\nexport const lorentzianFwhmToWidth = (fwhm) => {\n    return fwhm / ROOT_THREE;\n};\nexport const getLorentzianFactor = (area = 0.9999) => {\n    if (area >= 1) {\n        throw new Error('area should be (0 - 1)');\n    }\n    const halfResidual = (1 - area) * 0.5;\n    const quantileFunction = (p) => Math.tan(Math.PI * (p - 0.5));\n    return ((quantileFunction(1 - halfResidual) - quantileFunction(halfResidual)) / 2);\n};\nexport const getLorentzianData = (shape = {}, options = {}) => {\n    let { fwhm = 500 } = shape;\n    let { length, factor = getLorentzianFactor(), height = calculateLorentzianHeight({ fwhm, area: 1 }), } = options;\n    if (!length) {\n        length = Math.min(Math.ceil(fwhm * factor), Math.pow(2, 25) - 1);\n        if (length % 2 === 0)\n            length++;\n    }\n    const center = (length - 1) / 2;\n    const data = new Float64Array(length);\n    for (let i = 0; i <= center; i++) {\n        data[i] = lorentzianFct(i - center, fwhm) * height;\n        data[length - 1 - i] = data[i];\n    }\n    return data;\n};\n//# sourceMappingURL=Lorentzian.js.map","import { GAUSSIAN_EXP_FACTOR, ROOT_2LN2_MINUS_ONE, ROOT_PI_OVER_LN2, } from '../../../util/constants';\nimport { gaussianFct, getGaussianFactor } from '../gaussian/Gaussian';\nimport { lorentzianFct, getLorentzianFactor } from '../lorentzian/Lorentzian';\nexport class PseudoVoigt {\n    constructor(options = {}) {\n        const { fwhm = 500, mu = 0.5 } = options;\n        this.mu = mu;\n        this.fwhm = fwhm;\n    }\n    fwhmToWidth(fwhm = this.fwhm, mu = this.mu) {\n        return pseudoVoigtFwhmToWidth(fwhm, mu);\n    }\n    widthToFWHM(width, mu = this.mu) {\n        return pseudoVoigtWidthToFWHM(width, mu);\n    }\n    fct(x) {\n        return pseudoVoigtFct(x, this.fwhm, this.mu);\n    }\n    getArea(height = 1) {\n        return getPseudoVoigtArea({ fwhm: this.fwhm, height, mu: this.mu });\n    }\n    getFactor(area) {\n        return getPseudoVoigtFactor(area);\n    }\n    getData(options = {}) {\n        const { length, factor, height = calculatePseudoVoigtHeight({\n            fwhm: this.fwhm,\n            mu: this.mu,\n            area: 1,\n        }), } = options;\n        return getPseudoVoigtData(this, { factor, length, height });\n    }\n    calculateHeight(area = 1) {\n        return calculatePseudoVoigtHeight({ fwhm: this.fwhm, mu: this.mu, area });\n    }\n    getParameters() {\n        return ['fwhm', 'mu'];\n    }\n}\nexport const calculatePseudoVoigtHeight = (options = {}) => {\n    let { fwhm = 1, mu = 0.5, area = 1 } = options;\n    return (2 * area) / (fwhm * (mu * ROOT_PI_OVER_LN2 + (1 - mu) * Math.PI));\n};\nexport const pseudoVoigtFct = (x, fwhm, mu) => {\n    return (1 - mu) * lorentzianFct(x, fwhm) + mu * gaussianFct(x, fwhm);\n};\nexport const pseudoVoigtWidthToFWHM = (width, mu = 0.5) => {\n    return width * (mu * ROOT_2LN2_MINUS_ONE + 1);\n};\nexport const pseudoVoigtFwhmToWidth = (fwhm, mu = 0.5) => {\n    return fwhm / (mu * ROOT_2LN2_MINUS_ONE + 1);\n};\nexport const getPseudoVoigtArea = (options) => {\n    const { fwhm = 500, height = 1, mu = 0.5 } = options;\n    return (fwhm * height * (mu * ROOT_PI_OVER_LN2 + (1 - mu) * Math.PI)) / 2;\n};\nexport const getPseudoVoigtFactor = (area = 0.9999, mu = 0.5) => {\n    return mu < 1 ? getLorentzianFactor(area) : getGaussianFactor(area);\n};\nexport const getPseudoVoigtData = (shape = {}, options = {}) => {\n    let { fwhm = 500, mu = 0.5 } = shape;\n    let { length, factor = getPseudoVoigtFactor(0.999, mu), height = calculatePseudoVoigtHeight({ fwhm, mu, area: 1 }), } = options;\n    if (!height) {\n        height =\n            1 /\n                ((mu / Math.sqrt(-GAUSSIAN_EXP_FACTOR / Math.PI)) * fwhm +\n                    ((1 - mu) * fwhm * Math.PI) / 2);\n    }\n    if (!length) {\n        length = Math.min(Math.ceil(fwhm * factor), Math.pow(2, 25) - 1);\n        if (length % 2 === 0)\n            length++;\n    }\n    const center = (length - 1) / 2;\n    const data = new Float64Array(length);\n    for (let i = 0; i <= center; i++) {\n        data[i] = pseudoVoigtFct(i - center, fwhm, mu) * height;\n        data[length - 1 - i] = data[i];\n    }\n    return data;\n};\n//# sourceMappingURL=PseudoVoigt.js.map","import { ROOT_THREE } from '../../../util/constants';\n/**\n * This shape is a linear combination of rational function (n|n+2), for n = 0 (lorentzian function) and n = 2\n * the parameter that combines those two functions is `gamma` and it is called the kurtosis parameter, it is an\n * implementation of generalized lorentzian shape published by Stanislav Sykora in the SMASH 2010. DOI:10.3247/SL3nmr10.006\n * @link http://www.ebyte.it/stan/Talk_ML_UserMeeting_SMASH_2010_GeneralizedLorentzian.html\n */\nexport class GeneralizedLorentzian {\n    constructor(options = {}) {\n        const { fwhm = 500, gamma = 0.5 } = options;\n        this.fwhm = fwhm;\n        this.gamma = gamma;\n    }\n    fwhmToWidth(fwhm = this.fwhm) {\n        return generalizedLorentzianFwhmToWidth(fwhm);\n    }\n    widthToFWHM(width) {\n        return generalizedLorentzianWidthToFWHM(width);\n    }\n    fct(x) {\n        return generalizedLorentzianFct(x, this.fwhm, this.gamma);\n    }\n    getArea(height = 1) {\n        return getGeneralizedLorentzianArea({\n            fwhm: this.fwhm,\n            height,\n            gamma: this.gamma,\n        });\n    }\n    getFactor(area) {\n        return getGeneralizedLorentzianFactor(area);\n    }\n    getData(options = {}) {\n        return getGeneralizedLorentzianData(this, options);\n    }\n    calculateHeight(area = 1) {\n        const { gamma, fwhm } = this;\n        return calculateGeneralizedLorentzianHeight({ fwhm, area, gamma });\n    }\n    getParameters() {\n        return ['fwhm', 'gamma'];\n    }\n}\nexport const calculateGeneralizedLorentzianHeight = ({ fwhm = 1, gamma = 1, area = 1, }) => {\n    return (area / fwhm / (3.14159 - 0.420894 * gamma)) * 2;\n};\n/**\n * expression of integral generated by Mathematica of the function\n */\nexport const getGeneralizedLorentzianArea = (options) => {\n    const { fwhm = 500, height = 1, gamma = 1 } = options;\n    return (height * fwhm * (3.14159 - 0.420894 * gamma)) / 2;\n};\nexport const generalizedLorentzianFct = (x, fwhm, gamma) => {\n    const u = ((2 * x) / fwhm) ** 2;\n    return (1 - gamma) / (1 + u) + (gamma * (1 + u / 2)) / (1 + u + u ** 2);\n};\nexport const generalizedLorentzianWidthToFWHM = (width) => {\n    return width * ROOT_THREE;\n};\nexport const generalizedLorentzianFwhmToWidth = (fwhm) => {\n    return fwhm / ROOT_THREE;\n};\nexport const getGeneralizedLorentzianFactor = (area = 0.9999) => {\n    if (area >= 1) {\n        throw new Error('area should be (0 - 1)');\n    }\n    const halfResidual = (1 - area) * 0.5;\n    const quantileFunction = (p) => Math.tan(Math.PI * (p - 0.5));\n    return ((quantileFunction(1 - halfResidual) - quantileFunction(halfResidual)) / 2);\n};\nexport const getGeneralizedLorentzianData = (shape = {}, options = {}) => {\n    let { fwhm = 500, gamma = 1 } = shape;\n    let { length, factor = getGeneralizedLorentzianFactor(), height = calculateGeneralizedLorentzianHeight({ fwhm, area: 1, gamma }), } = options;\n    if (!length) {\n        length = Math.min(Math.ceil(fwhm * factor), Math.pow(2, 25) - 1);\n        if (length % 2 === 0)\n            length++;\n    }\n    const center = (length - 1) / 2;\n    const data = new Float64Array(length);\n    for (let i = 0; i <= center; i++) {\n        data[i] = generalizedLorentzianFct(i - center, fwhm, gamma) * height;\n        data[length - 1 - i] = data[i];\n    }\n    return data;\n};\n//# sourceMappingURL=GeneralizedLorentzian.js.map","import { GAUSSIAN_EXP_FACTOR } from '../../../util/constants';\nimport { getGaussianFactor, gaussianFwhmToWidth, gaussianWidthToFWHM, } from '../../1d/gaussian/Gaussian';\nexport class Gaussian2D {\n    constructor(options = {}) {\n        let { fwhm = 20, sd } = options;\n        fwhm = ensureFWHM2D(fwhm, sd);\n        this.fwhmX = fwhm.x;\n        this.fwhmY = fwhm.y;\n    }\n    fct(x, y) {\n        return gaussian2DFct(x, y, this.fwhmX, this.fwhmY);\n    }\n    getData(options = {}) {\n        return getGaussian2DData({\n            fwhm: { x: this.fwhmX, y: this.fwhmY },\n        }, options);\n    }\n    getFactor(volume = 1) {\n        return getGaussianFactor(volume);\n    }\n    getVolume(height = calculateGaussian2DHeight({\n        fwhm: { x: this.fwhmX, y: this.fwhmY },\n        volume: 1,\n    })) {\n        return getGaussian2DVolume({\n            fwhm: { x: this.fwhmX, y: this.fwhmY },\n            height,\n        });\n    }\n    widthToFWHM(width) {\n        return gaussianWidthToFWHM(width);\n    }\n    fwhmToWidth(fwhm) {\n        return gaussianFwhmToWidth(fwhm);\n    }\n    calculateHeight(volume = 1) {\n        return calculateGaussian2DHeight({\n            volume,\n            fwhm: { x: this.fwhmX, y: this.fwhmY },\n        });\n    }\n    set fwhm(fwhm) {\n        fwhm = ensureXYNumber(fwhm);\n        this.fwhmX = fwhm.x;\n        this.fwhmY = fwhm.y;\n    }\n}\nexport const gaussian2DFct = (x, y, xFWHM, yFWHM) => {\n    return Math.exp(GAUSSIAN_EXP_FACTOR * (Math.pow(x / xFWHM, 2) + Math.pow(y / yFWHM, 2)));\n};\nexport const getGaussian2DData = (shape, options = {}) => {\n    let { fwhm = 50, sd } = shape;\n    fwhm = ensureFWHM2D(fwhm, sd);\n    let { factor = getGaussianFactor(), length = { x: 0, y: 0 }, height = calculateGaussian2DHeight({ fwhm, volume: 1 }), } = options;\n    factor = ensureXYNumber(factor);\n    length = ensureXYNumber(length);\n    for (const axis of ['x', 'y']) {\n        if (!length[axis]) {\n            length[axis] = Math.min(Math.ceil(fwhm[axis] * factor[axis]), Math.pow(2, 25) - 1);\n            if (length[axis] % 2 === 0)\n                length[axis]++;\n        }\n    }\n    const xCenter = (length.x - 1) / 2;\n    const yCenter = (length.y - 1) / 2;\n    const data = new Array(length.x);\n    for (let i = 0; i < length.x; i++) {\n        data[i] = new Float64Array(length.y);\n    }\n    for (let i = 0; i < length.x; i++) {\n        for (let j = 0; j < length.y; j++) {\n            data[i][j] =\n                gaussian2DFct(i - xCenter, j - yCenter, fwhm.x, fwhm.y) * height;\n        }\n    }\n    return data;\n};\nexport const calculateGaussian2DHeight = (options = {}) => {\n    let { volume = 1, fwhm = 50, sd } = options;\n    fwhm = ensureFWHM2D(fwhm, sd);\n    return (volume * Math.LN2 * 4) / (Math.PI * fwhm.y * fwhm.x);\n};\nexport const getGaussian2DVolume = (options = {}) => {\n    let { fwhm = 50, height = 1, sd } = options;\n    fwhm = ensureFWHM2D(fwhm, sd);\n    return (height * Math.PI * fwhm.y * fwhm.x) / Math.LN2 / 4;\n};\nfunction ensureXYNumber(input) {\n    return typeof input !== 'object' ? { x: input, y: input } : { ...input };\n}\nfunction ensureFWHM2D(fwhm, sd) {\n    if (sd !== undefined) {\n        let sdObject = ensureXYNumber(sd);\n        return {\n            x: gaussianWidthToFWHM(2 * sdObject.x),\n            y: gaussianWidthToFWHM(2 * sdObject.y),\n        };\n    }\n    else if (fwhm !== undefined) {\n        return ensureXYNumber(fwhm);\n    }\n    else {\n        throw new Error('ensureFWHM2D must have either fwhm or sd defined');\n    }\n}\n//# sourceMappingURL=Gaussian2D.js.map","import { Gaussian } from './gaussian/Gaussian';\nimport { GeneralizedLorentzian } from './generalizedLorentzian/GeneralizedLorentzian';\nimport { Lorentzian } from './lorentzian/Lorentzian';\nimport { PseudoVoigt } from './pseudoVoigt/PseudoVoigt';\n/**\n * Generate a instance of a specific kind of shape.\n */\nexport function getShape1D(shape) {\n    const { kind } = shape;\n    switch (kind) {\n        case 'gaussian':\n            return new Gaussian(shape);\n        case 'lorentzian':\n            return new Lorentzian(shape);\n        case 'pseudoVoigt':\n            return new PseudoVoigt(shape);\n        case 'generalizedLorentzian':\n            return new GeneralizedLorentzian(shape);\n        default: {\n            throw Error(`Unknown distribution ${kind}`);\n        }\n    }\n}\n//# sourceMappingURL=getShape1D.js.map","import { Gaussian2D } from './gaussian2D/Gaussian2D';\n/**\n * Generate a instance of a specific kind of shape.\n */\nexport function getShape2D(shape) {\n    const { kind } = shape;\n    switch (kind) {\n        case 'gaussian':\n            return new Gaussian2D(shape);\n        default: {\n            const unHandled = kind;\n            // eslint-disable-next-line @typescript-eslint/restrict-template-expressions\n            throw Error(`Unknown distribution ${unHandled}`);\n        }\n    }\n}\n//# sourceMappingURL=getShape2D.js.map","// eslint-disable-next-line @typescript-eslint/unbound-method\nconst toString = Object.prototype.toString;\n/**\n * Checks if an object is an instance of an Array (array or typed array, except those that contain bigint values).\n *\n * @param value - Object to check.\n * @returns True if the object is an array or a typed array.\n */\nexport function isAnyArray(value) {\n    const tag = toString.call(value);\n    return tag.endsWith('Array]') && !tag.includes('Big');\n}\n//# sourceMappingURL=index.js.map","/**\n * Returns the closest index of a `target`\n * @param array - array of numbers\n * @param target - target\n * @param options\n * @returns - closest index\n */\nexport function xFindClosestIndex(array, target, options = {}) {\n    const { sorted = true } = options;\n    if (sorted) {\n        let low = 0;\n        let high = array.length - 1;\n        let middle = 0;\n        while (high - low > 1) {\n            middle = low + ((high - low) >> 1);\n            if (array[middle] < target) {\n                low = middle;\n            }\n            else if (array[middle] > target) {\n                high = middle;\n            }\n            else {\n                return middle;\n            }\n        }\n        if (low < array.length - 1) {\n            if (Math.abs(target - array[low]) < Math.abs(array[low + 1] - target)) {\n                return low;\n            }\n            else {\n                return low + 1;\n            }\n        }\n        else {\n            return low;\n        }\n    }\n    else {\n        let index = 0;\n        let diff = Number.POSITIVE_INFINITY;\n        for (let i = 0; i < array.length; i++) {\n            const currentDiff = Math.abs(array[i] - target);\n            if (currentDiff < diff) {\n                diff = currentDiff;\n                index = i;\n            }\n        }\n        return index;\n    }\n}\n//# sourceMappingURL=xFindClosestIndex.js.map","export function createNumberArray(ArrayConstructor, length) {\n    if (ArrayConstructor === Array) {\n        return new ArrayConstructor(length).fill(0);\n    }\n    else {\n        return new ArrayConstructor(length);\n    }\n}\nexport function createDoubleArray(ArrayConstructor, length) {\n    if (ArrayConstructor === Array) {\n        return new ArrayConstructor(length).fill(0);\n    }\n    else {\n        return new ArrayConstructor(length);\n    }\n}\n//# sourceMappingURL=createArray.js.map","const LOOP = 8;\nconst FLOAT_MUL = 1 / 16777216;\nconst sh1 = 15;\nconst sh2 = 18;\nconst sh3 = 11;\nfunction multiplyUint32(n, m) {\n    n >>>= 0;\n    m >>>= 0;\n    const nlo = n & 0xffff;\n    const nhi = n - nlo;\n    return (((nhi * m) >>> 0) + nlo * m) >>> 0;\n}\nexport class XSadd {\n    random;\n    state;\n    /**\n     * create an instance of XSadd with the specified seed\n     * @param [seed=Date.now()]\n     */\n    constructor(seed = Date.now()) {\n        this.state = new Uint32Array(4);\n        this.init(seed);\n        this.random = this.getFloat.bind(this);\n    }\n    /**\n     * Returns a 32-bit integer r (0 <= r < 2^32)\n     */\n    getUint32() {\n        this.nextState();\n        return (this.state[3] + this.state[2]) >>> 0;\n    }\n    /**\n     * Returns a floating point number r (0.0 <= r < 1.0)\n     */\n    getFloat() {\n        return (this.getUint32() >>> 8) * FLOAT_MUL;\n    }\n    init(seed) {\n        if (!Number.isInteger(seed)) {\n            throw new TypeError(\"seed must be an integer\");\n        }\n        this.state[0] = seed;\n        this.state[1] = 0;\n        this.state[2] = 0;\n        this.state[3] = 0;\n        for (let i = 1; i < LOOP; i++) {\n            this.state[i & 3] ^=\n                (i +\n                    multiplyUint32(1812433253, this.state[(i - 1) & 3] ^ ((this.state[(i - 1) & 3] >>> 30) >>> 0))) >>>\n                    0;\n        }\n        this.periodCertification();\n        for (let i = 0; i < LOOP; i++) {\n            this.nextState();\n        }\n    }\n    periodCertification() {\n        if (this.state[0] === 0 &&\n            this.state[1] === 0 &&\n            this.state[2] === 0 &&\n            this.state[3] === 0) {\n            this.state[0] = 88; // X\n            this.state[1] = 83; // S\n            this.state[2] = 65; // A\n            this.state[3] = 68; // D\n        }\n    }\n    nextState() {\n        let t = this.state[0];\n        t ^= t << sh1;\n        t ^= t >>> sh2;\n        t ^= this.state[3] << sh3;\n        this.state[0] = this.state[1];\n        this.state[1] = this.state[2];\n        this.state[2] = this.state[3];\n        this.state[3] = t;\n    }\n}\n//# sourceMappingURL=XSadd.js.map","import { XSadd } from 'ml-xsadd';\n/**\n * Create a random array of numbers of a specific length.\n * @param options\n * @returns - array of random floats normally distributed\n */\nexport function createRandomArray(options = {}) {\n    const { mean = 0, standardDeviation = 1, length = 1000, range = 1, seed, distribution = 'normal', } = options;\n    const generator = new XSadd(seed);\n    const returnArray = new Float64Array(length);\n    switch (distribution) {\n        case 'normal': {\n            const gaussianGenerator = new GaussianGenerator(mean, standardDeviation, generator);\n            for (let i = 0; i < length; i++) {\n                returnArray[i] = gaussianGenerator.generateGaussian();\n            }\n            break;\n        }\n        case 'uniform': {\n            for (let i = 0; i < length; i++) {\n                returnArray[i] = (generator.random() - 0.5) * range + mean;\n            }\n            break;\n        }\n        default: {\n            throw new Error(`unknown distribution: ${String(distribution)}`);\n        }\n    }\n    return returnArray;\n}\nclass GaussianGenerator {\n    #spare = 0;\n    #hasSpare = false;\n    #mean;\n    #standardDeviation;\n    #generator;\n    constructor(mean, standardDeviation, generator) {\n        this.#mean = mean;\n        this.#standardDeviation = standardDeviation;\n        this.#generator = generator;\n    }\n    generateGaussian() {\n        let val, u, v, s;\n        if (this.#hasSpare) {\n            this.#hasSpare = false;\n            val = this.#spare * this.#standardDeviation + this.#mean;\n        }\n        else {\n            do {\n                u = this.#generator.random() * 2 - 1;\n                v = this.#generator.random() * 2 - 1;\n                s = u * u + v * v;\n            } while (s >= 1 || s === 0);\n            s = Math.sqrt((-2 * Math.log(s)) / s);\n            this.#spare = v * s;\n            this.#hasSpare = true;\n            val = this.#mean + this.#standardDeviation * u * s;\n        }\n        return val;\n    }\n}\n//# sourceMappingURL=createRandomArray.js.map","import { xCheck } from './xCheck';\nimport { xGetFromToIndex } from './xGetFromToIndex';\n/**\n * Computes the maximal value of an array of values\n * @param array - array of numbers\n * @param options - options\n */\nexport function xMaxValue(array, options = {}) {\n    xCheck(array);\n    const { fromIndex, toIndex } = xGetFromToIndex(array, options);\n    let maxValue = array[fromIndex];\n    for (let i = fromIndex + 1; i <= toIndex; i++) {\n        if (array[i] > maxValue) {\n            maxValue = array[i];\n        }\n    }\n    return maxValue;\n}\n//# sourceMappingURL=xMaxValue.js.map","import { isAnyArray } from 'is-any-array';\n/**\n * Checks if the input is a non-empty array of numbers.\n * Only checks the first element.\n * @param input - Array to check.\n * @param options - Additional checks.\n */\nexport function xCheck(input, options = {}) {\n    const { minLength } = options;\n    if (!isAnyArray(input)) {\n        throw new TypeError('input must be an array');\n    }\n    if (input.length === 0) {\n        throw new TypeError('input must not be empty');\n    }\n    if (typeof input[0] !== 'number') {\n        throw new TypeError('input must contain numbers');\n    }\n    if (minLength && input.length < minLength) {\n        throw new Error(`input must have a length of at least ${minLength}`);\n    }\n}\n//# sourceMappingURL=xCheck.js.map","import { xFindClosestIndex } from './xFindClosestIndex';\n/**\n * Returns an object with {fromIndex, toIndex} for a specific from / to\n * @param x - array of numbers\n * @param options - Options\n */\nexport function xGetFromToIndex(x, options = {}) {\n    let { fromIndex, toIndex } = options;\n    const { from, to } = options;\n    if (fromIndex === undefined) {\n        if (from !== undefined) {\n            fromIndex = xFindClosestIndex(x, from);\n        }\n        else {\n            fromIndex = 0;\n        }\n    }\n    if (toIndex === undefined) {\n        if (to !== undefined) {\n            toIndex = xFindClosestIndex(x, to);\n        }\n        else {\n            toIndex = x.length - 1;\n        }\n    }\n    if (fromIndex < 0)\n        fromIndex = 0;\n    if (toIndex < 0)\n        toIndex = 0;\n    if (fromIndex >= x.length)\n        fromIndex = x.length - 1;\n    if (toIndex >= x.length)\n        toIndex = x.length - 1;\n    if (fromIndex > toIndex)\n        [fromIndex, toIndex] = [toIndex, fromIndex];\n    return { fromIndex, toIndex };\n}\n//# sourceMappingURL=xGetFromToIndex.js.map","import { createNumberArray } from '../utils';\n/**\n * Create a new matrix based on the size of the current one or by using specific dimensions.\n * @param options\n */\nexport function matrixCreateEmpty(options) {\n    const { matrix, nbRows = matrix?.length || 1, nbColumns = matrix?.[0].length || 1, ArrayConstructor = Float64Array, } = options;\n    const newMatrix = [];\n    for (let row = 0; row < nbRows; row++) {\n        newMatrix.push(createNumberArray(ArrayConstructor, nbColumns));\n    }\n    return newMatrix;\n}\n//# sourceMappingURL=matrixCreateEmpty.js.map","import { matrixCheck } from './matrixCheck';\n/**\n * Get min and max Z.\n * @param matrix - matrix [rows][cols].\n */\nexport function matrixMinMaxZ(matrix) {\n    matrixCheck(matrix);\n    const nbRows = matrix.length;\n    const nbColumns = matrix[0].length;\n    let min = matrix[0][0];\n    let max = matrix[0][0];\n    for (let column = 0; column < nbColumns; column++) {\n        for (let row = 0; row < nbRows; row++) {\n            if (matrix[row][column] < min)\n                min = matrix[row][column];\n            if (matrix[row][column] > max)\n                max = matrix[row][column];\n        }\n    }\n    return { min, max };\n}\n//# sourceMappingURL=matrixMinMaxZ.js.map","export function matrixCheck(data) {\n    if (data.length === 0 || data[0].length === 0) {\n        throw new RangeError('matrix must contain data');\n    }\n    const firstLength = data[0].length;\n    for (let i = 1; i < data.length; i++) {\n        if (data[i].length !== firstLength) {\n            throw new RangeError('all rows must has the same length');\n        }\n    }\n}\n//# sourceMappingURL=matrixCheck.js.map","import { xMaxValue, xAdd, createRandomArray } from 'ml-spectra-processing';\n/**\n * Adds noise to the spectrum data.\n * @param data - The spectrum data to modify.\n * @param options - Configuration for noise generation.\n * @returns The modified spectrum data.\n */\nexport default function addNoise(data, options = {}) {\n    const { seed = 0, distribution = 'normal', percent = 1 } = options;\n    const range = (xMaxValue(data.y) * percent) / 100;\n    const noise = createRandomArray({\n        distribution,\n        seed,\n        mean: 0,\n        standardDeviation: range,\n        range,\n        length: data.x.length,\n    });\n    data.y = xAdd(data.y, noise);\n    return data;\n}\n//# sourceMappingURL=addNoise.js.map","import { isAnyArray } from 'is-any-array';\n/**\n * This function xAdd the first array by the second array or a constant value to each element of the first array\n * @param array1 - the first array\n * @param array2 - the second array or number\n * @returns the result of the addition\n */\nexport function xAdd(array1, array2) {\n    let isConstant = false;\n    let constant = 0;\n    if (isAnyArray(array2)) {\n        if (array1.length !== array2.length) {\n            throw new Error('size of array1 and array2 must be identical');\n        }\n    }\n    else {\n        isConstant = true;\n        constant = array2;\n    }\n    const array3 = new Float64Array(array1.length);\n    if (isConstant) {\n        for (let i = 0; i < array1.length; i++) {\n            array3[i] = array1[i] + constant;\n        }\n    }\n    else {\n        for (let i = 0; i < array1.length; i++) {\n            array3[i] = array1[i] + array2[i];\n        }\n    }\n    return array3;\n}\n//# sourceMappingURL=xAdd.js.map","import { getShape1D } from 'ml-peak-shape-generator';\nimport addBaseline from './util/addBaseline';\nimport addNoise from './util/addNoise';\nexport class SpectrumGenerator {\n    constructor(options = {}) {\n        const { from = 0, to = 1000, nbPoints = 10001, peakWidthFct, shape = { kind: 'gaussian', fwhm: 5 }, } = options;\n        this.from = from;\n        this.to = to;\n        this.nbPoints = nbPoints;\n        this.interval = (this.to - this.from) / (this.nbPoints - 1);\n        this.peakWidthFct = peakWidthFct;\n        this.maxPeakHeight = Number.MIN_SAFE_INTEGER;\n        this.data = {\n            x: new Float64Array(this.nbPoints),\n            y: new Float64Array(this.nbPoints),\n        };\n        const shapeGenerator = getShape1D(shape);\n        this.shape = shapeGenerator;\n        assertNumber(this.from, 'from');\n        assertNumber(this.to, 'to');\n        assertInteger(this.nbPoints, 'nbPoints');\n        if (this.to <= this.from) {\n            throw new RangeError('to option must be larger than from');\n        }\n        if (this.peakWidthFct && typeof this.peakWidthFct !== 'function') {\n            throw new TypeError('peakWidthFct option must be a function');\n        }\n        this.reset();\n    }\n    /**\n     * Add a series of peaks to the spectrum.\n     * @param peaks - Peaks to add.\n     * @param options - Options for adding peaks.\n     * @returns The generator instance.\n     */\n    addPeaks(peaks, options) {\n        if (!Array.isArray(peaks) &&\n            (typeof peaks !== 'object' ||\n                peaks.x === undefined ||\n                peaks.y === undefined ||\n                !Array.isArray(peaks.x) ||\n                !Array.isArray(peaks.y) ||\n                peaks.x.length !== peaks.y.length)) {\n            throw new TypeError('peaks must be an array or an object containing x[] and y[]');\n        }\n        if (Array.isArray(peaks)) {\n            for (const peak of peaks) {\n                this.addPeak(peak, options);\n            }\n        }\n        else {\n            for (let i = 0; i < peaks.x.length; i++) {\n                this.addPeak([peaks.x[i], peaks.y[i]], options);\n            }\n        }\n        return this;\n    }\n    /**\n     * Add a single peak to the spectrum.\n     * A peak may be either defined as [x,y,fwhm,...] or as {x, y, shape}\n     * @param peak - The peak to add, defined as array or object.\n     * @param options - Options for adding the peak.\n     * @returns The generator instance.\n     */\n    addPeak(peak, options = {}) {\n        if (Array.isArray(peak) && peak.length < 2) {\n            throw new Error('peak must be an array with two (or three) values or an object with {x,y,width?}');\n        }\n        if (!Array.isArray(peak) &&\n            (peak.x === undefined || peak.y === undefined)) {\n            throw new Error('peak must be an array with two (or three) values or an object with {x,y,width?}');\n        }\n        let xPosition;\n        let intensity;\n        let peakFWHM;\n        let peakWidth;\n        let peakShapeOptions;\n        if (Array.isArray(peak)) {\n            [xPosition, intensity, peakFWHM, peakShapeOptions] = peak;\n        }\n        else {\n            xPosition = peak.x;\n            intensity = peak.y;\n            peakWidth = peak.width;\n            peakShapeOptions = peak.shape;\n        }\n        if (intensity > this.maxPeakHeight)\n            this.maxPeakHeight = intensity;\n        let { shape: shapeOptions } = options;\n        if (peakShapeOptions) {\n            shapeOptions = shapeOptions\n                ? { ...shapeOptions, ...peakShapeOptions }\n                : peakShapeOptions;\n        }\n        const shape = shapeOptions\n            ? getShape1D(shapeOptions)\n            : Object.assign(Object.create(Object.getPrototypeOf(this.shape)), this.shape);\n        const { width } = options;\n        let { widthLeft, widthRight } = options;\n        /*\n         if we don't force the fwhm we just take the one from the shape\n         however we have many way to force it:\n         - use [x,y,fwhm]\n         - define `width` that will be converted to fwhm\n         - define `widthLeft` and `widthRight` to define asymmetric peaks\n         - have a callback `peakWidthFct`\n         This should evolve in the future because we will not always have `fwhm`\n         */\n        const fwhm = peakFWHM !== undefined\n            ? peakFWHM\n            : peakWidth\n                ? shape.widthToFWHM(peakWidth)\n                : this.peakWidthFct\n                    ? this.peakWidthFct(xPosition)\n                    : width !== undefined\n                        ? width\n                        : shape.fwhm;\n        if (!widthLeft)\n            widthLeft = fwhm;\n        if (!widthRight)\n            widthRight = fwhm;\n        if (!widthLeft || !widthRight) {\n            throw new Error('Width left or right is undefined or zero');\n        }\n        const factor = options.factor === undefined ? shape.getFactor() : options.factor;\n        const firstValue = xPosition - (widthLeft / 2) * factor;\n        const lastValue = xPosition + (widthRight / 2) * factor;\n        const firstPoint = Math.max(0, Math.floor((firstValue - this.from) / this.interval));\n        const lastPoint = Math.min(this.nbPoints - 1, Math.ceil((lastValue - this.from) / this.interval));\n        const middlePoint = Math.round((xPosition - this.from) / this.interval);\n        // PEAK SHAPE MAY BE ASYMMETRC (widthLeft and widthRight) !\n        // we calculate the left part of the shape\n        shape.fwhm = widthLeft;\n        for (let index = firstPoint; index < Math.max(middlePoint, 0); index++) {\n            this.data.y[index] +=\n                intensity * shape.fct(this.data.x[index] - xPosition);\n        }\n        // we calculate the right part of the gaussian\n        shape.fwhm = widthRight;\n        for (let index = Math.min(middlePoint, lastPoint); index <= lastPoint; index++) {\n            this.data.y[index] +=\n                intensity * shape.fct(this.data.x[index] - xPosition);\n        }\n        return this;\n    }\n    /**\n     * Add a baseline to the spectrum.\n     * @param baselineFct - Mathematical function producing the baseline you want.\n     * @returns The generator instance.\n     */\n    addBaseline(baselineFct) {\n        addBaseline(this.data, baselineFct);\n        return this;\n    }\n    /**\n     * Add noise to the spectrum.\n     * @param options - Configuration for noise generation.\n     * @returns The generator instance.\n     */\n    addNoise(options) {\n        addNoise(this.data, options);\n        return this;\n    }\n    /**\n     * Get the generated spectrum.\n     * @param options - Options for getting the spectrum.\n     * @returns The generated spectrum data.\n     */\n    getSpectrum(options = {}) {\n        if (typeof options === 'boolean') {\n            options = { copy: options };\n        }\n        const { copy = true, threshold = 0 } = options;\n        if (threshold) {\n            const minPeakHeight = this.maxPeakHeight * threshold;\n            const x = [];\n            const y = [];\n            for (let i = 0; i < this.data.x.length; i++) {\n                if (this.data.y[i] >= minPeakHeight) {\n                    x.push(this.data.x[i]);\n                    y.push(this.data.y[i]);\n                }\n            }\n            return { x: Float64Array.from(x), y: Float64Array.from(y) };\n        }\n        if (copy) {\n            return {\n                x: this.data.x.slice(),\n                y: this.data.y.slice(),\n            };\n        }\n        else {\n            return this.data;\n        }\n    }\n    /**\n     * Resets the generator with an empty spectrum.\n     * @returns The generator instance.\n     */\n    reset() {\n        if (this.nbPoints === 1) {\n            this.data.x[0] = (this.from + this.to) / 2;\n        }\n        else {\n            for (let i = 0; i < this.nbPoints; i++) {\n                this.data.x[i] = this.from + i * this.interval;\n            }\n        }\n        return this;\n    }\n}\nfunction assertInteger(value, name) {\n    if (!Number.isInteger(value)) {\n        throw new TypeError(`${name} option must be an integer`);\n    }\n}\nfunction assertNumber(value, name) {\n    if (!Number.isFinite(value)) {\n        throw new TypeError(`${name} option must be a number`);\n    }\n}\n/**\n * Generates a spectrum and returns it.\n * @param peaks - List of peaks to put in the spectrum.\n * @param options - Configuration for spectrum generation.\n * @returns The generated spectrum data.\n */\nexport function generateSpectrum(peaks, options = {}) {\n    const { generator: generatorOptions, noise, baseline, threshold, peakOptions, } = options;\n    const generator = new SpectrumGenerator(generatorOptions);\n    generator.addPeaks(peaks, peakOptions);\n    if (baseline)\n        generator.addBaseline(baseline);\n    if (noise) {\n        generator.addNoise(noise);\n    }\n    return generator.getSpectrum({\n        threshold,\n    });\n}\n//# sourceMappingURL=SpectrumGenerator.js.map","/**\n * Adds a baseline to the spectrum data.\n * @param data - The spectrum data to modify.\n * @param baselineFct - Function that generates the baseline value for a given x.\n * @returns The modified spectrum data.\n */\nexport default function addBaseline(data, baselineFct) {\n    if (!baselineFct)\n        return data;\n    const xs = data.x;\n    const ys = data.y;\n    for (let i = 0; i < xs.length; i++) {\n        ys[i] += baselineFct(xs[i]);\n    }\n    return data;\n}\n//# sourceMappingURL=addBaseline.js.map","import { getShape2D } from 'ml-peak-shape-generator';\nimport { matrixCreateEmpty, matrixMinMaxZ } from 'ml-spectra-processing';\nconst axis2D = ['x', 'y'];\nconst peakCoordinates = ['x', 'y', 'z'];\n/**\n * Converts a width value to full width at half maximum (FWHM).\n * @param shape - The 2D shape to use for conversion.\n * @param width - The width value to convert.\n * @returns The converted FWHM value as an XYNumber object.\n */\nconst convertWidthToFWHM = (shape, width) => {\n    const widthData = ensureXYNumber(width);\n    for (const key of axis2D) {\n        widthData[key] = shape.widthToFWHM(widthData[key]);\n    }\n    return widthData;\n};\nexport class Spectrum2DGenerator {\n    constructor(options = {}) {\n        const { peakWidthFct = () => 5, shape = {\n            kind: 'gaussian',\n        }, } = options;\n        let { from = 0, to = 100, nbPoints = 1001 } = options;\n        from = ensureXYNumber(from);\n        to = ensureXYNumber(to);\n        nbPoints = ensureXYNumber(nbPoints);\n        for (const axis of axis2D) {\n            assertNumber(from[axis], `from-${axis}`);\n            assertNumber(to[axis], `to-${axis}`);\n            assertInteger(nbPoints[axis], `nbPoints-${axis}`);\n        }\n        this.from = from;\n        this.to = to;\n        this.nbPoints = nbPoints;\n        this.interval = calculeIntervals(from, to, nbPoints);\n        this.peakWidthFct = peakWidthFct;\n        this.maxPeakHeight = Number.MIN_SAFE_INTEGER;\n        const shapeGenerator = getShape2D(shape);\n        this.shape = shapeGenerator;\n        this.data = {\n            x: new Float64Array(nbPoints.x),\n            y: new Float64Array(nbPoints.y),\n            z: matrixCreateEmpty({\n                nbRows: this.nbPoints.y,\n                nbColumns: this.nbPoints.x,\n            }),\n        };\n        for (const axis of axis2D) {\n            if (this.to[axis] <= this.from[axis]) {\n                throw new RangeError('to option must be larger than from');\n            }\n        }\n        if (typeof this.peakWidthFct !== 'function') {\n            throw new TypeError('peakWidthFct option must be a function');\n        }\n        this.reset();\n    }\n    /**\n     * Adds multiple peaks to the 2D spectrum.\n     * @param peaks - Array of peaks or peak series to add.\n     * @param options - Options for adding peaks.\n     * @returns The generator instance.\n     */\n    addPeaks(peaks, options) {\n        if (!Array.isArray(peaks) &&\n            (typeof peaks !== 'object' ||\n                peaks.x === undefined ||\n                peaks.y === undefined ||\n                !Array.isArray(peaks.x) ||\n                !Array.isArray(peaks.y) ||\n                peaks.x.length !== peaks.y.length)) {\n            throw new TypeError('peaks must be an array or an object containing x[] and y[]');\n        }\n        if (Array.isArray(peaks)) {\n            for (const peak of peaks) {\n                this.addPeak(peak, options);\n            }\n        }\n        else {\n            const nbPeaks = peaks.x.length;\n            for (const c of peakCoordinates) {\n                if (peaks[c] &&\n                    Array.isArray(peaks[c]) &&\n                    nbPeaks !== peaks[c].length) {\n                    throw new Error('x, y, z should have the same length');\n                }\n            }\n            for (let i = 0; i < peaks.x.length; i++) {\n                this.addPeak([peaks.x[i], peaks.y[i], peaks.z[i]], options);\n            }\n        }\n        return this;\n    }\n    /**\n     * Adds a single peak to the 2D spectrum.\n     * @param peak - Peak to add, can be array or object format.\n     * @param options - Options for adding the peak.\n     * @returns The generator instance.\n     */\n    addPeak(peak, options = {}) {\n        if (Array.isArray(peak) && peak.length < 3) {\n            throw new Error('peak must be an array with three (or four) values or an object with {x,y,z,width?}');\n        }\n        if (!Array.isArray(peak) &&\n            peakCoordinates.some((e) => peak[e] === undefined)) {\n            throw new Error('peak must be an array with three (or four) values or an object with {x,y,z,width?}');\n        }\n        let xPosition;\n        let yPosition;\n        let intensity;\n        let peakFWHM;\n        let peakWidth;\n        let peakShapeOptions;\n        if (Array.isArray(peak)) {\n            [xPosition, yPosition, intensity, peakFWHM, peakShapeOptions] = peak;\n        }\n        else {\n            xPosition = peak.x;\n            yPosition = peak.y;\n            intensity = peak.z;\n            peakFWHM = peak.fwhm;\n            peakWidth = peak.width;\n            peakShapeOptions = peak.shape;\n        }\n        const position = { x: xPosition, y: yPosition };\n        if (intensity > this.maxPeakHeight)\n            this.maxPeakHeight = intensity;\n        const { width } = options;\n        let { shape: shapeOptions } = options;\n        if (peakShapeOptions) {\n            shapeOptions = shapeOptions\n                ? { ...shapeOptions, ...peakShapeOptions }\n                : peakShapeOptions;\n        }\n        const shape = shapeOptions\n            ? getShape2D(shapeOptions)\n            : Object.assign(Object.create(Object.getPrototypeOf(this.shape)), structuredClone(this.shape));\n        let { fwhm = peakFWHM !== undefined\n            ? peakFWHM\n            : peakWidth\n                ? convertWidthToFWHM(shape, peakWidth)\n                : width\n                    ? convertWidthToFWHM(shape, width)\n                    : this.peakWidthFct(xPosition, yPosition), } = options;\n        fwhm = ensureXYNumber(fwhm);\n        let factor = options.factor === undefined ? shape.getFactor() : options.factor;\n        factor = ensureXYNumber(factor);\n        const firstPoint = { x: 0, y: 0 };\n        const lastPoint = { x: 0, y: 0 };\n        for (const axis of axis2D) {\n            const first = position[axis] - (fwhm[axis] / 2) * factor[axis];\n            const last = position[axis] + (fwhm[axis] / 2) * factor[axis];\n            firstPoint[axis] = Math.max(0, Math.floor((first - this.from[axis]) / this.interval[axis]));\n            lastPoint[axis] = Math.min(this.nbPoints[axis], Math.ceil((last - this.from[axis]) / this.interval[axis]));\n        }\n        shape.fwhm = fwhm;\n        for (let xIndex = firstPoint.x; xIndex < lastPoint.x; xIndex++) {\n            for (let yIndex = firstPoint.y; yIndex < lastPoint.y; yIndex++) {\n                const value = intensity *\n                    shape.fct(this.data.x[xIndex] - position.x, this.data.y[yIndex] - position.y);\n                if (Math.abs(value) > 1e-6) {\n                    this.data.z[yIndex][xIndex] += value;\n                }\n            }\n        }\n        return this;\n    }\n    /**\n     * Gets the generated 2D spectrum data.\n     * @param options - Options for getting the spectrum.\n     * @returns The spectrum data object.\n     */\n    getSpectrum(options = {}) {\n        if (typeof options === 'boolean') {\n            options = { copy: options };\n        }\n        const { copy = true } = options;\n        const minMaxZ = matrixMinMaxZ(this.data.z);\n        return {\n            minX: this.from.x,\n            maxX: this.to.x,\n            maxY: this.to.y,\n            minY: this.from.y,\n            minZ: minMaxZ.min,\n            maxZ: minMaxZ.max,\n            z: copy ? this.data.z.slice() : this.data.z,\n        };\n    }\n    /**\n     * Resets the generator to initial state.\n     * @returns The generator instance.\n     */\n    reset() {\n        const spectrum = this.data;\n        for (const axis of axis2D) {\n            for (let i = 0; i < this.nbPoints[axis]; i++) {\n                spectrum[axis][i] = this.from[axis] + i * this.interval[axis];\n            }\n        }\n        for (const row of spectrum.z) {\n            for (let j = 0; j < row.length; j++) {\n                row[j] = 0;\n            }\n        }\n        return this;\n    }\n}\n/**\n * Generates a 2D spectrum with the given peaks.\n * @param peaks - Peaks to include in the spectrum.\n * @param options - Options for spectrum generation.\n * @returns The generated spectrum data.\n */\nexport function generateSpectrum2D(peaks, options = {}) {\n    const { generator: generatorOptions, peaks: addPeaksOptions } = options;\n    const generator = new Spectrum2DGenerator(generatorOptions);\n    generator.addPeaks(peaks, addPeaksOptions);\n    return generator.getSpectrum();\n}\n/**\n * Ensures the input is an XYNumber object.\n * @param input - Number or XYNumber to process.\n * @returns An XYNumber object.\n */\nfunction ensureXYNumber(input) {\n    return typeof input !== 'object' ? { x: input, y: input } : { ...input };\n}\n/**\n * Calculates the intervals between points for both x and y axes.\n * @param from - Starting point coordinates.\n * @param to - Ending point coordinates.\n * @param nbPoints - Number of points in each dimension.\n * @returns The calculated intervals as an XYNumber object.\n */\nfunction calculeIntervals(from, to, nbPoints) {\n    return {\n        x: (to.x - from.x) / (nbPoints.x - 1),\n        y: (to.y - from.y) / (nbPoints.y - 1),\n    };\n}\n/**\n * Validates if a value is an integer.\n * @param value - Number to validate.\n * @param name - Name of the parameter for error message.\n * @throws {TypeError} If value is not an integer.\n */\nfunction assertInteger(value, name) {\n    if (!Number.isInteger(value)) {\n        throw new TypeError(`${name} option must be an integer`);\n    }\n}\n/**\n * Validates if a value is a finite number.\n * @param value - Number to validate.\n * @param name - Name of the parameter for error message.\n * @throws {TypeError} If value is not a finite number.\n */\nfunction assertNumber(value, name) {\n    if (!Number.isFinite(value)) {\n        throw new TypeError(`${name} option must be a number`);\n    }\n}\n//# sourceMappingURL=Spectrum2DGenerator.js.map"],"names":["GAUSSIAN_EXP_FACTOR","Math","LN2","ROOT_PI_OVER_LN2","sqrt","PI","ROOT_THREE","ROOT_2LN2","ROOT_2LN2_MINUS_ONE","Gaussian","constructor","options","fwhm","sd","this","gaussianWidthToFWHM","fwhmToWidth","gaussianFwhmToWidth","widthToFWHM","width","fct","x","gaussianFct","getArea","height","calculateGaussianHeight","getGaussianArea","getFactor","area","getGaussianFactor","getData","shape","length","factor","min","ceil","pow","center","data","Float64Array","i","getGaussianData","calculateHeight","getParameters","exp","ln1MinusXSqrd","log","lnEtcBy2Plus2","firstSqrt","erfinv","Lorentzian","lorentzianFwhmToWidth","lorentzianWidthToFWHM","lorentzianFct","getLorentzianArea","getLorentzianFactor","getLorentzianData","calculateLorentzianHeight","Error","halfResidual","quantileFunction","p","tan","PseudoVoigt","mu","pseudoVoigtFwhmToWidth","pseudoVoigtWidthToFWHM","pseudoVoigtFct","getPseudoVoigtArea","getPseudoVoigtFactor","calculatePseudoVoigtHeight","getPseudoVoigtData","GeneralizedLorentzian","gamma","generalizedLorentzianFwhmToWidth","generalizedLorentzianWidthToFWHM","generalizedLorentzianFct","getGeneralizedLorentzianArea","getGeneralizedLorentzianFactor","getGeneralizedLorentzianData","calculateGeneralizedLorentzianHeight","u","Gaussian2D","ensureFWHM2D","fwhmX","fwhmY","y","gaussian2DFct","getGaussian2DData","volume","getVolume","calculateGaussian2DHeight","getGaussian2DVolume","ensureXYNumber","xFWHM","yFWHM","axis","xCenter","yCenter","Array","j","input","undefined","sdObject","getShape1D","kind","getShape2D","toString","Object","prototype","isAnyArray","value","tag","call","endsWith","includes","xFindClosestIndex","array","target","sorted","low","high","middle","abs","index","diff","Number","POSITIVE_INFINITY","currentDiff","createNumberArray","ArrayConstructor","fill","FLOAT_MUL","multiplyUint32","n","m","nlo","XSadd","random","state","seed","Date","now","Uint32Array","init","getFloat","bind","getUint32","nextState","isInteger","TypeError","periodCertification","t","GaussianGenerator","spare","hasSpare","mean","standardDeviation","generator","generateGaussian","val","v","s","xMaxValue","minLength","xCheck","fromIndex","toIndex","from","to","xGetFromToIndex","maxValue","matrixCreateEmpty","matrix","nbRows","nbColumns","newMatrix","row","push","matrixMinMaxZ","RangeError","firstLength","matrixCheck","max","column","addNoise","distribution","percent","range","noise","returnArray","gaussianGenerator","String","createRandomArray","array1","array2","isConstant","constant","array3","xAdd","SpectrumGenerator","nbPoints","peakWidthFct","interval","maxPeakHeight","MIN_SAFE_INTEGER","shapeGenerator","assertNumber","name","assertInteger","reset","addPeaks","peaks","isArray","peak","addPeak","xPosition","intensity","peakFWHM","peakWidth","peakShapeOptions","shapeOptions","assign","create","getPrototypeOf","widthLeft","widthRight","firstValue","lastValue","firstPoint","floor","lastPoint","middlePoint","round","addBaseline","baselineFct","xs","ys","getSpectrum","copy","threshold","minPeakHeight","slice","isFinite","axis2D","peakCoordinates","convertWidthToFWHM","widthData","key","Spectrum2DGenerator","calculeIntervals","z","nbPeaks","c","some","e","yPosition","position","structuredClone","first","last","xIndex","yIndex","minMaxZ","minX","maxX","maxY","minY","minZ","maxZ","spectrum","generatorOptions","baseline","peakOptions","addPeaksOptions"],"mappings":"yPAAO,MAAMA,GAAsB,EAAKC,KAAKC,IAChCC,EAAmBF,KAAKG,KAAKH,KAAKI,GAAKJ,KAAKC,KAC5CI,EAAaL,KAAKG,KAAK,GACvBG,EAAYN,KAAKG,KAAK,EAAIH,KAAKC,KAC/BM,EAAsBP,KAAKG,KAAK,EAAIH,KAAKC,KAAO,ECiDvD,MAAOO,EAOXC,WAAAA,CAAmBC,EAAgC,IACjD,MAAMC,KAAEA,EAAO,IAAGC,GAAEA,GAAOF,EAE3BG,KAAKF,KAAOC,EAAKE,EAAoB,EAAIF,GAAMD,CACjD,CAEOI,WAAAA,CAAYJ,EAAOE,KAAKF,MAC7B,OAAOK,EAAoBL,EAC7B,CAEOM,WAAAA,CAAYC,GACjB,OAAOJ,EAAoBI,EAC7B,CAEOC,GAAAA,CAAIC,GACT,OAAOC,EAAYD,EAAGP,KAAKF,KAC7B,CAEOW,OAAAA,CAAQC,EAASC,EAAwB,CAAEb,KAAME,KAAKF,QAC3D,OAiDE,SAA0BD,GAC9B,IAAIC,KAAEA,EAAO,IAAGC,GAAEA,EAAEW,OAAEA,EAAS,GAAMb,EAEjCE,IAAID,EAAOG,EAAoB,EAAIF,IAEvC,OAAQW,EAASrB,EAAmBS,EAAQ,CAC9C,CAvDWc,CAAgB,CAAEd,KAAME,KAAKF,KAAMY,UAC5C,CAEOG,SAAAA,CAAUC,GACf,OAAOC,EAAkBD,EAC3B,CAEOE,OAAAA,CAAQnB,EAA4B,IACzC,OAqDE,SACJoB,EAA8B,GAC9BpB,EAA4B,CAAA,GAE5B,IAAIC,KAAEA,EAAO,IAAGC,GAAEA,GAAOkB,EACrBlB,IAAID,EAAOG,EAAoB,EAAIF,IAEvC,IAAImB,OACFA,EAAMC,OACNA,EAASJ,IAAmBL,OAC5BA,EAASC,EAAwB,CAAEb,UACjCD,EAECqB,IACHA,EAAS/B,KAAKiC,IAAIjC,KAAKkC,KAAKvB,EAAOqB,GAAShC,KAAKmC,IAAI,EAAG,IAAM,GAC1DJ,EAAS,GAAM,GAAGA,KAGxB,MAAMK,GAAUL,EAAS,GAAK,EACxBM,EAAO,IAAIC,aAAaP,GAC9B,IAAK,IAAIQ,EAAI,EAAGA,GAAKH,EAAQG,IAC3BF,EAAKE,GAAKlB,EAAYkB,EAAIH,EAAQzB,GAAQY,EAC1Cc,EAAKN,EAAS,EAAIQ,GAAKF,EAAKE,GAG9B,OAAOF,CACT,CA/EWG,CAAgB3B,KAAMH,EAC/B,CAEO+B,eAAAA,CAAgBd,EAAO,GAC5B,OAAOH,EAAwB,CAAEb,KAAME,KAAKF,KAAMgB,QACpD,CAEOe,aAAAA,GACL,MAAO,CAAC,OACV,EAGI,SAAUlB,EACdd,GAEA,IAAIC,KAAEA,EAAO,IAAGgB,KAAEA,EAAO,EAACf,GAAEA,GAAOF,EAInC,OAFIE,IAAID,EAAOG,EAAoB,EAAIF,IAE/B,EAAIe,EAAQzB,EAAmBS,CACzC,CASM,SAAUU,EAAYD,EAAWT,GACrC,OAAOX,KAAK2C,IAAI5C,EAAsBC,KAAKmC,IAAIf,EAAIT,EAAM,GAC3D,CAEM,SAAUG,EAAoBI,GAClC,OAAOA,EAAQZ,CACjB,CAEM,SAAUU,EAAoBL,GAClC,OAAOA,EAAOL,CAChB,CAUM,SAAUsB,EAAkBD,EAAO,OACvC,OAAO3B,KAAKG,KAAK,GCpIL,SAAiBiB,GAE7B,GAAU,IAANA,EAAS,OAAO,EACpB,IAAIwB,EAAgB5C,KAAK6C,IAAI,EAAIzB,EAAIA,GACjC0B,EAAgBF,EAAgB,EAAI,GAHhC,KAGqC5C,KAAKI,IAC9C2C,EAAY/C,KAAKG,KAAK2C,GAAiB,EAAIF,EAJvC,MAMR,OADiB5C,KAAKG,KAAK4C,EAAYD,IAClB1B,EAAI,EAAI,GAAI,EACnC,CD4HwB4B,CAAOrB,EAC/B,CEjHM,MAAOsB,EAOXxC,WAAAA,CAAmBC,EAAkC,IACnD,MAAMC,KAAEA,EAAO,KAAQD,EAEvBG,KAAKF,KAAOA,CACd,CAEOI,WAAAA,CAAYJ,EAAOE,KAAKF,MAC7B,OAAOuC,EAAsBvC,EAC/B,CAEOM,WAAAA,CAAYC,GACjB,OAAOiC,EAAsBjC,EAC/B,CAEOC,GAAAA,CAAIC,GACT,OAAOgC,EAAchC,EAAGP,KAAKF,KAC/B,CAEOW,OAAAA,CAAQC,EAAS,GACtB,OAAO8B,EAAkB,CAAE1C,KAAME,KAAKF,KAAMY,UAC9C,CAEOG,SAAAA,CAAUC,GACf,OAAO2B,EAAoB3B,EAC7B,CAEOE,OAAAA,CAAQnB,EAA4B,IACzC,OAAO6C,EAAkB1C,KAAMH,EACjC,CAEO+B,eAAAA,CAAgBd,EAAO,GAC5B,OAAO6B,EAA0B,CAAE7C,KAAME,KAAKF,KAAMgB,QACtD,CAEOe,aAAAA,GACL,MAAO,CAAC,OACV,EAGK,MAAMc,EAA4BA,EAAG7C,OAAO,EAAGgB,OAAO,KACnD,EAAIA,EAAQ3B,KAAKI,GAAKO,EAGnB0C,EAAqB3C,IAChC,MAAMC,KAAEA,EAAO,IAAGY,OAAEA,EAAS,GAAMb,EACnC,OAAQa,EAASvB,KAAKI,GAAKO,EAAQ,CAAC,EAGzByC,EAAgBA,CAAChC,EAAWT,IAChCA,GAAQ,GAAK,EAAIS,GAAK,EAAIT,GAAQ,GAG9BwC,EAAyBjC,GAC7BA,EAAQb,EAGJ6C,EAAyBvC,GAC7BA,EAAON,EAGHiD,EAAsBA,CAAC3B,EAAO,SACzC,GAAIA,GAAQ,EACV,MAAM,IAAI8B,MAAM,0BAElB,MAAMC,EAA4B,IAAZ,EAAI/B,GACpBgC,EAAoBC,GAAc5D,KAAK6D,IAAI7D,KAAKI,IAAMwD,EAAI,KAChE,OACGD,EAAiB,EAAID,GAAgBC,EAAiBD,IAAiB,CAAC,EAIhEH,EAAoBA,CAC/BzB,EAAgC,GAChCpB,EAA4B,CAAA,KAE5B,IAAIC,KAAEA,EAAO,KAAQmB,GACjBC,OACFA,EAAMC,OACNA,EAASsB,IAAqB/B,OAC9BA,EAASiC,EAA0B,CAAE7C,OAAMgB,KAAM,KAC/CjB,EAECqB,IACHA,EAAS/B,KAAKiC,IAAIjC,KAAKkC,KAAKvB,EAAOqB,GAAShC,KAAKmC,IAAI,EAAG,IAAM,GAC1DJ,EAAS,GAAM,GAAGA,KAGxB,MAAMK,GAAUL,EAAS,GAAK,EACxBM,EAAO,IAAIC,aAAaP,GAC9B,IAAK,IAAIQ,EAAI,EAAGA,GAAKH,EAAQG,IAC3BF,EAAKE,GAAKa,EAAcb,EAAIH,EAAQzB,GAAQY,EAC5Cc,EAAKN,EAAS,EAAIQ,GAAKF,EAAKE,GAG9B,OAAOF,CAAI,ECtEP,MAAOyB,EAQXrD,WAAAA,CAAmBC,EAAmC,IACpD,MAAMC,KAAEA,EAAO,IAAGoD,GAAEA,EAAK,IAAQrD,EAEjCG,KAAKkD,GAAKA,EACVlD,KAAKF,KAAOA,CACd,CAEOI,WAAAA,CAAYJ,EAAOE,KAAKF,KAAMoD,EAAKlD,KAAKkD,IAC7C,OAAOC,EAAuBrD,EAAMoD,EACtC,CAEO9C,WAAAA,CAAYC,EAAe6C,EAAalD,KAAKkD,IAClD,OAAOE,EAAuB/C,EAAO6C,EACvC,CAEO5C,GAAAA,CAAIC,GACT,OAAO8C,EAAe9C,EAAGP,KAAKF,KAAME,KAAKkD,GAC3C,CAEOzC,OAAAA,CAAQC,EAAS,GACtB,OAAO4C,EAAmB,CAAExD,KAAME,KAAKF,KAAMY,SAAQwC,GAAIlD,KAAKkD,IAChE,CAEOrC,SAAAA,CAAUC,GACf,OAAOyC,EAAqBzC,EAC9B,CAEOE,OAAAA,CAAQnB,EAA4B,IACzC,MAAMqB,OACJA,EAAMC,OACNA,EAAMT,OACNA,EAAS8C,EAA2B,CAClC1D,KAAME,KAAKF,KACXoD,GAAIlD,KAAKkD,GACTpC,KAAM,KAENjB,EACJ,OAAO4D,EAAmBzD,KAAM,CAAEmB,SAAQD,SAAQR,UACpD,CAEOkB,eAAAA,CAAgBd,EAAO,GAC5B,OAAO0C,EAA2B,CAAE1D,KAAME,KAAKF,KAAMoD,GAAIlD,KAAKkD,GAAIpC,QACpE,CAEOe,aAAAA,GACL,MAAO,CAAC,OAAQ,KAClB,EAGK,MAAM2B,EAA6BA,CACxC3D,EAA8C,MAE9C,IAAIC,KAAEA,EAAO,EAACoD,GAAEA,EAAK,GAAGpC,KAAEA,EAAO,GAAMjB,EACvC,OAAQ,EAAIiB,GAAShB,GAAQoD,EAAK7D,GAAoB,EAAI6D,GAAM/D,KAAKI,IAAI,EAG9D8D,EAAiBA,CAAC9C,EAAWT,EAAcoD,KAC9C,EAAIA,GAAMX,EAAchC,EAAGT,GAAQoD,EAAK1C,EAAYD,EAAGT,GAGpDsD,EAAyBA,CAAC/C,EAAe6C,EAAK,KAClD7C,GAAS6C,EAAKxD,EAAsB,GAGhCyD,EAAyBA,CAACrD,EAAcoD,EAAK,KACjDpD,GAAQoD,EAAKxD,EAAsB,GAG/B4D,EAAsBzD,IACjC,MAAMC,KAAEA,EAAO,IAAGY,OAAEA,EAAS,EAACwC,GAAEA,EAAK,IAAQrD,EAC7C,OAAQC,EAAOY,GAAUwC,EAAK7D,GAAoB,EAAI6D,GAAM/D,KAAKI,IAAO,CAAC,EAG9DgE,EAAuBA,CAACzC,EAAO,MAAQoC,EAAK,KAChDA,EAAK,EAAIT,EAAoB3B,GAAQC,EAAkBD,GAGnD2C,EAAqBA,CAChCxC,EAAiC,GACjCpB,EAA4B,CAAA,KAE5B,IAAIC,KAAEA,EAAO,IAAGoD,GAAEA,EAAK,IAAQjC,GAC3BC,OACFA,EAAMC,OACNA,EAASoC,EAAqB,KAAOL,GAAGxC,OACxCA,EAAS8C,EAA2B,CAAE1D,OAAMoD,KAAIpC,KAAM,KACpDjB,EAECa,IACHA,EACE,GACEwC,EAAK/D,KAAKG,MAAMJ,EAAsBC,KAAKI,IAAOO,GAChD,EAAIoD,GAAMpD,EAAOX,KAAKI,GAAM,IAG/B2B,IACHA,EAAS/B,KAAKiC,IAAIjC,KAAKkC,KAAKvB,EAAOqB,GAAShC,KAAKmC,IAAI,EAAG,IAAM,GAC1DJ,EAAS,GAAM,GAAGA,KAGxB,MAAMK,GAAUL,EAAS,GAAK,EACxBM,EAAO,IAAIC,aAAaP,GAC9B,IAAK,IAAIQ,EAAI,EAAGA,GAAKH,EAAQG,IAC3BF,EAAKE,GAAK2B,EAAe3B,EAAIH,EAAQzB,EAAMoD,GAAMxC,EACjDc,EAAKN,EAAS,EAAIQ,GAAKF,EAAKE,GAG9B,OAAOF,CAAI,ECrIP,MAAOkC,EAYX9D,WAAAA,CAAmBC,EAA6C,IAC9D,MAAMC,KAAEA,EAAO,IAAG6D,MAAEA,EAAQ,IAAQ9D,EAEpCG,KAAKF,KAAOA,EACZE,KAAK2D,MAAQA,CACf,CAEOzD,WAAAA,CAAYJ,EAAOE,KAAKF,MAC7B,OAAO8D,EAAiC9D,EAC1C,CAEOM,WAAAA,CAAYC,GACjB,OAAOwD,EAAiCxD,EAC1C,CAEOC,GAAAA,CAAIC,GACT,OAAOuD,EAAyBvD,EAAGP,KAAKF,KAAME,KAAK2D,MACrD,CAEOlD,OAAAA,CAAQC,EAAS,GACtB,OAAOqD,EAA6B,CAClCjE,KAAME,KAAKF,KACXY,SACAiD,MAAO3D,KAAK2D,OAEhB,CAEO9C,SAAAA,CAAUC,GACf,OAAOkD,EAA+BlD,EACxC,CAEOE,OAAAA,CAAQnB,EAA4B,IACzC,OAAOoE,EAA6BjE,KAAMH,EAC5C,CAEO+B,eAAAA,CAAgBd,EAAO,GAC5B,MAAM6C,MAAEA,EAAK7D,KAAEA,GAASE,KACxB,OAAOkE,EAAqC,CAAEpE,OAAMgB,OAAM6C,SAC5D,CAEO9B,aAAAA,GACL,MAAO,CAAC,OAAQ,QAClB,EAGK,MAAMqC,EAAuCA,EAClDpE,OAAO,EACP6D,QAAQ,EACR7C,OAAO,KAECA,EAAOhB,GAAQ,QAAU,QAAW6D,GAAU,EAM3CI,EACXlE,IAEA,MAAMC,KAAEA,EAAO,IAAGY,OAAEA,EAAS,EAACiD,MAAEA,EAAQ,GAAM9D,EAC9C,OAAQa,EAASZ,GAAQ,QAAU,QAAW6D,GAAU,CAAC,EAG9CG,EAA2BA,CACtCvD,EACAT,EACA6D,KAEA,MAAMQ,GAAM,EAAI5D,EAAKT,IAAS,EAC9B,OAAQ,EAAI6D,IAAU,EAAIQ,GAAMR,GAAS,EAAIQ,EAAI,IAAO,EAAIA,EAAIA,GAAK,EAAE,EAG5DN,EAAoCxD,GACxCA,EAAQb,EAGJoE,EAAoC9D,GACxCA,EAAON,EAGHwE,EAAiCA,CAAClD,EAAO,SACpD,GAAIA,GAAQ,EACV,MAAM,IAAI8B,MAAM,0BAElB,MAAMC,EAA4B,IAAZ,EAAI/B,GACpBgC,EAAoBC,GAAc5D,KAAK6D,IAAI7D,KAAKI,IAAMwD,EAAI,KAChE,OACGD,EAAiB,EAAID,GAAgBC,EAAiBD,IAAiB,CAAC,EAQhEoB,EAA+BA,CAC1ChD,EAA2C,GAC3CpB,EAAwC,CAAA,KAExC,IAAIC,KAAEA,EAAO,IAAG6D,MAAEA,EAAQ,GAAM1C,GAC5BC,OACFA,EAAMC,OACNA,EAAS6C,IAAgCtD,OACzCA,EAASwD,EAAqC,CAAEpE,OAAMgB,KAAM,EAAG6C,WAC7D9D,EAECqB,IACHA,EAAS/B,KAAKiC,IAAIjC,KAAKkC,KAAKvB,EAAOqB,GAAShC,KAAKmC,IAAI,EAAG,IAAM,GAC1DJ,EAAS,GAAM,GAAGA,KAGxB,MAAMK,GAAUL,EAAS,GAAK,EACxBM,EAAO,IAAIC,aAAaP,GAC9B,IAAK,IAAIQ,EAAI,EAAGA,GAAKH,EAAQG,IAC3BF,EAAKE,GAAKoC,EAAyBpC,EAAIH,EAAQzB,EAAM6D,GAASjD,EAC9Dc,EAAKN,EAAS,EAAIQ,GAAKF,EAAKE,GAG9B,OAAOF,CAAI,ECzGP,MAAO4C,EAIXxE,WAAAA,CAAmBC,EAAkC,IACnD,IAAIC,KAAEA,EAAO,GAAEC,GAAEA,GAAOF,EAExBC,EAAOuE,EAAavE,EAAMC,GAE1BC,KAAKsE,MAAQxE,EAAKS,EAClBP,KAAKuE,MAAQzE,EAAK0E,CACpB,CAEOlE,GAAAA,CAAIC,EAAWiE,GACpB,OAAOC,EAAclE,EAAGiE,EAAGxE,KAAKsE,MAAOtE,KAAKuE,MAC9C,CAEOvD,OAAAA,CAAQnB,EAA4B,IACzC,OAAO6E,EACL,CACE5E,KAAM,CAAES,EAAGP,KAAKsE,MAAOE,EAAGxE,KAAKuE,QAEjC1E,EAEJ,CAEOgB,SAAAA,CAAU8D,EAAS,GACxB,OAAO5D,EAAkB4D,EAC3B,CAEOC,SAAAA,CACLlE,EAASmE,EAA0B,CACjC/E,KAAM,CAAES,EAAGP,KAAKsE,MAAOE,EAAGxE,KAAKuE,OAC/BI,OAAQ,KAGV,OAAOG,EAAoB,CACzBhF,KAAM,CAAES,EAAGP,KAAKsE,MAAOE,EAAGxE,KAAKuE,OAC/B7D,UAEJ,CAEON,WAAAA,CAAYC,GACjB,OAAOJ,EAAoBI,EAC7B,CAEOH,WAAAA,CAAYJ,GACjB,OAAOK,EAAoBL,EAC7B,CAEO8B,eAAAA,CAAgB+C,EAAS,GAC9B,OAAOE,EAA0B,CAC/BF,SACA7E,KAAM,CAAES,EAAGP,KAAKsE,MAAOE,EAAGxE,KAAKuE,QAEnC,CAEA,QAAWzE,CAAKA,GACdA,EAAOiF,EAAejF,GACtBE,KAAKsE,MAAQxE,EAAKS,EAClBP,KAAKuE,MAAQzE,EAAK0E,CACpB,EAGK,MAAMC,EAAgBA,CAC3BlE,EACAiE,EACAQ,EACAC,IAEO9F,KAAK2C,IACV5C,GAAuBC,KAAKmC,IAAIf,EAAIyE,EAAO,GAAK7F,KAAKmC,IAAIkD,EAAIS,EAAO,KAI3DP,EAAoBA,CAC/BzD,EACApB,EAA4B,MAE5B,IAAIC,KAAEA,EAAO,GAAEC,GAAEA,GAAOkB,EAExBnB,EAAOuE,EAAavE,EAAMC,GAE1B,IAAIoB,OACFA,EAASJ,IAAmBG,OAC5BA,EAAS,CAAEX,EAAG,EAAGiE,EAAG,GAAG9D,OACvBA,EAASmE,EAA0B,CAAE/E,OAAM6E,OAAQ,KACjD9E,EAEJsB,EAAS4D,EAAe5D,GAExBD,EAAS6D,EAAe7D,GAExB,IAAK,MAAMgE,IAAQ,CAAC,IAAK,KAClBhE,EAAOgE,KACVhE,EAAOgE,GAAQ/F,KAAKiC,IAClBjC,KAAKkC,KAAKvB,EAAKoF,GAAQ/D,EAAO+D,IAC9B/F,KAAKmC,IAAI,EAAG,IAAM,GAEhBJ,EAAOgE,GAAQ,GAAM,GAAGhE,EAAOgE,MAIvC,MAAMC,GAAWjE,EAAOX,EAAI,GAAK,EAC3B6E,GAAWlE,EAAOsD,EAAI,GAAK,EAC3BhD,EAAO,IAAI6D,MAAMnE,EAAOX,GAC9B,IAAK,IAAImB,EAAI,EAAGA,EAAIR,EAAOX,EAAGmB,IAC5BF,EAAKE,GAAK,IAAID,aAAaP,EAAOsD,GAEpC,IAAK,IAAI9C,EAAI,EAAGA,EAAIR,EAAOX,EAAGmB,IAC5B,IAAK,IAAI4D,EAAI,EAAGA,EAAIpE,EAAOsD,EAAGc,IAC5B9D,EAAKE,GAAG4D,GACNb,EAAc/C,EAAIyD,EAASG,EAAIF,EAAStF,EAAKS,EAAGT,EAAK0E,GAAK9D,EAGhE,OAAOc,CAAI,EAGAqD,EAA4BA,CACvChF,EAA4C,MAE5C,IAAI8E,OAAEA,EAAS,EAAC7E,KAAEA,EAAO,GAAEC,GAAEA,GAAOF,EAEpC,OADAC,EAAOuE,EAAavE,EAAMC,GAClB4E,EAASxF,KAAKC,IAAM,GAAMD,KAAKI,GAAKO,EAAK0E,EAAI1E,EAAKS,EAAE,EAGjDuE,EAAsBA,CACjCjF,EAAsC,MAEtC,IAAIC,KAAEA,EAAO,GAAEY,OAAEA,EAAS,EAACX,GAAEA,GAAOF,EAIpC,OAFAC,EAAOuE,EAAavE,EAAMC,GAElBW,EAASvB,KAAKI,GAAKO,EAAK0E,EAAI1E,EAAKS,EAAKpB,KAAKC,IAAM,CAAC,EAG5D,SAAS2F,EAAeQ,GACtB,MAAwB,iBAAVA,EAAqB,CAAEhF,EAAGgF,EAAOf,EAAGe,GAAU,IAAKA,EACnE,CAEA,SAASlB,EAAavE,EAA0BC,GAC9C,QAAWyF,IAAPzF,EAAkB,CACpB,IAAI0F,EAAWV,EAAehF,GAC9B,MAAO,CACLQ,EAAGN,EAAoB,EAAIwF,EAASlF,GACpCiE,EAAGvE,EAAoB,EAAIwF,EAASjB,GAEvC,CAAM,QAAagB,IAAT1F,EACT,OAAOiF,EAAejF,GAEtB,MAAM,IAAI8C,MAAM,mDAEpB,CC7MM,SAAU8C,EAAWzE,GACzB,MAAM0E,KAAEA,GAAS1E,EAEjB,OAAQ0E,GACN,IAAK,WACH,OAAO,IAAIhG,EAASsB,GACtB,IAAK,aACH,OAAO,IAAImB,EAAWnB,GACxB,IAAK,cACH,OAAO,IAAIgC,EAAYhC,GACzB,IAAK,wBACH,OAAO,IAAIyC,EAAsBzC,GACnC,QACE,MAAM2B,MAAM,wBAAwB+C,KAG1C,CCpBM,SAAUC,EAAW3E,GACzB,MAAM0E,KAAEA,GAAS1E,EACjB,GACO,aADC0E,EAEJ,OAAO,IAAIvB,EAAWnD,GAItB,MAAM2B,MAAM,wBAFa+C,IAK/B,CChBA,MAAME,EAAWC,OAAOC,UAAUF,SAoB5B,SAAUG,EAAWC,GACzB,MAAMC,EAAML,EAASM,KAAKF,GAC1B,OAAOC,EAAIE,SAAS,YAAcF,EAAIG,SAAS,MACjD,CCPM,SAAUC,EACdC,EACAC,EACA3G,EAAoC,CAAA,GAEpC,MAAM4G,OAAEA,GAAS,GAAS5G,EAC1B,GAAI4G,EAAQ,CACV,IAAIC,EAAM,EACNC,EAAOJ,EAAMrF,OAAS,EACtB0F,EAAS,EACb,KAAOD,EAAOD,EAAM,GAElB,GADAE,EAASF,GAAQC,EAAOD,GAAQ,GAC5BH,EAAMK,GAAUJ,EAClBE,EAAME,MACD,MAAIL,EAAMK,GAAUJ,GAGzB,OAAOI,EAFPD,EAAOC,CAGT,CAGF,OAAIF,EAAMH,EAAMrF,OAAS,EACnB/B,KAAK0H,IAAIL,EAASD,EAAMG,IAAQvH,KAAK0H,IAAIN,EAAMG,EAAM,GAAKF,GACrDE,EAEAA,EAAM,EAGRA,CAEX,CAAO,CACL,IAAII,EAAQ,EACRC,EAAOC,OAAOC,kBAClB,IAAK,IAAIvF,EAAI,EAAGA,EAAI6E,EAAMrF,OAAQQ,IAAK,CACrC,MAAMwF,EAAc/H,KAAK0H,IAAIN,EAAM7E,GAAK8E,GACpCU,EAAcH,IAChBA,EAAOG,EACPJ,EAAQpF,EAEZ,CACA,OAAOoF,CACT,CACF,CClBM,SAAUK,EAGdC,EACAlG,GAEA,OAAIkG,IAAqB/B,MAChB,IAAI+B,EAAiBlG,GAAQmG,KAClC,GAGK,IAAID,EACTlG,EAGN,CCxDA,MACMoG,EAAY,EAAI,SAMtB,SAASC,EAAeC,EAAWC,GAGjC,MAAMC,EAAU,OAFhBF,KAAO,GAIP,QADYA,EAAIE,IAFhBD,KAAO,KAGgB,GAAKC,EAAMD,IAAO,CAC3C,CAEM,MAAOE,EACJC,OAECC,MAMRjI,WAAAA,CAAYkI,EAAOC,KAAKC,OACtBhI,KAAK6H,MAAQ,IAAII,YAAY,GAC7BjI,KAAKkI,KAAKJ,GACV9H,KAAK4H,OAAS5H,KAAKmI,SAASC,KAAKpI,KACnC,CAKOqI,SAAAA,GAEL,OADArI,KAAKsI,YACGtI,KAAK6H,MAAM,GAAK7H,KAAK6H,MAAM,KAAQ,CAC7C,CAKOM,QAAAA,GACL,OAAQnI,KAAKqI,cAAgB,GAAKf,CACpC,CAEOY,IAAAA,CAAKJ,GACV,IAAKd,OAAOuB,UAAUT,GACpB,MAAM,IAAIU,UAAU,2BAEtBxI,KAAK6H,MAAM,GAAKC,EAChB9H,KAAK6H,MAAM,GAAK,EAChB7H,KAAK6H,MAAM,GAAK,EAChB7H,KAAK6H,MAAM,GAAK,EAChB,IAAK,IAAInG,EAAI,EAAGA,EArDP,EAqDiBA,IACxB1B,KAAK6H,MAAU,EAAJnG,IACRA,EACC6F,EACE,WACAvH,KAAK6H,MAAOnG,EAAI,EAAK,GAAO1B,KAAK6H,MAAOnG,EAAI,EAAK,KAAO,KAAQ,KAEpE,EAEJ1B,KAAKyI,sBACL,IAAK,IAAI/G,EAAI,EAAGA,EA/DP,EA+DiBA,IACxB1B,KAAKsI,WAET,CAEQG,mBAAAA,GAEc,IAAlBzI,KAAK6H,MAAM,IACO,IAAlB7H,KAAK6H,MAAM,IACO,IAAlB7H,KAAK6H,MAAM,IACO,IAAlB7H,KAAK6H,MAAM,KAEX7H,KAAK6H,MAAM,GAAK,GAChB7H,KAAK6H,MAAM,GAAK,GAChB7H,KAAK6H,MAAM,GAAK,GAChB7H,KAAK6H,MAAM,GAAK,GAEpB,CAEQS,SAAAA,GACN,IAAII,EAAI1I,KAAK6H,MAAM,GACnBa,GAAKA,GAjFG,GAkFRA,GAAKA,IAjFG,GAkFRA,GAAK1I,KAAK6H,MAAM,IAjFR,GAkFR7H,KAAK6H,MAAM,GAAK7H,KAAK6H,MAAM,GAC3B7H,KAAK6H,MAAM,GAAK7H,KAAK6H,MAAM,GAC3B7H,KAAK6H,MAAM,GAAK7H,KAAK6H,MAAM,GAC3B7H,KAAK6H,MAAM,GAAKa,CAClB,ECJF,MAAMC,EACJC,GAAS,EACTC,IAAY,EAEZC,GACAC,GACAC,GAEApJ,WAAAA,CAAYkJ,EAAcC,EAA2BC,GACnDhJ,MAAK8I,EAAQA,EACb9I,MAAK+I,EAAqBA,EAC1B/I,MAAKgJ,EAAaA,CACpB,CAEAC,gBAAAA,GACE,IAAIC,EAAK/E,EAAGgF,EAAGC,EAEf,GAAIpJ,MAAK6I,EACP7I,MAAK6I,GAAY,EACjBK,EAAMlJ,MAAK4I,EAAS5I,MAAK+I,EAAqB/I,MAAK8I,MAC9C,CACL,GACE3E,EAA+B,EAA3BnE,MAAKgJ,EAAWpB,SAAe,EACnCuB,EAA+B,EAA3BnJ,MAAKgJ,EAAWpB,SAAe,EAEnCwB,EAAIjF,EAAIA,EAAIgF,EAAIA,QACTC,GAAK,GAAW,IAANA,GAEnBA,EAAIjK,KAAKG,MAAQ,EAAGH,KAAK6C,IAAIoH,GAAMA,GAEnCpJ,MAAK4I,EAASO,EAAIC,EAClBpJ,MAAK6I,GAAY,EACjBK,EAAMlJ,MAAK8I,EAAQ9I,MAAK+I,EAAqB5E,EAAIiF,CACnD,CACA,OAAOF,CACT,EC/GI,SAAUG,EACd9C,EACA1G,EAAkC,KCI9B,SACJ0F,EACA1F,EAAyB,IAEzB,MAAMyJ,UAAEA,GAAczJ,EACtB,IAAKmG,EAAWT,GACd,MAAM,IAAIiD,UAAU,0BAEtB,GAAqB,IAAjBjD,EAAMrE,OACR,MAAM,IAAIsH,UAAU,2BAEtB,GAAwB,iBAAbjD,EAAM,GACf,MAAM,IAAIiD,UAAU,8BAEtB,GAAIc,GAAa/D,EAAMrE,OAASoI,EAC9B,MAAM,IAAI1G,MAAM,wCAAwC0G,IAE5D,CDnBEC,CAAOhD,GACP,MAAMiD,UAAEA,EAASC,QAAEA,GEiBf,SACJlJ,EACAV,EAAkC,IAElC,IAAI2J,UAAEA,EAASC,QAAEA,GAAY5J,EAC7B,MAAM6J,KAAEA,EAAIC,GAAEA,GAAO9J,EAsBrB,YApBkB2F,IAAdgE,IAEAA,OADWhE,IAATkE,EACUpD,EAAkB/F,EAAGmJ,GAErB,QAGAlE,IAAZiE,IAEAA,OADSjE,IAAPmE,EACQrD,EAAkB/F,EAAGoJ,GAErBpJ,EAAEW,OAAS,GAGrBsI,EAAY,IAAGA,EAAY,GAC3BC,EAAU,IAAGA,EAAU,GACvBD,GAAajJ,EAAEW,SAAQsI,EAAYjJ,EAAEW,OAAS,GAC9CuI,GAAWlJ,EAAEW,SAAQuI,EAAUlJ,EAAEW,OAAS,GAE1CsI,EAAYC,KAAUD,EAAWC,GAAW,CAACA,EAASD,IACnD,CAAEA,YAAWC,UACtB,CF7CiCG,CAAgBrD,EAAO1G,GACtD,IAAIgK,EAAWtD,EAAMiD,GAErB,IAAK,IAAI9H,EAAI8H,EAAY,EAAG9H,GAAK+H,EAAS/H,IACpC6E,EAAM7E,GAAKmI,IACbA,EAAWtD,EAAM7E,IAGrB,OAAOmI,CACT,CGUM,SAAUC,GAGdjK,GAEA,MAAMkK,OACJA,EAAMC,OACNA,EAASD,GAAQ7I,QAAU,EAAC+I,UAC5BA,EAAYF,IAAS,GAAG7I,QAAU,EAACkG,iBACnCA,EAAmB3F,cACjB5B,EAEEqK,EAA0D,GAChE,IAAK,IAAIC,EAAM,EAAGA,EAAMH,EAAQG,IAC9BD,EAAUE,KAAKjD,EAAkBC,EAAkB6C,IAErD,OAAOC,CACT,CC5CM,SAAUG,GAAcN,ICNxB,SAAsBvI,GAC1B,GAAoB,IAAhBA,EAAKN,QAAmC,IAAnBM,EAAK,GAAGN,OAC/B,MAAM,IAAIoJ,WAAW,4BAGvB,MAAMC,EAAc/I,EAAK,GAAGN,OAC5B,IAAK,IAAIQ,EAAI,EAAGA,EAAIF,EAAKN,OAAQQ,IAC/B,GAAIF,EAAKE,GAAGR,SAAWqJ,EACrB,MAAM,IAAID,WAAW,oCAG3B,CDDEE,CAAYT,GACZ,MAAMC,EAASD,EAAO7I,OAChB+I,EAAYF,EAAO,GAAG7I,OAE5B,IAAIE,EAAM2I,EAAO,GAAG,GAChBU,EAAMV,EAAO,GAAG,GAEpB,IAAK,IAAIW,EAAS,EAAGA,EAAST,EAAWS,IACvC,IAAK,IAAIP,EAAM,EAAGA,EAAMH,EAAQG,IAC1BJ,EAAOI,GAAKO,GAAUtJ,IAAKA,EAAM2I,EAAOI,GAAKO,IAC7CX,EAAOI,GAAKO,GAAUD,IAAKA,EAAMV,EAAOI,GAAKO,IAIrD,MAAO,CAAEtJ,MAAKqJ,MAChB,CEGc,SAAUE,GAASnJ,EAAc3B,EAAwB,IACrE,MAAMiI,KAAEA,EAAO,EAAC8C,aAAEA,EAAe,SAAQC,QAAEA,EAAU,GAAMhL,EACrDiL,EAASzB,EAAU7H,EAAKgD,GAAKqG,EAAW,IACxCE,EPaF,SACJlL,EAAoC,IAEpC,MAAMiJ,KACJA,EAAO,EAACC,kBACRA,EAAoB,EAAC7H,OACrBA,EAAS,IAAI4J,MACbA,EAAQ,EAAChD,KACTA,EAAI8C,aACJA,EAAe,UACb/K,EAEEmJ,EAAY,IAAIrB,EAAMG,GACtBkD,EAAc,IAAIvJ,aAAaP,GAErC,OAAQ0J,GACN,IAAK,SAAU,CACb,MAAMK,EAAoB,IAAItC,EAC5BG,EACAC,EACAC,GAEF,IAAK,IAAItH,EAAI,EAAGA,EAAIR,EAAQQ,IAC1BsJ,EAAYtJ,GAAKuJ,EAAkBhC,mBAErC,KACF,CACA,IAAK,UACH,IAAK,IAAIvH,EAAI,EAAGA,EAAIR,EAAQQ,IAC1BsJ,EAAYtJ,IAAMsH,EAAUpB,SAAW,IAAOkD,EAAQhC,EAExD,MAEF,QACE,MAAM,IAAIlG,MAAM,yBAAyBsI,OAAON,MAIpD,OAAOI,CACT,COpDgBG,CAAkB,CAC9BP,eACA9C,OACAgB,KAAM,EACNC,kBAAmB+B,EACnBA,QACA5J,OAAQM,EAAKjB,EAAEW,SAGjB,OADAM,EAAKgD,EChCD,SACJ4G,EACAC,GAEA,IAAIC,GAAa,EACbC,EAAW,EACf,GAAIvF,EAAWqF,IACb,GAAID,EAAOlK,SAAWmK,EAAOnK,OAC3B,MAAM,IAAI0B,MAAM,oDAGlB0I,GAAa,EACbC,EAAWF,EAGb,MAAMG,EAAS,IAAI/J,aAAa2J,EAAOlK,QACvC,GAAIoK,EACF,IAAK,IAAI5J,EAAI,EAAGA,EAAI0J,EAAOlK,OAAQQ,IACjC8J,EAAO9J,GAAK0J,EAAO1J,GAAK6J,OAG1B,IAAK,IAAI7J,EAAI,EAAGA,EAAI0J,EAAOlK,OAAQQ,IACjC8J,EAAO9J,GAAK0J,EAAO1J,GAAM2J,EAAuB3J,GAIpD,OAAO8J,CACT,CDKWC,CAAKjK,EAAKgD,EAAGuG,GACfvJ,CACT,CE4DM,MAAOkK,GASX9L,WAAAA,CAAmBC,EAAuB,IACxC,MAAM6J,KACJA,EAAO,EAACC,GACRA,EAAK,IAAIgC,SACTA,EAAW,MAAKC,aAChBA,EAAY3K,MACZA,EAAQ,CAAE0E,KAAM,WAAY7F,KAAM,IAChCD,EAEJG,KAAK0J,KAAOA,EACZ1J,KAAK2J,GAAKA,EACV3J,KAAK2L,SAAWA,EAChB3L,KAAK6L,UAAY7L,KAAK2J,GAAK3J,KAAK0J,OAAS1J,KAAK2L,SAAW,GACzD3L,KAAK4L,aAAeA,EACpB5L,KAAK8L,cAAgB9E,OAAO+E,iBAE5B/L,KAAKwB,KAAO,CACVjB,EAAG,IAAIkB,aAAazB,KAAK2L,UACzBnH,EAAG,IAAI/C,aAAazB,KAAK2L,WAG3B,MAAMK,EAAiBtG,EAAWzE,GAOlC,GANAjB,KAAKiB,MAAQ+K,EAEbC,GAAajM,KAAK0J,KAAM,QACxBuC,GAAajM,KAAK2J,GAAI,MA4O1B,SAAuB1D,EAAeiG,GACpC,IAAKlF,OAAOuB,UAAUtC,GACpB,MAAM,IAAIuC,UAAU,GAAG0D,8BAE3B,CA/OIC,CAAcnM,KAAK2L,SAAU,YAEzB3L,KAAK2J,IAAM3J,KAAK0J,KAClB,MAAM,IAAIY,WAAW,sCAGvB,GAAItK,KAAK4L,cAA6C,mBAAtB5L,KAAK4L,aACnC,MAAM,IAAIpD,UAAU,0CAGtBxI,KAAKoM,OACP,CAQOC,QAAAA,CAASC,EAA8BzM,GAC5C,KACGwF,MAAMkH,QAAQD,IACG,iBAAVA,QACM9G,IAAZ8G,EAAM/L,QACMiF,IAAZ8G,EAAM9H,GACLa,MAAMkH,QAAQD,EAAM/L,IACpB8E,MAAMkH,QAAQD,EAAM9H,IACrB8H,EAAM/L,EAAEW,SAAWoL,EAAM9H,EAAEtD,QAE7B,MAAM,IAAIsH,UACR,8DAGJ,GAAInD,MAAMkH,QAAQD,GAChB,IAAK,MAAME,KAAQF,EACjBtM,KAAKyM,QAAQD,EAAM3M,QAGrB,IAAK,IAAI6B,EAAI,EAAGA,EAAI4K,EAAM/L,EAAEW,OAAQQ,IAClC1B,KAAKyM,QAAQ,CAACH,EAAM/L,EAAEmB,GAAI4K,EAAM9H,EAAE9C,IAAK7B,GAG3C,OAAOG,IACT,CASOyM,OAAAA,CAAQD,EAAc3M,EAAuB,IAClD,GAAIwF,MAAMkH,QAAQC,IAASA,EAAKtL,OAAS,EACvC,MAAM,IAAI0B,MACR,mFAIJ,IACGyC,MAAMkH,QAAQC,UACHhH,IAAXgH,EAAKjM,QAA8BiF,IAAXgH,EAAKhI,GAE9B,MAAM,IAAI5B,MACR,mFAIJ,IAAI8J,EACAC,EACAC,EACAC,EACAC,EACAzH,MAAMkH,QAAQC,IACfE,EAAWC,EAAWC,EAAUE,GAAoBN,GAErDE,EAAYF,EAAKjM,EACjBoM,EAAYH,EAAKhI,EACjBqI,EAAYL,EAAKnM,MACjByM,EAAmBN,EAAKvL,OAEtB0L,EAAY3M,KAAK8L,gBAAe9L,KAAK8L,cAAgBa,GAEzD,IAAM1L,MAAO8L,GAAiBlN,EAE1BiN,IACFC,EAAeA,EACX,IAAKA,KAAiBD,GACtBA,GAGN,MAAM7L,EAAQ8L,EACVrH,EAAWqH,GACVjH,OAAOkH,OACNlH,OAAOmH,OAAOnH,OAAOoH,eAAelN,KAAKiB,QACzCjB,KAAKiB,QAGLZ,MAAEA,GAAUR,EAClB,IAAIsN,UAAEA,EAASC,WAAEA,GAAevN,EAWhC,MAAMC,OACS0F,IAAboH,EACIA,EACAC,EACE5L,EAAMb,YAAYyM,GAClB7M,KAAK4L,aACH5L,KAAK4L,aAAac,QACRlH,IAAVnF,EACEA,EACAY,EAAMnB,KAKlB,GAHKqN,IAAWA,EAAYrN,GACvBsN,IAAYA,EAAatN,IAEzBqN,IAAcC,EACjB,MAAM,IAAIxK,MAAM,4CAGlB,MAAMzB,OACeqE,IAAnB3F,EAAQsB,OAAuBF,EAAMJ,YAAchB,EAAQsB,OAEvDkM,EAAaX,EAAaS,EAAY,EAAKhM,EAC3CmM,EAAYZ,EAAaU,EAAa,EAAKjM,EAE3CoM,EAAapO,KAAKsL,IACtB,EACAtL,KAAKqO,OAAOH,EAAarN,KAAK0J,MAAQ1J,KAAK6L,WAEvC4B,EAAYtO,KAAKiC,IACrBpB,KAAK2L,SAAW,EAChBxM,KAAKkC,MAAMiM,EAAYtN,KAAK0J,MAAQ1J,KAAK6L,WAErC6B,EAAcvO,KAAKwO,OAAOjB,EAAY1M,KAAK0J,MAAQ1J,KAAK6L,UAK9D5K,EAAMnB,KAAOqN,EACb,IAAK,IAAIrG,EAAQyG,EAAYzG,EAAQ3H,KAAKsL,IAAIiD,EAAa,GAAI5G,IAC7D9G,KAAKwB,KAAKgD,EAAEsC,IACV6F,EAAY1L,EAAMX,IAAIN,KAAKwB,KAAKjB,EAAEuG,GAAS4F,GAI/CzL,EAAMnB,KAAOsN,EACb,IACE,IAAItG,EAAQ3H,KAAKiC,IAAIsM,EAAaD,GAClC3G,GAAS2G,EACT3G,IAEA9G,KAAKwB,KAAKgD,EAAEsC,IACV6F,EAAY1L,EAAMX,IAAIN,KAAKwB,KAAKjB,EAAEuG,GAAS4F,GAG/C,OAAO1M,IACT,CAOO4N,WAAAA,CAAYC,GAEjB,OC9SU,SAAsBrM,EAAcqM,GAChD,IAAKA,EAAa,OAAOrM,EACzB,MAAMsM,EAAKtM,EAAKjB,EACVwN,EAAKvM,EAAKgD,EAChB,IAAK,IAAI9C,EAAI,EAAGA,EAAIoM,EAAG5M,OAAQQ,IAC7BqM,EAAGrM,IAAMmM,EAAYC,EAAGpM,GAG5B,CDqSIkM,CAAY5N,KAAKwB,KAAMqM,GAChB7N,IACT,CAOO2K,QAAAA,CAAS9K,GAEd,OADA8K,GAAS3K,KAAKwB,KAAM3B,GACbG,IACT,CAOOgO,WAAAA,CAAYnO,EAAwC,IAClC,kBAAZA,IACTA,EAAU,CAAEoO,KAAMpO,IAEpB,MAAMoO,KAAEA,GAAO,EAAIC,UAAEA,EAAY,GAAMrO,EACvC,GAAIqO,EAAW,CACb,MAAMC,EAAgBnO,KAAK8L,cAAgBoC,EACrC3N,EAAI,GACJiE,EAAI,GACV,IAAK,IAAI9C,EAAI,EAAGA,EAAI1B,KAAKwB,KAAKjB,EAAEW,OAAQQ,IAClC1B,KAAKwB,KAAKgD,EAAE9C,IAAMyM,IACpB5N,EAAE6J,KAAKpK,KAAKwB,KAAKjB,EAAEmB,IACnB8C,EAAE4F,KAAKpK,KAAKwB,KAAKgD,EAAE9C,KAGvB,MAAO,CAAEnB,EAAGkB,aAAaiI,KAAKnJ,GAAIiE,EAAG/C,aAAaiI,KAAKlF,GACzD,CACA,OAAIyJ,EACK,CACL1N,EAAGP,KAAKwB,KAAKjB,EAAE6N,QACf5J,EAAGxE,KAAKwB,KAAKgD,EAAE4J,SAGVpO,KAAKwB,IAEhB,CAMO4K,KAAAA,GACL,GAAsB,IAAlBpM,KAAK2L,SACP3L,KAAKwB,KAAKjB,EAAE,IAAMP,KAAK0J,KAAO1J,KAAK2J,IAAM,OAEzC,IAAK,IAAIjI,EAAI,EAAGA,EAAI1B,KAAK2L,SAAUjK,IACjC1B,KAAKwB,KAAKjB,EAAEmB,GAAK1B,KAAK0J,KAAOhI,EAAI1B,KAAK6L,SAG1C,OAAO7L,IACT,EASF,SAASiM,GAAahG,EAAeiG,GACnC,IAAKlF,OAAOqH,SAASpI,GACnB,MAAM,IAAIuC,UAAU,GAAG0D,4BAE3B,CE/WA,MAAMoC,GAAmB,CAAC,IAAK,KAGzBC,GAAqC,CAAC,IAAK,IAAK,KAQhDC,GAAqBA,CAACvN,EAAqBZ,KAC/C,MAAMoO,EAAY1J,GAAe1E,GACjC,IAAK,MAAMqO,KAAOJ,GAChBG,EAAUC,GAAOzN,EAAMb,YAAYqO,EAAUC,IAE/C,OAAOD,CAAS,EAgFZ,MAAOE,GAUX/O,WAAAA,CAAmBC,EAAuB,IACxC,MAAM+L,aACJA,EAAeA,IAAM,EAAC3K,MACtBA,EAAQ,CACN0E,KAAM,aAEN9F,EACJ,IAAI6J,KAAEA,EAAO,EAACC,GAAEA,EAAK,IAAGgC,SAAEA,EAAW,MAAS9L,EAE9C6J,EAAO3E,GAAe2E,GACtBC,EAAK5E,GAAe4E,GACpBgC,EAAW5G,GAAe4G,GAE1B,IAAK,MAAMzG,KAAQoJ,GACjBrC,GAAavC,EAAKxE,GAAO,QAAQA,KACjC+G,GAAatC,EAAGzE,GAAO,MAAMA,KAC7BiH,GAAcR,EAASzG,GAAO,YAAYA,KAG5ClF,KAAK0J,KAAOA,EACZ1J,KAAK2J,GAAKA,EACV3J,KAAK2L,SAAWA,EAChB3L,KAAK6L,SAmQT,SAA0BnC,EAAgBC,EAAcgC,GACtD,MAAO,CACLpL,GAAIoJ,EAAGpJ,EAAImJ,EAAKnJ,IAAMoL,EAASpL,EAAI,GACnCiE,GAAImF,EAAGnF,EAAIkF,EAAKlF,IAAMmH,EAASnH,EAAI,GAEvC,CAxQoBoK,CAAiBlF,EAAMC,EAAIgC,GAE3C3L,KAAK4L,aAAeA,EACpB5L,KAAK8L,cAAgB9E,OAAO+E,iBAE5B,MAAMC,EAAiBpG,EAAW3E,GAClCjB,KAAKiB,MAAQ+K,EAEbhM,KAAKwB,KAAO,CACVjB,EAAG,IAAIkB,aAAakK,EAASpL,GAC7BiE,EAAG,IAAI/C,aAAakK,EAASnH,GAC7BqK,EAAG/E,GAAkB,CACnBE,OAAQhK,KAAK2L,SAASnH,EACtByF,UAAWjK,KAAK2L,SAASpL,KAI7B,IAAK,MAAM2E,KAAQoJ,GACjB,GAAItO,KAAK2J,GAAGzE,IAASlF,KAAK0J,KAAKxE,GAC7B,MAAM,IAAIoF,WAAW,sCAIzB,GAAiC,mBAAtBtK,KAAK4L,aACd,MAAM,IAAIpD,UAAU,0CAGtBxI,KAAKoM,OACP,CAQOC,QAAAA,CAASC,EAAgCzM,GAC9C,KACGwF,MAAMkH,QAAQD,IACG,iBAAVA,QACM9G,IAAZ8G,EAAM/L,QACMiF,IAAZ8G,EAAM9H,GACLa,MAAMkH,QAAQD,EAAM/L,IACpB8E,MAAMkH,QAAQD,EAAM9H,IACrB8H,EAAM/L,EAAEW,SAAWoL,EAAM9H,EAAEtD,QAE7B,MAAM,IAAIsH,UACR,8DAGJ,GAAInD,MAAMkH,QAAQD,GAChB,IAAK,MAAME,KAAQF,EACjBtM,KAAKyM,QAAQD,EAAM3M,OAEhB,CACL,MAAMiP,EAAUxC,EAAM/L,EAAEW,OACxB,IAAK,MAAM6N,KAAKR,GACd,GACEjC,EAAMyC,IACN1J,MAAMkH,QAAQD,EAAMyC,KACpBD,IAAYxC,EAAMyC,GAAG7N,OAErB,MAAM,IAAI0B,MAAM,uCAGpB,IAAK,IAAIlB,EAAI,EAAGA,EAAI4K,EAAM/L,EAAEW,OAAQQ,IAClC1B,KAAKyM,QAAQ,CAACH,EAAM/L,EAAEmB,GAAI4K,EAAM9H,EAAE9C,GAAI4K,EAAMuC,EAAEnN,IAAK7B,EAEvD,CAEA,OAAOG,IACT,CAQOyM,OAAAA,CAAQD,EAAc3M,EAA4B,IACvD,GAAIwF,MAAMkH,QAAQC,IAASA,EAAKtL,OAAS,EACvC,MAAM,IAAI0B,MACR,sFAIJ,IACGyC,MAAMkH,QAAQC,IACf+B,GAAgBS,MAAMC,QAAkBzJ,IAAZgH,EAAKyC,KAEjC,MAAM,IAAIrM,MACR,sFAIJ,IAAI8J,EACAwC,EACAvC,EACAC,EACAC,EACAC,EACAzH,MAAMkH,QAAQC,IACfE,EAAWwC,EAAWvC,EAAWC,EAAUE,GAAoBN,GAEhEE,EAAYF,EAAKjM,EACjB2O,EAAY1C,EAAKhI,EACjBmI,EAAYH,EAAKqC,EACjBjC,EAAWJ,EAAK1M,KAChB+M,EAAYL,EAAKnM,MACjByM,EAAmBN,EAAKvL,OAG1B,MAAMkO,EAAqB,CAAE5O,EAAGmM,EAAWlI,EAAG0K,GAC1CvC,EAAY3M,KAAK8L,gBAAe9L,KAAK8L,cAAgBa,GAEzD,MAAMtM,MAAEA,GAAUR,EAClB,IAAMoB,MAAO8L,GAAiBlN,EAE1BiN,IACFC,EAAeA,EACX,IAAKA,KAAiBD,GACtBA,GAGN,MAAM7L,EAAQ8L,EACVnH,EAAWmH,GACVjH,OAAOkH,OACNlH,OAAOmH,OAAOnH,OAAOoH,eAAelN,KAAKiB,QACzCmO,gBAAgBpP,KAAKiB,QAG3B,IAAInB,KACFA,QAAoB0F,IAAboH,EACHA,EACAC,EACE2B,GAAmBvN,EAAO4L,GAC1BxM,EACEmO,GAAmBvN,EAAOZ,GAC1BL,KAAK4L,aAAac,EAAWwC,KACnCrP,EAEJC,EAAOiF,GAAejF,GAEtB,IAAIqB,OACiBqE,IAAnB3F,EAAQsB,OAAuBF,EAAMJ,YAAchB,EAAQsB,OAE7DA,EAAS4D,GAAe5D,GAExB,MAAMoM,EAAuB,CAAEhN,EAAG,EAAGiE,EAAG,GAClCiJ,EAAsB,CAAElN,EAAG,EAAGiE,EAAG,GACvC,IAAK,MAAMU,KAAQoJ,GAAQ,CACzB,MAAMe,EAAQF,EAASjK,GAASpF,EAAKoF,GAAQ,EAAK/D,EAAO+D,GACnDoK,EAAOH,EAASjK,GAASpF,EAAKoF,GAAQ,EAAK/D,EAAO+D,GACxDqI,EAAWrI,GAAQ/F,KAAKsL,IACtB,EACAtL,KAAKqO,OAAO6B,EAAQrP,KAAK0J,KAAKxE,IAASlF,KAAK6L,SAAS3G,KAEvDuI,EAAUvI,GAAQ/F,KAAKiC,IACrBpB,KAAK2L,SAASzG,GACd/F,KAAKkC,MAAMiO,EAAOtP,KAAK0J,KAAKxE,IAASlF,KAAK6L,SAAS3G,IAEvD,CAEAjE,EAAMnB,KAAOA,EACb,IAAK,IAAIyP,EAAShC,EAAWhN,EAAGgP,EAAS9B,EAAUlN,EAAGgP,IACpD,IAAK,IAAIC,EAASjC,EAAW/I,EAAGgL,EAAS/B,EAAUjJ,EAAGgL,IAAU,CAC9D,MAAMvJ,EACJ0G,EACA1L,EAAMX,IACJN,KAAKwB,KAAKjB,EAAEgP,GAAUJ,EAAS5O,EAC/BP,KAAKwB,KAAKgD,EAAEgL,GAAUL,EAAS3K,GAE/BrF,KAAK0H,IAAIZ,GAAS,OACpBjG,KAAKwB,KAAKqN,EAAEW,GAAQD,IAAWtJ,EAEnC,CAGF,OAAOjG,IACT,CAOOgO,WAAAA,CAAYnO,EAA0C,IACpC,kBAAZA,IACTA,EAAU,CAAEoO,KAAMpO,IAEpB,MAAMoO,KAAEA,GAAO,GAASpO,EAClB4P,EAAUpF,GAAcrK,KAAKwB,KAAKqN,GAExC,MAAO,CACLa,KAAM1P,KAAK0J,KAAKnJ,EAChBoP,KAAM3P,KAAK2J,GAAGpJ,EACdqP,KAAM5P,KAAK2J,GAAGnF,EACdqL,KAAM7P,KAAK0J,KAAKlF,EAChBsL,KAAML,EAAQrO,IACd2O,KAAMN,EAAQhF,IACdoE,EAAGZ,EAAOjO,KAAKwB,KAAKqN,EAAET,QAAUpO,KAAKwB,KAAKqN,EAE9C,CAMOzC,KAAAA,GACL,MAAM4D,EAAmBhQ,KAAKwB,KAE9B,IAAK,MAAM0D,KAAQoJ,GACjB,IAAK,IAAI5M,EAAI,EAAGA,EAAI1B,KAAK2L,SAASzG,GAAOxD,IACvCsO,EAAS9K,GAAMxD,GAAK1B,KAAK0J,KAAKxE,GAAQxD,EAAI1B,KAAK6L,SAAS3G,GAG5D,IAAK,MAAMiF,KAAO6F,EAASnB,EACzB,IAAK,IAAIvJ,EAAI,EAAGA,EAAI6E,EAAIjJ,OAAQoE,IAC9B6E,EAAI7E,GAAK,EAGb,OAAOtF,IACT,EA0BF,SAAS+E,GAAeQ,GACtB,MAAwB,iBAAVA,EAAqB,CAAEhF,EAAGgF,EAAOf,EAAGe,GAAU,IAAKA,EACnE,CAsBA,SAAS4G,GAAclG,EAAeiG,GACpC,IAAKlF,OAAOuB,UAAUtC,GACpB,MAAM,IAAIuC,UAAU,GAAG0D,8BAE3B,CAQA,SAASD,GAAahG,EAAeiG,GACnC,IAAKlF,OAAOqH,SAASpI,GACnB,MAAM,IAAIuC,UAAU,GAAG0D,4BAE3B,oEFzCM,SACJI,EACAzM,EAAmC,IAEnC,MACEmJ,UAAWiH,EAAgBlF,MAC3BA,EAAKmF,SACLA,EAAQhC,UACRA,EAASiC,YACTA,GACEtQ,EAEEmJ,EAAY,IAAI0C,GAAkBuE,GAOxC,OALAjH,EAAUqD,SAASC,EAAO6D,GACtBD,GAAUlH,EAAU4E,YAAYsC,GAChCnF,GACF/B,EAAU2B,SAASI,GAEd/B,EAAUgF,YAAY,CAC3BE,aAEJ,uBEtCM,SACJ5B,EACAzM,EAAqC,IAErC,MAAQmJ,UAAWiH,EAAkB3D,MAAO8D,GAAoBvQ,EAE1DmJ,EAAY,IAAI2F,GAAoBsB,GAG1C,OADAjH,EAAUqD,SAASC,EAAO8D,GACnBpH,EAAUgF,aACnB","x_google_ignoreList":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,21]}