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

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

@@ -0,0 +1,285 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from barMark.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', 'fragment'] },
+  { 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 = 12
+
+export const UNIFORMS_SIZE_BYTES = 288
+
+// Word indices into a Float32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_F32 = {
+  bpRangeX: 0,
+  canvasHeight: 3,
+  domainMin: 4,
+  domainMax: 5,
+  valueSymlogConstant: 7,
+  rampMin: 9,
+  rampMax: 10,
+  rampMidNorm: 11,
+  originYPx: 60,
+  rowHeight: 61,
+  rowBandPx: 62,
+  rowOffsetPx: 63,
+  zero: 64,
+  minCellDenomPx: 65,
+  minWidthPx: 66,
+  devicePixelRatio: 67,
+} as const
+
+// Word indices into a Int32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_I32 = {
+  valueScaleType: 6,
+  rampMode: 8,
+  colorCutCount: 12,
+  rowTableKeys: 68,
+} 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
+  domainMin: number
+  domainMax: number
+  valueScaleType: number
+  valueSymlogConstant: 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]]
+  originYPx: number
+  rowHeight: number
+  rowBandPx: number
+  rowOffsetPx: number
+  zero: number
+  minCellDenomPx: number
+  minWidthPx: number
+  devicePixelRatio: number
+  rowTableKeys: 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.domainMin
+  f32[5] = uniforms.domainMax
+  i32[6] = uniforms.valueScaleType
+  f32[7] = uniforms.valueSymlogConstant
+  i32[8] = uniforms.rampMode
+  f32[9] = uniforms.rampMin
+  f32[10] = uniforms.rampMax
+  f32[11] = uniforms.rampMidNorm
+  i32[12] = 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]
+  f32[60] = uniforms.originYPx
+  f32[61] = uniforms.rowHeight
+  f32[62] = uniforms.rowBandPx
+  f32[63] = uniforms.rowOffsetPx
+  f32[64] = uniforms.zero
+  f32[65] = uniforms.minCellDenomPx
+  f32[66] = uniforms.minWidthPx
+  f32[67] = uniforms.devicePixelRatio
+  i32[68] = uniforms.rowTableKeys
+}
+
+export const INSTANCE_STRIDE_BYTES = 20
+export const INSTANCE_STRIDE_WORDS = 5
+
+// Word indices into a Float32Array view over the instance buffer.
+export const INSTANCE_OFFSET_F32 = {
+  y: 2,
+} as const
+
+// Word indices into a Uint32Array view over the instance buffer.
+export const INSTANCE_OFFSET_U32 = {
+  x: 0,
+  x2: 1,
+  color: 3,
+  row: 4,
+} 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_y', components: 1, type: 'float', offsetBytes: 8, integer: false },
+  { name: 'a_color', components: 1, type: 'uint', offsetBytes: 12, integer: true },
+  { name: 'a_row', components: 1, type: 'uint', offsetBytes: 16, integer: true },
+]
+
+export interface InstanceArrays {
+  x: ArrayLike<number>
+  x2: ArrayLike<number>
+  y: ArrayLike<number>
+  color: ArrayLike<number>
+  row: ArrayLike<number>
+}
+
+export function packInstances(
+  arrays: InstanceArrays,
+  numInstances: number,
+  buf: ArrayBuffer = new ArrayBuffer(numInstances * INSTANCE_STRIDE_BYTES),
+) {
+  const f32 = new Float32Array(buf)
+  const u32 = new Uint32Array(buf)
+  const { x, x2, y, color, row } = arrays
+  for (let i = 0; i < numInstances; i++) {
+    const o = i * INSTANCE_STRIDE_WORDS
+    u32[o + 0] = x[i]!
+    u32[o + 1] = x2[i]!
+    f32[o + 2] = y[i]!
+    u32[o + 3] = color[i]!
+    u32[o + 4] = row[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 `y`.
+export function getInstanceY(f32: Float32Array, i: number) {
+  return f32[i * INSTANCE_STRIDE_WORDS + 2]!
+}
+
+export function setInstanceY(f32: Float32Array, i: number, v: number) {
+  f32[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
+}
+
+// Instance `i`'s `row`.
+export function getInstanceRow(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS + 4]!
+}
+
+export function setInstanceRow(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS + 4] = v
+}
+
+// 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/barMark.slang

@@ -0,0 +1,167 @@
+// The `bar` mark shape: one instance is a rectangle spanning [x, x2) in
+// absolute genomic coords, standing between the `origin` value and the
+// instance's own `y` on the shared `valueScale` — the quantitative bar a BED
+// score column, a CNV segment ratio or a bedGraph-shaped interval draws as.
+//
+// The scale is `valueScale`'s and the anchor is this shape's own (ADR-097): the
+// point shape centres a glyph on the value, this one grows from the baseline to
+// it, and both read one `valueToYPxScaled`.
+//
+// The band is the instance's `row`: `rowHeight` px per row, the value scale
+// ruling each band on its own over `rowBandPx` of it, `rowOffsetPx` below the
+// row's top. A caller with one row sends the canvas height for both and no
+// offset, so its bars stand where they always did (ADR-126). The row is read
+// through the bound row table, as `spanMark` reads it; a pass binding none
+// (`rowTableKeys` −1) draws `row` as the band.
+//
+// The two HORIZONTAL cuts are antialiased analytically, as wiggle's xyplot bar
+// does: the top cut IS the datum, so it must not quantise to the output grid,
+// and the band is the difference of the two cuts' ramps so a bar thinner than a
+// device pixel covers exactly its area. The VERTICAL cuts stay hard on purpose
+// — intervals that tile (segments, bins) split a pixel column with their
+// neighbours, and per-fragment alpha on both sides of that shared edge
+// composites to less than full ink. That fractional edge is why this shader
+// declares no `//! coverage: analytic`: the multisampled target resolves it.
+//
+// The min-width floor grows off the START edge (`extendToMinWidthX`), which is
+// what the painter's `spanLeft` does too, so a reversed block's bars stay exact
+// mirrors of a forward block's.
+//! targets: wgsl, glsl
+
+import antialias;
+import hpmath;
+import colorPack;
+import clipStrip;
+import markColor;
+import rowTable;
+import valueScale;
+
+// Two quads: the bar, then its clip strip, which every vertex of a bar the
+// domain did not cut sends off clip.
+public static const uint VERTS_PER_INSTANCE = 12u;
+
+struct BarMarkInstance {
+  uint  x     : ATTR0;
+  uint  x2    : ATTR1;
+  float y     : ATTR2;
+  uint  color : ATTR3;
+  uint  row   : ATTR4;
+};
+
+struct Uniforms {
+  float3 bpRangeX;
+  float  canvasHeight;
+  float  domainMin;
+  float  domainMax;
+  // `scoreScale`'s scale-type vocabulary, for the declared y scale.
+  int    valueScaleType;
+  // Already resolved from the domain (`resolveSymlogConstant`); read by symlog alone.
+  float  valueSymlogConstant;
+  // `markColor`'s: 0 packs the instance's colour, 1 and 2 read its slot as a
+  // value through the ramp below.
+  int    rampMode;
+  float  rampMin;
+  float  rampMax;
+  float  rampMidNorm;
+  // `markColor`'s threshold: how many cuts are in use, the cuts four to a
+  // vector, and the colour of each interval between them.
+  int    colorCutCount;
+  float4 colorCuts[2];
+  float4 colorBands[MAX_COLOR_BANDS];
+  // The baseline, in CSS px below the band's top: `origin` through the same
+  // scale `y` reads, resolved on the CPU because the whole draw shares it.
+  float  originYPx;
+  // CSS px per row band; the canvas height where every instance is on row 0.
+  float  rowHeight;
+  // CSS px the value scale runs over inside a row, and how far below the
+  // row's top it starts: an inset, less a scroll.
+  float  rowBandPx;
+  float  rowOffsetPx;
+  float  zero;
+  // The CSS width `minWidthPx` is measured against; see rowRect's
+  // `minCellDenomPx` for why it has to be a CSS width.
+  float  minCellDenomPx;
+  float  minWidthPx;
+  float  devicePixelRatio;
+  int    rowTableKeys;
+};
+[[vk::binding(1, 0)]] ConstantBuffer<Uniforms> u;
+[[vk::binding(2, 0)]] Sampler2D<float4> colorRamp;
+[[vk::binding(4, 0)]] Sampler2D<float4> rowTable;
+
+struct VsOut {
+  float4 position : SV_Position;
+  float4 color    : COLOR0;
+  // CSS px below the bar's top cut (.x) and below its bottom cut (.y), both
+  // affine in screen y with unit slope, so the fragment's ramp is one output
+  // pixel wide by construction and takes no derivative.
+  float2 barInkPx : TEXCOORD0;
+};
+
+float bpToClipX(uint bp, Uniforms u) {
+  return hpToClipX(hpSplitUint(bp), u.bpRangeX, u.zero);
+}
+
+// One device pixel in CSS px: the ramp's full width and the geometry pad, so
+// every fragment the ramp gives ink to is fully covered by the quad.
+float barAaPx() {
+  return 2.0 * aaHalfPx(u.devicePixelRatio);
+}
+
+[shader("vertex")]
+VsOut vs_main(BarMarkInstance inst, uint vid : SV_VertexID) {
+  VsOut o;
+  RowPlacement place = 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));
+  if (place.hidden) {
+    o.position = float4(2.0, 2.0, 0.0, 1.0);
+    o.color = float4(0.0);
+    o.barInkPx = float2(0.0);
+    return o;
+  }
+  float sx1 = bpToClipX(inst.x, u);
+  float sx2 = bpToClipX(inst.x2, u);
+  sx2 = extendToMinWidthX(sx1, sx2, u.minWidthPx, u.minCellDenomPx);
+  float2 q = quadLocal(vid);
+
+  float bandTop = rowBandTopPx(u.rowOffsetPx, u.rowHeight, place.slot);
+  float padPx = barAaPx();
+
+  if (vid >= 6u) {
+    // The clip strip: CLIP_STRIP_PX across the band edge that cut the bar,
+    // antialiased at its two horizontal cuts as the bar is.
+    int side = clipSide(inst.y, u.domainMin, u.domainMax, u.valueScaleType);
+    if (side == 0) {
+      o.position = float4(2.0, 2.0, 0.0, 1.0);
+      o.color = float4(0.0);
+      o.barInkPx = float2(0.0);
+      return o;
+    }
+    float stripTop = side > 0 ? bandTop : bandTop + u.rowBandPx - CLIP_STRIP_PX;
+    float stripBot = stripTop + CLIP_STRIP_PX;
+    float sy = lerp(stripBot + padPx, stripTop - padPx, q.y);
+    o.position = float4(lerp(sx1, sx2, q.x), yPxToClipY(sy, u.canvasHeight), 0.0, 1.0);
+    o.color = float4(unpackRGBA(CLIP_STRIP_COLOR).rgb, 1.0);
+    o.barInkPx = float2(sy - stripTop, sy - stripBot);
+    return o;
+  }
+
+  float valueY = bandTop + valueToYPxScaled(inst.y, u.domainMin, u.domainMax, u.rowBandPx, u.valueScaleType, u.valueSymlogConstant);
+  float originY = bandTop + u.originYPx;
+  float topPx = min(valueY, originY);
+  float botPx = max(valueY, originY);
+  float vy = lerp(botPx + padPx, topPx - padPx, q.y);
+
+  o.position = float4(lerp(sx1, sx2, q.x), yPxToClipY(vy, u.canvasHeight), 0.0, 1.0);
+  o.color = place.color;
+  o.barInkPx = float2(vy - topPx, vy - botPx);
+  return o;
+}
+
+[shader("fragment")]
+float4 fs_main(VsOut v) : SV_Target {
+  float aaPx = barAaPx();
+  float alpha = aaRamp(v.barInkPx.x, aaPx) - aaRamp(v.barInkPx.y, aaPx);
+  if (alpha <= 0.0) discard;
+  return float4(v.color.rgb, v.color.a * alpha);
+}

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

