@cortex-js/compute-engine 0.55.3 → 0.55.4

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 import type { CompileTarget, CompiledOperators, CompiledFunctions, LanguageTarget, CompilationOptions, CompilationResult } from './types';
 /**
  * GPU shader operators shared by GLSL and WGSL.
@@ -108,6 +108,38 @@ export declare const GPU_BESSELJ_PREAMBLE_GLSL = "\nfloat _gpu_factorial(int 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";
+/**
+ * Double-single (DS) emulated double-precision arithmetic preamble (GLSL).
+ *
+ * A "double-single" number is stored as vec2(hi, lo) where value = hi + lo,
+ * giving ~48 bits of mantissa from two 32-bit floats.
+ * Algorithms: Dekker (1971) / Knuth (1997).
+ */
+export declare const GPU_DS_ARITHMETIC_PREAMBLE_GLSL = "\n// Split a float into high and low parts for exact multiplication\nvec2 ds_split(float a) {\n  const float SPLIT = 4097.0; // 2^12 + 1\n  float t = SPLIT * a;\n  float hi = t - (t - a);\n  float lo = a - hi;\n  return vec2(hi, lo);\n}\n\n// Create a double-single from a single float\nvec2 ds_from(float a) {\n  return vec2(a, 0.0);\n}\n\n// Error-free addition (Knuth TwoSum)\nvec2 ds_add(vec2 a, vec2 b) {\n  float s = a.x + b.x;\n  float v = s - a.x;\n  float e = (a.x - (s - v)) + (b.x - v);\n  float lo = (a.y + b.y) + e;\n  float hi = s + lo;\n  lo = lo - (hi - s);\n  return vec2(hi, lo);\n}\n\n// Double-single subtraction\nvec2 ds_sub(vec2 a, vec2 b) {\n  return ds_add(a, vec2(-b.x, -b.y));\n}\n\n// Error-free multiplication (Dekker TwoProduct)\nvec2 ds_mul(vec2 a, vec2 b) {\n  float p = a.x * b.x;\n  vec2 sa = ds_split(a.x);\n  vec2 sb = ds_split(b.x);\n  float err = ((sa.x * sb.x - p) + sa.x * sb.y + sa.y * sb.x) + sa.y * sb.y;\n  err += a.x * b.y + a.y * b.x;\n  float hi = p + err;\n  float lo = err - (hi - p);\n  return vec2(hi, lo);\n}\n\n// Optimized self-multiply\nvec2 ds_sqr(vec2 a) {\n  float p = a.x * a.x;\n  vec2 sa = ds_split(a.x);\n  float err = ((sa.x * sa.x - p) + 2.0 * sa.x * sa.y) + sa.y * sa.y;\n  err += 2.0 * a.x * a.y;\n  float hi = p + err;\n  float lo = err - (hi - p);\n  return vec2(hi, lo);\n}\n\n// Compare magnitude: returns -1, 0, or 1\nfloat ds_cmp(vec2 a, vec2 b) {\n  float d = a.x - b.x;\n  if (d != 0.0) return sign(d);\n  return sign(a.y - b.y);\n}\n";
+/**
+ * Double-single (DS) emulated double-precision arithmetic preamble (WGSL).
+ *
+ * A "double-single" number is stored as vec2f(hi, lo) where value = hi + lo,
+ * giving ~48 bits of mantissa from two 32-bit floats.
+ * Algorithms: Dekker (1971) / Knuth (1997).
+ */
+export declare const GPU_DS_ARITHMETIC_PREAMBLE_WGSL = "\nfn ds_split(a: f32) -> vec2f {\n  const SPLIT: f32 = 4097.0;\n  let t = SPLIT * a;\n  let hi = t - (t - a);\n  let lo = a - hi;\n  return vec2f(hi, lo);\n}\n\nfn ds_from(a: f32) -> vec2f {\n  return vec2f(a, 0.0);\n}\n\nfn ds_add(a: vec2f, b: vec2f) -> vec2f {\n  let s = a.x + b.x;\n  let v = s - a.x;\n  let e = (a.x - (s - v)) + (b.x - v);\n  let lo_t = (a.y + b.y) + e;\n  let hi = s + lo_t;\n  let lo = lo_t - (hi - s);\n  return vec2f(hi, lo);\n}\n\nfn ds_sub(a: vec2f, b: vec2f) -> vec2f {\n  return ds_add(a, vec2f(-b.x, -b.y));\n}\n\nfn ds_mul(a: vec2f, b: vec2f) -> vec2f {\n  let p = a.x * b.x;\n  let sa = ds_split(a.x);\n  let sb = ds_split(b.x);\n  var err = ((sa.x * sb.x - p) + sa.x * sb.y + sa.y * sb.x) + sa.y * sb.y;\n  err += a.x * b.y + a.y * b.x;\n  let hi = p + err;\n  let lo = err - (hi - p);\n  return vec2f(hi, lo);\n}\n\nfn ds_sqr(a: vec2f) -> vec2f {\n  let p = a.x * a.x;\n  let sa = ds_split(a.x);\n  var err = ((sa.x * sa.x - p) + 2.0 * sa.x * sa.y) + sa.y * sa.y;\n  err += 2.0 * a.x * a.y;\n  let hi = p + err;\n  let lo = err - (hi - p);\n  return vec2f(hi, lo);\n}\n\nfn ds_cmp(a: vec2f, b: vec2f) -> f32 {\n  let d = a.x - b.x;\n  if (d != 0.0) { return sign(d); }\n  return sign(a.y - b.y);\n}\n";
+/**
+ * Double-precision fractal preamble (GLSL syntax).
+ *
+ * Emulated double-precision (double-single) Mandelbrot and Julia iterations.
+ * Complex numbers are packed as vec4(re_hi, im_hi, re_lo, im_lo).
+ * Requires GPU_DS_ARITHMETIC_PREAMBLE_GLSL to be included first.
+ */
+export declare const GPU_FRACTAL_DP_PREAMBLE_GLSL = "\nfloat _fractal_mandelbrot_dp(vec4 c, int maxIter) {\n  // c = (re_hi, im_hi, re_lo, im_lo)\n  vec2 cr = vec2(c.x, c.z);  // real part as ds\n  vec2 ci = vec2(c.y, c.w);  // imag part as ds\n  vec2 zr = vec2(0.0, 0.0);\n  vec2 zi = vec2(0.0, 0.0);\n  for (int i = 0; i < maxIter; i++) {\n    vec2 zr2 = ds_sqr(zr);\n    vec2 zi2 = ds_sqr(zi);\n    // |z|^2 > 4.0 ?\n    vec2 mag2 = ds_add(zr2, zi2);\n    if (mag2.x > 4.0)\n      return clamp((float(i) - log2(log2(mag2.x)) + 4.0) / float(maxIter), 0.0, 1.0);\n    // z = z^2 + c\n    vec2 new_zi = ds_add(ds_mul(ds_add(zr, zr), zi), ci); // 2*zr*zi + ci\n    zr = ds_add(ds_sub(zr2, zi2), cr);                    // zr^2 - zi^2 + cr\n    zi = new_zi;\n  }\n  return 1.0;\n}\n\nfloat _fractal_julia_dp(vec4 z_in, vec4 c, int maxIter) {\n  vec2 zr = vec2(z_in.x, z_in.z);\n  vec2 zi = vec2(z_in.y, z_in.w);\n  vec2 cr = vec2(c.x, c.z);\n  vec2 ci = vec2(c.y, c.w);\n  for (int i = 0; i < maxIter; i++) {\n    vec2 zr2 = ds_sqr(zr);\n    vec2 zi2 = ds_sqr(zi);\n    vec2 mag2 = ds_add(zr2, zi2);\n    if (mag2.x > 4.0)\n      return clamp((float(i) - log2(log2(mag2.x)) + 4.0) / float(maxIter), 0.0, 1.0);\n    vec2 new_zi = ds_add(ds_mul(ds_add(zr, zr), zi), ci);\n    zr = ds_add(ds_sub(zr2, zi2), cr);\n    zi = new_zi;\n  }\n  return 1.0;\n}\n";
+/**
+ * Double-precision fractal preamble (WGSL syntax).
+ *
+ * Emulated double-precision (double-single) Mandelbrot and Julia iterations.
+ * Complex numbers are packed as vec4f(re_hi, im_hi, re_lo, im_lo).
+ * Requires GPU_DS_ARITHMETIC_PREAMBLE_WGSL to be included first.
+ */
+export declare const GPU_FRACTAL_DP_PREAMBLE_WGSL = "\nfn _fractal_mandelbrot_dp(c: vec4f, maxIter: i32) -> f32 {\n  let cr = vec2f(c.x, c.z);\n  let ci = vec2f(c.y, c.w);\n  var zr = vec2f(0.0, 0.0);\n  var zi = vec2f(0.0, 0.0);\n  for (var i: i32 = 0; i < maxIter; i++) {\n    let zr2 = ds_sqr(zr);\n    let zi2 = ds_sqr(zi);\n    let mag2 = ds_add(zr2, zi2);\n    if (mag2.x > 4.0) {\n      return clamp((f32(i) - log2(log2(mag2.x)) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n    let new_zi = ds_add(ds_mul(ds_add(zr, zr), zi), ci);\n    zr = ds_add(ds_sub(zr2, zi2), cr);\n    zi = new_zi;\n  }\n  return 1.0;\n}\n\nfn _fractal_julia_dp(z_in: vec4f, c: vec4f, maxIter: i32) -> f32 {\n  var zr = vec2f(z_in.x, z_in.z);\n  var zi = vec2f(z_in.y, z_in.w);\n  let cr = vec2f(c.x, c.z);\n  let ci = vec2f(c.y, c.w);\n  for (var i: i32 = 0; i < maxIter; i++) {\n    let zr2 = ds_sqr(zr);\n    let zi2 = ds_sqr(zi);\n    let mag2 = ds_add(zr2, zi2);\n    if (mag2.x > 4.0) {\n      return clamp((f32(i) - log2(log2(mag2.x)) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n    let new_zi = ds_add(ds_mul(ds_add(zr, zr), zi), ci);\n    zr = ds_add(ds_sub(zr2, zi2), cr);\n    zi = new_zi;\n  }\n  return 1.0;\n}\n";
 /**
  * Fractal preamble (GLSL syntax).
  *
@@ -120,6 +152,20 @@ export declare const GPU_FRACTAL_PREAMBLE_GLSL = "\nfloat _fractal_mandelbrot(ve
  * Fractal preamble (WGSL syntax).
  */
 export declare const GPU_FRACTAL_PREAMBLE_WGSL = "\nfn _fractal_mandelbrot(c: vec2f, maxIter: i32) -> f32 {\n  var z = vec2f(0.0, 0.0);\n  for (var i: i32 = 0; i < maxIter; i++) {\n    z = vec2f(z.x*z.x - z.y*z.y + c.x, 2.0*z.x*z.y + c.y);\n    if (dot(z, z) > 4.0) {\n      return clamp((f32(i) - log2(log2(dot(z, z))) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n  }\n  return 1.0;\n}\n\nfn _fractal_julia(z_in: vec2f, c: vec2f, maxIter: i32) -> f32 {\n  var z = z_in;\n  for (var i: i32 = 0; i < maxIter; i++) {\n    z = vec2f(z.x*z.x - z.y*z.y + c.x, 2.0*z.x*z.y + c.y);\n    if (dot(z, z) > 4.0) {\n      return clamp((f32(i) - log2(log2(dot(z, z))) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n  }\n  return 1.0;\n}\n";
