@jbrowse/render-core 5.0.0-beta.1 → 5.0.0-beta.10

package/src/shaders/sdEllipse.test.ts

@@ -0,0 +1,211 @@
+// `sdEllipse` (curveDistance.slang) measures a fragment's distance to the arc
+// band's and the link mark's half-ellipse. The shader runs it in float32, and
+// float32 is the whole subject here: the same source in float64 is exact to
+// 1e-9 at every shape below, so nothing a float64 test can reach says whether
+// the picture is right.
+//
+// What goes wrong when it is not: the stroke reads the distance through a ramp
+// one output pixel wide, so at dpr 2 an error of half a CSS px is the whole of
+// the ink and the curve breaks into scattered dots. That is what a dome at
+// 380x249 did on WebGL while Canvas2D drew it solid.
+//
+// The shader isn't unit-testable, so this models it — same approach and same
+// reason as glyphEdgeAlpha.test.ts. The reference is a brute-force minimum over
+// the parametric angle, as in marks/ellipseDistance.test.ts: the DEFINITION of
+// the distance rather than a second evaluation of the same formula.
+// SYNC: keep in step with curveDistance.slang's sdEllipse.
+import { edgeCoverage } from './antialias.js.generated.ts'
+
+const f = Math.fround
+const len = (x: number, y: number) => f(Math.sqrt(f(f(x * x) + f(y * y))))
+
+const ELLIPSE_CIRCLE_TOL = 0.05
+
+function nearCircleDistancePx(qx: number, qy: number, ax: number, ay: number) {
+  const k = len(f(qx / ax), f(qy / ay))
+  const g = len(f(qx / f(ax * ax)), f(qy / f(ay * ay)))
+  return Math.abs(f(f(k * f(k - 1)) / Math.max(g, 1e-20)))
+}
+
+function sdEllipse(px: number, py: number, radiusX: number, radiusY: number) {
+  let qx = Math.abs(px)
+  let qy = Math.abs(py)
+  let ax = radiusX
+  let ay = radiusY
+  if (qx > qy) {
+    ;[qx, qy] = [qy, qx]
+    ;[ax, ay] = [ay, ax]
+  }
+  if (f(f(ax * ax) + f(ay * ay)) <= 0) {
+    return len(qx, qy)
+  }
+  const l = f(f(ay * ay) - f(ax * ax))
+  if (Math.abs(l) < f(ELLIPSE_CIRCLE_TOL * f(ay * ay))) {
+    return nearCircleDistancePx(qx, qy, ax, ay)
+  }
+  const m = f(f(ax * qx) / l)
+  const m2 = f(m * m)
+  const n = f(f(ay * qy) / l)
+  const n2 = f(n * n)
+  const c = f(f(f(m2 + n2) - 1) / 3)
+  const c3 = f(f(c * c) * c)
+  const d = f(c3 + f(m2 * n2))
+  const g = f(m + f(m * n2))
+  let co: number
+  if (d < 0) {
+    const h = f(
+      Math.acos(Math.max(-1, Math.min(1, f(f(c3 + f(f(m2 * n2) * 2)) / c3)))) /
+        3,
+    )
+    const s = f(Math.cos(h))
+    const t = f(f(Math.sin(h)) * f(Math.sqrt(3)))
+    const rx = f(Math.sqrt(Math.max(f(f(-c * f(f(s + t) + 2)) + m2), 0)))
+    const ry = f(Math.sqrt(Math.max(f(f(-c * f(f(s - t) + 2)) + m2), 0)))
+    co = f(
+      f(
+        f(
+          f(ry + f(Math.sign(l) * rx)) +
+            f(Math.abs(g) / Math.max(f(rx * ry), 1e-20)),
+        ) - m,
+      ) / 2,
+    )
+  } else {
+    const qq = f(c3 + f(f(m2 * n2) * 2))
+    const h = f(f(f(2 * m) * n) * f(Math.sqrt(d)))
+    const w = f(qq + h)
+    const aw = Math.abs(w)
+    const s = f(Math.sign(w) * f(Math.pow(aw, f(1 / 3))))
+    const t =
+      aw > 0 ? f(Math.sign(w) * f(Math.pow(f(c3 * f(c3 / aw)), f(1 / 3)))) : 0
+    const rx = f(f(f(f(-s - t) - f(c * 4)) + f(2 * m2)))
+    const ry = f(f(s - t) * f(Math.sqrt(3)))
+    const rm = len(rx, ry)
+    co = f(
+      f(
+        f(
+          f(ry / f(Math.sqrt(Math.max(f(rm - rx), 1e-20)))) +
+            f(f(2 * g) / Math.max(rm, 1e-20)),
+        ) - m,
+      ) / 2,
+    )
+  }
+  co = Math.max(-1, Math.min(1, co))
+  const ex = f(ax * co)
+  const ey = f(ay * f(Math.sqrt(Math.max(f(1 - f(co * co)), 0))))
+  return len(f(ex - qx), f(ey - qy))
+}
+
+function bruteForce(px: number, py: number, rx: number, ry: number) {
+  const d = (t: number) =>
+    Math.hypot(rx * Math.cos(t) - Math.abs(px), ry * Math.sin(t) - Math.abs(py))
+  let bestT = 0
+  let best = Number.POSITIVE_INFINITY
+  const N = 1000
+  for (let i = 0; i <= N; i++) {
+    const t = (i / N) * (Math.PI / 2)
+    const v = d(t)
+    if (v < best) {
+      best = v
+      bestT = t
+    }
+  }
+  const cell = Math.PI / 2 / N
+  let lo = bestT - 2 * cell
+  let hi = bestT + 2 * cell
+  for (let i = 0; i < 60; i++) {
+    const m1 = lo + (hi - lo) / 3
+    const m2 = hi - (hi - lo) / 3
+    if (d(m1) < d(m2)) {
+      hi = m2
+    } else {
+      lo = m1
+    }
+  }
+  return d((lo + hi) / 2)
+}
+
+const HALF_PX = 1
+const DPR = 2
+
+/** Worst distance error, and the worst change in ink it causes, over a stroke. */
+function scoreStroke(rx: number, ry: number) {
+  let worst = 0
+  let worstInk = 0
+  for (let k = 0; k <= 40; k++) {
+    const t = (Math.PI * k) / 40
+    const ex = rx * Math.cos(t)
+    const ey = -ry * Math.sin(t)
+    // Along the outward normal, so the samples straddle the stroke the way a
+    // fragment inside the hull does rather than sitting on a grid.
+    const nx = ry * Math.cos(t)
+    const ny = -rx * Math.sin(t)
+    const nl = Math.hypot(nx, ny) || 1
+    for (let i = -7; i <= 7; i++) {
+      const off = (i / 7) * (HALF_PX + 1 / DPR)
+      const px = ex + (nx / nl) * off
+      const py = ey + (ny / nl) * off
+      const ref = bruteForce(px, py, rx, ry)
+      const got = sdEllipse(px, py, rx, ry)
+      worst = Math.max(worst, Math.abs(got - ref))
+      worstInk = Math.max(
+        worstInk,
+        Math.abs(
+          edgeCoverage(HALF_PX - got, DPR) - edgeCoverage(HALF_PX - ref, DPR),
+        ),
+      )
+    }
+  }
+  return { worst, worstInk }
+}
+
+// Every shape here broke before the two float32 guards went in, except the last
+// three, which are the regression side: the flat end the closed form was chosen
+// for, and the exact circle. The figures in the comments are what the COMPILED
+// shader did, measured through slangc's C++ target over the same samples.
+//
+// The rows are not interchangeable — between them they reach all three paths,
+// which is what keeps a change to one from being scored on another. Of each
+// row's 615 samples: the two domes and the 260x250 take the cube-root branch
+// (600, 600, 615), the two bands take the trigonometric one (589, 597) and are
+// the only cover it has, and the four near-circles take the new branch whole.
+test.each([
+  ['the captured speckled dome', 380.25, 249], // was 2.9px out
+  ['a dome at 3:2, 1000px tall', 1500, 1000], // was 10.9px out
+  ['radii a millionth apart', 250.00025, 250], // was 49px out, 70% of the stroke
+  ['radii 1e-5 apart', 250.0025, 250], // was 3.6px out
+  ['radii 1e-4 apart', 250.025, 250],
+  ['near-circle at 2000px', 2000.2, 2000], // was 2.3px out
+  ['just outside the near-circle branch', 260, 250],
+  ['a 44:1 band', 1100, 25],
+  ['a 88:1 band', 2200, 25],
+  ['an exact circle', 250, 250],
+])('float32 holds the stroke: %s', (_name, rx, ry) => {
+  const { worst, worstInk } = scoreStroke(rx, ry)
+  // A tenth of the ramp. The margin is what matters rather than the figure —
+  // the failures this pins were one to two hundred times the ramp's whole
+  // width, not a rounding away from it.
+  expect(worst).toBeLessThan(0.05)
+  expect(worstInk).toBeLessThan(0.1)
+})
+
+// A link pair whose two ends land on the same pixel hands `linkRadiiPx` a zero
+// half-width, and both of its shapes then take the apex to zero as well. The
+// near-circle test cannot catch that pair: its threshold is `ab.y * ab.y`, so
+// at zero radii it reads `abs(l) < 0`, and the solve divides by `l == 0`.
+test.each([
+  [3, 4, 5],
+  [0, 2, 2],
+  [-1.5, 0, 1.5],
+  [0, 0, 0],
+])('radii (0, 0) measure to the origin: (%p, %p)', (px, py, want) => {
+  const got = sdEllipse(px, py, 0, 0)
+  expect(Number.isNaN(got)).toBe(false)
+  expect(got).toBeCloseTo(want, 6)
+})
+
+test('the near-circle branch answers a true circle exactly', () => {
+  // With equal radii the implicit form reduces to |q| - r algebraically, which
+  // is why it can replace the dedicated circle branch it grew out of.
+  expect(nearCircleDistancePx(3, 4, 1, 1)).toBeCloseTo(4, 6)
+  expect(nearCircleDistancePx(30, 40, 100, 100)).toBeCloseTo(50, 4)
+})