@@ -0,0 +1,4 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from barMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const WGSL_SOURCE = "struct Uniforms_std140_0\n{\n    @align(16) bpRangeX_0 : vec3<f32>,\n    @align(4) canvasHeight_0 : f32,\n    @align(16) domainMin_0 : f32,\n    @align(4) domainMax_0 : f32,\n    @align(8) valueScaleType_0 : i32,\n    @align(4) valueSymlogConstant_0 : f32,\n    @align(16) rampMode_0 : i32,\n    @align(4) rampMin_0 : f32,\n    @align(8) rampMax_0 : f32,\n    @align(4) rampMidNorm_0 : f32,\n    @align(16) colorCutCount_0 : i32,\n    @align(16) colorCuts_0 : array<vec4<f32>, i32(2)>,\n    @align(16) colorBands_0 : array<vec4<f32>, i32(9)>,\n    @align(16) originYPx_0 : f32,\n    @align(4) rowHeight_0 : f32,\n    @align(8) rowBandPx_0 : f32,\n    @align(4) rowOffsetPx_0 : f32,\n    @align(16) zero_0 : f32,\n    @align(4) minCellDenomPx_0 : f32,\n    @align(8) minWidthPx_0 : f32,\n    @align(4) devicePixelRatio_0 : f32,\n    @align(16) rowTableKeys_0 : i32,\n};\n\n@binding(1) @group(0) var<uniform> u_0 : Uniforms_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_1 : f32,  domainMax_1 : 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_1, symlogConstant_0);\n        tMax_0 = symlogTransform_0(domainMax_1, symlogConstant_0);\n        tMin_0 = _S5;\n        t_2 = _S4;\n    }\n    else\n    {\n        if(scaleType_0 == i32(1))\n        {\n            if(domainMin_1 > 0.0f)\n            {\n                tMax_0 = domainMin_1;\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_1, tMax_0));\n            tMin_0 = _S7;\n            t_2 = _S6;\n        }\n        else\n        {\n            tMax_0 = domainMax_1;\n            tMin_0 = domainMin_1;\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_2 : f32,  domainMax_2 : f32,  scaleType_1 : i32,  symlogConstant_1 : f32) -> f32\n{\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_2, domainMax_2, 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 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 _S15 : f32;\n    if((abs(dx_0)) < minWidth_0)\n    {\n        if(dx_0 < 0.0f)\n        {\n            _S15 = - minWidth_0;\n        }\n        else\n        {\n            _S15 = minWidth_0;\n        }\n        _S15 = x1_0 + _S15;\n    }\n    else\n    {\n        _S15 = x2_0;\n    }\n    return _S15;\n}\n\nfn extendToMinWidthX_0( x1_1 : f32,  x2_1 : f32,  minWidthPx_1 : f32,  canvasWidth_0 : f32) -> f32\n{\n    return extendToMinWidthPx_0(x1_1, x2_1, pxToClipLen_0(minWidthPx_1, canvasWidth_0));\n}\n\nfn quadLocal_0( vid_0 : u32) -> vec2<f32>\n{\n    var v_0 : u32 = vid_0 % u32(6);\n    var _S16 : bool = v_0 == u32(0);\n    var _S17 : bool;\n    if(_S16)\n    {\n        _S17 = true;\n    }\n    else\n    {\n        _S17 = v_0 == u32(2);\n    }\n    if(_S17)\n    {\n        _S17 = true;\n    }\n    else\n    {\n        _S17 = v_0 == u32(3);\n    }\n    var _S18 : f32;\n    if(_S17)\n    {\n        _S18 = 0.0f;\n    }\n    else\n    {\n        _S18 = 1.0f;\n    }\n    if(_S16)\n    {\n        _S17 = true;\n    }\n    else\n    {\n        _S17 = v_0 == u32(1);\n    }\n    if(_S17)\n    {\n        _S17 = true;\n    }\n    else\n    {\n        _S17 = v_0 == u32(4);\n    }\n    var _S19 : f32;\n    if(_S17)\n    {\n        _S19 = 0.0f;\n    }\n    else\n    {\n        _S19 = 1.0f;\n    }\n    return vec2<f32>(_S18, _S19);\n}\n\nfn rowBandTopPx_0( rowOffsetPx_1 : f32,  rowHeight_1 : f32,  slot_1 : u32) -> f32\n{\n    return rowOffsetPx_1 + rowHeight_1 * f32(slot_1);\n}\n\nfn aaPx_0( dpr_0 : f32) -> f32\n{\n    return 1.0f / max(dpr_0, 1.0f);\n}\n\nfn aaHalfPx_0( dpr_1 : f32) -> f32\n{\n    return 0.5f * aaPx_0(dpr_1);\n}\n\nfn barAaPx_0() -> f32\n{\n    return 2.0f * aaHalfPx_0(u_0.devicePixelRatio_0);\n}\n\nfn clipSide_0( value_2 : f32,  domainMin_3 : f32,  domainMax_3 : f32,  scaleType_2 : i32) -> i32\n{\n    var tolerance_0 : f32 = (domainMax_3 - domainMin_3) * 0.00000999999974738f;\n    if(value_2 > (domainMax_3 + tolerance_0))\n    {\n        return i32(1);\n    }\n    var _S20 : bool;\n    if(value_2 < (domainMin_3 - tolerance_0))\n    {\n        _S20 = scaleType_2 != i32(1);\n    }\n    else\n    {\n        _S20 = false;\n    }\n    if(_S20)\n    {\n        return i32(-1);\n    }\n    return i32(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\nfn valueToYPxScaled_0( value_3 : f32,  domainMin_4 : f32,  domainMax_4 : f32,  h_1 : f32,  scaleType_3 : i32,  symlogConstant_2 : f32) -> f32\n{\n    return (1.0f - normalizeScore_0(value_3, domainMin_4, domainMax_4, scaleType_3, symlogConstant_2)) * h_1;\n}\n\nstruct VsOut_0\n{\n    @builtin(position) position_0 : vec4<f32>,\n    @location(0) color_1 : vec4<f32>,\n    @location(1) barInkPx_0 : vec2<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) x_2 : u32,\n    @location(1) x2_2 : u32,\n    @location(2) y_1 : f32,\n    @location(3) color_2 : u32,\n    @location(4) row_0 : u32,\n};\n\nfn markScaleColor_0( _S21 : texture_2d<f32>,  _S22 : sampler,  _S23 : u32,  _S24 : i32,  _S25 : f32,  _S26 : f32,  _S27 : f32,  _S28 : i32,  _S29 : vec4<f32>,  _S30 : vec4<f32>) -> vec4<f32>\n{\n    if(_S24 == i32(0))\n    {\n        return unpackRGBA_0(_S23);\n    }\n    if(_S23 == u32(2143289345))\n    {\n        return unpackRGBA_0(u32(4289703855));\n    }\n    var scaleType_4 : i32;\n    if(((_S23 & (u32(2139095040)))) == u32(2139095040))\n    {\n        if(((_S23 & (u32(8388607)))) != u32(0))\n        {\n            return unpackRGBA_0(u32(4286611584));\n        }\n        var negative_0 : bool = ((_S23 & (u32(2147483648)))) != u32(0);\n        if(_S24 == i32(3))\n        {\n            if(negative_0)\n            {\n                scaleType_4 = i32(0);\n            }\n            else\n            {\n                scaleType_4 = _S28;\n            }\n            return u_0.colorBands_0[scaleType_4];\n        }\n        var _S31 : f32;\n        if(negative_0)\n        {\n            _S31 = 0.0f;\n        }\n        else\n        {\n            _S31 = 1.0f;\n        }\n        return rampColor_0(_S21, _S22, rampMidT_0(_S31, _S27));\n    }\n    var value_4 : f32 = (bitcast<f32>((_S23)));\n    if(_S24 == i32(3))\n    {\n        return u_0.colorBands_0[thresholdBand_0(value_4, _S28, _S29, _S30)];\n    }\n    if(_S24 == i32(2))\n    {\n        scaleType_4 = i32(1);\n    }\n    else\n    {\n        scaleType_4 = i32(0);\n    }\n    return rampColor_0(_S21, _S22, rampMidT_0(normalizeScore_0(value_4, _S25, _S26, scaleType_4, 1.0f), _S27));\n}\n\nfn bpToClipX_0( _S32 : u32) -> f32\n{\n    return hpToClipX_0(hpSplitUint_0(_S32), u_0.bpRangeX_0, u_0.zero_0);\n}\n\n@vertex\nfn vs_main( _S33 : vertexInput_0, @builtin(vertex_index) vid_1 : u32) -> VsOut_0\n{\n    var o_0 : VsOut_0;\n    var _S34 : RowPlacement_0 = rowTableLookup_0(rowTable_texture_0, rowTable_sampler_0, _S33.row_0, u_0.rowTableKeys_0, markScaleColor_0(colorRamp_texture_0, colorRamp_sampler_0, _S33.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    if(_S34.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        o_0.barInkPx_0 = vec2<f32>(0.0f);\n        return o_0;\n    }\n    var _S35 : f32 = bpToClipX_0(_S33.x_2);\n    var sx2_0 : f32 = extendToMinWidthX_0(_S35, bpToClipX_0(_S33.x2_2), u_0.minWidthPx_0, u_0.minCellDenomPx_0);\n    var q_0 : vec2<f32> = quadLocal_0(vid_1);\n    var bandTop_0 : f32 = rowBandTopPx_0(u_0.rowOffsetPx_0, u_0.rowHeight_0, _S34.slot_0);\n    var padPx_0 : f32 = barAaPx_0();\n    if(vid_1 >= u32(6))\n    {\n        var side_0 : i32 = clipSide_0(_S33.y_1, u_0.domainMin_0, u_0.domainMax_0, u_0.valueScaleType_0);\n        if(side_0 == i32(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            o_0.barInkPx_0 = vec2<f32>(0.0f);\n            return o_0;\n        }\n        var stripTop_0 : f32;\n        if(side_0 > i32(0))\n        {\n            stripTop_0 = bandTop_0;\n        }\n        else\n        {\n            stripTop_0 = bandTop_0 + u_0.rowBandPx_0 - 2.0f;\n        }\n        var stripBot_0 : f32 = stripTop_0 + 2.0f;\n        var sy_0 : f32 = mix(stripBot_0 + padPx_0, stripTop_0 - padPx_0, q_0.y);\n        o_0.position_0 = vec4<f32>(mix(_S35, sx2_0, q_0.x), yPxToClipY_0(sy_0, u_0.canvasHeight_0), 0.0f, 1.0f);\n        o_0.color_1 = vec4<f32>(unpackRGBA_0(u32(4280821958)).xyz, 1.0f);\n        o_0.barInkPx_0 = vec2<f32>(sy_0 - stripTop_0, sy_0 - stripBot_0);\n        return o_0;\n    }\n    var valueY_0 : f32 = bandTop_0 + valueToYPxScaled_0(_S33.y_1, u_0.domainMin_0, u_0.domainMax_0, u_0.rowBandPx_0, u_0.valueScaleType_0, u_0.valueSymlogConstant_0);\n    var originY_0 : f32 = bandTop_0 + u_0.originYPx_0;\n    var _S36 : f32 = min(valueY_0, originY_0);\n    var _S37 : f32 = max(valueY_0, originY_0);\n    var vy_0 : f32 = mix(_S37 + padPx_0, _S36 - padPx_0, q_0.y);\n    o_0.position_0 = vec4<f32>(mix(_S35, sx2_0, q_0.x), yPxToClipY_0(vy_0, u_0.canvasHeight_0), 0.0f, 1.0f);\n    o_0.color_1 = _S34.color_0;\n    o_0.barInkPx_0 = vec2<f32>(vy_0 - _S36, vy_0 - _S37);\n    return o_0;\n}\n\nfn aaRamp_0( signedInk_0 : f32,  widthPx_0 : f32) -> f32\n{\n    return clamp(signedInk_0 / widthPx_0 + 0.5f, 0.0f, 1.0f);\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    @location(1) barInkPx_1 : vec2<f32>,\n};\n\n@fragment\nfn fs_main( _S38 : pixelInput_0, @builtin(position) position_1 : vec4<f32>) -> pixelOutput_0\n{\n    var aaPx_1 : f32 = barAaPx_0();\n    var alpha_0 : f32 = aaRamp_0(_S38.barInkPx_1.x, aaPx_1) - aaRamp_0(_S38.barInkPx_1.y, aaPx_1);\n    if(alpha_0 <= 0.0f)\n    {\n        discard;\n    }\n    var _S39 : pixelOutput_0 = pixelOutput_0( vec4<f32>(_S38.color_3.xyz, _S38.color_3.w * alpha_0) );\n    return _S39;\n}\n\n"

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