+/**
+ * Perturbation-theory fractal preamble (GLSL syntax).
+ *
+ * Provides `_fractal_mandelbrot_pt` and `_fractal_julia_pt` which iterate
+ * a perturbation delta relative to a precomputed reference orbit stored in
+ * a RG32F texture. Includes glitch detection with single-float fallback.
+ */
+export declare const GPU_FRACTAL_PT_PREAMBLE_GLSL = "\nuniform sampler2D _refOrbit;\nuniform int _refOrbitLen;\nuniform int _refOrbitTexWidth;\n\nvec2 _pt_fetch_orbit(int i) {\n  int y = i / _refOrbitTexWidth;\n  int x = i - y * _refOrbitTexWidth;\n  return texelFetch(_refOrbit, ivec2(x, y), 0).rg;\n}\n\nfloat _fractal_mandelbrot_pt(vec2 delta_c, int maxIter) {\n  float dr = 0.0;\n  float di = 0.0;\n  int orbitLen = min(maxIter, _refOrbitLen);\n  for (int i = 0; i < orbitLen; i++) {\n    vec2 Zn = _pt_fetch_orbit(i);\n    // delta_{n+1} = 2*Z_n*delta_n + delta_n^2 + delta_c\n    float new_dr = 2.0 * (Zn.x * dr - Zn.y * di) + dr * dr - di * di + delta_c.x;\n    float new_di = 2.0 * (Zn.x * di + Zn.y * dr) + 2.0 * dr * di + delta_c.y;\n    dr = new_dr;\n    di = new_di;\n    // Full z = Z_{n+1} + delta for escape check\n    vec2 Zn1 = (i + 1 < orbitLen) ? _pt_fetch_orbit(i + 1) : vec2(0.0);\n    float zr = Zn1.x + dr;\n    float zi = Zn1.y + di;\n    float mag2 = zr * zr + zi * zi;\n    if (mag2 > 4.0)\n      return clamp((float(i) - log2(log2(mag2)) + 4.0) / float(maxIter), 0.0, 1.0);\n    // Glitch detection: |delta|^2 > |Z|^2\n    float dmag2 = dr * dr + di * di;\n    float Zmag2 = Zn.x * Zn.x + Zn.y * Zn.y;\n    if (dmag2 > Zmag2 && Zmag2 > 0.0) {\n      // Rebase to absolute coordinates and continue with single-float\n      float abs_zr = Zn1.x + dr;\n      float abs_zi = Zn1.y + di;\n      // Reconstruct absolute c from reference + delta\n      // (Use ds_from for the concept, but single-float suffices for fallback)\n      float cx = abs_zr - dr + delta_c.x;\n      float cy = abs_zi - di + delta_c.y;\n      for (int j = i + 1; j < maxIter; j++) {\n        float new_zr = abs_zr * abs_zr - abs_zi * abs_zi + cx;\n        abs_zi = 2.0 * abs_zr * abs_zi + cy;\n        abs_zr = new_zr;\n        mag2 = abs_zr * abs_zr + abs_zi * abs_zi;\n        if (mag2 > 4.0)\n          return clamp((float(j) - log2(log2(mag2)) + 4.0) / float(maxIter), 0.0, 1.0);\n      }\n      return 1.0;\n    }\n  }\n  return 1.0;\n}\n\nfloat _fractal_julia_pt(vec2 z_delta, vec2 delta_c, int maxIter) {\n  float dr = z_delta.x;\n  float di = z_delta.y;\n  int orbitLen = min(maxIter, _refOrbitLen);\n  for (int i = 0; i < orbitLen; i++) {\n    vec2 Zn = _pt_fetch_orbit(i);\n    float new_dr = 2.0 * (Zn.x * dr - Zn.y * di) + dr * dr - di * di + delta_c.x;\n    float new_di = 2.0 * (Zn.x * di + Zn.y * dr) + 2.0 * dr * di + delta_c.y;\n    dr = new_dr;\n    di = new_di;\n    vec2 Zn1 = (i + 1 < orbitLen) ? _pt_fetch_orbit(i + 1) : vec2(0.0);\n    float zr = Zn1.x + dr;\n    float zi = Zn1.y + di;\n    float mag2 = zr * zr + zi * zi;\n    if (mag2 > 4.0)\n      return clamp((float(i) - log2(log2(mag2)) + 4.0) / float(maxIter), 0.0, 1.0);\n    float dmag2 = dr * dr + di * di;\n    float Zmag2 = Zn.x * Zn.x + Zn.y * Zn.y;\n    if (dmag2 > Zmag2 && Zmag2 > 0.0) {\n      float abs_zr = Zn1.x + dr;\n      float abs_zi = Zn1.y + di;\n      float cx = delta_c.x;\n      float cy = delta_c.y;\n      for (int j = i + 1; j < maxIter; j++) {\n        float new_zr = abs_zr * abs_zr - abs_zi * abs_zi + cx;\n        abs_zi = 2.0 * abs_zr * abs_zi + cy;\n        abs_zr = new_zr;\n        mag2 = abs_zr * abs_zr + abs_zi * abs_zi;\n        if (mag2 > 4.0)\n          return clamp((float(j) - log2(log2(mag2)) + 4.0) / float(maxIter), 0.0, 1.0);\n      }\n      return 1.0;\n    }\n  }\n  return 1.0;\n}\n";
+/**
+ * Perturbation-theory fractal preamble (WGSL syntax).
+ *
+ * Same logic as the GLSL version, adapted for WGSL syntax.
+ */
+export declare const GPU_FRACTAL_PT_PREAMBLE_WGSL = "\n@group(0) @binding(1) var _refOrbit: texture_2d<f32>;\nvar<uniform> _refOrbitLen: i32;\nvar<uniform> _refOrbitTexWidth: i32;\n\nfn _pt_fetch_orbit(i: i32) -> vec2f {\n  let y = i / _refOrbitTexWidth;\n  let x = i - y * _refOrbitTexWidth;\n  return textureLoad(_refOrbit, vec2i(x, y), 0).rg;\n}\n\nfn _fractal_mandelbrot_pt(delta_c: vec2f, maxIter: i32) -> f32 {\n  var dr: f32 = 0.0;\n  var di: f32 = 0.0;\n  let orbitLen = min(maxIter, _refOrbitLen);\n  for (var i: i32 = 0; i < orbitLen; i++) {\n    let Zn = _pt_fetch_orbit(i);\n    let new_dr = 2.0 * (Zn.x * dr - Zn.y * di) + dr * dr - di * di + delta_c.x;\n    let new_di = 2.0 * (Zn.x * di + Zn.y * dr) + 2.0 * dr * di + delta_c.y;\n    dr = new_dr;\n    di = new_di;\n    var Zn1 = vec2f(0.0);\n    if (i + 1 < orbitLen) { Zn1 = _pt_fetch_orbit(i + 1); }\n    let zr = Zn1.x + dr;\n    let zi = Zn1.y + di;\n    var mag2 = zr * zr + zi * zi;\n    if (mag2 > 4.0) {\n      return clamp((f32(i) - log2(log2(mag2)) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n    let dmag2 = dr * dr + di * di;\n    let Zmag2 = Zn.x * Zn.x + Zn.y * Zn.y;\n    if (dmag2 > Zmag2 && Zmag2 > 0.0) {\n      var f_zr = Zn1.x + dr;\n      var f_zi = Zn1.y + di;\n      let cx = delta_c.x;\n      let cy = delta_c.y;\n      for (var j: i32 = i + 1; j < maxIter; j++) {\n        let t_zr = f_zr * f_zr - f_zi * f_zi + cx;\n        f_zi = 2.0 * f_zr * f_zi + cy;\n        f_zr = t_zr;\n        mag2 = f_zr * f_zr + f_zi * f_zi;\n        if (mag2 > 4.0) {\n          return clamp((f32(j) - log2(log2(mag2)) + 4.0) / f32(maxIter), 0.0, 1.0);\n        }\n      }\n      return 1.0;\n    }\n  }\n  return 1.0;\n}\n\nfn _fractal_julia_pt(z_delta: vec2f, delta_c: vec2f, maxIter: i32) -> f32 {\n  var dr = z_delta.x;\n  var di = z_delta.y;\n  let orbitLen = min(maxIter, _refOrbitLen);\n  for (var i: i32 = 0; i < orbitLen; i++) {\n    let Zn = _pt_fetch_orbit(i);\n    let new_dr = 2.0 * (Zn.x * dr - Zn.y * di) + dr * dr - di * di + delta_c.x;\n    let new_di = 2.0 * (Zn.x * di + Zn.y * dr) + 2.0 * dr * di + delta_c.y;\n    dr = new_dr;\n    di = new_di;\n    var Zn1 = vec2f(0.0);\n    if (i + 1 < orbitLen) { Zn1 = _pt_fetch_orbit(i + 1); }\n    let zr = Zn1.x + dr;\n    let zi = Zn1.y + di;\n    var mag2 = zr * zr + zi * zi;\n    if (mag2 > 4.0) {\n      return clamp((f32(i) - log2(log2(mag2)) + 4.0) / f32(maxIter), 0.0, 1.0);\n    }\n    let dmag2 = dr * dr + di * di;\n    let Zmag2 = Zn.x * Zn.x + Zn.y * Zn.y;\n    if (dmag2 > Zmag2 && Zmag2 > 0.0) {\n      var f_zr = Zn1.x + dr;\n      var f_zi = Zn1.y + di;\n      let cx = delta_c.x;\n      let cy = delta_c.y;\n      for (var j: i32 = i + 1; j < maxIter; j++) {\n        let t_zr = f_zr * f_zr - f_zi * f_zi + cx;\n        f_zi = 2.0 * f_zr * f_zi + cy;\n        f_zr = t_zr;\n        mag2 = f_zr * f_zr + f_zi * f_zi;\n        if (mag2 > 4.0) {\n          return clamp((f32(j) - log2(log2(mag2)) + 4.0) / f32(maxIter), 0.0, 1.0);\n        }\n      }\n      return 1.0;\n    }\n  }\n  return 1.0;\n}\n";
 /**
  * GPU color space conversion preamble (GLSL syntax).
  *

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 /**
  * Source code in the target language
  */
