@carbonenginejs/runtime-resource 0.1.0

package/dist/formats/ogg/core/imdct.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"imdct.js","sources":["../../../../../src/formats/ogg/core/imdct.js"],"sourcesContent":["/**\n * Fast IMDCT for Vorbis synthesis.\n *\n * The IMDCT is computed as a DCT-IV of half the block size followed by the\n * standard three-segment rearrangement; the DCT-IV itself runs on a complex\n * FFT of a quarter block with pre/post twiddles. The factorization was\n * derived analytically and is verified against the naive transforms to float\n * precision in the test suite.\n */\n\nconst tableCache = new Map();\n\nfunction getTables(m)\n{\n    let tables = tableCache.get(m);\n    if (tables) return tables;\n\n    const h = m >> 1;\n    const preCos = new Float32Array(h);\n    const preSin = new Float32Array(h);\n    const postCos = new Float32Array(h);\n    const postSin = new Float32Array(h);\n    for (let i = 0; i < h; i++)\n    {\n        const pre = (Math.PI * i) / m;\n        preCos[i] = Math.cos(pre);\n        preSin[i] = Math.sin(pre);\n        const post = (Math.PI * (4 * i + 1)) / (4 * m);\n        postCos[i] = Math.cos(post);\n        postSin[i] = Math.sin(post);\n    }\n\n    const bitrev = new Uint32Array(h);\n    for (let i = 1, j = 0; i < h; i++)\n    {\n        let bit = h >> 1;\n        for (; j & bit; bit >>= 1) j ^= bit;\n        j ^= bit;\n        bitrev[i] = j;\n    }\n\n    const twiddles = [];\n    for (let len = 2; len <= h; len <<= 1)\n    {\n        const half = len >> 1;\n        const cos = new Float32Array(half);\n        const sin = new Float32Array(half);\n        for (let j = 0; j < half; j++)\n        {\n            const angle = (-2 * Math.PI * j) / len;\n            cos[j] = Math.cos(angle);\n            sin[j] = Math.sin(angle);\n        }\n        twiddles.push({ len, cos, sin });\n    }\n\n    tables = {\n        preCos,\n        preSin,\n        postCos,\n        postSin,\n        bitrev,\n        twiddles,\n        workRe: new Float32Array(h),\n        workIm: new Float32Array(h),\n        dct: new Float32Array(m)\n    };\n    tableCache.set(m, tables);\n    return tables;\n}\n\nfunction fftInPlace(re, im, tables)\n{\n    const n = re.length;\n    const bitrev = tables.bitrev;\n    for (let i = 1; i < n; i++)\n    {\n        const j = bitrev[i];\n        if (i < j)\n        {\n            const tr = re[i]; re[i] = re[j]; re[j] = tr;\n            const ti = im[i]; im[i] = im[j]; im[j] = ti;\n        }\n    }\n    for (const { len, cos, sin } of tables.twiddles)\n    {\n        const half = len >> 1;\n        for (let i = 0; i < n; i += len)\n        {\n            for (let j = 0; j < half; j++)\n            {\n                const wr = cos[j];\n                const wi = sin[j];\n                const a = i + j;\n                const b = a + half;\n                const vr = re[b] * wr - im[b] * wi;\n                const vi = re[b] * wi + im[b] * wr;\n                re[b] = re[a] - vr;\n                im[b] = im[a] - vi;\n                re[a] += vr;\n                im[a] += vi;\n            }\n        }\n    }\n}\n\n/**\n * In-place unnormalized DCT-IV of `input` (length must be a power of two ≥ 4)\n * into `output`.\n *\n * @param {Float32Array} input Input samples (not modified).\n * @param {Float32Array} output Output buffer of the same length.\n */\nexport function dctIv(input, output)\n{\n    const m = input.length;\n    const h = m >> 1;\n    const tables = getTables(m);\n    const { preCos, preSin, postCos, postSin, workRe, workIm } = tables;\n\n    for (let i = 0; i < h; i++)\n    {\n        const xr = input[2 * i];\n        const xi = input[m - 1 - 2 * i];\n        workRe[i] = xr * preCos[i] + xi * preSin[i];\n        workIm[i] = xi * preCos[i] - xr * preSin[i];\n    }\n    fftInPlace(workRe, workIm, tables);\n    for (let k = 0; k < h; k++)\n    {\n        const wr = workRe[k] * postCos[k] + workIm[k] * postSin[k];\n        const wi = workIm[k] * postCos[k] - workRe[k] * postSin[k];\n        output[2 * k] = wr;\n        output[m - 1 - 2 * k] = -wi;\n    }\n}\n\n/**\n * IMDCT: `spectrum` (n/2 coefficients) to `output` (n samples).\n *\n * Matches y[i] = Σ x[j]·cos(π/(2n)·(2i+1+n/2)·(2j+1)).\n *\n * @param {Float32Array} spectrum Input coefficients (length n/2, not modified).\n * @param {Float32Array} output Output buffer (length n).\n * @param {number} n Block size (power of two).\n */\nexport function imdct(spectrum, output, n)\n{\n    const n2 = n >> 1;\n    const n4 = n >> 2;\n    const n34 = n - n4;\n    const tables = getTables(n2);\n    const t = tables.dct;\n    dctIv(spectrum, t);\n    for (let i = 0; i < n4; i++) output[i] = t[i + n4];\n    for (let i = n4; i < n34; i++) output[i] = -t[n34 - i - 1];\n    for (let i = n34; i < n; i++) output[i] = -t[i - n34];\n}\n\nconst windowCache = new Map();\n\n/**\n * The Vorbis synthesis window slope of `length` samples:\n * sin(π/2 · sin²(π·(i+0.5)/(2·length))).\n *\n * @param {number} length Slope length in samples.\n * @returns {Float32Array} Cached window slope values.\n */\nexport function vorbisWindowSlope(length)\n{\n    let slope = windowCache.get(length);\n    if (slope) return slope;\n    slope = new Float32Array(length);\n    for (let i = 0; i < length; i++)\n    {\n        const s = Math.sin((Math.PI * (i + 0.5)) / (2 * length));\n        slope[i] = Math.sin((Math.PI / 2) * s * s);\n    }\n    windowCache.set(length, slope);\n    return slope;\n}\n"],"names":["tableCache","Map","getTables","m","tables","get","h","preCos","Float32Array","preSin","postCos","postSin","i","pre","Math","PI","cos","sin","post","bitrev","Uint32Array","j","bit","twiddles","len","half","angle","push","workRe","workIm","dct","set","fftInPlace","re","im","n","length","tr","ti","wr","wi","a","b","vr","vi","dctIv","input","output","xr","xi","k","imdct","spectrum","n2","n4","n34","t","windowCache","vorbisWindowSlope","slope","s"],"mappings":"AAAA;AACA;AACA;AACA;AACA;AACA;AACA;AACA;AACA;;AAEA,MAAMA,UAAU,GAAG,IAAIC,GAAG,EAAE;AAE5B,SAASC,SAASA,CAACC,CAAC,EACpB;AACI,EAAA,IAAIC,MAAM,GAAGJ,UAAU,CAACK,GAAG,CAACF,CAAC,CAAC;EAC9B,IAAIC,MAAM,EAAE,OAAOA,MAAM;AAEzB,EAAA,MAAME,CAAC,GAAGH,CAAC,IAAI,CAAC;AAChB,EAAA,MAAMI,MAAM,GAAG,IAAIC,YAAY,CAACF,CAAC,CAAC;AAClC,EAAA,MAAMG,MAAM,GAAG,IAAID,YAAY,CAACF,CAAC,CAAC;AAClC,EAAA,MAAMI,OAAO,GAAG,IAAIF,YAAY,CAACF,CAAC,CAAC;AACnC,EAAA,MAAMK,OAAO,GAAG,IAAIH,YAAY,CAACF,CAAC,CAAC;EACnC,KAAK,IAAIM,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGN,CAAC,EAAEM,CAAC,EAAE,EAC1B;IACI,MAAMC,GAAG,GAAIC,IAAI,CAACC,EAAE,GAAGH,CAAC,GAAIT,CAAC;IAC7BI,MAAM,CAACK,CAAC,CAAC,GAAGE,IAAI,CAACE,GAAG,CAACH,GAAG,CAAC;IACzBJ,MAAM,CAACG,CAAC,CAAC,GAAGE,IAAI,CAACG,GAAG,CAACJ,GAAG,CAAC;AACzB,IAAA,MAAMK,IAAI,GAAIJ,IAAI,CAACC,EAAE,IAAI,CAAC,GAAGH,CAAC,GAAG,CAAC,CAAC,IAAK,CAAC,GAAGT,CAAC,CAAC;IAC9CO,OAAO,CAACE,CAAC,CAAC,GAAGE,IAAI,CAACE,GAAG,CAACE,IAAI,CAAC;IAC3BP,OAAO,CAACC,CAAC,CAAC,GAAGE,IAAI,CAACG,GAAG,CAACC,IAAI,CAAC;AAC/B,EAAA;AAEA,EAAA,MAAMC,MAAM,GAAG,IAAIC,WAAW,CAACd,CAAC,CAAC;AACjC,EAAA,KAAK,IAAIM,CAAC,GAAG,CAAC,EAAES,CAAC,GAAG,CAAC,EAAET,CAAC,GAAGN,CAAC,EAAEM,CAAC,EAAE,EACjC;AACI,IAAA,IAAIU,GAAG,GAAGhB,CAAC,IAAI,CAAC;IAChB,OAAOe,CAAC,GAAGC,GAAG,EAAEA,GAAG,KAAK,CAAC,EAAED,CAAC,IAAIC,GAAG;AACnCD,IAAAA,CAAC,IAAIC,GAAG;AACRH,IAAAA,MAAM,CAACP,CAAC,CAAC,GAAGS,CAAC;AACjB,EAAA;EAEA,MAAME,QAAQ,GAAG,EAAE;AACnB,EAAA,KAAK,IAAIC,GAAG,GAAG,CAAC,EAAEA,GAAG,IAAIlB,CAAC,EAAEkB,GAAG,KAAK,CAAC,EACrC;AACI,IAAA,MAAMC,IAAI,GAAGD,GAAG,IAAI,CAAC;AACrB,IAAA,MAAMR,GAAG,GAAG,IAAIR,YAAY,CAACiB,IAAI,CAAC;AAClC,IAAA,MAAMR,GAAG,GAAG,IAAIT,YAAY,CAACiB,IAAI,CAAC;IAClC,KAAK,IAAIJ,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGI,IAAI,EAAEJ,CAAC,EAAE,EAC7B;MACI,MAAMK,KAAK,GAAI,EAAE,GAAGZ,IAAI,CAACC,EAAE,GAAGM,CAAC,GAAIG,GAAG;MACtCR,GAAG,CAACK,CAAC,CAAC,GAAGP,IAAI,CAACE,GAAG,CAACU,KAAK,CAAC;MACxBT,GAAG,CAACI,CAAC,CAAC,GAAGP,IAAI,CAACG,GAAG,CAACS,KAAK,CAAC;AAC5B,IAAA;IACAH,QAAQ,CAACI,IAAI,CAAC;MAAEH,GAAG;MAAER,GAAG;AAAEC,MAAAA;AAAI,KAAC,CAAC;AACpC,EAAA;AAEAb,EAAAA,MAAM,GAAG;IACLG,MAAM;IACNE,MAAM;IACNC,OAAO;IACPC,OAAO;IACPQ,MAAM;IACNI,QAAQ;AACRK,IAAAA,MAAM,EAAE,IAAIpB,YAAY,CAACF,CAAC,CAAC;AAC3BuB,IAAAA,MAAM,EAAE,IAAIrB,YAAY,CAACF,CAAC,CAAC;AAC3BwB,IAAAA,GAAG,EAAE,IAAItB,YAAY,CAACL,CAAC;GAC1B;AACDH,EAAAA,UAAU,CAAC+B,GAAG,CAAC5B,CAAC,EAAEC,MAAM,CAAC;AACzB,EAAA,OAAOA,MAAM;AACjB;AAEA,SAAS4B,UAAUA,CAACC,EAAE,EAAEC,EAAE,EAAE9B,MAAM,EAClC;AACI,EAAA,MAAM+B,CAAC,GAAGF,EAAE,CAACG,MAAM;AACnB,EAAA,MAAMjB,MAAM,GAAGf,MAAM,CAACe,MAAM;EAC5B,KAAK,IAAIP,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGuB,CAAC,EAAEvB,CAAC,EAAE,EAC1B;AACI,IAAA,MAAMS,CAAC,GAAGF,MAAM,CAACP,CAAC,CAAC;IACnB,IAAIA,CAAC,GAAGS,CAAC,EACT;AACI,MAAA,MAAMgB,EAAE,GAAGJ,EAAE,CAACrB,CAAC,CAAC;AAAEqB,MAAAA,EAAE,CAACrB,CAAC,CAAC,GAAGqB,EAAE,CAACZ,CAAC,CAAC;AAAEY,MAAAA,EAAE,CAACZ,CAAC,CAAC,GAAGgB,EAAE;AAC3C,MAAA,MAAMC,EAAE,GAAGJ,EAAE,CAACtB,CAAC,CAAC;AAAEsB,MAAAA,EAAE,CAACtB,CAAC,CAAC,GAAGsB,EAAE,CAACb,CAAC,CAAC;AAAEa,MAAAA,EAAE,CAACb,CAAC,CAAC,GAAGiB,EAAE;AAC/C,IAAA;AACJ,EAAA;AACA,EAAA,KAAK,MAAM;IAAEd,GAAG;IAAER,GAAG;AAAEC,IAAAA;AAAI,GAAC,IAAIb,MAAM,CAACmB,QAAQ,EAC/C;AACI,IAAA,MAAME,IAAI,GAAGD,GAAG,IAAI,CAAC;AACrB,IAAA,KAAK,IAAIZ,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGuB,CAAC,EAAEvB,CAAC,IAAIY,GAAG,EAC/B;MACI,KAAK,IAAIH,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGI,IAAI,EAAEJ,CAAC,EAAE,EAC7B;AACI,QAAA,MAAMkB,EAAE,GAAGvB,GAAG,CAACK,CAAC,CAAC;AACjB,QAAA,MAAMmB,EAAE,GAAGvB,GAAG,CAACI,CAAC,CAAC;AACjB,QAAA,MAAMoB,CAAC,GAAG7B,CAAC,GAAGS,CAAC;AACf,QAAA,MAAMqB,CAAC,GAAGD,CAAC,GAAGhB,IAAI;AAClB,QAAA,MAAMkB,EAAE,GAAGV,EAAE,CAACS,CAAC,CAAC,GAAGH,EAAE,GAAGL,EAAE,CAACQ,CAAC,CAAC,GAAGF,EAAE;AAClC,QAAA,MAAMI,EAAE,GAAGX,EAAE,CAACS,CAAC,CAAC,GAAGF,EAAE,GAAGN,EAAE,CAACQ,CAAC,CAAC,GAAGH,EAAE;QAClCN,EAAE,CAACS,CAAC,CAAC,GAAGT,EAAE,CAACQ,CAAC,CAAC,GAAGE,EAAE;QAClBT,EAAE,CAACQ,CAAC,CAAC,GAAGR,EAAE,CAACO,CAAC,CAAC,GAAGG,EAAE;AAClBX,QAAAA,EAAE,CAACQ,CAAC,CAAC,IAAIE,EAAE;AACXT,QAAAA,EAAE,CAACO,CAAC,CAAC,IAAIG,EAAE;AACf,MAAA;AACJ,IAAA;AACJ,EAAA;AACJ;;AAEA;AACA;AACA;AACA;AACA;AACA;AACA;AACO,SAASC,KAAKA,CAACC,KAAK,EAAEC,MAAM,EACnC;AACI,EAAA,MAAM5C,CAAC,GAAG2C,KAAK,CAACV,MAAM;AACtB,EAAA,MAAM9B,CAAC,GAAGH,CAAC,IAAI,CAAC;AAChB,EAAA,MAAMC,MAAM,GAAGF,SAAS,CAACC,CAAC,CAAC;EAC3B,MAAM;IAAEI,MAAM;IAAEE,MAAM;IAAEC,OAAO;IAAEC,OAAO;IAAEiB,MAAM;AAAEC,IAAAA;AAAO,GAAC,GAAGzB,MAAM;EAEnE,KAAK,IAAIQ,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGN,CAAC,EAAEM,CAAC,EAAE,EAC1B;AACI,IAAA,MAAMoC,EAAE,GAAGF,KAAK,CAAC,CAAC,GAAGlC,CAAC,CAAC;IACvB,MAAMqC,EAAE,GAAGH,KAAK,CAAC3C,CAAC,GAAG,CAAC,GAAG,CAAC,GAAGS,CAAC,CAAC;AAC/BgB,IAAAA,MAAM,CAAChB,CAAC,CAAC,GAAGoC,EAAE,GAAGzC,MAAM,CAACK,CAAC,CAAC,GAAGqC,EAAE,GAAGxC,MAAM,CAACG,CAAC,CAAC;AAC3CiB,IAAAA,MAAM,CAACjB,CAAC,CAAC,GAAGqC,EAAE,GAAG1C,MAAM,CAACK,CAAC,CAAC,GAAGoC,EAAE,GAAGvC,MAAM,CAACG,CAAC,CAAC;AAC/C,EAAA;AACAoB,EAAAA,UAAU,CAACJ,MAAM,EAAEC,MAAM,EAAEzB,MAAM,CAAC;EAClC,KAAK,IAAI8C,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAG5C,CAAC,EAAE4C,CAAC,EAAE,EAC1B;IACI,MAAMX,EAAE,GAAGX,MAAM,CAACsB,CAAC,CAAC,GAAGxC,OAAO,CAACwC,CAAC,CAAC,GAAGrB,MAAM,CAACqB,CAAC,CAAC,GAAGvC,OAAO,CAACuC,CAAC,CAAC;IAC1D,MAAMV,EAAE,GAAGX,MAAM,CAACqB,CAAC,CAAC,GAAGxC,OAAO,CAACwC,CAAC,CAAC,GAAGtB,MAAM,CAACsB,CAAC,CAAC,GAAGvC,OAAO,CAACuC,CAAC,CAAC;AAC1DH,IAAAA,MAAM,CAAC,CAAC,GAAGG,CAAC,CAAC,GAAGX,EAAE;IAClBQ,MAAM,CAAC5C,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG+C,CAAC,CAAC,GAAG,CAACV,EAAE;AAC/B,EAAA;AACJ;;AAEA;AACA;AACA;AACA;AACA;AACA;AACA;AACA;AACA;AACO,SAASW,KAAKA,CAACC,QAAQ,EAAEL,MAAM,EAAEZ,CAAC,EACzC;AACI,EAAA,MAAMkB,EAAE,GAAGlB,CAAC,IAAI,CAAC;AACjB,EAAA,MAAMmB,EAAE,GAAGnB,CAAC,IAAI,CAAC;AACjB,EAAA,MAAMoB,GAAG,GAAGpB,CAAC,GAAGmB,EAAE;AAClB,EAAA,MAAMlD,MAAM,GAAGF,SAAS,CAACmD,EAAE,CAAC;AAC5B,EAAA,MAAMG,CAAC,GAAGpD,MAAM,CAAC0B,GAAG;AACpBe,EAAAA,KAAK,CAACO,QAAQ,EAAEI,CAAC,CAAC;EAClB,KAAK,IAAI5C,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAG0C,EAAE,EAAE1C,CAAC,EAAE,EAAEmC,MAAM,CAACnC,CAAC,CAAC,GAAG4C,CAAC,CAAC5C,CAAC,GAAG0C,EAAE,CAAC;EAClD,KAAK,IAAI1C,CAAC,GAAG0C,EAAE,EAAE1C,CAAC,GAAG2C,GAAG,EAAE3C,CAAC,EAAE,EAAEmC,MAAM,CAACnC,CAAC,CAAC,GAAG,CAAC4C,CAAC,CAACD,GAAG,GAAG3C,CAAC,GAAG,CAAC,CAAC;EAC1D,KAAK,IAAIA,CAAC,GAAG2C,GAAG,EAAE3C,CAAC,GAAGuB,CAAC,EAAEvB,CAAC,EAAE,EAAEmC,MAAM,CAACnC,CAAC,CAAC,GAAG,CAAC4C,CAAC,CAAC5C,CAAC,GAAG2C,GAAG,CAAC;AACzD;AAEA,MAAME,WAAW,GAAG,IAAIxD,GAAG,EAAE;;AAE7B;AACA;AACA;AACA;AACA;AACA;AACA;AACO,SAASyD,iBAAiBA,CAACtB,MAAM,EACxC;AACI,EAAA,IAAIuB,KAAK,GAAGF,WAAW,CAACpD,GAAG,CAAC+B,MAAM,CAAC;EACnC,IAAIuB,KAAK,EAAE,OAAOA,KAAK;AACvBA,EAAAA,KAAK,GAAG,IAAInD,YAAY,CAAC4B,MAAM,CAAC;EAChC,KAAK,IAAIxB,CAAC,GAAG,CAAC,EAAEA,CAAC,GAAGwB,MAAM,EAAExB,CAAC,EAAE,EAC/B;AACI,IAAA,MAAMgD,CAAC,GAAG9C,IAAI,CAACG,GAAG,CAAEH,IAAI,CAACC,EAAE,IAAIH,CAAC,GAAG,GAAG,CAAC,IAAK,CAAC,GAAGwB,MAAM,CAAC,CAAC;AACxDuB,IAAAA,KAAK,CAAC/C,CAAC,CAAC,GAAGE,IAAI,CAACG,GAAG,CAAEH,IAAI,CAACC,EAAE,GAAG,CAAC,GAAI6C,CAAC,GAAGA,CAAC,CAAC;AAC9C,EAAA;AACAH,EAAAA,WAAW,CAAC1B,GAAG,CAACK,MAAM,EAAEuB,KAAK,CAAC;AAC9B,EAAA,OAAOA,KAAK;AAChB;;;;"}