@@ -0,0 +1,14 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from capsule.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+// Scalar twins of capsule.slang, transliterated from slangc's WGSL so
+// the Canvas2D and SVG paths run the shader's own math. See adr-051.
+
+function _max(a: number, b: number) {
+  return b > a || Number.isNaN(a) ? b : a
+}
+
+export function capsuleDistPx(alongPx: number, acrossPx: number, halfLenPx: number): number {
+  let _t0 = _max((Math.abs(alongPx) - halfLenPx), 0.0)
+  return Math.sqrt(((_t0 * _t0) + (acrossPx * acrossPx)))
+}

package/src/shaders/capsule.slang

@@ -15,6 +15,7 @@
 // comes off a clamp and the 2x2 quads straddling the cap/body junction would be
 // differentiating a discontinuity.
 //! export-consts: CAPSULE_MIN_LEN_PX
+//! js-export: capsuleDistPx
 module capsule;
 
 import antialias;
@@ -54,9 +55,16 @@ public float2 capsuleQuadLocal(
   return float2(q.x * (halfLenPx + ext), q.y * ext);
 }
 
-// CSS-px distance from `local` (the segment's frame, center-anchored) to the
-// segment {(-halfLen,0)..(halfLen,0)}. Collapses to distance-from-origin for a
-// degenerate segment — the circular dot of radius halfWidth.
+// CSS-px distance from a point in the segment's frame (center-anchored, x
+// along, y across) to the segment {(-halfLen,0)..(halfLen,0)}. Collapses to
+// distance-from-origin for a degenerate segment — the circular dot of radius
+// halfWidth. Scalar so `//! js-export` can lift it: the dotplot's pick measures
+// the cursor against the same capsule the fragment inks, end caps included.
+public float capsuleDistPx(float alongPx, float acrossPx, float halfLenPx) {
+  float beyond = max(abs(alongPx) - halfLenPx, 0.0);
+  return sqrt(beyond * beyond + acrossPx * acrossPx);
+}
+
 public float capsuleDist(float2 local, float halfLenPx) {
-  return length(float2(max(abs(local.x) - halfLenPx, 0.0), local.y));
+  return capsuleDistPx(local.x, local.y, halfLenPx);
 }