@@ -53,6 +53,8 @@ export interface CompileTarget<Expr = unknown> {
     block?: (statements: string[]) => string;
     /** Target language identifier (for debugging/logging) */
     language?: string;
+    /** Compilation hints (viewport, etc.) passed through from options. */
+    hints?: CompilationOptions['hints'];
 }
 /**
  * Base interface for language-specific compilation targets
@@ -147,6 +149,22 @@ export interface CompilationOptions<Expr = unknown> {
      * results (e.g., plotting).
      */
     realOnly?: boolean;
+    /**
+     * Compilation hints for precision-adaptive strategies.
+     *
+     * The compiler may use these hints to select different code generation
+     * strategies (e.g., emulated double precision for deep-zoom fractals).
+     * Non-fractal functions ignore hints.
+     */
+    hints?: {
+        /** Current viewport for precision-adaptive compilation. */
+        viewport?: {
+            /** Center of the viewport as [re, im]. */
+            center: [number, number];
+            /** Viewport radius (half-width in complex plane units). */
+            radius: number;
+        };
+    };
 }
 /**
  * Built-in targets that produce an executable `run` function.
@@ -247,6 +265,35 @@ export type CompilationResult<T extends string = string, R = number | ComplexRes
     calling?: 'expression' | 'lambda';
     /** Executable function (present for JS-executable targets only). */
     run?: CompiledRunner<R>;