package/src/shaders/spanMark.generated.ts

@@ -0,0 +1,11 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from spanMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+import type { ShaderSource } from '@jbrowse/render-core/hal'
+
+export * from './spanMark.iface.generated.ts'
+
+export const SOURCE: ShaderSource = {
+  wgsl: () => import('./spanMark.wgsl.generated.ts'),
+  glsl: () => import('./spanMark.glsl.generated.ts'),
+}

package/src/shaders/spanMark.glsl.generated.ts

@@ -0,0 +1,6 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from spanMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const GLSL_VERTEX = "#version 300 es\nprecision highp float;\nprecision highp int;\nprecision highp sampler2D;\n\n\n\nlayout(std140) uniform Uniforms\n{\n    vec3 bpRangeX_0;\n    float canvasHeight_0;\n    float minCellDenomPx_0;\n    float minCellPx_0;\n    float zero_0;\n    float rowHeight_0;\n    float rowProportion_0;\n    float scrollTop_0;\n    int rowTableKeys_0;\n    int rampMode_0;\n    float rampMin_0;\n    float rampMax_0;\n    float rampMidNorm_0;\n    int colorCutCount_0;\n    vec4  colorCuts_0[2];\n    vec4  colorBands_0[9];\n}u_0;\n\nuniform sampler2D u_colorRamp;\n\n\nuniform sampler2D u_rowTable;\n\n\nvec4 unpackRGBA_0(uint c_0)\n{\n\n    return vec4(float((c_0 >> 0U) & 255U), float((c_0 >> 8U) & 255U), float((c_0 >> 16U) & 255U), float((c_0 >> 24U) & 255U)) / 255.0;\n}\n\n\nfloat rampMidT_0(float t_0, float mid_0)\n{\n\n    float _S1 = max(mid_0, 1.0 - mid_0);\n    return clamp((t_0 + (_S1 - mid_0)) / (2.0 * _S1), 0.0, 1.0);\n}\n\n\nfloat colorCutAt_0(int i_0, vec4 cutsLo_0, vec4 cutsHi_0)\n{\n\n    float _S2;\n    if(i_0 < 4)\n    {\n\n        _S2 = cutsLo_0[i_0];\n\n    }\n    else\n    {\n\n        _S2 = cutsHi_0[i_0 - 4];\n\n    }\n\n    return _S2;\n}\n\n\nint thresholdBand_0(float value_0, int cutCount_0, vec4 cutsLo_1, vec4 cutsHi_1)\n{\n\n    int band_0 = 0;\n\n    int i_1 = 0;\n\n    for(;;)\n    {\n\n        if(i_1 < cutCount_0)\n        {\n        }\n        else\n        {\n\n            break;\n        }\n\n        if(value_0 >= (colorCutAt_0(i_1, cutsLo_1, cutsHi_1)))\n        {\n\n            band_0 = i_1 + 1;\n\n        }\n\n        i_1 = i_1 + 1;\n\n    }\n\n    return band_0;\n}\n\n\nfloat log1pf_0(float x_0)\n{\n\n    float u_1 = 1.0 + x_0;\n\n    float _S3;\n    if(u_1 == 1.0)\n    {\n\n        _S3 = x_0;\n\n    }\n    else\n    {\n\n        _S3 = log(u_1) * (x_0 / (u_1 - 1.0));\n\n    }\n\n    return _S3;\n}\n\nfloat symlogTransform_0(float x_1, float c_1)\n{\n\n    return float((int(sign((x_1))))) * log1pf_0(abs(x_1 / c_1));\n}\n\n\nfloat normalizeScoreUnclamped_0(float score_0, float domainMin_0, float domainMax_0, int scaleType_0, float symlogConstant_0)\n{\n\n    float tMax_0;\n\n    float tMin_0;\n\n    float t_1;\n\n\n\n    if(scaleType_0 == 2)\n    {\n\n        float _S4 = symlogTransform_0(score_0, symlogConstant_0);\n        float _S5 = symlogTransform_0(domainMin_0, symlogConstant_0);\n\n        tMax_0 = symlogTransform_0(domainMax_0, symlogConstant_0);\n\n        tMin_0 = _S5;\n\n        t_1 = _S4;\n\n    }\n    else\n    {\n\n        if(scaleType_0 == 1)\n        {\n            if(domainMin_0 > 0.0)\n            {\n\n                tMax_0 = domainMin_0;\n\n            }\n            else\n            {\n\n                tMax_0 = 1.0;\n\n            }\n            float _S6 = log2(max(score_0, tMax_0));\n            float _S7 = log2(tMax_0);\n\n            tMax_0 = log2(max(domainMax_0, tMax_0));\n\n            tMin_0 = _S7;\n\n            t_1 = _S6;\n\n        }\n        else\n        {\n\n            tMax_0 = domainMax_0;\n\n            tMin_0 = domainMin_0;\n\n            t_1 = score_0;\n\n        }\n\n    }\n\n    float tRange_0 = tMax_0 - tMin_0;\n    if(tRange_0 <= 0.0)\n    {\n\n        if(t_1 > tMin_0)\n        {\n\n            tMax_0 = 1.0;\n\n        }\n        else\n        {\n\n            tMax_0 = 0.0;\n\n        }\n\n        return tMax_0;\n    }\n    return (t_1 - tMin_0) / tRange_0;\n}\n\nfloat normalizeScore_0(float score_1, float domainMin_1, float domainMax_1, int scaleType_1, float symlogConstant_1)\n{\n\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1), 0.0, 1.0);\n}\n\n\nuint rowTableWidth_0(uint keys_0)\n{\n\n    return max(min(keys_0, 2048U), 1U);\n}\n\nuint rowTablePlaneHeight_0(uint keys_1)\n{\n\n    uint w_0 = rowTableWidth_0(keys_1);\n    uint _S8 = (keys_1 + w_0 - 1U) / w_0;\n\n    return max(_S8, 1U);\n}\n\nuint rowTableTexelX_0(uint key_0, uint keys_2)\n{\n\n    uint _S9 = key_0 % rowTableWidth_0(keys_2);\n\n    return _S9;\n}\n\n\nuint rowTableTexelY_0(uint key_1, uint keys_3)\n{\n\n    uint _S10 = key_1 / rowTableWidth_0(keys_3);\n\n    return _S10;\n}\n\n\nuint byteOf_0(float channel_0)\n{\n\n    return uint(channel_0 * 255.0 + 0.5);\n}\n\n\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S11 = v_0 == 0U;\n\n    bool _S12;\n\n    if(_S11)\n    {\n\n        _S12 = true;\n\n    }\n    else\n    {\n\n        _S12 = v_0 == 2U;\n\n    }\n\n    if(_S12)\n    {\n\n        _S12 = true;\n\n    }\n    else\n    {\n\n        _S12 = v_0 == 3U;\n\n    }\n\n    float _S13;\n\n    if(_S12)\n    {\n\n        _S13 = 0.0;\n\n    }\n    else\n    {\n\n        _S13 = 1.0;\n\n    }\n    if(_S11)\n    {\n\n        _S12 = true;\n\n    }\n    else\n    {\n\n        _S12 = v_0 == 1U;\n\n    }\n\n    if(_S12)\n    {\n\n        _S12 = true;\n\n    }\n    else\n    {\n\n        _S12 = v_0 == 4U;\n\n    }\n\n    float _S14;\n\n    if(_S12)\n    {\n\n        _S14 = 0.0;\n\n    }\n    else\n    {\n\n        _S14 = 1.0;\n\n    }\n\n    return vec2(_S13, _S14);\n}\n\n\nfloat drawnRowHeightPx_0(float rowHeight_1, float rowProportion_1)\n{\n\n    return max(rowHeight_1 * rowProportion_1, 1.0);\n}\n\nfloat rowBandOffsetPx_0(float rowHeight_2, float rowProportion_2)\n{\n\n    return (rowHeight_2 - drawnRowHeightPx_0(rowHeight_2, rowProportion_2)) * 0.5;\n}\n\n\nstruct RowBand_0\n{\n    float height_0;\n    float offset_0;\n};\n\n\nRowBand_0 rowBandPx_0(float rowHeight_3, float rowProportion_3)\n{\n\n    RowBand_0 band_1;\n    band_1.height_0 = drawnRowHeightPx_0(rowHeight_3, rowProportion_3);\n    band_1.offset_0 = rowBandOffsetPx_0(rowHeight_3, rowProportion_3);\n    return band_1;\n}\n\n\nvec2 hpSplitUint_0(uint value_1)\n{\n\n    uint lo_0 = value_1 & 4095U;\n\n    return vec2(float(value_1 - lo_0), float(lo_0));\n}\n\n\nfloat hpToClipX_0(vec2 splitPos_0, vec3 bpRange_0, float hpZero_0)\n{\n\n    float floorBp_0 = -3.00000000549775576e+38 + hpZero_0;\n    float step_0 = 2.0 / bpRange_0.z;\n\n\n    return dot(vec3(-1.0, max(splitPos_0.x - bpRange_0.x, floorBp_0), max(splitPos_0.y - bpRange_0.y, floorBp_0)), vec3(1.0, step_0, step_0));\n}\n\n\nfloat pxToClipLen_0(float px_0, float size_0)\n{\n\n    return px_0 * 2.0 / size_0;\n}\n\n\nfloat extendToMinWidthPx_0(float x1_0, float x2_0, float minWidth_0)\n{\n\n    float dx_0 = x2_0 - x1_0;\n\n    float _S15;\n    if((abs(dx_0)) < minWidth_0)\n    {\n\n        if(dx_0 < 0.0)\n        {\n\n            _S15 = - minWidth_0;\n\n        }\n        else\n        {\n\n            _S15 = minWidth_0;\n\n        }\n\n        _S15 = x1_0 + _S15;\n\n    }\n    else\n    {\n\n        _S15 = x2_0;\n\n    }\n\n    return _S15;\n}\n\nfloat extendToMinWidthX_0(float x1_1, float x2_1, float minWidthPx_0, float canvasWidth_0)\n{\n\n    return extendToMinWidthPx_0(x1_1, x2_1, pxToClipLen_0(minWidthPx_0, canvasWidth_0));\n}\n\n\nfloat yPxToClipY_0(float yPx_0, float canvasHeight_1)\n{\n\n    return 1.0 - yPx_0 / canvasHeight_1 * 2.0;\n}\n\n\nstruct RowRectUniforms_0\n{\n    vec3 bpRangeX_1;\n    float canvasHeight_2;\n    float minCellDenomPx_1;\n    float minCellPx_1;\n    float zero_1;\n    float rowHeight_4;\n    float rowProportion_4;\n    float scrollTop_1;\n};\n\n\nvec4 rowRectClipPos_0(uint startBp_0, uint endBp_0, float rowIndex_0, uint vid_1, RowRectUniforms_0 u_2)\n{\n\n    vec2 local_0 = quadLocal_0(vid_1);\n\n    RowBand_0 band_2 = rowBandPx_0(u_2.rowHeight_4, u_2.rowProportion_4);\n    float rowY_0 = band_2.offset_0 + u_2.rowHeight_4 * rowIndex_0 - u_2.scrollTop_1;\n\n    float sx1_0 = hpToClipX_0(hpSplitUint_0(startBp_0), u_2.bpRangeX_1, u_2.zero_1);\n\n    return vec4(mix(sx1_0, extendToMinWidthX_0(sx1_0, hpToClipX_0(hpSplitUint_0(endBp_0), u_2.bpRangeX_1, u_2.zero_1), u_2.minCellPx_1, u_2.minCellDenomPx_1), local_0.x), yPxToClipY_0(mix(rowY_0, rowY_0 + band_2.height_0, local_0.y), u_2.canvasHeight_2), 0.0, 1.0);\n}\n\n\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_x;\n\n\nlayout(location = 1)\nin uint a_x2;\n\n\nlayout(location = 2)\nin uint a_row;\n\n\nlayout(location = 3)\nin uint a_color;\n\n\nvec4 rampColor_0(float _S16)\n{\n\n    return (textureLod((u_colorRamp), (vec2((clamp(_S16, 0.0, 1.0) * 255.0 + 0.5) / 256.0, 0.5)), (0.0)));\n}\n\n\nvec4 markScaleColor_0(uint _S17, int _S18, float _S19, float _S20, float _S21, int _S22, vec4 _S23, vec4 _S24)\n{\n\n    if(_S18 == 0)\n    {\n\n        return unpackRGBA_0(_S17);\n    }\n    if(_S17 == 2143289345U)\n    {\n\n        return unpackRGBA_0(4289703855U);\n    }\n\n    int scaleType_2;\n\n    if((_S17 & 2139095040U) == 2139095040U)\n    {\n\n        if((_S17 & 8388607U) != 0U)\n        {\n\n            return unpackRGBA_0(4286611584U);\n        }\n\n        bool negative_0 = (_S17 & 2147483648U) != 0U;\n        if(_S18 == 3)\n        {\n\n            if(negative_0)\n            {\n\n                scaleType_2 = 0;\n\n            }\n            else\n            {\n\n                scaleType_2 = _S22;\n\n            }\n\n            return u_0.colorBands_0[scaleType_2];\n        }\n\n        float _S25;\n\n        if(negative_0)\n        {\n\n            _S25 = 0.0;\n\n        }\n        else\n        {\n\n            _S25 = 1.0;\n\n        }\n\n        return rampColor_0(rampMidT_0(_S25, _S21));\n    }\n    float value_2 = uintBitsToFloat(_S17);\n    if(_S18 == 3)\n    {\n\n        return u_0.colorBands_0[thresholdBand_0(value_2, _S22, _S23, _S24)];\n    }\n    if(_S18 == 2)\n    {\n\n        scaleType_2 = 1;\n\n    }\n    else\n    {\n\n        scaleType_2 = 0;\n\n    }\n    return rampColor_0(rampMidT_0(normalizeScore_0(value_2, _S19, _S20, scaleType_2, 1.0), _S21));\n}\n\n\nstruct RowPlacement_0\n{\n    bool hidden_0;\n    uint slot_0;\n    vec4 color_0;\n};\n\n\nRowPlacement_0 rowTableLookup_0(uint _S26, int _S27, vec4 _S28)\n{\n\n    RowPlacement_0 p_0;\n    p_0.hidden_0 = false;\n    p_0.slot_0 = _S26;\n    p_0.color_0 = _S28;\n    if(_S27 < 0)\n    {\n\n        return p_0;\n    }\n    uint n_0 = uint(_S27);\n    if(_S26 >= n_0)\n    {\n\n        p_0.hidden_0 = true;\n        return p_0;\n    }\n    uint h_0 = rowTablePlaneHeight_0(n_0);\n    uint _S29 = rowTableTexelX_0(_S26, n_0);\n\n    float u_3 = (float(_S29) + 0.5) / float(rowTableWidth_0(n_0));\n    float rows_0 = float(2U * h_0);\n    uint y_0 = rowTableTexelY_0(_S26, n_0);\n    vec4 s_0 = (textureLod((u_rowTable), (vec2(u_3, (float(y_0) + 0.5) / rows_0)), (0.0)));\n    if((s_0.w) < 0.5)\n    {\n\n        p_0.hidden_0 = true;\n        return p_0;\n    }\n    p_0.slot_0 = ((byteOf_0(s_0.x)) | ((byteOf_0(s_0.y)) << 8)) | ((byteOf_0(s_0.z)) << 16);\n    vec4 c_2 = (textureLod((u_rowTable), (vec2(u_3, (float(y_0 + h_0) + 0.5) / rows_0)), (0.0)));\n\n    vec4 _S30;\n    if((c_2.w) > 0.0)\n    {\n\n        _S30 = c_2;\n\n    }\n    else\n    {\n\n        _S30 = _S28;\n\n    }\n\n    p_0.color_0 = _S30;\n    return p_0;\n}\n\n\nstruct RowRectVsOut_0\n{\n    vec4 position_0;\n    vec4 color_1;\n};\n\n\nvoid main()\n{\n\n    RowPlacement_0 _S31 = rowTableLookup_0(a_row, u_0.rowTableKeys_0, markScaleColor_0(a_color, u_0.rampMode_0, u_0.rampMin_0, u_0.rampMax_0, u_0.rampMidNorm_0, u_0.colorCutCount_0, u_0.colorCuts_0[0], u_0.colorCuts_0[1]));\n\n\n    RowRectVsOut_0 o_0;\n    if(_S31.hidden_0)\n    {\n\n        o_0.position_0 = vec4(2.0, 2.0, 0.0, 1.0);\n        o_0.color_1 = vec4(0.0);\n        RowRectVsOut_0 _S32 = o_0;\n\n        gl_Position = o_0.position_0;\n\n        v_color = _S32.color_1;\n\n        return;\n    }\n    RowRectUniforms_0 r_0;\n    r_0.bpRangeX_1 = u_0.bpRangeX_0;\n    r_0.canvasHeight_2 = u_0.canvasHeight_0;\n    r_0.minCellDenomPx_1 = u_0.minCellDenomPx_0;\n    r_0.minCellPx_1 = u_0.minCellPx_0;\n    r_0.zero_1 = u_0.zero_0;\n    r_0.rowHeight_4 = u_0.rowHeight_0;\n    r_0.rowProportion_4 = u_0.rowProportion_0;\n    r_0.scrollTop_1 = u_0.scrollTop_0;\n    o_0.position_0 = rowRectClipPos_0(a_x, a_x2, float(_S31.slot_0), uint(gl_VertexID), r_0);\n    o_0.color_1 = _S31.color_0;\n    RowRectVsOut_0 _S33 = o_0;\n\n    gl_Position = o_0.position_0;\n\n    v_color = _S33.color_1;\n\n    return;\n}\n\n"
+
+export const GLSL_FRAGMENT = "#version 300 es\nprecision highp float;\nprecision highp int;\nprecision highp sampler2D;\nlayout(location = 0)\nout vec4 entryPointParam_fs_main_0;\n\n\nin vec4 v_color;\n\n\nvoid main()\n{\n\n    entryPointParam_fs_main_0 = v_color;\n\n    return;\n}\n\n"

