@cortex-js/compute-engine 0.53.0 → 0.53.1

package/dist/types/compute-engine/compilation/gpu-target.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { CompileTarget, CompiledOperators, CompiledFunctions, LanguageTarget, CompilationOptions, CompilationResult } from './types';
 /**
  * GPU shader operators shared by GLSL and WGSL.
@@ -41,6 +41,73 @@ export declare const GPU_GAMMA_PREAMBLE = "\nfloat _gpu_gamma(float z) {\n  cons
  * Maximum error: |epsilon(x)| <= 1.5e-7.
  */
 export declare const GPU_ERF_PREAMBLE = "\nfloat _gpu_erf(float x) {\n  float ax = abs(x);\n  float t = 1.0 / (1.0 + 0.3275911 * ax);\n  float y = ((((1.061405429 * t - 1.453152027) * t + 1.421413741) * t - 0.284496736) * t + 0.254829592) * t;\n  float result = 1.0 - y * exp(-ax * ax);\n  return x < 0.0 ? -result : result;\n}\n\nfloat _gpu_erfinv(float x) {\n  float pi = 3.14159265358979;\n  float x2 = x * x;\n  float x3 = x * x2;\n  float x5 = x3 * x2;\n  float x7 = x5 * x2;\n  float x9 = x7 * x2;\n  return sqrt(pi) * 0.5 * (x + (pi / 12.0) * x3 + (7.0 * pi * pi / 480.0) * x5 + (127.0 * pi * pi * pi / 40320.0) * x7 + (4369.0 * pi * pi * pi * pi / 5806080.0) * x9);\n}\n";