+    /**
+     * Cheap staleness check for precision-adaptive compilation.
+     *
+     * The plot engine checks these thresholds on each viewport change
+     * (a few number comparisons). When any condition is met, the expression
+     * should be recompiled with updated hints.
+     */
+    staleWhen?: {
+        /** Recompile when viewport radius drops below this value. */
+        radiusBelow?: number;
+        /** Recompile when viewport radius rises above this value. */
+        radiusAbove?: number;
+        /** Recompile when center moves more than this distance. */
+        centerDistance?: number;
+    };
+    /** Scalar uniform values the shader needs. */
+    uniforms?: Record<string, number>;
+    /**
+     * Texture data the shader needs (e.g., reference orbit).
+     *
+     * Separated from `uniforms` because the GPU upload path is fundamentally
+     * different (createTexture + sampler vs uniform1f).
+     */
+    textures?: Record<string, {
+        data: Float32Array;
+        width: number;
+        height: number;
+        format: 'r32f' | 'rg32f' | 'rgba32f';
+    }>;
 } & (T extends ExecutableTarget ? {
     calling: 'expression' | 'lambda';
     run: CompiledRunner<R>;

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

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

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 /**
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
     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.55.3 */
+/* 0.55.4 */
 export type CommonSymbolTable = {
     [symbol: string]: null | Expression;
 };

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 type ResetHooks = {
     refreshNumericConstants: () => void;
     resetCommonSymbols: () => void;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 import { BigDecimal } from '../big-decimal';
 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.55.3 */
+/* 0.55.4 */
 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-library-bootstrap.d.ts

@@ -1,3 +1,3 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 import type { AngularUnit } from './types-definitions';
 export declare class EngineNumericConfiguration {
     private _precision;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
     private _timeLimit;
     private _iterationLimit;
     private _recursionLimit;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 /**
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 type ResolverScope = {
     parent: ResolverScope | null;
     types?: Record<string, TypeReference>;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 import type { Type, TypeString } from '../common/type/types';
 import type { LatexString, ParseLatexOptions } 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 import { BigDecimal } from '../big-decimal';
 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.55.3 */
+/* 0.55.4 */
  * Basic interval arithmetic operations
  *
  * @module interval/arithmetic

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * Comparison and conditional interval operations
  *
  * @module interval/comparison

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
  * Trigonometric interval functions
  *
  * @module interval/trigonometric

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

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

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * Utility functions for interval arithmetic
  *
  * @module interval/util

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * Default LaTeX dictionary assembly.
  *
  * This file imports all individual domain dictionaries and assembles them

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_ALGEBRA: LatexDictionary;

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_ARITHMETIC: LatexDictionary;

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_CALCULUS: LatexDictionary;

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_COMPLEX: LatexDictionary;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_LINEAR_ALGEBRA: LatexDictionary;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 import type { MathJsonExpression } from '../../../math-json';
 export declare const DEFINITIONS_LOGIC: LatexDictionary;
 export declare function parseQuantifier(kind: 'NotForAll' | 'NotExists' | 'ForAll' | 'Exists' | 'ExistsUnique'): (parser: Parser, terminator: Readonly<Terminator>) => MathJsonExpression | null;

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_OTHERS: LatexDictionary;

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

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

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_SETS: LatexDictionary;

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

@@ -1,2 +1,2 @@
-/* 0.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_STATISTICS: LatexDictionary;

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

@@ -1,3 +1,3 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 export declare const DEFINITIONS_TRIGONOMETRY: LatexDictionary;

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
  * 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/latex-syntax.d.ts

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * Standalone LaTeX <-> MathJSON parsing and serialization.
  *
  * This module provides a `LatexSyntax` class and free-standing `parse()` /

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
  * 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
  * ## 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 export declare const CALCULUS_LIBRARY: SymbolDefinitions[];

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

@@ -1,4 +1,4 @@
-/* 0.55.3 */
+/* 0.55.4 */
 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.55.3 */
+/* 0.55.4 */
 export declare const COLORS_LIBRARY: SymbolDefinitions;