package/src/shaders/clipStrip.generated.ts

@@ -0,0 +1,6 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from clipStrip.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const CLIP_STRIP_PX = 2
+
+export const CLIP_STRIP_COLOR = 4280821958

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

@@ -0,0 +1,22 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from clipStrip.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+// Scalar twins of clipStrip.slang, transliterated from slangc's WGSL so
+// the Canvas2D and SVG paths run the shader's own math. See adr-051.
+
+export function clipSide(value: number, domainMin: number, domainMax: number, scaleType: number): number {
+  let tolerance = ((domainMax - domainMin) * 0.00000999999974738)
+  if ((value > (domainMax + tolerance))) {
+    return 1
+  }
+  let _t0: boolean
+  if ((value < (domainMin - tolerance))) {
+    _t0 = (scaleType != 1)
+  } else {
+    _t0 = false
+  }
+  if (_t0) {
+    return ((-1) | 0)
+  }
+  return 0
+}

package/src/shaders/clipStrip.slang

@@ -0,0 +1,34 @@
+// The clip strip: 2 CSS px of red across the cut edge of a bar whose value
+// the domain excludes, so a bar clamped to full height does not read as its
+// value (ADR-183). Wiggle's xyplot, the coverage band and the `bar` shape each
+// draw it as a second quad, or a pass, off the same instance, and their
+// Canvas2D painters draw it off `clipSide`'s twin.
+//! export-consts: CLIP_STRIP_PX, CLIP_STRIP_COLOR
+//! js-export: clipSide
+module clipStrip;
+
+import scoreScale;
+
+public static const float CLIP_STRIP_PX = 2.0;
+// #c62828, packed ABGR as every instance colour is: a red a shade darker than
+// the coverage band's hard-clip triangles, which stand just above the strip.
+public static const uint CLIP_STRIP_COLOR = 0xFF2828C6u;
+
+// Which edge of the axis cut the bar: 1 the top, -1 the bottom, 0 neither. A
+// log scale has no bottom to cut, since a value under its floor stands on the
+// baseline as a bar of no height, which reads as what it is.
+//
+// The cut is measured with a tolerance of a hundred-thousandth of the domain:
+// the coverage band rebuilds a depth as float32 relDepth times the region's
+// peak, and a quantile that lands on an observed depth would otherwise strip
+// on the backend whose rounding came out high and not on the other.
+public int clipSide(float value, float domainMin, float domainMax, int scaleType) {
+  float tolerance = (domainMax - domainMin) * 1e-5;
+  if (value > domainMax + tolerance) {
+    return 1;
+  }
+  if (value < domainMin - tolerance && scaleType != SCALE_TYPE_LOG) {
+    return -1;
+  }
+  return 0;
+}