+/**
+ * GPU Heaviside step function preamble (GLSL syntax).
+ * Returns 0 for x<0, 0.5 at x=0, 1 for x>0.
+ */
+export declare const GPU_HEAVISIDE_PREAMBLE_GLSL = "\nfloat _gpu_heaviside(float x) {\n  if (x < 0.0) return 0.0;\n  if (x > 0.0) return 1.0;\n  return 0.5;\n}\n";
+/**
+ * GPU Heaviside step function preamble (WGSL syntax).
+ */
+export declare const GPU_HEAVISIDE_PREAMBLE_WGSL = "\nfn _gpu_heaviside(x: f32) -> f32 {\n  if (x < 0.0) { return 0.0; }\n  if (x > 0.0) { return 1.0; }\n  return 0.5;\n}\n";
+/**
+ * GPU sinc function preamble (GLSL syntax).
+ * sinc(x) = sin(x)/x, sinc(0) = 1.
+ */
+export declare const GPU_SINC_PREAMBLE_GLSL = "\nfloat _gpu_sinc(float x) {\n  if (abs(x) < 1e-10) return 1.0;\n  return sin(x) / x;\n}\n";
+/**
+ * GPU sinc function preamble (WGSL syntax).
+ */
+export declare const GPU_SINC_PREAMBLE_WGSL = "\nfn _gpu_sinc(x: f32) -> f32 {\n  if (abs(x) < 1e-10) { return 1.0; }\n  return sin(x) / x;\n}\n";
+/**
+ * GPU Horner polynomial evaluation helper (GLSL syntax).
+ * Shared by FresnelC and FresnelS preambles.
+ */
+export declare const GPU_POLEVL_PREAMBLE_GLSL = "\nfloat _gpu_polevl(float x, float c[12], int n) {\n  float ans = c[0];\n  for (int i = 1; i < n; i++) ans = ans * x + c[i];\n  return ans;\n}\n";
+/**
+ * GPU Horner polynomial evaluation helper (WGSL syntax).
+ */
+export declare const GPU_POLEVL_PREAMBLE_WGSL = "\nfn _gpu_polevl(x: f32, c: array<f32, 12>, n: i32) -> f32 {\n  var ans = c[0];\n  for (var i: i32 = 1; i < n; i++) { ans = ans * x + c[i]; }\n  return ans;\n}\n";
+/**
+ * GPU Fresnel cosine integral preamble (GLSL syntax).
+ *
+ * C(x) = integral from 0 to x of cos(pi*t^2/2) dt.
+ * Uses rational Chebyshev approximation (Cephes/scipy) with three regions:
+ * |x|<1.6, 1.6<=|x|<36, |x|>=36.
+ * Requires _gpu_polevl preamble.
+ */
+export declare const GPU_FRESNELC_PREAMBLE_GLSL = "\nfloat _gpu_fresnelC(float x_in) {\n  float sgn = x_in < 0.0 ? -1.0 : 1.0;\n  float x = abs(x_in);\n\n  if (x < 1.6) {\n    float x2 = x * x;\n    float t = x2 * x2;\n    float cn[6] = float[6](\n      -4.98843114573573548651e-8, 9.50428062829859605134e-6,\n      -6.45191435683965050962e-4, 1.88843319396703850064e-2,\n      -2.05525900955013891793e-1, 9.99999999999999998822e-1\n    );\n    float cd[7] = float[7](\n      3.99982968972495980367e-12, 9.15439215774657478799e-10,\n      1.25001862479598821474e-7, 1.22262789024179030997e-5,\n      8.68029542941784300606e-4, 4.12142090722199792936e-2, 1.0\n    );\n    return sgn * x * _gpu_polevl(t, cn, 6) / _gpu_polevl(t, cd, 7);\n  }\n\n  if (x < 36.0) {\n    float x2 = x * x;\n    float t = 3.14159265358979 * x2;\n    float u = 1.0 / (t * t);\n    float fn[10] = float[10](\n      4.21543555043677546506e-1, 1.43407919780758885261e-1,\n      1.15220955073585758835e-2, 3.450179397825740279e-4,\n      4.63613749287867322088e-6, 3.05568983790257605827e-8,\n      1.02304514164907233465e-10, 1.72010743268161828879e-13,\n      1.34283276233062758925e-16, 3.76329711269987889006e-20\n    );\n    float fd[11] = float[11](\n      1.0, 7.51586398353378947175e-1,\n      1.16888925859191382142e-1, 6.44051526508858611005e-3,\n      1.55934409164153020873e-4, 1.8462756734893054587e-6,\n      1.12699224763999035261e-8, 3.60140029589371370404e-11,\n      5.8875453362157841001e-14, 4.52001434074129701496e-17,\n      1.25443237090011264384e-20\n    );\n    float gn[11] = float[11](\n      5.04442073643383265887e-1, 1.97102833525523411709e-1,\n      1.87648584092575249293e-2, 6.84079380915393090172e-4,\n      1.15138826111884280931e-5, 9.82852443688422223854e-8,\n      4.45344415861750144738e-10, 1.08268041139020870318e-12,\n      1.37555460633261799868e-15, 8.36354435630677421531e-19,\n      1.86958710162783235106e-22\n    );\n    float gd[12] = float[12](\n      1.0, 1.47495759925128324529,\n      3.37748989120019970451e-1, 2.53603741420338795122e-2,\n      8.14679107184306179049e-4, 1.27545075667729118702e-5,\n      1.04314589657571990585e-7, 4.60680728515232032307e-10,\n      1.10273215066240270757e-12, 1.38796531259578871258e-15,\n      8.39158816283118707363e-19, 1.86958710162783236342e-22\n    );\n    float f = 1.0 - u * _gpu_polevl(u, fn, 10) / _gpu_polevl(u, fd, 11);\n    float g = (1.0 / t) * _gpu_polevl(u, gn, 11) / _gpu_polevl(u, gd, 12);\n    float z = 1.5707963267948966 * x2;\n    float c = cos(z);\n    float s = sin(z);\n    return sgn * (0.5 + (f * s - g * c) / (3.14159265358979 * x));\n  }\n\n  return sgn * 0.5;\n}\n";
+/**
+ * GPU Fresnel cosine integral preamble (WGSL syntax).
+ * Requires _gpu_polevl preamble.
+ */
+export declare const GPU_FRESNELC_PREAMBLE_WGSL = "\nfn _gpu_fresnelC(x_in: f32) -> f32 {\n  let sgn: f32 = select(1.0, -1.0, x_in < 0.0);\n  let x = abs(x_in);\n\n  if (x < 1.6) {\n    let x2 = x * x;\n    let t = x2 * x2;\n    var cn = array<f32, 12>(\n      -4.98843114573573548651e-8, 9.50428062829859605134e-6,\n      -6.45191435683965050962e-4, 1.88843319396703850064e-2,\n      -2.05525900955013891793e-1, 9.99999999999999998822e-1,\n      0.0, 0.0, 0.0, 0.0, 0.0, 0.0\n    );\n    var cd = array<f32, 12>(\n      3.99982968972495980367e-12, 9.15439215774657478799e-10,\n      1.25001862479598821474e-7, 1.22262789024179030997e-5,\n      8.68029542941784300606e-4, 4.12142090722199792936e-2, 1.0,\n      0.0, 0.0, 0.0, 0.0, 0.0\n    );\n    return sgn * x * _gpu_polevl(t, cn, 6) / _gpu_polevl(t, cd, 7);\n  }\n\n  if (x < 36.0) {\n    let x2 = x * x;\n    let t = 3.14159265358979 * x2;\n    let u = 1.0 / (t * t);\n    var fn = array<f32, 12>(\n      4.21543555043677546506e-1, 1.43407919780758885261e-1,\n      1.15220955073585758835e-2, 3.450179397825740279e-4,\n      4.63613749287867322088e-6, 3.05568983790257605827e-8,\n      1.02304514164907233465e-10, 1.72010743268161828879e-13,\n      1.34283276233062758925e-16, 3.76329711269987889006e-20,\n      0.0, 0.0\n    );\n    var fd = array<f32, 12>(\n      1.0, 7.51586398353378947175e-1,\n      1.16888925859191382142e-1, 6.44051526508858611005e-3,\n      1.55934409164153020873e-4, 1.8462756734893054587e-6,\n      1.12699224763999035261e-8, 3.60140029589371370404e-11,\n      5.8875453362157841001e-14, 4.52001434074129701496e-17,\n      1.25443237090011264384e-20, 0.0\n    );\n    var gn = array<f32, 12>(\n      5.04442073643383265887e-1, 1.97102833525523411709e-1,\n      1.87648584092575249293e-2, 6.84079380915393090172e-4,\n      1.15138826111884280931e-5, 9.82852443688422223854e-8,\n      4.45344415861750144738e-10, 1.08268041139020870318e-12,\n      1.37555460633261799868e-15, 8.36354435630677421531e-19,\n      1.86958710162783235106e-22, 0.0\n    );\n    var gd = array<f32, 12>(\n      1.0, 1.47495759925128324529,\n      3.37748989120019970451e-1, 2.53603741420338795122e-2,\n      8.14679107184306179049e-4, 1.27545075667729118702e-5,\n      1.04314589657571990585e-7, 4.60680728515232032307e-10,\n      1.10273215066240270757e-12, 1.38796531259578871258e-15,\n      8.39158816283118707363e-19, 1.86958710162783236342e-22\n    );\n    let f = 1.0 - u * _gpu_polevl(u, fn, 10) / _gpu_polevl(u, fd, 11);\n    let g = (1.0 / t) * _gpu_polevl(u, gn, 11) / _gpu_polevl(u, gd, 12);\n    let z = 1.5707963267948966 * x2;\n    let c = cos(z);\n    let s = sin(z);\n    return sgn * (0.5 + (f * s - g * c) / (3.14159265358979 * x));\n  }\n\n  return sgn * 0.5;\n}\n";
+/**
+ * GPU Fresnel sine integral preamble (GLSL syntax).
+ *
+ * S(x) = integral from 0 to x of sin(pi*t^2/2) dt.
+ * Uses rational Chebyshev approximation (Cephes/scipy) with three regions.
+ * Requires _gpu_polevl preamble.
+ */
+export declare const GPU_FRESNELS_PREAMBLE_GLSL = "\nfloat _gpu_fresnelS(float x_in) {\n  float sgn = x_in < 0.0 ? -1.0 : 1.0;\n  float x = abs(x_in);\n\n  if (x < 1.6) {\n    float x2 = x * x;\n    float t = x2 * x2;\n    float sn[6] = float[6](\n      -2.99181919401019853726e3, 7.08840045257738576863e5,\n      -6.29741486205862506537e7, 2.54890880573376359104e9,\n      -4.42979518059697779103e10, 3.18016297876567817986e11\n    );\n    float sd[7] = float[7](\n      1.0, 2.81376268889994315696e2, 4.55847810806532581675e4,\n      5.1734388877009640073e6, 4.19320245898111231129e8, 2.2441179564534092094e10,\n      6.07366389490084914091e11\n    );\n    return sgn * x * x2 * _gpu_polevl(t, sn, 6) / _gpu_polevl(t, sd, 7);\n  }\n\n  if (x < 36.0) {\n    float x2 = x * x;\n    float t = 3.14159265358979 * x2;\n    float u = 1.0 / (t * t);\n    float fn[10] = float[10](\n      4.21543555043677546506e-1, 1.43407919780758885261e-1,\n      1.15220955073585758835e-2, 3.450179397825740279e-4,\n      4.63613749287867322088e-6, 3.05568983790257605827e-8,\n      1.02304514164907233465e-10, 1.72010743268161828879e-13,\n      1.34283276233062758925e-16, 3.76329711269987889006e-20\n    );\n    float fd[11] = float[11](\n      1.0, 7.51586398353378947175e-1,\n      1.16888925859191382142e-1, 6.44051526508858611005e-3,\n      1.55934409164153020873e-4, 1.8462756734893054587e-6,\n      1.12699224763999035261e-8, 3.60140029589371370404e-11,\n      5.8875453362157841001e-14, 4.52001434074129701496e-17,\n      1.25443237090011264384e-20\n    );\n    float gn[11] = float[11](\n      5.04442073643383265887e-1, 1.97102833525523411709e-1,\n      1.87648584092575249293e-2, 6.84079380915393090172e-4,\n      1.15138826111884280931e-5, 9.82852443688422223854e-8,\n      4.45344415861750144738e-10, 1.08268041139020870318e-12,\n      1.37555460633261799868e-15, 8.36354435630677421531e-19,\n      1.86958710162783235106e-22\n    );\n    float gd[12] = float[12](\n      1.0, 1.47495759925128324529,\n      3.37748989120019970451e-1, 2.53603741420338795122e-2,\n      8.14679107184306179049e-4, 1.27545075667729118702e-5,\n      1.04314589657571990585e-7, 4.60680728515232032307e-10,\n      1.10273215066240270757e-12, 1.38796531259578871258e-15,\n      8.39158816283118707363e-19, 1.86958710162783236342e-22\n    );\n    float f = 1.0 - u * _gpu_polevl(u, fn, 10) / _gpu_polevl(u, fd, 11);\n    float g = (1.0 / t) * _gpu_polevl(u, gn, 11) / _gpu_polevl(u, gd, 12);\n    float z = 1.5707963267948966 * x2;\n    float c = cos(z);\n    float s = sin(z);\n    return sgn * (0.5 - (f * c + g * s) / (3.14159265358979 * x));\n  }\n\n  return sgn * 0.5;\n}\n";
+/**
+ * GPU Fresnel sine integral preamble (WGSL syntax).
+ * Requires _gpu_polevl preamble.
+ */
+export declare const GPU_FRESNELS_PREAMBLE_WGSL = "\nfn _gpu_fresnelS(x_in: f32) -> f32 {\n  let sgn: f32 = select(1.0, -1.0, x_in < 0.0);\n  let x = abs(x_in);\n\n  if (x < 1.6) {\n    let x2 = x * x;\n    let t = x2 * x2;\n    var sn = array<f32, 12>(\n      -2.99181919401019853726e3, 7.08840045257738576863e5,\n      -6.29741486205862506537e7, 2.54890880573376359104e9,\n      -4.42979518059697779103e10, 3.18016297876567817986e11,\n      0.0, 0.0, 0.0, 0.0, 0.0, 0.0\n    );\n    var sd = array<f32, 12>(\n      1.0, 2.81376268889994315696e2, 4.55847810806532581675e4,\n      5.1734388877009640073e6, 4.19320245898111231129e8, 2.2441179564534092094e10,\n      6.07366389490084914091e11,\n      0.0, 0.0, 0.0, 0.0, 0.0\n    );\n    return sgn * x * x2 * _gpu_polevl(t, sn, 6) / _gpu_polevl(t, sd, 7);\n  }\n\n  if (x < 36.0) {\n    let x2 = x * x;\n    let t = 3.14159265358979 * x2;\n    let u = 1.0 / (t * t);\n    var fn = array<f32, 12>(\n      4.21543555043677546506e-1, 1.43407919780758885261e-1,\n      1.15220955073585758835e-2, 3.450179397825740279e-4,\n      4.63613749287867322088e-6, 3.05568983790257605827e-8,\n      1.02304514164907233465e-10, 1.72010743268161828879e-13,\n      1.34283276233062758925e-16, 3.76329711269987889006e-20,\n      0.0, 0.0\n    );\n    var fd = array<f32, 12>(\n      1.0, 7.51586398353378947175e-1,\n      1.16888925859191382142e-1, 6.44051526508858611005e-3,\n      1.55934409164153020873e-4, 1.8462756734893054587e-6,\n      1.12699224763999035261e-8, 3.60140029589371370404e-11,\n      5.8875453362157841001e-14, 4.52001434074129701496e-17,\n      1.25443237090011264384e-20, 0.0\n    );\n    var gn = array<f32, 12>(\n      5.04442073643383265887e-1, 1.97102833525523411709e-1,\n      1.87648584092575249293e-2, 6.84079380915393090172e-4,\n      1.15138826111884280931e-5, 9.82852443688422223854e-8,\n      4.45344415861750144738e-10, 1.08268041139020870318e-12,\n      1.37555460633261799868e-15, 8.36354435630677421531e-19,\n      1.86958710162783235106e-22, 0.0\n    );\n    var gd = array<f32, 12>(\n      1.0, 1.47495759925128324529,\n      3.37748989120019970451e-1, 2.53603741420338795122e-2,\n      8.14679107184306179049e-4, 1.27545075667729118702e-5,\n      1.04314589657571990585e-7, 4.60680728515232032307e-10,\n      1.10273215066240270757e-12, 1.38796531259578871258e-15,\n      8.39158816283118707363e-19, 1.86958710162783236342e-22\n    );\n    let f = 1.0 - u * _gpu_polevl(u, fn, 10) / _gpu_polevl(u, fd, 11);\n    let g = (1.0 / t) * _gpu_polevl(u, gn, 11) / _gpu_polevl(u, gd, 12);\n    let z = 1.5707963267948966 * x2;\n    let c = cos(z);\n    let s = sin(z);\n    return sgn * (0.5 - (f * c + g * s) / (3.14159265358979 * x));\n  }\n\n  return sgn * 0.5;\n}\n";
+/**
+ * GPU Bessel J function preamble (GLSL syntax).
+ *
+ * J_n(x) for integer order n. Uses three algorithms:
+ * - Power series for small x (x < 5+n)
+ * - Hankel asymptotic for large x (x > 25+n^2/2)
+ * - Miller's backward recurrence for intermediate x
+ */
+export declare const GPU_BESSELJ_PREAMBLE_GLSL = "\nfloat _gpu_factorial(int n) {\n  float f = 1.0;\n  for (int i = 2; i <= n; i++) f *= float(i);\n  return f;\n}\n\nfloat _gpu_besselJ_series(int n, float x) {\n  float halfX = x / 2.0;\n  float negQ = -(x * x) / 4.0;\n  float term = 1.0;\n  for (int i = 1; i <= n; i++) term /= float(i);\n  float s = term;\n  for (int k = 1; k <= 60; k++) {\n    term *= negQ / (float(k) * float(n + k));\n    s += term;\n    if (abs(term) < abs(s) * 1e-7) break;\n  }\n  return s * pow(halfX, float(n));\n}\n\nfloat _gpu_besselJ_asymptotic(int n, float x) {\n  float mu = 4.0 * float(n) * float(n);\n  float P = 1.0;\n  float Q = 0.0;\n  float ak = 1.0;\n  float e8x = 8.0 * x;\n  for (int k = 1; k <= 12; k++) {\n    float twokm1 = float(2 * k - 1);\n    ak *= mu - twokm1 * twokm1;\n    float denom = _gpu_factorial(k) * pow(e8x, float(k));\n    float contrib = ak / denom;\n    if (k == 1 || k == 3 || k == 5 || k == 7 || k == 9 || k == 11) {\n      // odd k: contributes to Q\n      if (((k - 1) / 2) % 2 == 0) Q += contrib;\n      else Q -= contrib;\n    } else {\n      // even k: contributes to P\n      if ((k / 2) % 2 == 1) P -= contrib;\n      else P += contrib;\n    }\n    if (abs(contrib) < 1e-7) break;\n  }\n  float chi = x - (float(n) / 2.0 + 0.25) * 3.14159265358979;\n  return sqrt(2.0 / (3.14159265358979 * x)) * (P * cos(chi) - Q * sin(chi));\n}\n\nfloat _gpu_besselJ(int n, float x) {\n  if (x == 0.0) return n == 0 ? 1.0 : 0.0;\n  float sgn = 1.0;\n  if (n < 0) {\n    n = -n;\n    if (n % 2 != 0) sgn = -1.0;\n  }\n  if (x < 0.0) {\n    x = -x;\n    if (n % 2 != 0) sgn *= -1.0;\n  }\n  if (x > 25.0 + float(n * n) / 2.0) return sgn * _gpu_besselJ_asymptotic(n, x);\n  if (x < 5.0 + float(n)) return sgn * _gpu_besselJ_series(n, x);\n  // Miller's backward recurrence\n  int M = max(n + 20, int(ceil(x)) + 30);\n  if (M > 200) return sgn * _gpu_besselJ_series(n, x);\n  float vals[201];\n  float jp1 = 0.0;\n  float jk = 1.0;\n  vals[M] = jk;\n  for (int k = M; k >= 1; k--) {\n    float jm1 = (2.0 * float(k) / x) * jk - jp1;\n    jp1 = jk;\n    jk = jm1;\n    vals[k - 1] = jk;\n  }\n  float norm = vals[0];\n  for (int k = 2; k <= M; k += 2) norm += 2.0 * vals[k];\n  return sgn * vals[n] / norm;\n}\n";
+/**
+ * GPU Bessel J function preamble (WGSL syntax).
+ */
+export declare const GPU_BESSELJ_PREAMBLE_WGSL = "\nfn _gpu_factorial(n: i32) -> f32 {\n  var f: f32 = 1.0;\n  for (var i: i32 = 2; i <= n; i++) { f *= f32(i); }\n  return f;\n}\n\nfn _gpu_besselJ_series(n_in: i32, x: f32) -> f32 {\n  let halfX = x / 2.0;\n  let negQ = -(x * x) / 4.0;\n  var term: f32 = 1.0;\n  for (var i: i32 = 1; i <= n_in; i++) { term /= f32(i); }\n  var s = term;\n  for (var k: i32 = 1; k <= 60; k++) {\n    term *= negQ / (f32(k) * f32(n_in + k));\n    s += term;\n    if (abs(term) < abs(s) * 1e-7) { break; }\n  }\n  return s * pow(halfX, f32(n_in));\n}\n\nfn _gpu_besselJ_asymptotic(n_in: i32, x: f32) -> f32 {\n  let mu = 4.0 * f32(n_in) * f32(n_in);\n  var P: f32 = 1.0;\n  var Q: f32 = 0.0;\n  var ak: f32 = 1.0;\n  let e8x = 8.0 * x;\n  for (var k: i32 = 1; k <= 12; k++) {\n    let twokm1 = f32(2 * k - 1);\n    ak *= mu - twokm1 * twokm1;\n    let denom = _gpu_factorial(k) * pow(e8x, f32(k));\n    let contrib = ak / denom;\n    if (k == 1 || k == 3 || k == 5 || k == 7 || k == 9 || k == 11) {\n      if (((k - 1) / 2) % 2 == 0) { Q += contrib; }\n      else { Q -= contrib; }\n    } else {\n      if ((k / 2) % 2 == 1) { P -= contrib; }\n      else { P += contrib; }\n    }\n    if (abs(contrib) < 1e-7) { break; }\n  }\n  let chi = x - (f32(n_in) / 2.0 + 0.25) * 3.14159265358979;\n  return sqrt(2.0 / (3.14159265358979 * x)) * (P * cos(chi) - Q * sin(chi));\n}\n\nfn _gpu_besselJ(n_in: i32, x_in: f32) -> f32 {\n  var n = n_in;\n  var x = x_in;\n  if (x == 0.0) { return select(0.0, 1.0, n == 0); }\n  var sgn: f32 = 1.0;\n  if (n < 0) {\n    n = -n;\n    if (n % 2 != 0) { sgn = -1.0; }\n  }\n  if (x < 0.0) {\n    x = -x;\n    if (n % 2 != 0) { sgn *= -1.0; }\n  }\n  if (x > 25.0 + f32(n * n) / 2.0) { return sgn * _gpu_besselJ_asymptotic(n, x); }\n  if (x < 5.0 + f32(n)) { return sgn * _gpu_besselJ_series(n, x); }\n  // Miller's backward recurrence\n  var M = max(n + 20, i32(ceil(x)) + 30);\n  if (M > 200) { return sgn * _gpu_besselJ_series(n, x); }\n  var vals: array<f32, 201>;\n  var jp1: f32 = 0.0;\n  var jk: f32 = 1.0;\n  vals[M] = jk;\n  for (var k: i32 = M; k >= 1; k--) {\n    let jm1 = (2.0 * f32(k) / x) * jk - jp1;\n    jp1 = jk;\n    jk = jm1;\n    vals[k - 1] = jk;\n  }\n  var norm = vals[0];\n  for (var k2: i32 = 2; k2 <= M; k2 += 2) { norm += 2.0 * vals[k2]; }\n  return sgn * vals[n] / norm;\n}\n";
 /**
  * Fractal preamble (GLSL syntax).
  *

package/dist/types/compute-engine/compilation/interval-javascript-target.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * JavaScript interval arithmetic compilation target
  *
  * Compiles mathematical expressions to JavaScript code using interval arithmetic

package/dist/types/compute-engine/compilation/javascript-target.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { chop, factorial, factorial2, gcd, lcm, limit } from '../numerics/numeric';
 import { gamma, gammaln, erf, erfc, erfInv, beta, digamma, trigamma, polygamma, zeta, lambertW, besselJ, besselY, besselI, besselK, airyAi, airyBi, fresnelS, fresnelC, sinc } from '../numerics/special-functions';
 import { choose } from '../boxed-expression/expand';

package/dist/types/compute-engine/compilation/python-target.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { CompileTarget, CompiledOperators, CompiledFunctions, LanguageTarget, CompilationOptions, CompilationResult } from './types';
 /**
  * Python/NumPy language target implementation

package/dist/types/compute-engine/compilation/types.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 /**
  * Source code in the target language
  */