package/src/shaders/spanMark.iface.generated.ts

@@ -0,0 +1,289 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from spanMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+import type { VertexAttributeLayout, ShaderBinding, TextureBinding } from '@jbrowse/render-core/hal'
+
+export const BINDINGS: readonly ShaderBinding[] = [
+  { index: 1, kind: 'uniform', name: 'u', stages: ['vertex'] },
+  { index: 2, kind: 'texture', name: 'colorRamp', stages: ['vertex'] },
+  { index: 3, kind: 'sampler', name: 'colorRamp', stages: ['vertex'] },
+  { index: 4, kind: 'texture', name: 'rowTable', stages: ['vertex'] },
+  { index: 5, kind: 'sampler', name: 'rowTable', stages: ['vertex'] },
+]
+
+export const VERTS_PER_INSTANCE = 6
+
+export const UNIFORMS_SIZE_BYTES = 240
+
+// Word indices into a Float32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_F32 = {
+  bpRangeX: 0,
+  canvasHeight: 3,
+  minCellDenomPx: 4,
+  minCellPx: 5,
+  zero: 6,
+  rowHeight: 7,
+  rowProportion: 8,
+  scrollTop: 9,
+  rampMin: 12,
+  rampMax: 13,
+  rampMidNorm: 14,
+} as const
+
+// Word indices into a Int32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_I32 = {
+  rowTableKeys: 10,
+  rampMode: 11,
+  colorCutCount: 15,
+} as const
+
+
+// Word indices of each array field’s elements, into a 4-byte-word
+// view over the uniform buffer (Uint32Array or Float32Array — the
+// field’s scalar type picks, same as UNIFORM_OFFSET_*). NOT
+// consecutive: std140 pads every array element to 16 bytes.
+export const UNIFORM_SLOT_ARRAYS = {
+  colorCuts: [16, 20] as const,
+  colorBands: [24, 28, 32, 36, 40, 44, 48, 52, 56] as const,
+} as const
+
+// Element `i` of the `colorCuts` uniform array (4 components).
+export function setUniformColorCuts(
+  f32: Float32Array,
+  i: number,
+  v0: number,
+  v1: number,
+  v2: number,
+  v3: number,
+) {
+  const o = UNIFORM_SLOT_ARRAYS.colorCuts[i]!
+  f32[o] = v0
+  f32[o + 1] = v1
+  f32[o + 2] = v2
+  f32[o + 3] = v3
+}
+
+// Element `i` of the `colorBands` uniform array (4 components).
+export function setUniformColorBands(
+  f32: Float32Array,
+  i: number,
+  v0: number,
+  v1: number,
+  v2: number,
+  v3: number,
+) {
+  const o = UNIFORM_SLOT_ARRAYS.colorBands[i]!
+  f32[o] = v0
+  f32[o + 1] = v1
+  f32[o + 2] = v2
+  f32[o + 3] = v3
+}
+
+
+export interface Uniforms {
+  bpRangeX: [number, number, number]
+  canvasHeight: number
+  minCellDenomPx: number
+  minCellPx: number
+  zero: number
+  rowHeight: number
+  rowProportion: number
+  scrollTop: number
+  rowTableKeys: number
+  rampMode: number
+  rampMin: number
+  rampMax: number
+  rampMidNorm: number
+  colorCutCount: number
+  colorCuts: [[number, number, number, number], [number, number, number, number]]
+  colorBands: [[number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number], [number, number, number, number]]
+}
+
+export function writeUniforms(buf: ArrayBuffer, uniforms: Uniforms) {
+  const f32 = new Float32Array(buf)
+  const i32 = new Int32Array(buf)
+  f32[0] = uniforms.bpRangeX[0]
+  f32[1] = uniforms.bpRangeX[1]
+  f32[2] = uniforms.bpRangeX[2]
+  f32[3] = uniforms.canvasHeight
+  f32[4] = uniforms.minCellDenomPx
+  f32[5] = uniforms.minCellPx
+  f32[6] = uniforms.zero
+  f32[7] = uniforms.rowHeight
+  f32[8] = uniforms.rowProportion
+  f32[9] = uniforms.scrollTop
+  i32[10] = uniforms.rowTableKeys
+  i32[11] = uniforms.rampMode
+  f32[12] = uniforms.rampMin
+  f32[13] = uniforms.rampMax
+  f32[14] = uniforms.rampMidNorm
+  i32[15] = uniforms.colorCutCount
+  f32[16] = uniforms.colorCuts[0][0]
+  f32[17] = uniforms.colorCuts[0][1]
+  f32[18] = uniforms.colorCuts[0][2]
+  f32[19] = uniforms.colorCuts[0][3]
+  f32[20] = uniforms.colorCuts[1][0]
+  f32[21] = uniforms.colorCuts[1][1]
+  f32[22] = uniforms.colorCuts[1][2]
+  f32[23] = uniforms.colorCuts[1][3]
+  f32[24] = uniforms.colorBands[0][0]
+  f32[25] = uniforms.colorBands[0][1]
+  f32[26] = uniforms.colorBands[0][2]
+  f32[27] = uniforms.colorBands[0][3]
+  f32[28] = uniforms.colorBands[1][0]
+  f32[29] = uniforms.colorBands[1][1]
+  f32[30] = uniforms.colorBands[1][2]
+  f32[31] = uniforms.colorBands[1][3]
+  f32[32] = uniforms.colorBands[2][0]
+  f32[33] = uniforms.colorBands[2][1]
+  f32[34] = uniforms.colorBands[2][2]
+  f32[35] = uniforms.colorBands[2][3]
+  f32[36] = uniforms.colorBands[3][0]
+  f32[37] = uniforms.colorBands[3][1]
+  f32[38] = uniforms.colorBands[3][2]
+  f32[39] = uniforms.colorBands[3][3]
+  f32[40] = uniforms.colorBands[4][0]
+  f32[41] = uniforms.colorBands[4][1]
+  f32[42] = uniforms.colorBands[4][2]
+  f32[43] = uniforms.colorBands[4][3]
+  f32[44] = uniforms.colorBands[5][0]
+  f32[45] = uniforms.colorBands[5][1]
+  f32[46] = uniforms.colorBands[5][2]
+  f32[47] = uniforms.colorBands[5][3]
+  f32[48] = uniforms.colorBands[6][0]
+  f32[49] = uniforms.colorBands[6][1]
+  f32[50] = uniforms.colorBands[6][2]
+  f32[51] = uniforms.colorBands[6][3]
+  f32[52] = uniforms.colorBands[7][0]
+  f32[53] = uniforms.colorBands[7][1]
+  f32[54] = uniforms.colorBands[7][2]
+  f32[55] = uniforms.colorBands[7][3]
+  f32[56] = uniforms.colorBands[8][0]
+  f32[57] = uniforms.colorBands[8][1]
+  f32[58] = uniforms.colorBands[8][2]
+  f32[59] = uniforms.colorBands[8][3]
+}
+
+export const INSTANCE_STRIDE_BYTES = 16
+export const INSTANCE_STRIDE_WORDS = 4
+
+// Word indices into a Uint32Array view over the instance buffer.
+export const INSTANCE_OFFSET_U32 = {
+  x: 0,
+  x2: 1,
+  row: 2,
+  color: 3,
+} as const
+
+export const VERTEX_ATTRIBUTES: readonly VertexAttributeLayout[] = [
+  { name: 'a_x', components: 1, type: 'uint', offsetBytes: 0, integer: true },
+  { name: 'a_x2', components: 1, type: 'uint', offsetBytes: 4, integer: true },
+  { name: 'a_row', components: 1, type: 'uint', offsetBytes: 8, integer: true },
+  { name: 'a_color', components: 1, type: 'uint', offsetBytes: 12, integer: true },
+]
+
+export interface InstanceArrays {
+  x: ArrayLike<number>
+  x2: ArrayLike<number>
+  row: ArrayLike<number>
+  color: ArrayLike<number>
+}
+
+export function packInstances(
+  arrays: InstanceArrays,
+  numInstances: number,
+  buf: ArrayBuffer = new ArrayBuffer(numInstances * INSTANCE_STRIDE_BYTES),
+) {
+  const u32 = new Uint32Array(buf)
+  const { x, x2, row, color } = arrays
+  for (let i = 0; i < numInstances; i++) {
+    const o = i * INSTANCE_STRIDE_WORDS
+    u32[o + 0] = x[i]!
+    u32[o + 1] = x2[i]!
+    u32[o + 2] = row[i]!
+    u32[o + 3] = color[i]!
+  }
+  return buf
+}
+
+// Instance `i`'s `x`.
+export function getInstanceX(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS]!
+}
+
+export function setInstanceX(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS] = v
+}
+
+// Instance `i`'s `x2`.
+export function getInstanceX2(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS + 1]!
+}
+
+export function setInstanceX2(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS + 1] = v
+}
+
+// Instance `i`'s `row`.
+export function getInstanceRow(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS + 2]!
+}
+
+export function setInstanceRow(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS + 2] = v
+}
+
+// Instance `i`'s `color`.
+export function getInstanceColor(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS + 3]!
+}
+
+export function setInstanceColor(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS + 3] = v
+}
+
+// Appends instances one at a time, for an encoder that cannot say up front
+// how many it will emit. Seed the constructor with an upper bound and the
+// common path allocates exactly once; the doubling is a correctness
+// backstop, not the expected route.
+export class InstanceWriter {
+  private buf: ArrayBuffer
+  private u32: Uint32Array
+  private capacity: number
+  count = 0
+
+  constructor(capacity: number) {
+    this.capacity = Math.max(1, capacity)
+    this.buf = new ArrayBuffer(this.capacity * INSTANCE_STRIDE_BYTES)
+    this.u32 = new Uint32Array(this.buf)
+  }
+
+  push(x: number, x2: number, row: number, color: number) {
+    if (this.count === this.capacity) {
+      this.capacity *= 2
+      const grown = new ArrayBuffer(this.capacity * INSTANCE_STRIDE_BYTES)
+      new Uint8Array(grown).set(new Uint8Array(this.buf))
+      this.buf = grown
+      this.u32 = new Uint32Array(grown)
+    }
+    const o = this.count * INSTANCE_STRIDE_WORDS
+    this.u32[o] = x
+    this.u32[o + 1] = x2
+    this.u32[o + 2] = row
+    this.u32[o + 3] = color
+    this.count++
+  }
+
+  // A right-sized COPY, not a subarray view — see the class comment. Skipped
+  // entirely when the seed turned out to be exact, which is the common path.
+  finish() {
+    const used = this.count * INSTANCE_STRIDE_BYTES
+    return used === this.buf.byteLength ? this.buf : this.buf.slice(0, used)
+  }
+}
+
+// Combined `Sampler2D` bindings. Texture unit indices start at 0.
+export const TEXTURES: readonly [TextureBinding, ...TextureBinding[]] = [
+  { name: 'colorRamp', glTextureUnit: 0, glUniformName: 'u_colorRamp', filter: 'linear' },
+  { name: 'rowTable', glTextureUnit: 1, glUniformName: 'u_rowTable', filter: 'nearest' },
+]