package/src/shaders/clipStrip.test.ts

@@ -0,0 +1,39 @@
+import {
+  SCALE_TYPE_LINEAR,
+  SCALE_TYPE_LOG,
+  SCALE_TYPE_SYMLOG,
+} from '../scoreScale.ts'
+import { CLIP_STRIP_COLOR, CLIP_STRIP_PX } from './clipStrip.generated.ts'
+import { clipSide } from './clipStrip.js.generated.ts'
+
+// The twin every Canvas2D painter reads: which edge of the axis cut a bar. A
+// log scale has no bottom to cut, since a value under its floor is a bar of no
+// height on the baseline.
+test('clipSide names the cut edge, and a log scale has no bottom', () => {
+  for (const scaleType of [SCALE_TYPE_LINEAR, SCALE_TYPE_SYMLOG]) {
+    expect(clipSide(12, 0, 10, scaleType)).toBe(1)
+    expect(clipSide(10, 0, 10, scaleType)).toBe(0)
+    expect(clipSide(5, 0, 10, scaleType)).toBe(0)
+    expect(clipSide(0, 0, 10, scaleType)).toBe(0)
+    expect(clipSide(-3, 0, 10, scaleType)).toBe(-1)
+  }
+  expect(clipSide(200, 1, 100, SCALE_TYPE_LOG)).toBe(1)
+  expect(clipSide(0.5, 1, 100, SCALE_TYPE_LOG)).toBe(0)
+})
+
+// A depth rebuilt from float32 relDepth lands a few ulps off the quantile it
+// equals, on one backend and not the other.
+test('a value one float32 rounding past the domain is not cut', () => {
+  const max = 37
+  const rebuilt = Math.fround(max / 775) * 775
+  expect(rebuilt).not.toBe(max)
+  expect(clipSide(rebuilt, 0, max, SCALE_TYPE_LINEAR)).toBe(0)
+  expect(clipSide(max * 1.001, 0, max, SCALE_TYPE_LINEAR)).toBe(1)
+})
+
+test('the strip is 2 px of #c62828', () => {
+  expect(CLIP_STRIP_PX).toBe(2)
+  expect(
+    [0, 8, 16, 24].map(shift => (CLIP_STRIP_COLOR >>> shift) & 0xff),
+  ).toEqual([0xc6, 0x28, 0x28, 0xff])
+})

package/src/shaders/colorPack.slang

@@ -15,7 +15,11 @@ module colorPack;
 //
 // Declare such a palette as `float4[N]` instead — a vector element needs no
 // wrapper struct, lowers correctly, AND costs the same 16 bytes per slot that
-// std140 was padding the uint to anyway, so packing bought nothing here. Packing
+// std140 was padding the uint to anyway, so packing one colour per element
+// bought nothing here. FOUR to a `uint4` element, read as `p[i >> 2][i & 3]`,
+// does compile and does save — it is the only form that beats the padding — and
+// on alignments' 23-colour palette it was measured and declined:
+// agent-docs/reference/ARCHITECTURAL_LIMITS.md §"The uniform ring". Packing
 // still earns its keep in VERTEX ATTRIBUTES, which is what this is for.
 public float4 unpackRGBA(uint c) {
   return float4(

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

@@ -0,0 +1,22 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from colorRampLut.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+// Scalar twins of colorRampLut.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 rampMidT(t: number, mid: number): number {
+  let _t0 = _max(mid, (1.0 - mid))
+  return _clamp(((t + (_t0 - mid)) / (2.0 * _t0)), 0.0, 1.0)
+}