@jbrowse/render-core 5.0.0-beta.11 → 5.0.0-beta.14

package/src/shaders/curveDistance.js.generated.ts

@@ -4,15 +4,27 @@
 // Scalar twins of curveDistance.slang, transliterated from slangc's WGSL so
 // the Canvas2D and SVG paths run the shader's own math. See adr-051.
 
+function _clamp(x: number, lo: number, hi: number) {
+  return _min(_max(x, lo), hi)
+}
+
 function _max(a: number, b: number) {
   return b > a || Number.isNaN(a) ? b : a
 }
 
+function _min(a: number, b: number) {
+  return b < a || Number.isNaN(a) ? b : a
+}
+
 export function distToWideCirclePx(x: number, y: number, r: number): number {
   let rk = ((2.0 * x) + (((x * x) + (y * y)) / r))
   return Math.abs((rk / (Math.sqrt(_max((1.0 + (rk / r)), 0.0)) + 1.0)))
 }
 
+export function legSweepAngle(radiusPx: number, yReachPx: number): number {
+  return Math.asin(_clamp((yReachPx / _max(radiusPx, 9.99999997475242708e-07)), 0.0, 1.0))
+}
+
 export function wideCircleLeg(seg: number, segments: number, legSweep: number): [number, number] {
   let half = Math.trunc(segments / 2)
   if ((seg <= half)) {
@@ -89,3 +101,103 @@ export function wideCircleHullPoint(sinB: number, sinHalfB: number, cosHalf: num
   let sh2 = (sinHalfB * sinHalfB)
   return [((grow * (1.0 - (2.0 * sh2))) - ((2.0 * radiusPx) * sh2)), ((radiusPx + grow) * sinB)]
 }
+
+export function farDomeParam(dPx: number, rxPx: number): number {
+  return Math.sqrt(_clamp((dPx / (2.0 * _max(rxPx, 9.99999997475242708e-07))), 0.0, 0.8399999737739563))
+}
+
+export function farDomeHullPoint(u: number, stepU: number, rxPx: number, ryPx: number, padPx: number): [number, number] {
+  let _t0 = _max(ryPx, 9.99999997475242708e-07)
+  let _t1 = _max(u, (0.5 * stepU))
+  let nx = (-_t0 * (1.0 - ((2.0 * _t1) * _t1)))
+  let _t2 = (2.0 * rxPx)
+  let ny = ((_t2 * _t1) * Math.sqrt((1.0 - (_t1 * _t1))))
+  let _t3 = _max(((nx * nx) + (ny * ny)), 9.99999968265522539e-21)
+  let reach = (padPx + (((((((Math.sign(padPx)) | 0)) * 0.5) * rxPx) * stepU) * stepU))
+  let _t4 = Math.sqrt(_t3)
+  let along = (reach / _t4)
+  let d = ((_t2 * u) * u)
+  let _t5 = (u * u)
+  let h = (((2.0 * _t0) * u) * Math.sqrt((1.0 - _t5)))
+  let folded: boolean
+  if ((padPx < 0.0)) {
+    folded = (_t5 < 0.5)
+  } else {
+    folded = false
+  }
+  if (folded) {
+    folded = ((_t3 * _t4) < (((-2.0 * reach) * rxPx) * _t0))
+  } else {
+    folded = false
+  }
+  let _t6: number
+  if (folded) {
+    if ((u > 0.0)) {
+      _t6 = (d - reach)
+    } else {
+      _t6 = d
+    }
+  } else {
+    _t6 = (d + (nx * along))
+  }
+  let _t7: number
+  if (folded) {
+    _t7 = (h + reach)
+  } else {
+    _t7 = (h + (ny * along))
+  }
+  return [_t6, _t7]
+}
+
+function farDomeDistSq(u: number, dPx: number, hPx: number, rx: number, ry: number): number {
+  let ox = ((((2.0 * rx) * u) * u) - dPx)
+  let oy = ((((2.0 * ry) * u) * Math.sqrt(_max((1.0 - (u * u)), 0.0))) - hPx)
+  return ((ox * ox) + (oy * oy))
+}
+
+function farDomeStep(u: number, dPx: number, hPx: number, rx: number, ry: number): number {
+  let w = Math.sqrt(_max((1.0 - (u * u)), 9.99999997475242708e-07))
+  let ox = ((((2.0 * rx) * u) * u) - dPx)
+  let _t0 = (2.0 * ry)
+  let _t1 = (_t0 * u)
+  let oy = ((_t1 * w) - hPx)
+  let _t2 = (4.0 * rx)
+  let dx = (_t2 * u)
+  let _t3 = ((2.0 * u) * u)
+  let dy = ((_t0 * (1.0 - _t3)) / w)
+  let whole = _clamp((u - (((ox * dx) + (oy * dy)) / _max(((((dx * dx) + (dy * dy)) + (ox * _t2)) + (oy * ((_t1 * (_t3 - 3.0)) / ((w * w) * w)))), 9.999999960041972e-13))), 0.0, 0.99900001287460327)
+  let half = (0.5 * (u + whole))
+  let here = ((ox * ox) + (oy * oy))
+  let _t4: number
+  if ((farDomeDistSq(whole, dPx, hPx, rx, ry) < here)) {
+    _t4 = whole
+  } else {
+    if ((farDomeDistSq(half, dPx, hPx, rx, ry) < here)) {
+      _t4 = half
+    } else {
+      _t4 = u
+    }
+  }
+  return _t4
+}
+
+export function farDomeDistancePx(dPx: number, hPx: number, rxPx: number, ryPx: number): number {
+  let _t0 = _max(rxPx, 9.99999997475242708e-07)
+  let _t1 = _max(ryPx, 9.99999997475242708e-07)
+  let q = _clamp((hPx / _t1), 0.0, 1.0)
+  let c = Math.sqrt((1.0 - (q * q)))
+  let uH: number
+  if ((dPx > _t0)) {
+    uH = Math.sqrt((0.5 * (1.0 + c)))
+  } else {
+    uH = (q / Math.sqrt((2.0 * (1.0 + c))))
+  }
+  let uD = Math.sqrt(_clamp((dPx / (2.0 * _t0)), 0.0, 1.0))
+  let u: number
+  if ((farDomeDistSq(uH, dPx, hPx, _t0, _t1) < farDomeDistSq(uD, dPx, hPx, _t0, _t1))) {
+    u = uH
+  } else {
+    u = uD
+  }
+  return Math.sqrt(farDomeDistSq(farDomeStep(farDomeStep(farDomeStep(farDomeStep(u, dPx, hPx, _t0, _t1), dPx, hPx, _t0, _t1), dPx, hPx, _t0, _t1), dPx, hPx, _t0, _t1), dPx, hPx, _t0, _t1))
+}

package/src/shaders/curveDistance.slang

@@ -7,10 +7,14 @@
 // The fragment ramps its alpha to zero at `halfWidth + aaHalfPx`, so the hull
 // has to cover everything nearer the curve than that. Both flanks of a hull
 // sampled from the offset curve fall short of it at a clamped dome's feet.
-//! js-export: distToWideCirclePx, wideCircleLeg, wideCircleLegStep, ellipseHullPoint, wideCircleHullPoint
+//! js-export: distToWideCirclePx, legSweepAngle, wideCircleLeg, wideCircleLegStep, ellipseHullPoint, wideCircleHullPoint, farDomeParam, farDomeHullPoint, farDomeDistancePx
 //! js-skip: nearCircleDistancePx — a float32 remedy with no float64 caller. The hit tests' copy of this file's solve (marks/ellipseDistance.ts) holds the ellipse at every aspect, so nothing outside the shader has an approximation to want.
+//! js-skip: farDomeDistSq — a step of farDomeDistancePx's solve, which its twin carries; nothing outside it wants one step
+//! js-skip: farDomeStep — a step of farDomeDistancePx's solve, so it inherits farDomeDistSq's reason exactly
 //! oracle-skip: ellipseHullPoint — a miter between near-parallel tangent lines, so float32 and float64 part on the probe's arbitrary directions and radii; over the angles the hull steps through they agree to 0.004 px, and arcHull.test.ts referees the twin against every inked point instead
 //! oracle-skip: wideCircleHullPoint — the foot-relative frame cancels at half steps and angles the band never reaches; arcHull.test.ts referees it with the legs
+//! oracle-skip: farDomeDistancePx — at the probe's 65536 px ry, q = h / ry and 1 - q² round to the same float32, so the solve starts from the wrong foot and lands up to 18 px out; ry is the arc band, tens of px, where float32 and float64 agree, and farDome.test.ts referees the twin against an exact solve
+//! oracle-skip: farDomeStep — a step of farDomeDistancePx's solve, so it inherits that reason; the probe also feeds it negative radii and widths the solve never reaches
 
 module curveDistance;
 
@@ -163,6 +167,105 @@ public float distToWideCircle(float2 legLocal, float r) {
   return distToWideCirclePx(legLocal.x, legLocal.y, r);
 }
 
+// A far pair whose apex is low is not a circle: it is a dome thousands of px
+// wide and tens tall, which over one canvas is a nearly level line. The
+// functions below measure it from a foot, so nothing cancels at these radii.
+//
+// u = sin(b / 2) for the angle b from a foot at which the dome is `d` px inward
+// of that foot, by the half-angle identity d = 2 rx sin^2(b / 2). Held short of
+// the far foot, where the height's slope in u has no bound: the foot is the
+// nearer one, and a far pair's canvas reaches under a third of the way past
+// the apex.
+public float farDomeParam(float dPx, float rxPx) {
+  return sqrt(clamp(dPx / (2.0 * max(rxPx, 1e-6)), 0.0, 0.84));
+}
+
+// One vertex of a far dome's strip in the frame `farDomeDistancePx` measures
+// in, inward of a foot and up: the dome's point at `u`, `padPx` along its own
+// normal, positive on the outer flank. `stepU` is the strip's step in u. The
+// dome leaves its foot upright and is level within a px, so the foot's normal
+// is its first chord's, and a chord cuts inside the curve by its sagitta,
+// which the pad takes on.
+//
+// A flat dome turns tighter than the pad at its foot (radius ry^2 / rx), where
+// the inner offset folds and its chords cut the ink. There the inner flank is
+// the dome moved a pad inward and a pad down instead: the dome rises as it
+// runs inward, so that copy clears the pad around every point of it, and it is
+// as concave as the dome, so its chords clear it too. The foot's own vertex
+// stays under the foot, where the strip's end crosses it.
+public float2 farDomeHullPoint(float u, float stepU, float rxPx, float ryPx, float padPx) {
+  float ry = max(ryPx, 1e-6);
+  float un = max(u, 0.5 * stepU);
+  float nx = -ry * (1.0 - 2.0 * un * un);
+  float ny = 2.0 * rxPx * un * sqrt(1.0 - un * un);
+  float lenSq = max(nx * nx + ny * ny, 1e-20);
+  float reach = padPx + sign(padPx) * 0.5 * rxPx * stepU * stepU;
+  float along = reach / sqrt(lenSq);
+  float d = 2.0 * rxPx * u * u;
+  float h = 2.0 * ry * u * sqrt(1.0 - u * u);
+  // the radius of curvature, against twice the pad
+  bool folded = padPx < 0.0 && u * u < 0.5
+    && lenSq * sqrt(lenSq) < -2.0 * reach * rxPx * ry;
+  return float2(folded ? (u > 0.0 ? d - reach : d) : d + nx * along,
+                folded ? h + reach : h + ny * along);
+}
+
+// A far dome's point in the frame of one of its feet, by u = sin(b / 2) for the
+// angle b from that foot: 2 rx u^2 inward of the foot, 2 ry u sqrt(1 - u^2)
+// up. No trigonometry, so the twin the hit tests' referee runs is this code.
+float farDomeDistSq(float u, float dPx, float hPx, float rx, float ry) {
+  float ox = 2.0 * rx * u * u - dPx;
+  float oy = 2.0 * ry * u * sqrt(max(1.0 - u * u, 0.0)) - hPx;
+  return ox * ox + oy * oy;
+}
+
+// One Newton step toward the u whose dome point is nearest the fragment. Held
+// short of 1, the far foot, where the height's slope in u has no bound; the
+// frame is the nearer foot's, so the answer is never there.
+float farDomeStep(float u, float dPx, float hPx, float rx, float ry) {
+  float w = sqrt(max(1.0 - u * u, 1e-6));
+  float ox = 2.0 * rx * u * u - dPx;
+  float oy = 2.0 * ry * u * w - hPx;
+  float dx = 4.0 * rx * u;
+  float dy = 2.0 * ry * (1.0 - 2.0 * u * u) / w;
+  float ddx = 4.0 * rx;
+  float ddy = 2.0 * ry * u * (2.0 * u * u - 3.0) / (w * w * w);
+  float f = ox * dx + oy * dy;
+  float df = dx * dx + dy * dy + ox * ddx + oy * ddy;
+  // A foot is a tip sharper than a px, where the step overshoots: it is taken
+  // whole or by half, and only where it lands nearer.
+  float whole = clamp(u - f / max(df, 1e-12), 0.0, 0.999);
+  float half = 0.5 * (u + whole);
+  float here = ox * ox + oy * oy;
+  return farDomeDistSq(whole, dPx, hPx, rx, ry) < here
+    ? whole
+    : (farDomeDistSq(half, dPx, hPx, rx, ry) < here ? half : u);
+}
+
+// Distance from a fragment to a far dome, in the frame of one of its feet: the
+// fragment `dPx` inward of the foot and `hPx` above the baseline. Newton's
+// method from the nearer of the dome's point level with the fragment, which is
+// the near one where the dome is steep at its feet, and its point above or
+// below it, which is near everywhere else. Every term is assembled from the
+// foot, so a width in the millions of px cancels nothing.
+public float farDomeDistancePx(float dPx, float hPx, float rxPx, float ryPx) {
+  float rx = max(rxPx, 1e-6);
+  float ry = max(ryPx, 1e-6);
+  float q = clamp(hPx / ry, 0.0, 1.0);
+  float c = sqrt(1.0 - q * q);
+  // sin(b / 2) for sin(b) = q, on this foot's half or the far one's
+  float uH = dPx > rx ? sqrt(0.5 * (1.0 + c)) : q / sqrt(2.0 * (1.0 + c));
+  float uD = sqrt(clamp(dPx / (2.0 * rx), 0.0, 1.0));
+  float u = farDomeDistSq(uH, dPx, hPx, rx, ry) < farDomeDistSq(uD, dPx, hPx, rx, ry)
+    ? uH
+    : uD;
+  u = farDomeStep(u, dPx, hPx, rx, ry);
+  u = farDomeStep(u, dPx, hPx, rx, ry);
+  u = farDomeStep(u, dPx, hPx, rx, ry);
+  u = farDomeStep(u, dPx, hPx, rx, ry);
+  return sqrt(farDomeDistSq(u, dPx, hPx, rx, ry));
+}
+
 // How far around a far circle each leg's hull sweeps, in radians from its foot:
 // to the band's edge and no further. Sweeping the whole quarter spends nearly
 // every segment off-screen, and at r = 3e6 px the whole band then falls inside

package/src/shaders/farDome.test.ts

@@ -0,0 +1,117 @@
+import { ellipseNearest } from '../marks/ellipseDistance.ts'
+import {
+  farDomeDistancePx,
+  farDomeParam,
+} from './curveDistance.js.generated.ts'
+
+// A far dome is an ellipse thousands to millions of px wide and tens tall. The
+// fragment measures to it from one of its feet (`farDomeDistancePx`), where
+// nothing cancels, and these hold that to the exact ellipse distance at the
+// widths a far pair reaches and the strokes a link draws, over the half of the
+// dome nearer the foot it measures from, which is the half the shader asks of
+// it. Under the baseline the stroke is still the rising half's, where the exact
+// solve would answer with the mirrored half, so the referee there is the foot
+// or the level point, whichever is nearer.
+
+// A point `offPx` along the dome's outward normal from its point at angle `b`
+// from the left foot: inward of that foot by `d` and up by `h`.
+function probe(rx: number, ry: number, b: number, offPx: number) {
+  const nx = -Math.cos(b) * ry
+  const ny = Math.sin(b) * rx
+  const len = Math.hypot(nx, ny)
+  return {
+    d: 2 * rx * Math.sin(b / 2) ** 2 + (nx / len) * offPx,
+    h: ry * Math.sin(b) + (ny / len) * offPx,
+  }
+}
+
+describe('farDomeDistancePx', () => {
+  const domes: [number, number][] = [
+    [3000, 40],
+    [3000, 300],
+    [100_000, 40],
+    [1_000_000, 40],
+    [1_000_000, 90],
+    [10_000_000, 90],
+  ]
+  const angles = [
+    0.0005,
+    0.002,
+    0.01,
+    0.05,
+    0.2,
+    0.6,
+    1,
+    1.4,
+    Math.PI / 2,
+    1.7,
+    1.9,
+  ]
+  const offsets = [-8, -4, -1.5, -0.5, 0, 0.5, 1.5, 4, 8]
+
+  test.each(domes)('rx %d, ry %d agrees with the ellipse', (rx, ry) => {
+    let worst = 0
+    for (const b of angles) {
+      for (const off of offsets) {
+        const { d, h } = probe(rx, ry, b, off)
+        if (h >= 0) {
+          const exact = ellipseNearest(d - rx, -h, rx, ry).dist
+          worst = Math.max(
+            worst,
+            Math.abs(farDomeDistancePx(d, h, rx, ry) - exact),
+          )
+        }
+      }
+    }
+    // the worst is inside a foot, a tip sharper than a px, where no step
+    // lands nearer than the start
+    expect(worst).toBeLessThan(0.05)
+  })
+
+  test.each(domes)(
+    'rx %d, ry %d measures to the rising half from under the baseline',
+    (rx, ry) => {
+      for (const d of [-3, 0, 0.5, 3, 40, 900]) {
+        for (const h of [-0.5, -2, -5]) {
+          let exact = Math.hypot(d, h)
+          for (let k = 1; k <= 20_000; k++) {
+            const x = (k / 20_000) * Math.max(4 * d, 40)
+            const y = ry * Math.sqrt(Math.max(1 - (1 - x / rx) ** 2, 0))
+            exact = Math.min(exact, Math.hypot(d - x, h - y))
+          }
+          // the hull hands the fragment nothing further off than its pad
+          if (exact < 10) {
+            expect(
+              Math.abs(farDomeDistancePx(d, h, rx, ry) - exact),
+            ).toBeLessThan(0.05)
+          }
+        }
+      }
+    },
+  )
+
+  it('measures to the foot from outside the pair', () => {
+    expect(farDomeDistancePx(-5, 0, 1_000_000, 40)).toBeCloseTo(5)
+  })
+
+  it('is the height short of the apex in the middle of the pair', () => {
+    expect(farDomeDistancePx(1_000_000, 37, 1_000_000, 40)).toBeCloseTo(3)
+    expect(farDomeDistancePx(1_000_500, 44, 1_000_000, 40)).toBeCloseTo(4)
+  })
+})
+
+describe('farDomeParam', () => {
+  it('is sin(b / 2) where the dome is that far inward of its foot', () => {
+    for (const rx of [3000, 1_000_000]) {
+      for (const b of [0.001, 0.3, 1.2, Math.PI / 2, 2.2]) {
+        const d = 2 * rx * Math.sin(b / 2) ** 2
+        expect(farDomeParam(d, rx)).toBeCloseTo(Math.sin(b / 2), 6)
+      }
+    }
+  })
+
+  it('holds at the foot outside the pair and short of the far foot', () => {
+    expect(farDomeParam(-50, 3000)).toBe(0)
+    expect(farDomeParam(7000, 3000)).toBeCloseTo(Math.sqrt(0.84))
+  })
+})

package/src/shaders/hpClampFloor.test.ts

@@ -12,8 +12,6 @@ import { GLSL_VERTEX as barMarkGlsl } from './barMark.glsl.generated.ts'
 import { WGSL_SOURCE as barMarkWgsl } from './barMark.wgsl.generated.ts'
 import { GLSL_VERTEX as coverageBarGlsl } from './coverageBar.glsl.generated.ts'
 import { WGSL_SOURCE as coverageBarWgsl } from './coverageBar.wgsl.generated.ts'
-import { GLSL_VERTEX as coverageIndicatorGlsl } from './coverageIndicator.glsl.generated.ts'
-import { WGSL_SOURCE as coverageIndicatorWgsl } from './coverageIndicator.wgsl.generated.ts'
 import { GLSL_VERTEX as coverageInterbaseGlsl } from './coverageInterbase.glsl.generated.ts'
 import { WGSL_SOURCE as coverageInterbaseWgsl } from './coverageInterbase.wgsl.generated.ts'
 import { GLSL_VERTEX as coverageModGlsl } from './coverageMod.glsl.generated.ts'
@@ -30,8 +28,6 @@ const SOURCES = {
   'barMark glsl': barMarkGlsl,
   'coverageBar wgsl': coverageBarWgsl,
   'coverageBar glsl': coverageBarGlsl,
-  'coverageIndicator wgsl': coverageIndicatorWgsl,
-  'coverageIndicator glsl': coverageIndicatorGlsl,
   'coverageInterbase wgsl': coverageInterbaseWgsl,
   'coverageInterbase glsl': coverageInterbaseGlsl,
   'coverageMod wgsl': coverageModWgsl,

package/src/shaders/insertionGlyph.generated.ts

@@ -0,0 +1,14 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from insertionGlyph.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const LONG_INSERTION_MIN_LENGTH = 10
+
+export const LONG_INSERTION_TEXT_THRESHOLD_PX = 15
+
+export const MIN_HEIGHT_FOR_TEXT = 5
+
+export const INSERTION_SERIF_MIN_PX_PER_BP = 3
+
+export const SERIF_HALF_W_PX = 2
+
+export const SERIF_H_PX = 2

package/src/shaders/insertionGlyph.slang

@@ -0,0 +1,59 @@
+// The insertion glyph's geometry: the marker's width, its serif caps, and
+// whether it outgrows the reference span it stands on. The pileup's insertion
+// pass and alignments-core's `insertion` mark both import it, and the Canvas2D,
+// SVG and hit-test paths run its generated twins (adr-051).
+//
+// A marker is a 1px tick (small), a 1-5px bar (long), or a box sized for its
+// count label (large: long, zoomed in past LONG_INSERTION_TEXT_THRESHOLD_PX,
+// and in a row at least MIN_HEIGHT_FOR_TEXT tall). The long bar floors at 1px
+// so the width stays monotonic in length; a sub-pixel bar rasterizes to a
+// fraction of a column on either backend.
+//
+// Serif caps are wedges SERIF_HALF_W_PX either side of the bar, tapering
+// SERIF_H_PX into the row, on small insertions zoomed in past
+// INSERTION_SERIF_MIN_PX_PER_BP.
+//! export-consts: LONG_INSERTION_MIN_LENGTH, LONG_INSERTION_TEXT_THRESHOLD_PX, MIN_HEIGHT_FOR_TEXT, INSERTION_SERIF_MIN_PX_PER_BP, SERIF_HALF_W_PX, SERIF_H_PX
+//! consts-out: packages/alignments-core/src/insertionLabel.generated.ts
+//! js-export: insertionBarWidthPx, textWidth, insertionSerifsDraw, insertionMarkerDraws
+//! js-export-out: packages/alignments-core/src/insertionWidth.generated.ts
+module insertionGlyph;
+
+public static const uint  LONG_INSERTION_MIN_LENGTH = 10u;
+public static const float LONG_INSERTION_TEXT_THRESHOLD_PX = 15.0;
+public static const float INSERTION_SERIF_MIN_PX_PER_BP = 3.0;
+public static const float MIN_HEIGHT_FOR_TEXT = 5.0;
+public static const float SERIF_HALF_W_PX = 2.0;
+public static const float SERIF_H_PX = 2.0;
+
+static const float LABEL_CHAR_W = 6.0;
+static const float LABEL_PAD = 10.0;
+
+public float textWidth(uint num) {
+  if (num < 10u)    { return LABEL_CHAR_W + LABEL_PAD; }
+  if (num < 100u)   { return LABEL_CHAR_W * 2.0 + LABEL_PAD; }
+  if (num < 1000u)  { return LABEL_CHAR_W * 3.0 + LABEL_PAD; }
+  if (num < 10000u) { return LABEL_CHAR_W * 4.0 + LABEL_PAD; }
+  return LABEL_CHAR_W * 5.0 + LABEL_PAD;
+}
+
+public float insertionBarWidthPx(uint length, float pxPerBp, float featHeight) {
+  float insWPx = float(length) * pxPerBp;
+  bool isLong = length >= LONG_INSERTION_MIN_LENGTH;
+  bool isLarge = isLong
+    && insWPx >= LONG_INSERTION_TEXT_THRESHOLD_PX
+    && featHeight >= MIN_HEIGHT_FOR_TEXT;
+  if (isLarge) { return textWidth(length); }
+  if (isLong)  { return clamp(insWPx / 3.0, 1.0, 5.0); }
+  return 1.0;
+}
+
+public bool insertionSerifsDraw(uint length, float pxPerBp) {
+  return length < LONG_INSERTION_MIN_LENGTH && pxPerBp >= INSERTION_SERIF_MIN_PX_PER_BP;
+}
+
+// A marker standing on a reference span draws only where it is wider than
+// that span, floored at `spanFloorPx`, the narrowest the span itself paints.
+// An interbase insertion has no span and always draws.
+public bool insertionMarkerDraws(float barWidthPx, float spanPx, float spanFloorPx) {
+  return barWidthPx > max(spanFloorPx, abs(spanPx));
+}

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

@@ -1,6 +1,6 @@
 // AUTO-GENERATED by packages/shader-tools/src/shader-codegen from lineCenterMark.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 domainMin_0;\n    float domainMax_0;\n    int valueScaleType_0;\n    float valueSymlogConstant_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    float rowHeight_0;\n    float rowBandPx_0;\n    float rowOffsetPx_0;\n    float origin_0;\n    float lineWidth_0;\n    float zero_0;\n    float viewportWidth_0;\n    float devicePixelRatio_0;\n    int rowTableKeys_0;\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_1, float domainMax_1, 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_1, symlogConstant_0);\n\n        tMax_0 = symlogTransform_0(domainMax_1, 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_1 > 0.0)\n            {\n\n                tMax_0 = domainMin_1;\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_1, tMax_0));\n\n            tMin_0 = _S7;\n\n            t_1 = _S6;\n\n        }\n        else\n        {\n\n            tMax_0 = domainMax_1;\n\n            tMin_0 = domainMin_1;\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_2, float domainMax_2, int scaleType_1, float symlogConstant_1)\n{\n\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_2, domainMax_2, 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 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 rowBandTopPx_0(float rowOffsetPx_1, float rowHeight_1, uint slot_0)\n{\n\n    return rowOffsetPx_1 + rowHeight_1 * float(slot_0);\n}\n\n\nfloat valueToYPxScaled_0(float value_2, float domainMin_3, float domainMax_3, float h_0, int scaleType_2, float symlogConstant_2)\n{\n\n    return (1.0 - normalizeScore_0(value_2, domainMin_3, domainMax_3, scaleType_2, symlogConstant_2)) * h_0;\n}\n\n\nstruct CapsuleFrame_0\n{\n    vec2 tangent_0;\n    vec2 normal_0;\n};\n\n\nCapsuleFrame_0 capsuleFrame_0(vec2 delta_0, float len_0)\n{\n\n    CapsuleFrame_0 f_0;\n\n    vec2 _S11;\n    if(len_0 > 0.00009999999747379)\n    {\n\n        _S11 = delta_0 / len_0;\n\n    }\n    else\n    {\n\n        _S11 = vec2(1.0, 0.0);\n\n    }\n\n    f_0.tangent_0 = _S11;\n    f_0.normal_0 = vec2(- _S11.y, _S11.x);\n    return f_0;\n}\n\n\nfloat aaPx_0(float dpr_0)\n{\n\n    return 1.0 / max(dpr_0, 1.0);\n}\n\n\nfloat aaHalfPx_0(float dpr_1)\n{\n\n    return 0.5 * aaPx_0(dpr_1);\n}\n\n\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S12 = v_0 == 0U;\n\n    bool _S13;\n\n    if(_S12)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 2U;\n\n    }\n\n    if(_S13)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 3U;\n\n    }\n\n    float _S14;\n\n    if(_S13)\n    {\n\n        _S14 = 0.0;\n\n    }\n    else\n    {\n\n        _S14 = 1.0;\n\n    }\n    if(_S12)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 1U;\n\n    }\n\n    if(_S13)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 4U;\n\n    }\n\n    float _S15;\n\n    if(_S13)\n    {\n\n        _S15 = 0.0;\n\n    }\n    else\n    {\n\n        _S15 = 1.0;\n\n    }\n\n    return vec2(_S14, _S15);\n}\n\n\nvec2 capsuleQuadLocal_0(uint vid_1, float halfLenPx_0, float halfWidthPx_0, float dpr_2)\n{\n\n    float ext_0 = halfWidthPx_0 + aaHalfPx_0(dpr_2);\n    vec2 q_0 = quadLocal_0(vid_1) * 2.0 - 1.0;\n    return vec2(q_0.x * (halfLenPx_0 + ext_0), q_0.y * ext_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 yPxToClipY_0(float yPx_0, float canvasHeight_1)\n{\n\n    return 1.0 - yPx_0 / canvasHeight_1 * 2.0;\n}\n\n\nflat out vec4 v_color;\n\n\nout vec2 v_localPos;\n\n\nflat out float v_segHalfLen;\n\n\nflat out vec2 v_colorLine;\n\n\nflat out vec4 v_cutYsLo;\n\n\nflat out vec4 v_cutYsHi;\n\n\nlayout(location = 0)\nin uint a_x;\n\n\nlayout(location = 1)\nin uint a_x2;\n\n\nlayout(location = 2)\nin float a_y;\n\n\nlayout(location = 3)\nin uint a_color;\n\n\nlayout(location = 4)\nin uint a_row;\n\n\nlayout(location = 5)\nin uint a_prevX;\n\n\nlayout(location = 6)\nin uint a_prevX2;\n\n\nlayout(location = 7)\nin float a_prevY;\n\n\nstruct VsOut_0\n{\n    vec4 position_0;\n    vec4 color_0;\n    vec2 localPos_0;\n    float segHalfLen_0;\n    vec2 colorLine_0;\n    vec4 cutYsLo_0;\n    vec4 cutYsHi_0;\n};\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_3;\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_3 = 0;\n\n            }\n            else\n            {\n\n                scaleType_3 = _S22;\n\n            }\n\n            return u_0.colorBands_0[scaleType_3];\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_3 = uintBitsToFloat(_S17);\n    if(_S18 == 3)\n    {\n\n        return u_0.colorBands_0[thresholdBand_0(value_3, _S22, _S23, _S24)];\n    }\n    if(_S18 == 2)\n    {\n\n        scaleType_3 = 1;\n\n    }\n    else\n    {\n\n        scaleType_3 = 0;\n\n    }\n    return rampColor_0(rampMidT_0(normalizeScore_0(value_3, _S19, _S20, scaleType_3, 1.0), _S21));\n}\n\n\nstruct RowPlacement_0\n{\n    bool hidden_0;\n    uint slot_1;\n    vec4 color_1;\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_1 = _S26;\n    p_0.color_1 = _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_1 = rowTablePlaneHeight_0(n_0);\n    uint _S29 = rowTableTexelX_0(_S26, n_0);\n\n    float u_2 = (float(_S29) + 0.5) / float(rowTableWidth_0(n_0));\n    float rows_0 = float(2U * h_1);\n    uint y_0 = rowTableTexelY_0(_S26, n_0);\n    vec4 s_0 = (textureLod((u_rowTable), (vec2(u_2, (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_1 = ((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_2, (float(y_0 + h_1) + 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_1 = _S30;\n    return p_0;\n}\n\n\nfloat lineBpToClipX_0(uint _S31)\n{\n\n    return hpToClipX_0(hpSplitUint_0(_S31), u_0.bpRangeX_0, u_0.zero_0);\n}\n\n\nfloat lineValueYPx_0(float _S32, float _S33)\n{\n\n    return _S33 + valueToYPxScaled_0(_S32, u_0.domainMin_0, u_0.domainMax_0, u_0.rowBandPx_0, u_0.valueScaleType_0, u_0.valueSymlogConstant_0);\n}\n\n\nfloat lineColorYPx_0(float _S34, float _S35)\n{\n\n    return _S35 + (1.0 - normalizeScoreUnclamped_0(_S34, u_0.domainMin_0, u_0.domainMax_0, u_0.valueScaleType_0, u_0.valueSymlogConstant_0)) * u_0.rowBandPx_0;\n}\n\n\nvec4 lineCutYsPx_0(int _S36, float _S37)\n{\n\n    vec4 ys_0 = vec4(0.0);\n\n    int k_0 = 0;\n    for(;;)\n    {\n\n        if(k_0 < 4)\n        {\n        }\n        else\n        {\n\n            break;\n        }\n\n        int _S38 = _S36 + k_0;\n\n        if(_S38 < (u_0.colorCutCount_0))\n        {\n\n            ys_0[k_0] = lineColorYPx_0(colorCutAt_0(_S38, u_0.colorCuts_0[0], u_0.colorCuts_0[1]), _S37);\n\n        }\n\n        k_0 = k_0 + 1;\n\n    }\n\n    return ys_0;\n}\n\n\nvoid main()\n{\n\n    VsOut_0 o_0;\n\n    RowPlacement_0 _S39 = 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    if(_S39.hidden_0)\n    {\n\n        o_0.position_0 = vec4(2.0, 2.0, 0.0, 1.0);\n        const vec4 _S40 = vec4(0.0);\n\n        o_0.color_0 = _S40;\n        const vec2 _S41 = vec2(0.0);\n\n        o_0.localPos_0 = _S41;\n        o_0.segHalfLen_0 = 0.0;\n        o_0.colorLine_0 = _S41;\n        o_0.cutYsLo_0 = _S40;\n        o_0.cutYsHi_0 = _S40;\n        VsOut_0 _S42 = o_0;\n\n        gl_Position = o_0.position_0;\n\n        v_color = _S42.color_0;\n\n        v_localPos = _S42.localPos_0;\n\n        v_segHalfLen = _S42.segHalfLen_0;\n\n        v_colorLine = _S42.colorLine_0;\n\n        v_cutYsLo = _S42.cutYsLo_0;\n\n        v_cutYsHi = _S42.cutYsHi_0;\n\n        return;\n    }\n\n\n    float bandTop_0 = rowBandTopPx_0(u_0.rowOffsetPx_0, u_0.rowHeight_0, _S39.slot_1);\n    bool hasPrev_0 = a_prevX != 4294967295U;\n    float midCurrX_0 = (lineBpToClipX_0(a_x) + lineBpToClipX_0(a_x2)) * 0.5;\n\n    float midPrevX_0;\n\n\n\n    if(hasPrev_0)\n    {\n\n        midPrevX_0 = (lineBpToClipX_0(a_prevX) + lineBpToClipX_0(a_prevX2)) * 0.5;\n\n    }\n    else\n    {\n\n        midPrevX_0 = midCurrX_0;\n\n    }\n\n    float _S43 = lineValueYPx_0(a_y, bandTop_0);\n\n    float prevYPx_0;\n\n\n    if(hasPrev_0)\n    {\n\n        prevYPx_0 = lineValueYPx_0(a_prevY, bandTop_0);\n\n    }\n    else\n    {\n\n        prevYPx_0 = _S43;\n\n    }\n\n    vec2 deltaPx_0 = vec2((midCurrX_0 - midPrevX_0) * u_0.viewportWidth_0 * 0.5, _S43 - prevYPx_0);\n\n    float len_1 = length(deltaPx_0);\n    CapsuleFrame_0 f_1 = capsuleFrame_0(deltaPx_0, len_1);\n\n    float _S44 = len_1 * 0.5;\n\n    vec2 local_0 = capsuleQuadLocal_0(uint(gl_VertexID), _S44, u_0.lineWidth_0 * 0.5, u_0.devicePixelRatio_0);\n    vec2 off_0 = f_1.tangent_0 * local_0.x + f_1.normal_0 * local_0.y;\n    o_0.position_0 = vec4((midPrevX_0 + midCurrX_0) * 0.5 + pxToClipLen_0(off_0.x, u_0.viewportWidth_0), yPxToClipY_0((prevYPx_0 + _S43) * 0.5 + off_0.y, u_0.canvasHeight_0), 0.0, 1.0);\n\n    o_0.color_0 = _S39.color_1;\n    o_0.localPos_0 = local_0;\n    o_0.segHalfLen_0 = _S44;\n\n    float _S45 = lineColorYPx_0(a_y, bandTop_0);\n\n    float prevColorYPx_0;\n\n    if(hasPrev_0)\n    {\n\n        prevColorYPx_0 = lineColorYPx_0(a_prevY, bandTop_0);\n\n    }\n    else\n    {\n\n        prevColorYPx_0 = _S45;\n\n    }\n\n    float _S46 = (prevColorYPx_0 + _S45) * 0.5;\n    if(len_1 > 0.00009999999747379)\n    {\n\n        midPrevX_0 = (_S45 - prevColorYPx_0) / len_1;\n\n    }\n    else\n    {\n\n        midPrevX_0 = 0.0;\n\n    }\n\n    o_0.colorLine_0 = vec2(_S46, midPrevX_0);\n\n\n    o_0.cutYsLo_0 = lineCutYsPx_0(0, bandTop_0);\n    o_0.cutYsHi_0 = lineCutYsPx_0(4, bandTop_0);\n    VsOut_0 _S47 = o_0;\n\n    gl_Position = o_0.position_0;\n\n    v_color = _S47.color_0;\n\n    v_localPos = _S47.localPos_0;\n\n    v_segHalfLen = _S47.segHalfLen_0;\n\n    v_colorLine = _S47.colorLine_0;\n\n    v_cutYsLo = _S47.cutYsLo_0;\n\n    v_cutYsHi = _S47.cutYsHi_0;\n\n    return;\n}\n\n"
+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 domainMin_0;\n    float domainMax_0;\n    int valueScaleType_0;\n    float valueSymlogConstant_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    int colorFromY_0;\n    float rowHeight_0;\n    float rowBandPx_0;\n    float rowOffsetPx_0;\n    float origin_0;\n    float lineWidth_0;\n    float zero_0;\n    float viewportWidth_0;\n    float devicePixelRatio_0;\n    int rowTableKeys_0;\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_1, float domainMax_1, 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_1, symlogConstant_0);\n\n        tMax_0 = symlogTransform_0(domainMax_1, 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_1 > 0.0)\n            {\n\n                tMax_0 = domainMin_1;\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_1, tMax_0));\n\n            tMin_0 = _S7;\n\n            t_1 = _S6;\n\n        }\n        else\n        {\n\n            tMax_0 = domainMax_1;\n\n            tMin_0 = domainMin_1;\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_2, float domainMax_2, int scaleType_1, float symlogConstant_1)\n{\n\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_2, domainMax_2, 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 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 rowBandTopPx_0(float rowOffsetPx_1, float rowHeight_1, uint slot_0)\n{\n\n    return rowOffsetPx_1 + rowHeight_1 * float(slot_0);\n}\n\n\nfloat valueToYPxScaled_0(float value_2, float domainMin_3, float domainMax_3, float h_0, int scaleType_2, float symlogConstant_2)\n{\n\n    return (1.0 - normalizeScore_0(value_2, domainMin_3, domainMax_3, scaleType_2, symlogConstant_2)) * h_0;\n}\n\n\nfloat aaPx_0(float dpr_0)\n{\n\n    return 1.0 / max(dpr_0, 1.0);\n}\n\n\nfloat aaHalfPx_0(float dpr_1)\n{\n\n    return 0.5 * aaPx_0(dpr_1);\n}\n\n\nstruct CapsuleFrame_0\n{\n    vec2 tangent_0;\n    vec2 normal_0;\n};\n\n\nCapsuleFrame_0 capsuleFrame_0(vec2 delta_0, float len_0)\n{\n\n    CapsuleFrame_0 f_0;\n\n    vec2 _S11;\n    if(len_0 > 0.00009999999747379)\n    {\n\n        _S11 = delta_0 / len_0;\n\n    }\n    else\n    {\n\n        _S11 = vec2(1.0, 0.0);\n\n    }\n\n    f_0.tangent_0 = _S11;\n    f_0.normal_0 = vec2(- _S11.y, _S11.x);\n    return f_0;\n}\n\n\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S12 = v_0 == 0U;\n\n    bool _S13;\n\n    if(_S12)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 2U;\n\n    }\n\n    if(_S13)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 3U;\n\n    }\n\n    float _S14;\n\n    if(_S13)\n    {\n\n        _S14 = 0.0;\n\n    }\n    else\n    {\n\n        _S14 = 1.0;\n\n    }\n    if(_S12)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 1U;\n\n    }\n\n    if(_S13)\n    {\n\n        _S13 = true;\n\n    }\n    else\n    {\n\n        _S13 = v_0 == 4U;\n\n    }\n\n    float _S15;\n\n    if(_S13)\n    {\n\n        _S15 = 0.0;\n\n    }\n    else\n    {\n\n        _S15 = 1.0;\n\n    }\n\n    return vec2(_S14, _S15);\n}\n\n\nstruct SegmentCorner_0\n{\n    bool atEnd_0;\n    vec2 offsetPx_0;\n    vec3 local_0;\n};\n\n\nSegmentCorner_0 segmentCorner_0(uint vid_1, vec2 deltaPx_0, float alongPadPx_0, float acrossPadPx_0)\n{\n\n    float len_1 = length(deltaPx_0);\n    CapsuleFrame_0 f_1 = capsuleFrame_0(deltaPx_0, len_1);\n    vec2 q_0 = quadLocal_0(vid_1) * 2.0 - 1.0;\n    float across_0 = q_0.y * acrossPadPx_0;\n    SegmentCorner_0 c_2;\n    float _S16 = q_0.x;\n\n    bool _S17 = _S16 > 0.0;\n\n    c_2.atEnd_0 = _S17;\n    c_2.offsetPx_0 = f_1.tangent_0 * (_S16 * alongPadPx_0) + f_1.normal_0 * across_0;\n    float far_0 = len_1 + alongPadPx_0;\n\n    vec3 _S18;\n\n    if(_S17)\n    {\n\n        _S18 = vec3(far_0, - alongPadPx_0, across_0);\n\n    }\n    else\n    {\n\n        _S18 = vec3(- alongPadPx_0, far_0, across_0);\n\n    }\n\n    c_2.local_0 = _S18;\n\n\n    return c_2;\n}\n\n\n\nSegmentCorner_0 capsuleCorner_0(uint vid_2, vec2 deltaPx_1, float halfWidthPx_0, float dpr_2)\n{\n\n    float ext_0 = halfWidthPx_0 + aaHalfPx_0(dpr_2);\n    return segmentCorner_0(vid_2, deltaPx_1, ext_0, ext_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 yPxToClipY_0(float yPx_0, float canvasHeight_1)\n{\n\n    return 1.0 - yPx_0 / canvasHeight_1 * 2.0;\n}\n\n\nbool lineColorAlongY_0(int rampMode_1, int colorFromY_1, bool overridden_0)\n{\n\n    bool _S19;\n    if(rampMode_1 == 3)\n    {\n\n        _S19 = colorFromY_1 != 0;\n\n    }\n    else\n    {\n\n        _S19 = false;\n\n    }\n\n    if(_S19)\n    {\n\n        _S19 = !overridden_0;\n\n    }\n    else\n    {\n\n        _S19 = false;\n\n    }\n\n    return _S19;\n}\n\n\nflat out vec4 v_color;\n\n\nout vec3 v_localPos;\n\n\nflat out vec2 v_colorLine;\n\n\nflat out vec4 v_cutYsLo;\n\n\nflat out vec4 v_cutYsHi;\n\n\nflat out float v_colorAlongY;\n\n\nlayout(location = 0)\nin uint a_x;\n\n\nlayout(location = 1)\nin uint a_x2;\n\n\nlayout(location = 2)\nin float a_y;\n\n\nlayout(location = 3)\nin uint a_color;\n\n\nlayout(location = 4)\nin uint a_row;\n\n\nlayout(location = 5)\nin uint a_prevX;\n\n\nlayout(location = 6)\nin uint a_prevX2;\n\n\nlayout(location = 7)\nin float a_prevY;\n\n\nstruct VsOut_0\n{\n    vec4 position_0;\n    vec4 color_0;\n    vec3 localPos_0;\n    vec2 colorLine_0;\n    vec4 cutYsLo_0;\n    vec4 cutYsHi_0;\n    float colorAlongY_0;\n};\n\n\nvec4 rampColor_0(float _S20)\n{\n\n    return (textureLod((u_colorRamp), (vec2((clamp(_S20, 0.0, 1.0) * 255.0 + 0.5) / 256.0, 0.5)), (0.0)));\n}\n\n\nvec4 markScaleColor_0(uint _S21, int _S22, float _S23, float _S24, float _S25, int _S26, vec4 _S27, vec4 _S28)\n{\n\n    if(_S22 == 0)\n    {\n\n        return unpackRGBA_0(_S21);\n    }\n    if(_S21 == 2143289345U)\n    {\n\n        return unpackRGBA_0(4289703855U);\n    }\n\n    int scaleType_3;\n\n    if((_S21 & 2139095040U) == 2139095040U)\n    {\n\n        if((_S21 & 8388607U) != 0U)\n        {\n\n            return unpackRGBA_0(4286611584U);\n        }\n\n        bool negative_0 = (_S21 & 2147483648U) != 0U;\n        if(_S22 == 3)\n        {\n\n            if(negative_0)\n            {\n\n                scaleType_3 = 0;\n\n            }\n            else\n            {\n\n                scaleType_3 = _S26;\n\n            }\n\n            return u_0.colorBands_0[scaleType_3];\n        }\n\n        float _S29;\n\n        if(negative_0)\n        {\n\n            _S29 = 0.0;\n\n        }\n        else\n        {\n\n            _S29 = 1.0;\n\n        }\n\n        return rampColor_0(rampMidT_0(_S29, _S25));\n    }\n    float value_3 = uintBitsToFloat(_S21);\n    if(_S22 == 3)\n    {\n\n        return u_0.colorBands_0[thresholdBand_0(value_3, _S26, _S27, _S28)];\n    }\n    if(_S22 == 2)\n    {\n\n        scaleType_3 = 1;\n\n    }\n    else\n    {\n\n        scaleType_3 = 0;\n\n    }\n    return rampColor_0(rampMidT_0(normalizeScore_0(value_3, _S23, _S24, scaleType_3, 1.0), _S25));\n}\n\n\nstruct RowPlacement_0\n{\n    bool hidden_0;\n    uint slot_1;\n    vec4 color_1;\n    bool overridden_1;\n};\n\n\nRowPlacement_0 rowTableLookup_0(uint _S30, int _S31, vec4 _S32)\n{\n\n    RowPlacement_0 p_0;\n    p_0.hidden_0 = false;\n    p_0.slot_1 = _S30;\n    p_0.color_1 = _S32;\n    p_0.overridden_1 = false;\n    if(_S31 < 0)\n    {\n\n        return p_0;\n    }\n    uint n_0 = uint(_S31);\n    if(_S30 >= n_0)\n    {\n\n        p_0.hidden_0 = true;\n        return p_0;\n    }\n    uint h_1 = rowTablePlaneHeight_0(n_0);\n    uint _S33 = rowTableTexelX_0(_S30, n_0);\n\n    float u_2 = (float(_S33) + 0.5) / float(rowTableWidth_0(n_0));\n    float rows_0 = float(2U * h_1);\n    uint y_0 = rowTableTexelY_0(_S30, n_0);\n    vec4 s_0 = (textureLod((u_rowTable), (vec2(u_2, (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_1 = ((byteOf_0(s_0.x)) | ((byteOf_0(s_0.y)) << 8)) | ((byteOf_0(s_0.z)) << 16);\n    vec4 c_3 = (textureLod((u_rowTable), (vec2(u_2, (float(y_0 + h_1) + 0.5) / rows_0)), (0.0)));\n    bool _S34 = (c_3.w) > 0.0;\n\n    p_0.overridden_1 = _S34;\n\n    vec4 _S35;\n    if(_S34)\n    {\n\n        _S35 = c_3;\n\n    }\n    else\n    {\n\n        _S35 = _S32;\n\n    }\n\n    p_0.color_1 = _S35;\n    return p_0;\n}\n\n\nfloat lineBpToClipX_0(uint _S36)\n{\n\n    return hpToClipX_0(hpSplitUint_0(_S36), u_0.bpRangeX_0, u_0.zero_0);\n}\n\n\nfloat lineValueYPx_0(float _S37, float _S38)\n{\n\n    return _S38 + valueToYPxScaled_0(_S37, u_0.domainMin_0, u_0.domainMax_0, u_0.rowBandPx_0, u_0.valueScaleType_0, u_0.valueSymlogConstant_0);\n}\n\n\nfloat lineColorYPx_0(float _S39, float _S40)\n{\n\n    return _S40 + (1.0 - normalizeScoreUnclamped_0(_S39, u_0.domainMin_0, u_0.domainMax_0, u_0.valueScaleType_0, u_0.valueSymlogConstant_0)) * u_0.rowBandPx_0;\n}\n\n\nvec4 lineCutYsPx_0(int _S41, float _S42)\n{\n\n    vec4 ys_0 = vec4(0.0);\n\n    int k_0 = 0;\n    for(;;)\n    {\n\n        if(k_0 < 4)\n        {\n        }\n        else\n        {\n\n            break;\n        }\n\n        int _S43 = _S41 + k_0;\n\n        if(_S43 < (u_0.colorCutCount_0))\n        {\n\n            ys_0[k_0] = lineColorYPx_0(colorCutAt_0(_S43, u_0.colorCuts_0[0], u_0.colorCuts_0[1]), _S42);\n\n        }\n\n        k_0 = k_0 + 1;\n\n    }\n\n    return ys_0;\n}\n\n\nvoid main()\n{\n\n    VsOut_0 o_0;\n\n    RowPlacement_0 _S44 = 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    if(_S44.hidden_0)\n    {\n\n        o_0.position_0 = vec4(2.0, 2.0, 0.0, 1.0);\n        const vec4 _S45 = vec4(0.0);\n\n        o_0.color_0 = _S45;\n        o_0.localPos_0 = vec3(0.0);\n        o_0.colorLine_0 = vec2(0.0);\n        o_0.cutYsLo_0 = _S45;\n        o_0.cutYsHi_0 = _S45;\n        o_0.colorAlongY_0 = 0.0;\n        VsOut_0 _S46 = o_0;\n\n        gl_Position = o_0.position_0;\n\n        v_color = _S46.color_0;\n\n        v_localPos = _S46.localPos_0;\n\n        v_colorLine = _S46.colorLine_0;\n\n        v_cutYsLo = _S46.cutYsLo_0;\n\n        v_cutYsHi = _S46.cutYsHi_0;\n\n        v_colorAlongY = _S46.colorAlongY_0;\n\n        return;\n    }\n\n\n    float bandTop_0 = rowBandTopPx_0(u_0.rowOffsetPx_0, u_0.rowHeight_0, _S44.slot_1);\n    bool hasPrev_0 = a_prevX != 4294967295U;\n    float midCurrX_0 = (lineBpToClipX_0(a_x) + lineBpToClipX_0(a_x2)) * 0.5;\n\n    float midPrevX_0;\n\n\n\n    if(hasPrev_0)\n    {\n\n        midPrevX_0 = (lineBpToClipX_0(a_prevX) + lineBpToClipX_0(a_prevX2)) * 0.5;\n\n    }\n    else\n    {\n\n        midPrevX_0 = midCurrX_0;\n\n    }\n\n    float _S47 = lineValueYPx_0(a_y, bandTop_0);\n\n    float prevYPx_0;\n\n\n    if(hasPrev_0)\n    {\n\n        prevYPx_0 = lineValueYPx_0(a_prevY, bandTop_0);\n\n    }\n    else\n    {\n\n        prevYPx_0 = _S47;\n\n    }\n\n\n    vec2 deltaPx_2 = vec2((midCurrX_0 - midPrevX_0) * u_0.viewportWidth_0 * 0.5, _S47 - prevYPx_0);\n\n    float len_2 = length(deltaPx_2);\n\n    SegmentCorner_0 c_4 = capsuleCorner_0(uint(gl_VertexID), deltaPx_2, u_0.lineWidth_0 * 0.5, u_0.devicePixelRatio_0);\n\n    if(c_4.atEnd_0)\n    {\n\n        midPrevX_0 = midCurrX_0;\n\n    }\n\n    float _S48 = midPrevX_0 + pxToClipLen_0(c_4.offsetPx_0.x, u_0.viewportWidth_0);\n    if(c_4.atEnd_0)\n    {\n\n        midPrevX_0 = _S47;\n\n    }\n    else\n    {\n\n        midPrevX_0 = prevYPx_0;\n\n    }\n\n    o_0.position_0 = vec4(_S48, yPxToClipY_0(midPrevX_0 + c_4.offsetPx_0.y, u_0.canvasHeight_0), 0.0, 1.0);\n\n    o_0.color_0 = _S44.color_1;\n    o_0.localPos_0 = c_4.local_0;\n\n    float _S49 = lineColorYPx_0(a_y, bandTop_0);\n\n    float prevColorYPx_0;\n\n    if(hasPrev_0)\n    {\n\n        prevColorYPx_0 = lineColorYPx_0(a_prevY, bandTop_0);\n\n    }\n    else\n    {\n\n        prevColorYPx_0 = _S49;\n\n    }\n\n\n    if(len_2 > 0.00009999999747379)\n    {\n\n        midPrevX_0 = (_S49 - prevColorYPx_0) / len_2;\n\n    }\n    else\n    {\n\n        midPrevX_0 = 0.0;\n\n    }\n\n    o_0.colorLine_0 = vec2(prevColorYPx_0, midPrevX_0);\n\n\n    o_0.cutYsLo_0 = lineCutYsPx_0(0, bandTop_0);\n    o_0.cutYsHi_0 = lineCutYsPx_0(4, bandTop_0);\n\n    if(lineColorAlongY_0(u_0.rampMode_0, u_0.colorFromY_0, _S44.overridden_1))\n    {\n\n        midPrevX_0 = 1.0;\n\n    }\n    else\n    {\n\n        midPrevX_0 = 0.0;\n\n    }\n\n    o_0.colorAlongY_0 = midPrevX_0;\n\n    VsOut_0 _S50 = o_0;\n\n    gl_Position = o_0.position_0;\n\n    v_color = _S50.color_0;\n\n    v_localPos = _S50.localPos_0;\n\n    v_colorLine = _S50.colorLine_0;\n\n    v_cutYsLo = _S50.cutYsLo_0;\n\n    v_cutYsHi = _S50.cutYsHi_0;\n\n    v_colorAlongY = _S50.colorAlongY_0;\n\n    return;\n}\n\n"
 
-export const GLSL_FRAGMENT = "#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 domainMin_0;\n    float domainMax_0;\n    int valueScaleType_0;\n    float valueSymlogConstant_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    float rowHeight_0;\n    float rowBandPx_0;\n    float rowOffsetPx_0;\n    float origin_0;\n    float lineWidth_0;\n    float zero_0;\n    float viewportWidth_0;\n    float devicePixelRatio_0;\n    int rowTableKeys_0;\n}u_0;\n\nfloat capsuleDistPx_0(float alongPx_0, float acrossPx_0, float halfLenPx_0)\n{\n\n    float _S1 = max(abs(alongPx_0) - halfLenPx_0, 0.0);\n    return sqrt(_S1 * _S1 + acrossPx_0 * acrossPx_0);\n}\n\nfloat capsuleDist_0(vec2 local_0, float halfLenPx_1)\n{\n\n    return capsuleDistPx_0(local_0.x, local_0.y, halfLenPx_1);\n}\n\n\nfloat aaPx_0(float dpr_0)\n{\n\n    return 1.0 / max(dpr_0, 1.0);\n}\n\n\nfloat aaRamp_0(float signedInk_0, float widthPx_0)\n{\n\n    return clamp(signedInk_0 / widthPx_0 + 0.5, 0.0, 1.0);\n}\n\n\nfloat edgeCoverage_0(float signedInkCssPx_0, float dpr_1)\n{\n\n    return aaRamp_0(signedInkCssPx_0, aaPx_0(dpr_1));\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 thresholdBandAtYPx_0(float yPx_0, int cutCount_0, vec4 cutYsLo_0, vec4 cutYsHi_0)\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(yPx_0 <= (colorCutAt_0(i_1, cutYsLo_0, cutYsHi_0)))\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\nlayout(location = 0)\nout vec4 entryPointParam_fs_main_0;\n\n\nflat in vec4 v_color;\n\n\nin vec2 v_localPos;\n\n\nflat in float v_segHalfLen;\n\n\nflat in vec2 v_colorLine;\n\n\nflat in vec4 v_cutYsLo;\n\n\nflat in vec4 v_cutYsHi;\n\n\nvec4 lineFragmentColor_0(float _S3, vec4 _S4, vec4 _S5, vec4 _S6)\n{\n\n    if((u_0.rampMode_0) != 3)\n    {\n\n        return _S6;\n    }\n    return u_0.colorBands_0[thresholdBandAtYPx_0(_S3, u_0.colorCutCount_0, _S4, _S5)];\n}\n\n\nvoid main()\n{\n    float alpha_0 = edgeCoverage_0(u_0.lineWidth_0 * 0.5 - capsuleDist_0(v_localPos, v_segHalfLen), u_0.devicePixelRatio_0);\n    if(alpha_0 <= 0.0)\n    {\n\n        discard;\n\n    }\n\n    vec4 _S7 = lineFragmentColor_0(v_colorLine.x + v_colorLine.y * clamp(v_localPos.x, - v_segHalfLen, v_segHalfLen), v_cutYsLo, v_cutYsHi, v_color);\n\n    entryPointParam_fs_main_0 = vec4(_S7.xyz * alpha_0, _S7.w * alpha_0);\n\n    return;\n}\n\n"
+export const GLSL_FRAGMENT = "#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 domainMin_0;\n    float domainMax_0;\n    int valueScaleType_0;\n    float valueSymlogConstant_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    int colorFromY_0;\n    float rowHeight_0;\n    float rowBandPx_0;\n    float rowOffsetPx_0;\n    float origin_0;\n    float lineWidth_0;\n    float zero_0;\n    float viewportWidth_0;\n    float devicePixelRatio_0;\n    int rowTableKeys_0;\n}u_0;\n\nfloat capsuleDistPx_0(float fromStartPx_0, float fromEndPx_0, float acrossPx_0)\n{\n\n    float _S1 = max(max(- fromStartPx_0, - fromEndPx_0), 0.0);\n    return sqrt(_S1 * _S1 + acrossPx_0 * acrossPx_0);\n}\n\nfloat capsuleDist_0(vec3 local_0)\n{\n\n    return capsuleDistPx_0(local_0.x, local_0.y, local_0.z);\n}\n\n\nfloat aaPx_0(float dpr_0)\n{\n\n    return 1.0 / max(dpr_0, 1.0);\n}\n\n\nfloat aaRamp_0(float signedInk_0, float widthPx_0)\n{\n\n    return clamp(signedInk_0 / widthPx_0 + 0.5, 0.0, 1.0);\n}\n\n\nfloat edgeCoverage_0(float signedInkCssPx_0, float dpr_1)\n{\n\n    return aaRamp_0(signedInkCssPx_0, aaPx_0(dpr_1));\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 thresholdBandAtYPx_0(float yPx_0, int cutCount_0, vec4 cutYsLo_0, vec4 cutYsHi_0)\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(yPx_0 <= (colorCutAt_0(i_1, cutYsLo_0, cutYsHi_0)))\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\nlayout(location = 0)\nout vec4 entryPointParam_fs_main_0;\n\n\nflat in vec4 v_color;\n\n\nin vec3 v_localPos;\n\n\nflat in vec2 v_colorLine;\n\n\nflat in vec4 v_cutYsLo;\n\n\nflat in vec4 v_cutYsHi;\n\n\nflat in float v_colorAlongY;\n\n\nvec4 lineFragmentColor_0(float _S3, vec4 _S4, vec4 _S5, vec4 _S6, float _S7)\n{\n\n    if(_S7 < 0.5)\n    {\n\n        return _S6;\n    }\n    return u_0.colorBands_0[thresholdBandAtYPx_0(_S3, u_0.colorCutCount_0, _S4, _S5)];\n}\n\n\nvoid main()\n{\n    float alpha_0 = edgeCoverage_0(u_0.lineWidth_0 * 0.5 - capsuleDist_0(v_localPos), u_0.devicePixelRatio_0);\n    if(alpha_0 <= 0.0)\n    {\n\n        discard;\n\n    }\n\n    vec4 _S8 = lineFragmentColor_0(v_colorLine.x + v_colorLine.y * clamp(v_localPos.x, 0.0, v_localPos.x + v_localPos.y), v_cutYsLo, v_cutYsHi, v_color, v_colorAlongY);\n\n    entryPointParam_fs_main_0 = vec4(_S8.xyz * alpha_0, _S8.w * alpha_0);\n\n    return;\n}\n\n"

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

@@ -29,14 +29,14 @@ export const UNIFORM_OFFSET_F32 = {
   rampMin: 9,
   rampMax: 10,
   rampMidNorm: 11,
-  rowHeight: 60,
-  rowBandPx: 61,
-  rowOffsetPx: 62,
-  origin: 63,
-  lineWidth: 64,
-  zero: 65,
-  viewportWidth: 66,
-  devicePixelRatio: 67,
+  rowHeight: 61,
+  rowBandPx: 62,
+  rowOffsetPx: 63,
+  origin: 64,
+  lineWidth: 65,
+  zero: 66,
+  viewportWidth: 67,
+  devicePixelRatio: 68,
 } as const
 
 // Word indices into a Int32Array view over the uniform buffer.
@@ -44,7 +44,8 @@ export const UNIFORM_OFFSET_I32 = {
   valueScaleType: 6,
   rampMode: 8,
   colorCutCount: 12,
-  rowTableKeys: 68,
+  colorFromY: 60,
+  rowTableKeys: 69,
 } as const
 
 
@@ -104,6 +105,7 @@ export interface Uniforms {
   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]]
+  colorFromY: number
   rowHeight: number
   rowBandPx: number
   rowOffsetPx: number
@@ -175,15 +177,16 @@ export function writeUniforms(buf: ArrayBuffer, uniforms: Uniforms) {
   f32[57] = uniforms.colorBands[8][1]
   f32[58] = uniforms.colorBands[8][2]
   f32[59] = uniforms.colorBands[8][3]
-  f32[60] = uniforms.rowHeight
-  f32[61] = uniforms.rowBandPx
-  f32[62] = uniforms.rowOffsetPx
-  f32[63] = uniforms.origin
-  f32[64] = uniforms.lineWidth
-  f32[65] = uniforms.zero
-  f32[66] = uniforms.viewportWidth
-  f32[67] = uniforms.devicePixelRatio
-  i32[68] = uniforms.rowTableKeys
+  i32[60] = uniforms.colorFromY
+  f32[61] = uniforms.rowHeight
+  f32[62] = uniforms.rowBandPx
+  f32[63] = uniforms.rowOffsetPx
+  f32[64] = uniforms.origin
+  f32[65] = uniforms.lineWidth
+  f32[66] = uniforms.zero
+  f32[67] = uniforms.viewportWidth
+  f32[68] = uniforms.devicePixelRatio
+  i32[69] = uniforms.rowTableKeys
 }
 
 export const INSTANCE_STRIDE_BYTES = 32