package/src/shaders/spanMark.slang

@@ -0,0 +1,77 @@
+// The `span` mark shape: one instance is a coloured rectangle spanning
+// [x, x2) in absolute genomic coords, on the row band given by the `row`
+// channel. The geometry — row band, hp clip-x span, min-width widening — is
+// `rowRect`'s; this is the entry point that makes it a pass a mark can draw.
+//
+// `row` is read through the bound row table (`rowTable`): a key resolved to
+// its drawn slot and colour, or collapsed to a zero-area quad off clip space
+// where the table hides it. A pass binding no table (`rowTableKeys` −1) draws
+// `row` as the slot, which is what every consumer did before the table.
+//
+// The colour slot is `markColor`'s: a packed ABGR, or under a ramp or a
+// threshold the value's bits resolved here through the uniforms, which is
+// what makes a span coloured by a field a heatmap whose domain moves as a
+// uniform.
+//! targets: wgsl, glsl
+//! instance-writer
+
+import colorPack;
+import markColor;
+import rowRect;
+import rowTable;
+
+public static const uint VERTS_PER_INSTANCE = 6u;
+
+struct SpanMarkUniforms {
+  float3 bpRangeX;
+  float  canvasHeight;
+  float  minCellDenomPx;
+  float  minCellPx;
+  float  zero;
+  float  rowHeight;
+  float  rowProportion;
+  float  scrollTop;
+  int    rowTableKeys;
+  // `markColor`'s: 0 packs the instance's colour, 1 and 2 read its slot as a
+  // value through the ramp, 3 through the threshold below.
+  int    rampMode;
+  float  rampMin;
+  float  rampMax;
+  float  rampMidNorm;
+  int    colorCutCount;
+  float4 colorCuts[2];
+  float4 colorBands[MAX_COLOR_BANDS];
+};
+
+[[vk::binding(1, 0)]] ConstantBuffer<SpanMarkUniforms> u;
+[[vk::binding(2, 0)]] Sampler2D<float4> colorRamp;
+[[vk::binding(4, 0)]] Sampler2D<float4> rowTable;
+
+[shader("vertex")]
+RowRectVsOut vs_main(RowRectInstance inst, uint vid : SV_VertexID) {
+  RowPlacement p = rowTableLookup(rowTable, inst.row, u.rowTableKeys,
+    markScaleColor(colorRamp, inst.color, u.rampMode, u.rampMin, u.rampMax, u.rampMidNorm, u.colorCutCount, u.colorCuts[0], u.colorCuts[1], u.colorBands));
+  RowRectVsOut o;
+  if (p.hidden) {
+    o.position = float4(2.0, 2.0, 0.0, 1.0);
+    o.color = float4(0.0);
+    return o;
+  }
+  RowRectUniforms r;
+  r.bpRangeX       = u.bpRangeX;
+  r.canvasHeight   = u.canvasHeight;
+  r.minCellDenomPx = u.minCellDenomPx;
+  r.minCellPx      = u.minCellPx;
+  r.zero           = u.zero;
+  r.rowHeight      = u.rowHeight;
+  r.rowProportion  = u.rowProportion;
+  r.scrollTop      = u.scrollTop;
+  o.position = rowRectClipPos(inst.x, inst.x2, float(p.slot), vid, r);
+  o.color = p.color;
+  return o;
+}
+
+[shader("fragment")]
+float4 fs_main(RowRectVsOut v) : SV_Target {
+  return v.color;
+}