@@ -43,6 +43,14 @@ export interface CompileTarget<Expr = unknown> {
     preamble: string;
     /** Current indentation level */
     indent: number;
+    /** Format a variable declaration. Default: `let ${name}`.
+     *  `typeHint` is an optional target-language type string (e.g. `'vec2'`)
+     *  inferred from subsequent assignments. */
+    declare?: (name: string, typeHint?: string) => string;
+    /** Format a block expression. Receives compiled statements; the last
+     *  element is the block's return value (without `return` prefix).
+     *  Default: JavaScript IIFE. */
+    block?: (statements: string[]) => string;
     /** Target language identifier (for debugging/logging) */
     language?: string;
 }

package/dist/types/compute-engine/compilation/wgsl-target.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { CompiledFunctions } from './types';
 import { GPUShaderTarget } from './gpu-target';
 /**

package/dist/types/compute-engine/cost-function.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 /**
  * The default cost function, used to determine if a new expression is simpler
  * than the old one.

package/dist/types/compute-engine/engine-assumptions.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { MathJsonSymbol } from '../math-json/types';
 export declare function ask(ce: IComputeEngine, pattern: Expression): BoxedSubstitution[];
 export declare function verify(ce: IComputeEngine, query: Expression): boolean | undefined;

package/dist/types/compute-engine/engine-cache.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
     private _entries;
     getOrBuild<T>(cacheName: string, build: () => T, purge?: (t: T) => T | undefined): T;
     invalidate(cacheName: string): void;

package/dist/types/compute-engine/engine-common-symbols.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export type CommonSymbolTable = {
     [symbol: string]: null | Expression;
 };

package/dist/types/compute-engine/engine-compilation-targets.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { LanguageTarget } from './compilation/types';
 /**
  * Internal registry for compilation targets.

package/dist/types/compute-engine/engine-configuration-lifecycle.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 type ResetHooks = {
     refreshNumericConstants: () => void;
     resetCommonSymbols: () => void;

package/dist/types/compute-engine/engine-declarations.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { BoxedType } from '../common/type/boxed-type';
 import type { MathJsonSymbol } from '../math-json/types';
 import type { ValueDefinition, OperatorDefinition, AssignValue, Expression, BoxedDefinition, SymbolDefinition, IComputeEngine, Scope } from './global-types';

package/dist/types/compute-engine/engine-expression-entrypoints.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { Decimal } from 'decimal.js';
 import type { MathJsonNumberObject, MathJsonSymbol } from '../math-json/types';
 import { NumericValue } from './numeric-value/types';

package/dist/types/compute-engine/engine-extension-contracts.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { LanguageTarget, CompilationOptions } from './compilation/types';
 export declare function assertCompilationTargetName(name: unknown): string;
 export declare function assertLibraryName(name: unknown): string;

package/dist/types/compute-engine/engine-latex-dictionary-state.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { LatexDictionaryEntry } from './latex-syntax/types';
 type DefaultDictionaryProvider = () => Readonly<LatexDictionaryEntry[]>;
 export declare class EngineLatexDictionaryState {

package/dist/types/compute-engine/engine-library-bootstrap.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { LatexDictionaryEntry, LibraryCategory } from './latex-syntax/types';
 export declare function resolveBootstrapLibraries(libraries?: readonly (string | LibraryDefinition)[]): LibraryDefinition[];
 export declare function loadLibraryDefinitions(engine: ComputeEngine, libraries: readonly LibraryDefinition[]): void;

package/dist/types/compute-engine/engine-numeric-configuration.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { AngularUnit } from './types-definitions';
 export declare class EngineNumericConfiguration {
     private _bignum;

package/dist/types/compute-engine/engine-parse-entrypoint.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { LatexString, ParseLatexOptions } from './latex-syntax/types';
 import type { FormOption } from './types-serialization';
 export type ParseEntrypointOptions = Partial<ParseLatexOptions> & {

package/dist/types/compute-engine/engine-runtime-state.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
     private _timeLimit;
     private _iterationLimit;
     private _recursionLimit;

package/dist/types/compute-engine/engine-scope.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare function pushScope(ce: IComputeEngine, scope?: Scope, name?: string): void;
 export declare function popScope(ce: IComputeEngine): void;
 export declare function pushEvalContext(ce: IComputeEngine, scope: Scope, name?: string): void;

package/dist/types/compute-engine/engine-sequences.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare function declareSequence(ce: IComputeEngine, name: string, def: SequenceDefinition): IComputeEngine;
 export declare function getSequenceStatus(ce: IComputeEngine, name: string): SequenceStatus;
 export declare function getSequence(ce: IComputeEngine, name: string): SequenceInfo | undefined;

package/dist/types/compute-engine/engine-simplification-rules.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 /**
  * Internal holder for simplification rules and their cache-staleness marker.
  *

package/dist/types/compute-engine/engine-startup-coordinator.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { Expression, IComputeEngine as ComputeEngine, LibraryDefinition } from './global-types';
 export type CommonNumberBindings = {
     Zero: Expression;

package/dist/types/compute-engine/engine-type-resolver.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 type ResolverScope = {
     parent: ResolverScope | null;
     types?: Record<string, TypeReference>;

package/dist/types/compute-engine/engine-validation-entrypoints.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { BoxedType } from '../common/type/boxed-type';
 import type { Expression, Metadata, ExpressionInput } from './global-types';
 type ValidationHost = {

package/dist/types/compute-engine/free-functions.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { Type, TypeString } from '../common/type/types';
 import type { LatexString } from './latex-syntax/types';
 import { compile as compileExpr } from './compilation/compile-expression';

package/dist/types/compute-engine/function-utils.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { BoxedDefinition, Expression, IComputeEngine as ComputeEngine, Scope } from './global-types';
 /***
  * ### THEORY OF OPERATIONS

package/dist/types/compute-engine/global-types.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export type * from './types-serialization';
 export type { ValueDefinition, SequenceDefinition, SequenceStatus, SequenceInfo, OEISSequenceInfo, OEISOptions, OperatorDefinition, BaseDefinition, SimplifyOptions, SymbolDefinition, SymbolDefinitions, LibraryDefinition, AngularUnit, Sign, BaseCollectionHandlers, IndexedCollectionHandlers, CollectionHandlers, TaggedValueDefinition, TaggedOperatorDefinition, BoxedDefinition, BoxedBaseDefinition, BoxedValueDefinition, OperatorDefinitionFlags, BoxedOperatorDefinition, } from './types-definitions';
 export type * from './types-evaluation';

package/dist/types/compute-engine/index.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { Decimal } from 'decimal.js';
 import { Type, TypeResolver, TypeString } from '../common/type/types';
 import { BoxedType } from '../common/type/boxed-type';

package/dist/types/compute-engine/interval/arithmetic.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Basic interval arithmetic operations
  *
  * @module interval/arithmetic

package/dist/types/compute-engine/interval/comparison.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Comparison and conditional interval operations
  *
  * @module interval/comparison

package/dist/types/compute-engine/interval/elementary.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Elementary interval functions (sqrt, pow, exp, ln, abs, floor, ceil, min, max, mod)
  *
  * @module interval/elementary

package/dist/types/compute-engine/interval/index.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Interval arithmetic library for reliable function evaluation
  *
  * This module provides interval versions of mathematical operations

package/dist/types/compute-engine/interval/trigonometric.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Trigonometric interval functions
  *
  * @module interval/trigonometric

package/dist/types/compute-engine/interval/types.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Interval arithmetic types for reliable function evaluation
  *
  * @module interval/types

package/dist/types/compute-engine/interval/util.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Utility functions for interval arithmetic
  *
  * @module interval/util

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-algebra.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_ALGEBRA: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-arithmetic.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_ARITHMETIC: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-calculus.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_CALCULUS: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-complex.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_COMPLEX: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-core.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_CORE: LatexDictionary;
 export declare const DELIMITERS_SHORTHAND: {
     '(': string;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-linear-algebra.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_LINEAR_ALGEBRA: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-logic.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_LOGIC: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-other.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_OTHERS: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-relational-operators.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_INEQUALITIES: LatexDictionaryEntry[];

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-sets.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_SETS: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-statistics.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_STATISTICS: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-symbols.d.ts

@@ -1,3 +1,3 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const SYMBOLS: [string, string, number][];
 export declare const DEFINITIONS_SYMBOLS: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-trigonometry.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFINITIONS_TRIGONOMETRY: LatexDictionary;

package/dist/types/compute-engine/latex-syntax/dictionary/definitions-units.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * LaTeX dictionary entries for parsing and serializing physical quantities
  * with units.
  *

package/dist/types/compute-engine/latex-syntax/dictionary/definitions.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { LatexDictionaryEntry } from '../types';
 export type { CommonEntry, IndexedSymbolEntry, IndexedExpressionEntry, IndexedFunctionEntry, IndexedMatchfixEntry, IndexedInfixEntry, IndexedPrefixEntry, IndexedPostfixEntry, IndexedEnvironmentEntry, IndexedLatexDictionaryEntry, IndexedLatexDictionary, } from './indexed-types';
 import type { IndexedSymbolEntry, IndexedExpressionEntry, IndexedFunctionEntry, IndexedMatchfixEntry, IndexedInfixEntry, IndexedPostfixEntry, IndexedEnvironmentEntry, IndexedLatexDictionaryEntry, IndexedLatexDictionary } from './indexed-types';

package/dist/types/compute-engine/latex-syntax/dictionary/indexed-types.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Type definitions for the indexed LaTeX dictionary.
  *
  * These are separated from definitions.ts to break a circular dependency:

package/dist/types/compute-engine/latex-syntax/parse-number.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * Number parsing extracted from the _Parser class for modularity.
  *
  * All functions take a `Parser` interface and a `NumberFormatTokens` config

package/dist/types/compute-engine/latex-syntax/parse-symbol.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { Parser } from './types';
 /** For error handling, if we have a symbol prefix, assume
  * the symbol is invalid (it would have been captured by

package/dist/types/compute-engine/latex-syntax/parse.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { ParseLatexOptions, LatexToken, Terminator, Parser, SymbolTable } from './types';
 import type { IndexedLatexDictionary, IndexedLatexDictionaryEntry, IndexedInfixEntry, IndexedPostfixEntry, IndexedPrefixEntry, IndexedSymbolEntry, IndexedExpressionEntry, IndexedFunctionEntry } from './dictionary/definitions';
 import { BoxedType } from '../../common/type/boxed-type';

package/dist/types/compute-engine/latex-syntax/serialize-dms.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * DMS (Degrees-Minutes-Seconds) serialization utilities.
  */
 export interface DMSComponents {

package/dist/types/compute-engine/latex-syntax/serialize-number.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { NumberSerializationFormat } from './types';
 /**
  * @param expr - A number, can be represented as a string

package/dist/types/compute-engine/latex-syntax/serializer-style.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { DelimiterScale } from './types';
 export declare function getApplyFunctionStyle(_expr: MathJsonExpression, _level: number): DelimiterScale;
 export declare function getGroupStyle(_expr: MathJsonExpression, _level: number): DelimiterScale;

package/dist/types/compute-engine/latex-syntax/serializer.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import { LatexString, SerializeLatexOptions, DelimiterScale } from './types';
 import type { IndexedLatexDictionary, IndexedLatexDictionaryEntry } from './dictionary/definitions';
 export declare class Serializer {

package/dist/types/compute-engine/latex-syntax/tokenizer.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
  * ## Reference
  * TeX source code:
  * {@link  http://tug.org/texlive/devsrc/Build/source/texk/web2c/tex.web | Tex.web}

package/dist/types/compute-engine/latex-syntax/types.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 import type { MathJsonExpression, MathJsonSymbol } from '../../math-json/types';
 import type { TypeString } from '../../common/type/types';
 import { BoxedType } from '../../common/type/boxed-type';

package/dist/types/compute-engine/latex-syntax/utils.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare function asLatexString(s: unknown): string | null;
 export declare function isRelationalOperator(name: string | undefined): boolean;
 export declare function isInequalityOperator(operator: string | undefined): boolean;

package/dist/types/compute-engine/library/arithmetic.d.ts

@@ -1,3 +1,3 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export type CanonicalArithmeticOperators = 'Add' | 'Negate' | 'Multiply' | 'Divide' | 'Power' | 'Sqrt' | 'Root' | 'Ln';
 export declare const ARITHMETIC_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/calculus.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const CALCULUS_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/collections.d.ts

@@ -1,4 +1,4 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const DEFAULT_LINSPACE_COUNT = 50;
 export declare const COLLECTIONS_LIBRARY: SymbolDefinitions;
 /**

package/dist/types/compute-engine/library/colors.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const COLORS_LIBRARY: SymbolDefinitions;

package/dist/types/compute-engine/library/combinatorics.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const COMBINATORICS_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/complex.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const COMPLEX_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/control-structures.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const CONTROL_STRUCTURES_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/core.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const CORE_LIBRARY: SymbolDefinitions[];

package/dist/types/compute-engine/library/fractals.d.ts

@@ -1,2 +1,2 @@
-/* 0.53.0 */
+/* 0.53.1 */
 export declare const FRACTALS_LIBRARY: SymbolDefinitions[];