package/src/shaders/spanMark.wgsl.generated.ts

@@ -0,0 +1,4 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from spanMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const WGSL_SOURCE = "struct SpanMarkUniforms_std140_0\n{\n    @align(16) bpRangeX_0 : vec3<f32>,\n    @align(4) canvasHeight_0 : f32,\n    @align(16) minCellDenomPx_0 : f32,\n    @align(4) minCellPx_0 : f32,\n    @align(8) zero_0 : f32,\n    @align(4) rowHeight_0 : f32,\n    @align(16) rowProportion_0 : f32,\n    @align(4) scrollTop_0 : f32,\n    @align(8) rowTableKeys_0 : i32,\n    @align(4) rampMode_0 : i32,\n    @align(16) rampMin_0 : f32,\n    @align(4) rampMax_0 : f32,\n    @align(8) rampMidNorm_0 : f32,\n    @align(4) colorCutCount_0 : i32,\n    @align(16) colorCuts_0 : array<vec4<f32>, i32(2)>,\n    @align(16) colorBands_0 : array<vec4<f32>, i32(9)>,\n};\n\n@binding(1) @group(0) var<uniform> u_0 : SpanMarkUniforms_std140_0;\n@binding(2) @group(0) var colorRamp_texture_0 : texture_2d<f32>;\n\n@binding(3) @group(0) var colorRamp_sampler_0 : sampler;\n\n@binding(4) @group(0) var rowTable_texture_0 : texture_2d<f32>;\n\n@binding(5) @group(0) var rowTable_sampler_0 : sampler;\n\nfn unpackRGBA_0( c_0 : u32) -> vec4<f32>\n{\n    return vec4<f32>(f32((((c_0 >> (u32(0)))) & (u32(255)))), f32((((c_0 >> (u32(8)))) & (u32(255)))), f32((((c_0 >> (u32(16)))) & (u32(255)))), f32((((c_0 >> (u32(24)))) & (u32(255))))) / vec4<f32>(255.0f);\n}\n\nfn rampMidT_0( t_0 : f32,  mid_0 : f32) -> f32\n{\n    var _S1 : f32 = max(mid_0, 1.0f - mid_0);\n    return clamp((t_0 + (_S1 - mid_0)) / (2.0f * _S1), 0.0f, 1.0f);\n}\n\nfn rampColor_0( ramp_texture_0 : texture_2d<f32>,  ramp_sampler_0 : sampler,  t_1 : f32) -> vec4<f32>\n{\n    return (textureSampleLevel((ramp_texture_0), (ramp_sampler_0), (vec2<f32>((clamp(t_1, 0.0f, 1.0f) * 255.0f + 0.5f) / 256.0f, 0.5f)), (0.0f)));\n}\n\nfn colorCutAt_0( i_0 : i32,  cutsLo_0 : vec4<f32>,  cutsHi_0 : vec4<f32>) -> f32\n{\n    var _S2 : f32;\n    if(i_0 < i32(4))\n    {\n        _S2 = cutsLo_0[i_0];\n    }\n    else\n    {\n        _S2 = cutsHi_0[i_0 - i32(4)];\n    }\n    return _S2;\n}\n\nfn thresholdBand_0( value_0 : f32,  cutCount_0 : i32,  cutsLo_1 : vec4<f32>,  cutsHi_1 : vec4<f32>) -> i32\n{\n    var band_0 : i32 = i32(0);\n    var i_1 : i32 = i32(0);\n    for(;;)\n    {\n        if(i_1 < cutCount_0)\n        {\n        }\n        else\n        {\n            break;\n        }\n        if(value_0 >= (colorCutAt_0(i_1, cutsLo_1, cutsHi_1)))\n        {\n            band_0 = i_1 + i32(1);\n        }\n        i_1 = i_1 + i32(1);\n    }\n    return band_0;\n}\n\nfn log1pf_0( x_0 : f32) -> f32\n{\n    var u_1 : f32 = 1.0f + x_0;\n    var _S3 : f32;\n    if(u_1 == 1.0f)\n    {\n        _S3 = x_0;\n    }\n    else\n    {\n        _S3 = log(u_1) * (x_0 / (u_1 - 1.0f));\n    }\n    return _S3;\n}\n\nfn symlogTransform_0( x_1 : f32,  c_1 : f32) -> f32\n{\n    return f32((i32(sign((x_1))))) * log1pf_0(abs(x_1 / c_1));\n}\n\nfn normalizeScoreUnclamped_0( score_0 : f32,  domainMin_0 : f32,  domainMax_0 : f32,  scaleType_0 : i32,  symlogConstant_0 : f32) -> f32\n{\n    var tMax_0 : f32;\n    var tMin_0 : f32;\n    var t_2 : f32;\n    if(scaleType_0 == i32(2))\n    {\n        var _S4 : f32 = symlogTransform_0(score_0, symlogConstant_0);\n        var _S5 : f32 = symlogTransform_0(domainMin_0, symlogConstant_0);\n        tMax_0 = symlogTransform_0(domainMax_0, symlogConstant_0);\n        tMin_0 = _S5;\n        t_2 = _S4;\n    }\n    else\n    {\n        if(scaleType_0 == i32(1))\n        {\n            if(domainMin_0 > 0.0f)\n            {\n                tMax_0 = domainMin_0;\n            }\n            else\n            {\n                tMax_0 = 1.0f;\n            }\n            var _S6 : f32 = log2(max(score_0, tMax_0));\n            var _S7 : f32 = log2(tMax_0);\n            tMax_0 = log2(max(domainMax_0, tMax_0));\n            tMin_0 = _S7;\n            t_2 = _S6;\n        }\n        else\n        {\n            tMax_0 = domainMax_0;\n            tMin_0 = domainMin_0;\n            t_2 = score_0;\n        }\n    }\n    var tRange_0 : f32 = tMax_0 - tMin_0;\n    if(tRange_0 <= 0.0f)\n    {\n        if(t_2 > tMin_0)\n        {\n            tMax_0 = 1.0f;\n        }\n        else\n        {\n            tMax_0 = 0.0f;\n        }\n        return tMax_0;\n    }\n    return (t_2 - tMin_0) / tRange_0;\n}\n\nfn normalizeScore_0( score_1 : f32,  domainMin_1 : f32,  domainMax_1 : f32,  scaleType_1 : i32,  symlogConstant_1 : f32) -> f32\n{\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1), 0.0f, 1.0f);\n}\n\nfn rowTableWidth_0( keys_0 : u32) -> u32\n{\n    return max(min(keys_0, u32(2048)), u32(1));\n}\n\nfn rowTablePlaneHeight_0( keys_1 : u32) -> u32\n{\n    var w_0 : u32 = rowTableWidth_0(keys_1);\n    var _S8 : u32 = (keys_1 + w_0 - u32(1)) / w_0;\n    return max(_S8, u32(1));\n}\n\nfn rowTableTexelX_0( key_0 : u32,  keys_2 : u32) -> u32\n{\n    var _S9 : u32 = key_0 % rowTableWidth_0(keys_2);\n    return _S9;\n}\n\nfn rowTableTexelY_0( key_1 : u32,  keys_3 : u32) -> u32\n{\n    var _S10 : u32 = key_1 / rowTableWidth_0(keys_3);\n    return _S10;\n}\n\nfn byteOf_0( channel_0 : f32) -> u32\n{\n    return u32(channel_0 * 255.0f + 0.5f);\n}\n\nstruct RowPlacement_0\n{\n     hidden_0 : bool,\n     slot_0 : u32,\n     color_0 : vec4<f32>,\n};\n\nfn rowTableLookup_0( table_texture_0 : texture_2d<f32>,  table_sampler_0 : sampler,  key_2 : u32,  keys_4 : i32,  own_0 : vec4<f32>) -> RowPlacement_0\n{\n    var p_0 : RowPlacement_0;\n    p_0.hidden_0 = false;\n    p_0.slot_0 = key_2;\n    p_0.color_0 = own_0;\n    if(keys_4 < i32(0))\n    {\n        return p_0;\n    }\n    var n_0 : u32 = u32(keys_4);\n    if(key_2 >= n_0)\n    {\n        p_0.hidden_0 = true;\n        return p_0;\n    }\n    var h_0 : u32 = rowTablePlaneHeight_0(n_0);\n    var _S11 : u32 = rowTableTexelX_0(key_2, n_0);\n    var u_2 : f32 = (f32(_S11) + 0.5f) / f32(rowTableWidth_0(n_0));\n    var rows_0 : f32 = f32(u32(2) * h_0);\n    var y_0 : u32 = rowTableTexelY_0(key_2, n_0);\n    var _S12 : vec4<f32> = (textureSampleLevel((table_texture_0), (table_sampler_0), (vec2<f32>(u_2, (f32(y_0) + 0.5f) / rows_0)), (0.0f)));\n    if((_S12.w) < 0.5f)\n    {\n        p_0.hidden_0 = true;\n        return p_0;\n    }\n    p_0.slot_0 = ((((byteOf_0(_S12.x)) | ((((byteOf_0(_S12.y)) << (u32(8))))))) | ((((byteOf_0(_S12.z)) << (u32(16))))));\n    var _S13 : vec4<f32> = (textureSampleLevel((table_texture_0), (table_sampler_0), (vec2<f32>(u_2, (f32(y_0 + h_0) + 0.5f) / rows_0)), (0.0f)));\n    var _S14 : vec4<f32>;\n    if((_S13.w) > 0.0f)\n    {\n        _S14 = _S13;\n    }\n    else\n    {\n        _S14 = own_0;\n    }\n    p_0.color_0 = _S14;\n    return p_0;\n}\n\nfn quadLocal_0( vid_0 : u32) -> vec2<f32>\n{\n    var v_0 : u32 = vid_0 % u32(6);\n    var _S15 : bool = v_0 == u32(0);\n    var _S16 : bool;\n    if(_S15)\n    {\n        _S16 = true;\n    }\n    else\n    {\n        _S16 = v_0 == u32(2);\n    }\n    if(_S16)\n    {\n        _S16 = true;\n    }\n    else\n    {\n        _S16 = v_0 == u32(3);\n    }\n    var _S17 : f32;\n    if(_S16)\n    {\n        _S17 = 0.0f;\n    }\n    else\n    {\n        _S17 = 1.0f;\n    }\n    if(_S15)\n    {\n        _S16 = true;\n    }\n    else\n    {\n        _S16 = v_0 == u32(1);\n    }\n    if(_S16)\n    {\n        _S16 = true;\n    }\n    else\n    {\n        _S16 = v_0 == u32(4);\n    }\n    var _S18 : f32;\n    if(_S16)\n    {\n        _S18 = 0.0f;\n    }\n    else\n    {\n        _S18 = 1.0f;\n    }\n    return vec2<f32>(_S17, _S18);\n}\n\nfn drawnRowHeightPx_0( rowHeight_1 : f32,  rowProportion_1 : f32) -> f32\n{\n    return max(rowHeight_1 * rowProportion_1, 1.0f);\n}\n\nfn rowBandOffsetPx_0( rowHeight_2 : f32,  rowProportion_2 : f32) -> f32\n{\n    return (rowHeight_2 - drawnRowHeightPx_0(rowHeight_2, rowProportion_2)) * 0.5f;\n}\n\nstruct RowBand_0\n{\n     height_0 : f32,\n     offset_0 : f32,\n};\n\nfn rowBandPx_0( rowHeight_3 : f32,  rowProportion_3 : f32) -> RowBand_0\n{\n    var band_1 : RowBand_0;\n    band_1.height_0 = drawnRowHeightPx_0(rowHeight_3, rowProportion_3);\n    band_1.offset_0 = rowBandOffsetPx_0(rowHeight_3, rowProportion_3);\n    return band_1;\n}\n\nfn hpSplitUint_0( value_1 : u32) -> vec2<f32>\n{\n    var lo_0 : u32 = (value_1 & (u32(4095)));\n    return vec2<f32>(f32(value_1 - lo_0), f32(lo_0));\n}\n\nfn hpToClipX_0( splitPos_0 : vec2<f32>,  bpRange_0 : vec3<f32>,  hpZero_0 : f32) -> f32\n{\n    var floorBp_0 : f32 = -3.00000000549775576e+38f + hpZero_0;\n    var step_0 : f32 = 2.0f / bpRange_0.z;\n    return dot(vec3<f32>(-1.0f, max(splitPos_0.x - bpRange_0.x, floorBp_0), max(splitPos_0.y - bpRange_0.y, floorBp_0)), vec3<f32>(1.0f, step_0, step_0));\n}\n\nfn pxToClipLen_0( px_0 : f32,  size_0 : f32) -> f32\n{\n    return px_0 * 2.0f / size_0;\n}\n\nfn extendToMinWidthPx_0( x1_0 : f32,  x2_0 : f32,  minWidth_0 : f32) -> f32\n{\n    var dx_0 : f32 = x2_0 - x1_0;\n    var _S19 : f32;\n    if((abs(dx_0)) < minWidth_0)\n    {\n        if(dx_0 < 0.0f)\n        {\n            _S19 = - minWidth_0;\n        }\n        else\n        {\n            _S19 = minWidth_0;\n        }\n        _S19 = x1_0 + _S19;\n    }\n    else\n    {\n        _S19 = x2_0;\n    }\n    return _S19;\n}\n\nfn extendToMinWidthX_0( x1_1 : f32,  x2_1 : f32,  minWidthPx_0 : f32,  canvasWidth_0 : f32) -> f32\n{\n    return extendToMinWidthPx_0(x1_1, x2_1, pxToClipLen_0(minWidthPx_0, canvasWidth_0));\n}\n\nfn yPxToClipY_0( yPx_0 : f32,  canvasHeight_1 : f32) -> f32\n{\n    return 1.0f - yPx_0 / canvasHeight_1 * 2.0f;\n}\n\nstruct RowRectUniforms_0\n{\n     bpRangeX_1 : vec3<f32>,\n     canvasHeight_2 : f32,\n     minCellDenomPx_1 : f32,\n     minCellPx_1 : f32,\n     zero_1 : f32,\n     rowHeight_4 : f32,\n     rowProportion_4 : f32,\n     scrollTop_1 : f32,\n};\n\nfn rowRectClipPos_0( startBp_0 : u32,  endBp_0 : u32,  rowIndex_0 : f32,  vid_1 : u32,  u_3 : RowRectUniforms_0) -> vec4<f32>\n{\n    var local_0 : vec2<f32> = quadLocal_0(vid_1);\n    var band_2 : RowBand_0 = rowBandPx_0(u_3.rowHeight_4, u_3.rowProportion_4);\n    var rowY_0 : f32 = band_2.offset_0 + u_3.rowHeight_4 * rowIndex_0 - u_3.scrollTop_1;\n    var sx1_0 : f32 = hpToClipX_0(hpSplitUint_0(startBp_0), u_3.bpRangeX_1, u_3.zero_1);\n    return vec4<f32>(mix(sx1_0, extendToMinWidthX_0(sx1_0, hpToClipX_0(hpSplitUint_0(endBp_0), u_3.bpRangeX_1, u_3.zero_1), u_3.minCellPx_1, u_3.minCellDenomPx_1), local_0.x), yPxToClipY_0(mix(rowY_0, rowY_0 + band_2.height_0, local_0.y), u_3.canvasHeight_2), 0.0f, 1.0f);\n}\n\nstruct RowRectVsOut_0\n{\n    @builtin(position) position_0 : vec4<f32>,\n    @location(0) color_1 : vec4<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) x_2 : u32,\n    @location(1) x2_2 : u32,\n    @location(2) row_0 : u32,\n    @location(3) color_2 : u32,\n};\n\nfn markScaleColor_0( _S20 : texture_2d<f32>,  _S21 : sampler,  _S22 : u32,  _S23 : i32,  _S24 : f32,  _S25 : f32,  _S26 : f32,  _S27 : i32,  _S28 : vec4<f32>,  _S29 : vec4<f32>) -> vec4<f32>\n{\n    if(_S23 == i32(0))\n    {\n        return unpackRGBA_0(_S22);\n    }\n    if(_S22 == u32(2143289345))\n    {\n        return unpackRGBA_0(u32(4289703855));\n    }\n    var scaleType_2 : i32;\n    if(((_S22 & (u32(2139095040)))) == u32(2139095040))\n    {\n        if(((_S22 & (u32(8388607)))) != u32(0))\n        {\n            return unpackRGBA_0(u32(4286611584));\n        }\n        var negative_0 : bool = ((_S22 & (u32(2147483648)))) != u32(0);\n        if(_S23 == i32(3))\n        {\n            if(negative_0)\n            {\n                scaleType_2 = i32(0);\n            }\n            else\n            {\n                scaleType_2 = _S27;\n            }\n            return u_0.colorBands_0[scaleType_2];\n        }\n        var _S30 : f32;\n        if(negative_0)\n        {\n            _S30 = 0.0f;\n        }\n        else\n        {\n            _S30 = 1.0f;\n        }\n        return rampColor_0(_S20, _S21, rampMidT_0(_S30, _S26));\n    }\n    var value_2 : f32 = (bitcast<f32>((_S22)));\n    if(_S23 == i32(3))\n    {\n        return u_0.colorBands_0[thresholdBand_0(value_2, _S27, _S28, _S29)];\n    }\n    if(_S23 == i32(2))\n    {\n        scaleType_2 = i32(1);\n    }\n    else\n    {\n        scaleType_2 = i32(0);\n    }\n    return rampColor_0(_S20, _S21, rampMidT_0(normalizeScore_0(value_2, _S24, _S25, scaleType_2, 1.0f), _S26));\n}\n\n@vertex\nfn vs_main( _S31 : vertexInput_0, @builtin(vertex_index) vid_2 : u32) -> RowRectVsOut_0\n{\n    var _S32 : RowPlacement_0 = rowTableLookup_0(rowTable_texture_0, rowTable_sampler_0, _S31.row_0, u_0.rowTableKeys_0, markScaleColor_0(colorRamp_texture_0, colorRamp_sampler_0, _S31.color_2, u_0.rampMode_0, u_0.rampMin_0, u_0.rampMax_0, u_0.rampMidNorm_0, u_0.colorCutCount_0, u_0.colorCuts_0[i32(0)], u_0.colorCuts_0[i32(1)]));\n    var o_0 : RowRectVsOut_0;\n    if(_S32.hidden_0)\n    {\n        o_0.position_0 = vec4<f32>(2.0f, 2.0f, 0.0f, 1.0f);\n        o_0.color_1 = vec4<f32>(0.0f);\n        return o_0;\n    }\n    var r_0 : RowRectUniforms_0;\n    r_0.bpRangeX_1 = u_0.bpRangeX_0;\n    r_0.canvasHeight_2 = u_0.canvasHeight_0;\n    r_0.minCellDenomPx_1 = u_0.minCellDenomPx_0;\n    r_0.minCellPx_1 = u_0.minCellPx_0;\n    r_0.zero_1 = u_0.zero_0;\n    r_0.rowHeight_4 = u_0.rowHeight_0;\n    r_0.rowProportion_4 = u_0.rowProportion_0;\n    r_0.scrollTop_1 = u_0.scrollTop_0;\n    o_0.position_0 = rowRectClipPos_0(_S31.x_2, _S31.x2_2, f32(_S32.slot_0), vid_2, r_0);\n    o_0.color_1 = _S32.color_0;\n    return o_0;\n}\n\nstruct pixelOutput_0\n{\n    @location(0) output_0 : vec4<f32>,\n};\n\nstruct pixelInput_0\n{\n    @location(0) color_3 : vec4<f32>,\n};\n\n@fragment\nfn fs_main( _S33 : pixelInput_0, @builtin(position) position_1 : vec4<f32>) -> pixelOutput_0\n{\n    var _S34 : pixelOutput_0 = pixelOutput_0( _S33.color_3 );\n    return _S34;\n}\n\n"