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

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

@@ -0,0 +1,52 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from connectorMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+// Scalar twins of connectorMark.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
+}
+
+function _min(a: number, b: number) {
+  return b < a || Number.isNaN(a) ? b : a
+}
+
+export function connectorBudgetPx(dx: number, dy: number): number {
+  let _t0 = _min(Math.abs(dy), 150.0)
+  return (Math.sqrt(((dx * dx) + (_t0 * _t0))) * 0.30000001192092896)
+}
+
+export function connectorHandlePx(budget: number, maxHandlePx: number): number {
+  return _min(maxHandlePx, budget)
+}
+
+export function connectorLiftPx(budget: number, dy: number, bend: number): number {
+  let _t0: number
+  if ((bend < -1.00000001504746622e+29)) {
+    _t0 = _min(budget, _max(0.0, (30.0 - Math.abs(dy))))
+  } else {
+    _t0 = (-_min(bend, 110.0) / 0.75)
+  }
+  return _t0
+}
+
+export function connectorTangent(strand: number, leading: number, reversed: number): number {
+  let _t0: number
+  if ((leading != 0.0)) {
+    _t0 = -1.0
+  } else {
+    _t0 = 1.0
+  }
+  let _t1 = (strand * _t0)
+  if ((reversed != 0.0)) {
+    _t0 = -1.0
+  } else {
+    _t0 = 1.0
+  }
+  return (_t1 * _t0)
+}
+
+export function connectorStrokePx(width: number, dpr: number): number {
+  return _max(width, (1.5 / dpr))
+}

package/src/shaders/connectorMark.slang

@@ -0,0 +1,287 @@
+// The `connector` mark shape: one instance is a stroked cubic between two
+// placed points, `x` on the displayed region whose payload holds the instance
+// and `x2` on the region `x2Region` names, each at the centre of its own row.
+// Both place through a uniform table of the view's displayed regions, so a
+// pan or zoom writes the table and no buffer. The curve is the breakpoint and
+// read-pair connector both views draw: its handles leave each end along the
+// read's direction on screen, flipped at a split junction's leading end, and
+// the pair of controls lifts into a bow over a shared row or drops by `bend`
+// px into a dip, the depth a caller sizes to the room it has.
+//
+// The hull is a triangle strip along the curve inflated by the half stroke and
+// one antialiasing ramp, and the fragment fades on its distance across, which
+// the strip's straight sections carry exactly. An instance can end in an open
+// arrowhead, and a straight one can be dashed.
+//! targets: wgsl, glsl
+//! topology: triangle-strip
+//! coverage: analytic
+//! export-consts: CONNECTOR_CURVE_SEGMENTS, CONNECTOR_BOW, CONNECTOR_MAX_REACH_PX, CONNECTOR_APEX_RATIO, CONNECTOR_MAX_HANDLE_PX, CONNECTOR_MINUS_1, CONNECTOR_MINUS_2, CONNECTOR_LEADING_2, CONNECTOR_STRAIGHT, CONNECTOR_DASHED, CONNECTOR_ARROW
+//! js-export: connectorBudgetPx, connectorHandlePx, connectorLiftPx, connectorTangent, connectorStrokePx
+
+import antialias;
+import colorPack;
+import hpmath;
+import viewRegion;
+
+public static const uint CONNECTOR_CURVE_SEGMENTS = 48u;
+static const uint CURVE_VERTS = (CONNECTOR_CURVE_SEGMENTS + 1u) * 2u;
+// Each arrowhead arm's block of the strip: the join's repeated vertex, then its
+// quad with the first corner doubled, so both joins are degenerate triangles.
+static const uint ARM_VERTS = 6u;
+public static const uint VERTS_PER_INSTANCE = CURVE_VERTS + 2u * ARM_VERTS;
+
+// The `bend` lane's value for a connector that bows up over its row rather
+// than dipping by a set depth.
+public static const float CONNECTOR_BOW = -1.0e30;
+// Every shaping term stays under this fraction of the span the curve covers,
+// so a handle never folds the cubic back over itself.
+static const float CONNECTOR_SPAN_FACTOR = 0.3;
+public static const float CONNECTOR_MAX_HANDLE_PX = 200.0;
+// The handles are horizontal, so a tall vertical gap only widens the bulge;
+// past this it earns no more.
+static const float CONNECTOR_MAX_VERTICAL_SPAN_PX = 150.0;
+// A symmetric cubic's extreme sits 3/4 of the way from its ends to its
+// controls, so a dip asked to reach a depth puts the controls 4/3 past it.
+public static const float CONNECTOR_APEX_RATIO = 0.75;
+static const float CONNECTOR_MAX_BOW_PX = 30.0;
+// The deepest a dip may reach below its ends, a ceiling on `bend`.
+public static const float CONNECTOR_MAX_REACH_PX = 110.0;
+
+// The `bits` lane.
+public static const uint CONNECTOR_MINUS_1 = 1u;
+public static const uint CONNECTOR_MINUS_2 = 2u;
+public static const uint CONNECTOR_LEADING_2 = 4u;
+public static const uint CONNECTOR_STRAIGHT = 8u;
+public static const uint CONNECTOR_DASHED = 16u;
+public static const uint CONNECTOR_ARROW = 32u;
+
+struct ConnectorInstance {
+  uint  x        : ATTR0;
+  uint  x2       : ATTR1;
+  uint  x2Region : ATTR2;
+  uint  row      : ATTR3;
+  uint  row2     : ATTR4;
+  float bend     : ATTR5;
+  float width    : ATTR6;
+  uint  color    : ATTR7;
+  uint  bits     : ATTR8;
+};
+
+struct Uniforms {
+  float  canvasHeight;
+  float  canvasWidth;
+  float  devicePixelRatio;
+  // Canvas y of row 0's centre, and the px between two rows' centres.
+  float  rowOffsetPx;
+  float  rowPitchPx;
+  float  maxHandlePx;
+  float  dashPx;
+  float  gapPx;
+  float  arrowPx;
+  uint   regionCount;
+  float  zero;
+  // The table entry of the region whose payload is drawing, which every
+  // instance's `x` is on.
+  float4 ownEntry;
+  float4 regionTable[VIEW_REGION_MAX];
+};
+[[vk::binding(1, 0)]] ConstantBuffer<Uniforms> u;
+
+struct VsOut {
+  float4 position : SV_Position;
+  float4 color    : COLOR0;
+  // Signed px from the centreline, and px along it from the `x` end.
+  float  across   : COLOR1;
+  float  along    : COLOR2;
+  nointerpolation float halfWidthPx : COLOR3;
+  nointerpolation float dashed      : COLOR4;
+};
+
+// How much shaping a curve may spend: a share of the distance it covers.
+public float connectorBudgetPx(float dx, float dy) {
+  float h = min(abs(dy), CONNECTOR_MAX_VERTICAL_SPAN_PX);
+  return sqrt(dx * dx + h * h) * CONNECTOR_SPAN_FACTOR;
+}
+
+public float connectorHandlePx(float budget, float maxHandlePx) {
+  return min(maxHandlePx, budget);
+}
+
+// How far both controls rise above their ends, negative for a dip. A bow fades
+// as the rows separate, since a cross-row curve already has height; a dip is
+// not held to the budget, because a vertical offset cannot fold the curve.
+public float connectorLiftPx(float budget, float dy, float bend) {
+  return bend < -1.0e29
+    ? min(budget, max(0.0, CONNECTOR_MAX_BOW_PX - abs(dy)))
+    : -min(bend, CONNECTOR_MAX_REACH_PX) / CONNECTOR_APEX_RATIO;
+}
+
+// The screen sign of an end's handle: along the read's direction at a trailing
+// end, against it at a split junction's leading end, mirrored on a reversed
+// region.
+public float connectorTangent(float strand, float leading, float reversed) {
+  return strand * (leading != 0.0 ? -1.0 : 1.0) * (reversed != 0.0 ? -1.0 : 1.0);
+}
+
+// Floored at 1.5 device px, so the single antialiasing ramp has room.
+public float connectorStrokePx(float width, float dpr) {
+  return max(width, 1.5 / dpr);
+}
+
+struct Cubic {
+  float2 p0;
+  float2 p1;
+  float2 p2;
+  float2 p3;
+};
+
+float2 cubicAt(Cubic c, float t) {
+  float s = 1.0 - t;
+  return s * s * s * c.p0 + 3.0 * s * s * t * c.p1 + 3.0 * s * t * t * c.p2
+    + t * t * t * c.p3;
+}
+
+float2 cubicDirection(Cubic c, float t) {
+  float s = 1.0 - t;
+  float2 d = 3.0 * s * s * (c.p1 - c.p0) + 6.0 * s * t * (c.p2 - c.p1)
+    + 3.0 * t * t * (c.p3 - c.p2);
+  float len = length(d);
+  if (len > 1e-6) { return d / len; }
+  float2 chord = c.p3 - c.p0;
+  float chordLen = length(chord);
+  return chordLen > 1e-6 ? chord / chordLen : float2(1.0, 0.0);
+}
+
+float bitSet(uint bits, uint bit) {
+  return (bits & bit) != 0u ? 1.0 : 0.0;
+}
+
+Cubic placeCubic(ConnectorInstance inst, float4 entry2) {
+  Cubic c;
+  c.p0 = float2(viewRegionPx(inst.x, u.ownEntry, u.zero),
+                u.rowOffsetPx + float(inst.row) * u.rowPitchPx);
+  c.p3 = float2(viewRegionPx(inst.x2, entry2, u.zero),
+                u.rowOffsetPx + float(inst.row2) * u.rowPitchPx);
+  if ((inst.bits & CONNECTOR_STRAIGHT) != 0u) {
+    c.p1 = c.p0;
+    c.p2 = c.p3;
+    return c;
+  }
+  float2 d = c.p3 - c.p0;
+  float budget = connectorBudgetPx(d.x, d.y);
+  float handle = connectorHandlePx(budget, u.maxHandlePx);
+  float lift = connectorLiftPx(budget, d.y, inst.bend);
+  float t1 = connectorTangent(1.0 - 2.0 * bitSet(inst.bits, CONNECTOR_MINUS_1),
+                              0.0, u.ownEntry.w < 0.0 ? 1.0 : 0.0);
+  float t2 = connectorTangent(1.0 - 2.0 * bitSet(inst.bits, CONNECTOR_MINUS_2),
+                              bitSet(inst.bits, CONNECTOR_LEADING_2),
+                              entry2.w < 0.0 ? 1.0 : 0.0);
+  c.p1 = float2(c.p0.x + handle * t1, c.p0.y - lift);
+  c.p2 = float2(c.p3.x + handle * t2, c.p3.y - lift);
+  return c;
+}
+
+// The curve's length up to segment `seg`, as the strip's chords measure it.
+float lengthTo(Cubic c, uint seg) {
+  float total = 0.0;
+  float2 prev = c.p0;
+  for (uint k = 1u; k <= seg; k++) {
+    float2 next = cubicAt(c, float(k) / float(CONNECTOR_CURVE_SEGMENTS));
+    total += length(next - prev);
+    prev = next;
+  }
+  return total;
+}
+
+struct StripVertex {
+  float2 screen;
+  float across;
+  float along;
+};
+
+StripVertex curveVertex(Cubic c, uint seg, float side, float hw, bool dashed) {
+  float t = float(seg) / float(CONNECTOR_CURVE_SEGMENTS);
+  float2 dir = cubicDirection(c, t);
+  StripVertex v;
+  v.across = side * hw;
+  v.screen = cubicAt(c, t) + float2(-dir.y, dir.x) * v.across;
+  v.along = dashed ? lengthTo(c, seg) : 0.0;
+  return v;
+}
+
+// An open arrowhead's arm `a` at the `x2` end, pointing the way the read
+// continues: from the tip back and off the row by `arrowPx`.
+float2 armEnd(Cubic c, uint a, float strand2, bool reversed2) {
+  float back = (reversed2 ? -strand2 : strand2) * -u.arrowPx;
+  return c.p3 + float2(back, a == 0u ? -u.arrowPx : u.arrowPx);
+}
+
+StripVertex armCorner(Cubic c, uint a, uint corner, float hw, float strand2, bool reversed2) {
+  float2 end = armEnd(c, a, strand2, reversed2);
+  float2 axis = end - c.p3;
+  float len = max(length(axis), 1e-6);
+  float2 dir = axis / len;
+  StripVertex v;
+  v.across = (corner % 2u == 0u) ? -hw : hw;
+  v.along = 0.0;
+  float at = corner >= 2u ? len : 0.0;
+  v.screen = c.p3 + dir * at + float2(-dir.y, dir.x) * v.across;
+  return v;
+}
+
+[shader("vertex")]
+VsOut vs_main(ConnectorInstance inst, uint vid : SV_VertexID) {
+  bool placed = inst.x2Region < u.regionCount;
+  float4 entry2 = u.regionTable[min(inst.x2Region, VIEW_REGION_MAX - 1u)];
+  Cubic c = placeCubic(inst, entry2);
+  float halfWidth = 0.5 * connectorStrokePx(inst.width, u.devicePixelRatio);
+  float hw = halfWidth + aaHalfPx(u.devicePixelRatio);
+  bool dashed = (inst.bits & CONNECTOR_DASHED) != 0u;
+  bool arrow = (inst.bits & CONNECTOR_ARROW) != 0u;
+  float strand2 = 1.0 - 2.0 * bitSet(inst.bits, CONNECTOR_MINUS_2);
+  bool reversed2 = entry2.w < 0.0;
+
+  StripVertex v;
+  bool onArm = false;
+  if (vid < CURVE_VERTS) {
+    v = curveVertex(c, vid / 2u, (vid % 2u == 0u) ? -1.0 : 1.0, hw, dashed);
+  } else {
+    uint j = vid - CURVE_VERTS;
+    uint a = j / ARM_VERTS;
+    uint k = j % ARM_VERTS;
+    // The block starts from the vertex before it: the curve's last, or the
+    // first arm's last corner.
+    StripVertex start = a == 0u
+      ? curveVertex(c, CONNECTOR_CURVE_SEGMENTS, 1.0, hw, dashed)
+      : (arrow ? armCorner(c, 0u, 3u, hw, strand2, reversed2)
+               : curveVertex(c, CONNECTOR_CURVE_SEGMENTS, 1.0, hw, dashed));
+    if (k == 0u || !arrow) {
+      v = start;
+    } else {
+      v = armCorner(c, a, k <= 2u ? 0u : k - 2u, hw, strand2, reversed2);
+      onArm = true;
+    }
+  }
+
+  VsOut o;
+  o.position = placed
+    ? float4(v.screen.x / u.canvasWidth * 2.0 - 1.0,
+             yPxToClipY(v.screen.y, u.canvasHeight), 0.0, 1.0)
+    : float4(-2.0, -2.0, 0.0, 1.0);
+  o.color = unpackRGBA(inst.color);
+  o.across = v.across;
+  o.along = v.along;
+  o.halfWidthPx = halfWidth;
+  o.dashed = dashed && !onArm ? 1.0 : 0.0;
+  return o;
+}
+
+[shader("fragment")]
+float4 fs_main(VsOut fragIn) : SV_Target {
+  float dash = fragIn.dashed > 0.5 && u.dashPx > 0.0
+    ? dashCoverage(fragIn.along, u.dashPx, u.gapPx, u.devicePixelRatio)
+    : 1.0;
+  float alpha = edgeCoverage(fragIn.halfWidthPx - abs(fragIn.across), u.devicePixelRatio)
+    * dash;
+  return float4(fragIn.color.rgb, fragIn.color.a * alpha);
+}

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

@@ -0,0 +1,4 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from connectorMark.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+export const WGSL_SOURCE = "struct Uniforms_std140_0\n{\n    @align(16) canvasHeight_0 : f32,\n    @align(4) canvasWidth_0 : f32,\n    @align(8) devicePixelRatio_0 : f32,\n    @align(4) rowOffsetPx_0 : f32,\n    @align(16) rowPitchPx_0 : f32,\n    @align(4) maxHandlePx_0 : f32,\n    @align(8) dashPx_0 : f32,\n    @align(4) gapPx_0 : f32,\n    @align(16) arrowPx_0 : f32,\n    @align(4) regionCount_0 : u32,\n    @align(8) zero_0 : f32,\n    @align(16) ownEntry_0 : vec4<f32>,\n    @align(16) regionTable_0 : array<vec4<f32>, i32(256)>,\n};\n\n@binding(1) @group(0) var<uniform> u_0 : Uniforms_std140_0;\nfn hpSplitUint_0( value_0 : u32) -> vec2<f32>\n{\n    var lo_0 : u32 = (value_0 & (u32(4095)));\n    return vec2<f32>(f32(value_0 - lo_0), f32(lo_0));\n}\n\nfn hpViewportOffsetBp_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    return max(bpRange_0.x - splitPos_0.x, floorBp_0) + max(bpRange_0.y - splitPos_0.y, floorBp_0);\n}\n\nfn viewRegionPx_0( bp_0 : u32,  r_0 : vec4<f32>,  hpZero_1 : f32) -> f32\n{\n    return r_0.x - hpViewportOffsetBp_0(hpSplitUint_0(bp_0), vec3<f32>(r_0.y, r_0.z, 0.0f), hpZero_1) * r_0.w;\n}\n\nfn connectorBudgetPx_0( dx_0 : f32,  dy_0 : f32) -> f32\n{\n    var _S1 : f32 = min(abs(dy_0), 150.0f);\n    return sqrt(dx_0 * dx_0 + _S1 * _S1) * 0.30000001192092896f;\n}\n\nfn connectorHandlePx_0( budget_0 : f32,  maxHandlePx_1 : f32) -> f32\n{\n    return min(maxHandlePx_1, budget_0);\n}\n\nfn connectorLiftPx_0( budget_1 : f32,  dy_1 : f32,  bend_0 : f32) -> f32\n{\n    var _S2 : f32;\n    if(bend_0 < -1.00000001504746622e+29f)\n    {\n        _S2 = min(budget_1, max(0.0f, 30.0f - abs(dy_1)));\n    }\n    else\n    {\n        _S2 = - min(bend_0, 110.0f) / 0.75f;\n    }\n    return _S2;\n}\n\nfn bitSet_0( bits_0 : u32,  bit_0 : u32) -> f32\n{\n    var _S3 : f32;\n    if(((bits_0 & (bit_0))) != u32(0))\n    {\n        _S3 = 1.0f;\n    }\n    else\n    {\n        _S3 = 0.0f;\n    }\n    return _S3;\n}\n\nfn connectorTangent_0( strand_0 : f32,  leading_0 : f32,  reversed_0 : f32) -> f32\n{\n    var _S4 : f32;\n    if(leading_0 != 0.0f)\n    {\n        _S4 = -1.0f;\n    }\n    else\n    {\n        _S4 = 1.0f;\n    }\n    var _S5 : f32 = strand_0 * _S4;\n    if(reversed_0 != 0.0f)\n    {\n        _S4 = -1.0f;\n    }\n    else\n    {\n        _S4 = 1.0f;\n    }\n    return _S5 * _S4;\n}\n\nstruct Cubic_0\n{\n     p0_0 : vec2<f32>,\n     p1_0 : vec2<f32>,\n     p2_0 : vec2<f32>,\n     p3_0 : vec2<f32>,\n};\n\nstruct ConnectorInstance_0\n{\n     x_0 : u32,\n     x2_0 : u32,\n     x2Region_0 : u32,\n     row_0 : u32,\n     row2_0 : u32,\n     bend_1 : f32,\n     width_0 : f32,\n     color_0 : u32,\n     bits_1 : u32,\n};\n\nfn placeCubic_0( inst_0 : ConnectorInstance_0,  entry2_0 : vec4<f32>) -> Cubic_0\n{\n    var c_0 : Cubic_0;\n    c_0.p0_0 = vec2<f32>(viewRegionPx_0(inst_0.x_0, u_0.ownEntry_0, u_0.zero_0), u_0.rowOffsetPx_0 + f32(inst_0.row_0) * u_0.rowPitchPx_0);\n    c_0.p3_0 = vec2<f32>(viewRegionPx_0(inst_0.x2_0, entry2_0, u_0.zero_0), u_0.rowOffsetPx_0 + f32(inst_0.row2_0) * u_0.rowPitchPx_0);\n    if((((inst_0.bits_1) & (u32(8)))) != u32(0))\n    {\n        c_0.p1_0 = c_0.p0_0;\n        c_0.p2_0 = c_0.p3_0;\n        return c_0;\n    }\n    var d_0 : vec2<f32> = c_0.p3_0 - c_0.p0_0;\n    var _S6 : f32 = d_0.y;\n    var budget_2 : f32 = connectorBudgetPx_0(d_0.x, _S6);\n    var handle_0 : f32 = connectorHandlePx_0(budget_2, u_0.maxHandlePx_0);\n    var lift_0 : f32 = connectorLiftPx_0(budget_2, _S6, inst_0.bend_1);\n    var _S7 : f32 = 1.0f - 2.0f * bitSet_0(inst_0.bits_1, u32(1));\n    var _S8 : f32;\n    if((u_0.ownEntry_0.w) < 0.0f)\n    {\n        _S8 = 1.0f;\n    }\n    else\n    {\n        _S8 = 0.0f;\n    }\n    var t1_0 : f32 = connectorTangent_0(_S7, 0.0f, _S8);\n    var _S9 : f32 = 1.0f - 2.0f * bitSet_0(inst_0.bits_1, u32(2));\n    var _S10 : f32 = bitSet_0(inst_0.bits_1, u32(4));\n    if((entry2_0.w) < 0.0f)\n    {\n        _S8 = 1.0f;\n    }\n    else\n    {\n        _S8 = 0.0f;\n    }\n    var t2_0 : f32 = connectorTangent_0(_S9, _S10, _S8);\n    c_0.p1_0 = vec2<f32>(c_0.p0_0.x + handle_0 * t1_0, c_0.p0_0.y - lift_0);\n    c_0.p2_0 = vec2<f32>(c_0.p3_0.x + handle_0 * t2_0, c_0.p3_0.y - lift_0);\n    return c_0;\n}\n\nfn connectorStrokePx_0( width_1 : f32,  dpr_0 : f32) -> f32\n{\n    return max(width_1, 1.5f / dpr_0);\n}\n\nfn aaPx_0( dpr_1 : f32) -> f32\n{\n    return 1.0f / max(dpr_1, 1.0f);\n}\n\nfn aaHalfPx_0( dpr_2 : f32) -> f32\n{\n    return 0.5f * aaPx_0(dpr_2);\n}\n\nfn cubicDirection_0( c_1 : Cubic_0,  t_0 : f32) -> vec2<f32>\n{\n    var s_0 : f32 = 1.0f - t_0;\n    var d_1 : vec2<f32> = vec2<f32>((3.0f * s_0 * s_0)) * (c_1.p1_0 - c_1.p0_0) + vec2<f32>((6.0f * s_0 * t_0)) * (c_1.p2_0 - c_1.p1_0) + vec2<f32>((3.0f * t_0 * t_0)) * (c_1.p3_0 - c_1.p2_0);\n    var len_0 : f32 = length(d_1);\n    if(len_0 > 9.99999997475242708e-07f)\n    {\n        return d_1 / vec2<f32>(len_0);\n    }\n    var chord_0 : vec2<f32> = c_1.p3_0 - c_1.p0_0;\n    var chordLen_0 : f32 = length(chord_0);\n    var _S11 : vec2<f32>;\n    if(chordLen_0 > 9.99999997475242708e-07f)\n    {\n        _S11 = chord_0 / vec2<f32>(chordLen_0);\n    }\n    else\n    {\n        _S11 = vec2<f32>(1.0f, 0.0f);\n    }\n    return _S11;\n}\n\nfn cubicAt_0( c_2 : Cubic_0,  t_1 : f32) -> vec2<f32>\n{\n    var s_1 : f32 = 1.0f - t_1;\n    var _S12 : f32 = 3.0f * s_1;\n    return vec2<f32>((s_1 * s_1 * s_1)) * c_2.p0_0 + vec2<f32>((_S12 * s_1 * t_1)) * c_2.p1_0 + vec2<f32>((_S12 * t_1 * t_1)) * c_2.p2_0 + vec2<f32>((t_1 * t_1 * t_1)) * c_2.p3_0;\n}\n\nfn lengthTo_0( c_3 : Cubic_0,  seg_0 : u32) -> f32\n{\n    var prev_0 : vec2<f32> = c_3.p0_0;\n    var k_0 : u32 = u32(1);\n    var total_0 : f32 = 0.0f;\n    for(;;)\n    {\n        if(k_0 <= seg_0)\n        {\n        }\n        else\n        {\n            break;\n        }\n        var next_0 : vec2<f32> = cubicAt_0(c_3, f32(k_0) / 48.0f);\n        var total_1 : f32 = total_0 + length(next_0 - prev_0);\n        var _S13 : u32 = k_0 + u32(1);\n        prev_0 = next_0;\n        k_0 = _S13;\n        total_0 = total_1;\n    }\n    return total_0;\n}\n\nstruct StripVertex_0\n{\n     screen_0 : vec2<f32>,\n     across_0 : f32,\n     along_0 : f32,\n};\n\nfn curveVertex_0( c_4 : Cubic_0,  seg_1 : u32,  side_0 : f32,  hw_0 : f32,  dashed_0 : bool) -> StripVertex_0\n{\n    var t_2 : f32 = f32(seg_1) / 48.0f;\n    var dir_0 : vec2<f32> = cubicDirection_0(c_4, t_2);\n    var v_0 : StripVertex_0;\n    var _S14 : f32 = side_0 * hw_0;\n    v_0.across_0 = _S14;\n    v_0.screen_0 = cubicAt_0(c_4, t_2) + vec2<f32>(- dir_0.y, dir_0.x) * vec2<f32>(_S14);\n    var _S15 : f32;\n    if(dashed_0)\n    {\n        _S15 = lengthTo_0(c_4, seg_1);\n    }\n    else\n    {\n        _S15 = 0.0f;\n    }\n    v_0.along_0 = _S15;\n    return v_0;\n}\n\nfn armEnd_0( c_5 : Cubic_0,  a_0 : u32,  strand2_0 : f32,  reversed2_0 : bool) -> vec2<f32>\n{\n    var _S16 : f32;\n    if(reversed2_0)\n    {\n        _S16 = - strand2_0;\n    }\n    else\n    {\n        _S16 = strand2_0;\n    }\n    var back_0 : f32 = _S16 * - u_0.arrowPx_0;\n    if(a_0 == u32(0))\n    {\n        _S16 = - u_0.arrowPx_0;\n    }\n    else\n    {\n        _S16 = u_0.arrowPx_0;\n    }\n    return c_5.p3_0 + vec2<f32>(back_0, _S16);\n}\n\nfn armCorner_0( c_6 : Cubic_0,  a_1 : u32,  corner_0 : u32,  hw_1 : f32,  strand2_1 : f32,  reversed2_1 : bool) -> StripVertex_0\n{\n    var axis_0 : vec2<f32> = armEnd_0(c_6, a_1, strand2_1, reversed2_1) - c_6.p3_0;\n    var _S17 : f32 = max(length(axis_0), 9.99999997475242708e-07f);\n    var dir_1 : vec2<f32> = axis_0 / vec2<f32>(_S17);\n    var v_1 : StripVertex_0;\n    var at_0 : f32;\n    if((corner_0 % u32(2)) == u32(0))\n    {\n        at_0 = - hw_1;\n    }\n    else\n    {\n        at_0 = hw_1;\n    }\n    v_1.across_0 = at_0;\n    v_1.along_0 = 0.0f;\n    if(corner_0 >= u32(2))\n    {\n        at_0 = _S17;\n    }\n    else\n    {\n        at_0 = 0.0f;\n    }\n    v_1.screen_0 = c_6.p3_0 + dir_1 * vec2<f32>(at_0) + vec2<f32>(- dir_1.y, dir_1.x) * vec2<f32>(v_1.across_0);\n    return v_1;\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 unpackRGBA_0( c_7 : u32) -> vec4<f32>\n{\n    return vec4<f32>(f32((((c_7 >> (u32(0)))) & (u32(255)))), f32((((c_7 >> (u32(8)))) & (u32(255)))), f32((((c_7 >> (u32(16)))) & (u32(255)))), f32((((c_7 >> (u32(24)))) & (u32(255))))) / vec4<f32>(255.0f);\n}\n\nstruct VsOut_0\n{\n    @builtin(position) position_0 : vec4<f32>,\n    @location(0) color_1 : vec4<f32>,\n    @location(1) across_1 : f32,\n    @location(2) along_1 : f32,\n    @interpolate(flat) @location(3) halfWidthPx_0 : f32,\n    @interpolate(flat) @location(4) dashed_1 : f32,\n};\n\nstruct vertexInput_0\n{\n    @location(0) x_1 : u32,\n    @location(1) x2_1 : u32,\n    @location(2) x2Region_1 : u32,\n    @location(3) row_1 : u32,\n    @location(4) row2_1 : u32,\n    @location(5) bend_2 : f32,\n    @location(6) width_2 : f32,\n    @location(7) color_2 : u32,\n    @location(8) bits_2 : u32,\n};\n\n@vertex\nfn vs_main( _S18 : vertexInput_0, @builtin(vertex_index) vid_0 : u32) -> VsOut_0\n{\n    var _S19 : ConnectorInstance_0 = ConnectorInstance_0( _S18.x_1, _S18.x2_1, _S18.x2Region_1, _S18.row_1, _S18.row2_1, _S18.bend_2, _S18.width_2, _S18.color_2, _S18.bits_2 );\n    var placed_0 : bool = (_S18.x2Region_1) < (u_0.regionCount_0);\n    var c_8 : Cubic_0 = placeCubic_0(_S19, u_0.regionTable_0[min(_S18.x2Region_1, u32(255))]);\n    var halfWidth_0 : f32 = 0.5f * connectorStrokePx_0(_S18.width_2, u_0.devicePixelRatio_0);\n    var hw_2 : f32 = halfWidth_0 + aaHalfPx_0(u_0.devicePixelRatio_0);\n    var dashed_2 : bool = (((_S18.bits_2) & (u32(16)))) != u32(0);\n    var arrow_0 : bool = (((_S18.bits_2) & (u32(32)))) != u32(0);\n    var strand2_2 : f32 = 1.0f - 2.0f * bitSet_0(_S18.bits_2, u32(2));\n    var reversed2_2 : bool = (u_0.regionTable_0[min(_S18.x2Region_1, u32(255))].w) < 0.0f;\n    var v_2 : StripVertex_0;\n    var onArm_0 : bool;\n    var _S20 : f32;\n    if(vid_0 < u32(98))\n    {\n        var _S21 : u32 = vid_0 / u32(2);\n        if((vid_0 % u32(2)) == u32(0))\n        {\n            _S20 = -1.0f;\n        }\n        else\n        {\n            _S20 = 1.0f;\n        }\n        v_2 = curveVertex_0(c_8, _S21, _S20, hw_2, dashed_2);\n        onArm_0 = false;\n    }\n    else\n    {\n        var j_0 : u32 = vid_0 - u32(98);\n        var a_2 : u32 = j_0 / u32(6);\n        var k_1 : u32 = j_0 % u32(6);\n        var start_0 : StripVertex_0;\n        if(a_2 == u32(0))\n        {\n            start_0 = curveVertex_0(c_8, u32(48), 1.0f, hw_2, dashed_2);\n        }\n        else\n        {\n            if(arrow_0)\n            {\n                start_0 = armCorner_0(c_8, u32(0), u32(3), hw_2, strand2_2, reversed2_2);\n            }\n            else\n            {\n                start_0 = curveVertex_0(c_8, u32(48), 1.0f, hw_2, dashed_2);\n            }\n        }\n        if(k_1 == u32(0))\n        {\n            onArm_0 = true;\n        }\n        else\n        {\n            onArm_0 = !arrow_0;\n        }\n        if(onArm_0)\n        {\n            v_2 = start_0;\n            onArm_0 = false;\n        }\n        else\n        {\n            var _S22 : u32;\n            if(k_1 <= u32(2))\n            {\n                _S22 = u32(0);\n            }\n            else\n            {\n                _S22 = k_1 - u32(2);\n            }\n            v_2 = armCorner_0(c_8, a_2, _S22, hw_2, strand2_2, reversed2_2);\n            onArm_0 = true;\n        }\n    }\n    var o_0 : VsOut_0;\n    var _S23 : vec4<f32>;\n    if(placed_0)\n    {\n        _S23 = vec4<f32>(v_2.screen_0.x / u_0.canvasWidth_0 * 2.0f - 1.0f, yPxToClipY_0(v_2.screen_0.y, u_0.canvasHeight_0), 0.0f, 1.0f);\n    }\n    else\n    {\n        _S23 = vec4<f32>(-2.0f, -2.0f, 0.0f, 1.0f);\n    }\n    o_0.position_0 = _S23;\n    o_0.color_1 = unpackRGBA_0(_S18.color_2);\n    o_0.across_1 = v_2.across_0;\n    o_0.along_1 = v_2.along_0;\n    o_0.halfWidthPx_0 = halfWidth_0;\n    if(dashed_2)\n    {\n        onArm_0 = !onArm_0;\n    }\n    else\n    {\n        onArm_0 = false;\n    }\n    if(onArm_0)\n    {\n        _S20 = 1.0f;\n    }\n    else\n    {\n        _S20 = 0.0f;\n    }\n    o_0.dashed_1 = _S20;\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\nfn dashCoverageAt_0( coord_0 : f32,  dashPx_1 : f32,  gapPx_1 : f32,  aaPx_1 : f32) -> f32\n{\n    var period_0 : f32 = dashPx_1 + gapPx_1;\n    var shifted_0 : f32 = coord_0 + gapPx_1 * 0.5f;\n    return aaRamp_0(dashPx_1 * 0.5f - abs(shifted_0 - period_0 * floor(shifted_0 / period_0) - period_0 * 0.5f), aaPx_1);\n}\n\nfn dashCoverage_0( coord_1 : f32,  dashPx_2 : f32,  gapPx_2 : f32,  dpr_3 : f32) -> f32\n{\n    return dashCoverageAt_0(coord_1, dashPx_2, gapPx_2, aaPx_0(dpr_3));\n}\n\nfn edgeCoverage_0( signedInkCssPx_0 : f32,  dpr_4 : f32) -> f32\n{\n    return aaRamp_0(signedInkCssPx_0, aaPx_0(dpr_4));\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) across_2 : f32,\n    @location(2) along_2 : f32,\n    @interpolate(flat) @location(3) halfWidthPx_1 : f32,\n    @interpolate(flat) @location(4) dashed_3 : f32,\n};\n\n@fragment\nfn fs_main( _S24 : pixelInput_0, @builtin(position) position_1 : vec4<f32>) -> pixelOutput_0\n{\n    var _S25 : bool;\n    if((_S24.dashed_3) > 0.5f)\n    {\n        _S25 = (u_0.dashPx_0) > 0.0f;\n    }\n    else\n    {\n        _S25 = false;\n    }\n    var dash_0 : f32;\n    if(_S25)\n    {\n        dash_0 = dashCoverage_0(_S24.along_2, u_0.dashPx_0, u_0.gapPx_0, u_0.devicePixelRatio_0);\n    }\n    else\n    {\n        dash_0 = 1.0f;\n    }\n    var _S26 : pixelOutput_0 = pixelOutput_0( vec4<f32>(_S24.color_3.xyz, _S24.color_3.w * (edgeCoverage_0(_S24.halfWidthPx_1 - abs(_S24.across_2), u_0.devicePixelRatio_0) * dash_0)) );\n    return _S26;\n}\n\n"

package/src/shaders/coverageBand.slang

@@ -95,16 +95,14 @@ public struct CoverageBandUniforms {
 };
 
 // Coverage-band clip indicator triangle, in CSS px: half-width and height.
-// Four consumers depend on this one number: coverageIndicator.slang draws the
-// triangle, coverageInterbase.slang hangs its bars directly below one, and the
-// Canvas2D painter + hit test read it from `@jbrowse/alignments-core` — which
-// the `consts-out` directive above generates from these declarations, so the
-// cross-package copy can't drift either.
+// coverageInterbase.slang hangs its bars directly below one, and the point
+// mark that draws it reads the size from `@jbrowse/alignments-core`, which the
+// `consts-out` directive above generates from these declarations, so the
+// cross-package copy can't drift.
 public static const float TRIANGLE_HW = 3.5;
 public static const float TRIANGLE_H  = 4.5;
 
-// Shared vertex output for the band's four solid-fill passes. The indicator
-// pass carries a barycentric varying as well and declares its own.
+// Shared vertex output for the band's solid-fill passes.
 public struct CoverageVsOut {
   public float4 position : SV_Position;
   public float4 color    : COLOR0;
@@ -295,11 +293,11 @@ public float3 covSnpColor(float ct, CoverageBandUniforms u) {
   }
 }
 
-// Interbase mark fill, for both the histogram bars and the indicator triangles:
-// 1=insertion 2=softclip 3=hardclip, with anything else taking the insertion
-// color — the same fallback `drawIndicators`' lut takes for an out-of-range
-// kind. Nothing reaches that arm today (the interbase compute only ever emits
-// 1..3), which is exactly why the two spellings disagreeing would be invisible.
+// The interbase histogram's fill: 1=insertion 2=softclip 3=hardclip, with
+// anything else taking the insertion color — the same fallback
+// `clipKindColor` takes for an out-of-range kind. Nothing reaches that arm
+// today (the interbase compute only ever emits 1..3), which is exactly why the
+// two spellings disagreeing would be invisible.
 public float3 covClipKindColor(float ct, CoverageBandUniforms u) {
   int ci = int(ct);
   if (ci == 2) { return unpackRGBA(u.colorSoftclipIndicator).xyz; }

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

@@ -1,6 +1,6 @@
 // AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageInterbase.slang.
 // Do not edit. Run `pnpm gen:shaders` to regenerate.
 
-export const GLSL_VERTEX = "#version 300 es\nprecision highp float;\nprecision highp int;\nstruct CoverageBandUniforms_0\n{\n    float bpHi_0;\n    float bpLo_0;\n    float bpLen_0;\n    float hpZero_0;\n    float canvasW_0;\n    float canvasH_0;\n    float covHeight_0;\n    float covYOffset_0;\n    float covTop_0;\n    float regionMaxDepth_0;\n    float depthDomainMax_0;\n    float depthDomainMin_0;\n    float coverageSymlogConstant_0;\n    float binSize_0;\n    float interbaseHeight_0;\n    float snpMinFreq_0;\n    float reversed_0;\n    int coverageScaleType_0;\n    uint colorCoverage_0;\n    uint colorBaseA_0;\n    uint colorBaseC_0;\n    uint colorBaseG_0;\n    uint colorBaseT_0;\n    uint colorBaseN_0;\n    uint colorInsertionIndicator_0;\n    uint colorSoftclipIndicator_0;\n    uint colorHardclipIndicator_0;\n};\n\n\nlayout(std140) uniform Uniforms\n{\n    float bpHi_0;\n    float bpLo_0;\n    float bpLen_0;\n    float hpZero_0;\n    float canvasW_0;\n    float canvasH_0;\n    float covHeight_0;\n    float covYOffset_0;\n    float covTop_0;\n    float regionMaxDepth_0;\n    float depthDomainMax_0;\n    float depthDomainMin_0;\n    float coverageSymlogConstant_0;\n    float binSize_0;\n    float interbaseHeight_0;\n    float snpMinFreq_0;\n    float reversed_0;\n    int coverageScaleType_0;\n    uint colorCoverage_0;\n    uint colorBaseA_0;\n    uint colorBaseC_0;\n    uint colorBaseG_0;\n    uint colorBaseT_0;\n    uint colorBaseN_0;\n    uint colorInsertionIndicator_0;\n    uint colorSoftclipIndicator_0;\n    uint colorHardclipIndicator_0;\n}u_0;\n\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S1 = v_0 == 0U;\n\n    bool _S2;\n\n    if(_S1)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 2U;\n\n    }\n\n    if(_S2)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 3U;\n\n    }\n\n    float _S3;\n\n    if(_S2)\n    {\n\n        _S3 = 0.0;\n\n    }\n    else\n    {\n\n        _S3 = 1.0;\n\n    }\n    if(_S1)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 1U;\n\n    }\n\n    if(_S2)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 4U;\n\n    }\n\n    float _S4;\n\n    if(_S2)\n    {\n\n        _S4 = 0.0;\n\n    }\n    else\n    {\n\n        _S4 = 1.0;\n\n    }\n\n    return vec2(_S3, _S4);\n}\n\n\nvec2 hpSplitUint_0(uint value_0)\n{\n\n    uint lo_0 = value_0 & 4095U;\n\n    return vec2(float(value_0 - lo_0), float(lo_0));\n}\n\n\nfloat hpToClipX_0(vec2 splitPos_0, vec3 bpRange_0, float hpZero_1)\n{\n\n    float floorBp_0 = -3.00000000549775576e+38 + hpZero_1;\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 covBpToClipX_0(uint bp_0, CoverageBandUniforms_0 u_1)\n{\n\n    return hpToClipX_0(hpSplitUint_0(bp_0), vec3(u_1.bpHi_0, u_1.bpLo_0, u_1.bpLen_0), u_1.hpZero_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 interbaseEdgePx_0(float offset_0, float interbaseHeight_1)\n{\n\n    return floor(4.5 + offset_0 * interbaseHeight_1 + 0.5);\n}\n\n\nfloat yPxToClipY_0(float yPx_0, float canvasHeight_0)\n{\n\n    return 1.0 - yPx_0 / canvasHeight_0 * 2.0;\n}\n\n\nfloat covAreaTop_0(CoverageBandUniforms_0 u_2)\n{\n\n    return yPxToClipY_0(u_2.covTop_0, u_2.canvasH_0);\n}\n\n\nfloat covFlipX_0(float x_0, CoverageBandUniforms_0 u_3)\n{\n\n    return mix(x_0, - x_0, u_3.reversed_0);\n}\n\n\nvec4 covFlippedQuad_0(float x0_0, float x1_0, float yBot_0, float yTop_0, vec2 local_0, CoverageBandUniforms_0 u_4)\n{\n    return vec4(covFlipX_0(mix(x0_0, x1_0, local_0.x), u_4), mix(yBot_0, yTop_0, local_0.y), 0.0, 1.0);\n}\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\nvec3 covClipKindColor_0(float ct_0, CoverageBandUniforms_0 u_5)\n{\n\n    int ci_0 = int(ct_0);\n    if(ci_0 == 2)\n    {\n\n        return unpackRGBA_0(u_5.colorSoftclipIndicator_0).xyz;\n    }\n\n    if(ci_0 == 3)\n    {\n\n        return unpackRGBA_0(u_5.colorHardclipIndicator_0).xyz;\n    }\n\n    return unpackRGBA_0(u_5.colorInsertionIndicator_0).xyz;\n}\n\n\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_position;\n\n\nlayout(location = 1)\nin float a_yOffset;\n\n\nlayout(location = 2)\nin float a_segHeight;\n\n\nlayout(location = 3)\nin float a_colorType;\n\n\nstruct CoverageVsOut_0\n{\n    vec4 position_0;\n    vec4 color_0;\n};\n\n\nvoid main()\n{\n\n    CoverageBandUniforms_0 _S5 = CoverageBandUniforms_0(u_0.bpHi_0, u_0.bpLo_0, u_0.bpLen_0, u_0.hpZero_0, u_0.canvasW_0, u_0.canvasH_0, u_0.covHeight_0, u_0.covYOffset_0, u_0.covTop_0, u_0.regionMaxDepth_0, u_0.depthDomainMax_0, u_0.depthDomainMin_0, u_0.coverageSymlogConstant_0, u_0.binSize_0, u_0.interbaseHeight_0, u_0.snpMinFreq_0, u_0.reversed_0, u_0.coverageScaleType_0, u_0.colorCoverage_0, u_0.colorBaseA_0, u_0.colorBaseC_0, u_0.colorBaseG_0, u_0.colorBaseT_0, u_0.colorBaseN_0, u_0.colorInsertionIndicator_0, u_0.colorSoftclipIndicator_0, u_0.colorHardclipIndicator_0);\n\n\n\n    float cx_0 = covBpToClipX_0(a_position, _S5);\n\n    float _S6 = covAreaTop_0(_S5);\n\n    CoverageVsOut_0 o_0;\n\n    float _S7 = pxToClipLen_0(1.0, u_0.canvasW_0) * 0.5;\n\n    o_0.position_0 = covFlippedQuad_0(cx_0 - _S7, cx_0 + _S7, _S6 - pxToClipLen_0(interbaseEdgePx_0(a_yOffset + a_segHeight, u_0.interbaseHeight_0), u_0.canvasH_0), _S6 - pxToClipLen_0(interbaseEdgePx_0(a_yOffset, u_0.interbaseHeight_0), u_0.canvasH_0), quadLocal_0(uint(gl_VertexID)), _S5);\n    o_0.color_0 = vec4(covClipKindColor_0(a_colorType, _S5), 1.0);\n    CoverageVsOut_0 _S8 = o_0;\n\n    gl_Position = o_0.position_0;\n\n    v_color = _S8.color_0;\n\n    return;\n}\n\n"
+export const GLSL_VERTEX = "#version 300 es\nprecision highp float;\nprecision highp int;\nstruct CoverageBandUniforms_0\n{\n    float bpHi_0;\n    float bpLo_0;\n    float bpLen_0;\n    float hpZero_0;\n    float canvasW_0;\n    float canvasH_0;\n    float covHeight_0;\n    float covYOffset_0;\n    float covTop_0;\n    float regionMaxDepth_0;\n    float depthDomainMax_0;\n    float depthDomainMin_0;\n    float coverageSymlogConstant_0;\n    float binSize_0;\n    float interbaseHeight_0;\n    float snpMinFreq_0;\n    float reversed_0;\n    int coverageScaleType_0;\n    uint colorCoverage_0;\n    uint colorBaseA_0;\n    uint colorBaseC_0;\n    uint colorBaseG_0;\n    uint colorBaseT_0;\n    uint colorBaseN_0;\n    uint colorInsertionIndicator_0;\n    uint colorSoftclipIndicator_0;\n    uint colorHardclipIndicator_0;\n};\n\n\nlayout(std140) uniform Uniforms\n{\n    float bpHi_0;\n    float bpLo_0;\n    float bpLen_0;\n    float hpZero_0;\n    float canvasW_0;\n    float canvasH_0;\n    float covHeight_0;\n    float covYOffset_0;\n    float covTop_0;\n    float regionMaxDepth_0;\n    float depthDomainMax_0;\n    float depthDomainMin_0;\n    float coverageSymlogConstant_0;\n    float binSize_0;\n    float interbaseHeight_0;\n    float snpMinFreq_0;\n    float reversed_0;\n    int coverageScaleType_0;\n    uint colorCoverage_0;\n    uint colorBaseA_0;\n    uint colorBaseC_0;\n    uint colorBaseG_0;\n    uint colorBaseT_0;\n    uint colorBaseN_0;\n    uint colorInsertionIndicator_0;\n    uint colorSoftclipIndicator_0;\n    uint colorHardclipIndicator_0;\n}u_0;\n\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S1 = v_0 == 0U;\n\n    bool _S2;\n\n    if(_S1)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 2U;\n\n    }\n\n    if(_S2)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 3U;\n\n    }\n\n    float _S3;\n\n    if(_S2)\n    {\n\n        _S3 = 0.0;\n\n    }\n    else\n    {\n\n        _S3 = 1.0;\n\n    }\n    if(_S1)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 1U;\n\n    }\n\n    if(_S2)\n    {\n\n        _S2 = true;\n\n    }\n    else\n    {\n\n        _S2 = v_0 == 4U;\n\n    }\n\n    float _S4;\n\n    if(_S2)\n    {\n\n        _S4 = 0.0;\n\n    }\n    else\n    {\n\n        _S4 = 1.0;\n\n    }\n\n    return vec2(_S3, _S4);\n}\n\n\nvec2 hpSplitUint_0(uint value_0)\n{\n\n    uint lo_0 = value_0 & 4095U;\n\n    return vec2(float(value_0 - lo_0), float(lo_0));\n}\n\n\nfloat hpToClipX_0(vec2 splitPos_0, vec3 bpRange_0, float hpZero_1)\n{\n\n    float floorBp_0 = -3.00000000549775576e+38 + hpZero_1;\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 covBpToClipX_0(uint bp_0, CoverageBandUniforms_0 u_1)\n{\n\n    return hpToClipX_0(hpSplitUint_0(bp_0), vec3(u_1.bpHi_0, u_1.bpLo_0, u_1.bpLen_0), u_1.hpZero_0);\n}\n\nfloat covFlipX_0(float x_0, CoverageBandUniforms_0 u_2)\n{\n\n    return mix(x_0, - x_0, u_2.reversed_0);\n}\n\n\nfloat clipXToPx_0(float clipX_0, float canvasWidth_0)\n{\n\n    return (clipX_0 + 1.0) * 0.5 * canvasWidth_0;\n}\n\n\nfloat interbaseBarLeftPx_0(float centerPx_0)\n{\n\n    return floor(centerPx_0);\n}\n\n\nfloat pxToClipX_0(float xPx_0, float canvasWidth_1)\n{\n\n    return xPx_0 / canvasWidth_1 * 2.0 - 1.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 interbaseEdgePx_0(float offset_0, float interbaseHeight_1)\n{\n\n    return floor(4.5) + floor(offset_0 * interbaseHeight_1 + 0.5);\n}\n\n\nfloat yPxToClipY_0(float yPx_0, float canvasHeight_0)\n{\n\n    return 1.0 - yPx_0 / canvasHeight_0 * 2.0;\n}\n\n\nfloat covAreaTop_0(CoverageBandUniforms_0 u_3)\n{\n\n    return yPxToClipY_0(u_3.covTop_0, u_3.canvasH_0);\n}\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\nvec3 covClipKindColor_0(float ct_0, CoverageBandUniforms_0 u_4)\n{\n\n    int ci_0 = int(ct_0);\n    if(ci_0 == 2)\n    {\n\n        return unpackRGBA_0(u_4.colorSoftclipIndicator_0).xyz;\n    }\n\n    if(ci_0 == 3)\n    {\n\n        return unpackRGBA_0(u_4.colorHardclipIndicator_0).xyz;\n    }\n\n    return unpackRGBA_0(u_4.colorInsertionIndicator_0).xyz;\n}\n\n\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_position;\n\n\nlayout(location = 1)\nin float a_yOffset;\n\n\nlayout(location = 2)\nin float a_segHeight;\n\n\nlayout(location = 3)\nin float a_colorType;\n\n\nstruct CoverageVsOut_0\n{\n    vec4 position_0;\n    vec4 color_0;\n};\n\n\nvoid main()\n{\n    vec2 local_0 = quadLocal_0(uint(gl_VertexID));\n\n    CoverageBandUniforms_0 _S5 = CoverageBandUniforms_0(u_0.bpHi_0, u_0.bpLo_0, u_0.bpLen_0, u_0.hpZero_0, u_0.canvasW_0, u_0.canvasH_0, u_0.covHeight_0, u_0.covYOffset_0, u_0.covTop_0, u_0.regionMaxDepth_0, u_0.depthDomainMax_0, u_0.depthDomainMin_0, u_0.coverageSymlogConstant_0, u_0.binSize_0, u_0.interbaseHeight_0, u_0.snpMinFreq_0, u_0.reversed_0, u_0.coverageScaleType_0, u_0.colorCoverage_0, u_0.colorBaseA_0, u_0.colorBaseC_0, u_0.colorBaseG_0, u_0.colorBaseT_0, u_0.colorBaseN_0, u_0.colorInsertionIndicator_0, u_0.colorSoftclipIndicator_0, u_0.colorHardclipIndicator_0);\n\n\n    float x0_0 = pxToClipX_0(interbaseBarLeftPx_0(clipXToPx_0(covFlipX_0(covBpToClipX_0(a_position, _S5), _S5), u_0.canvasW_0)), u_0.canvasW_0);\n\n    float _S6 = covAreaTop_0(_S5);\n\n    CoverageVsOut_0 o_0;\n\n    o_0.position_0 = vec4(mix(x0_0, x0_0 + pxToClipLen_0(1.0, u_0.canvasW_0), local_0.x), mix(_S6 - pxToClipLen_0(interbaseEdgePx_0(a_yOffset + a_segHeight, u_0.interbaseHeight_0), u_0.canvasH_0), _S6 - pxToClipLen_0(interbaseEdgePx_0(a_yOffset, u_0.interbaseHeight_0), u_0.canvasH_0), local_0.y), 0.0, 1.0);\n    o_0.color_0 = vec4(covClipKindColor_0(a_colorType, _S5), 1.0);\n    CoverageVsOut_0 _S7 = o_0;\n\n    gl_Position = o_0.position_0;\n\n    v_color = _S7.color_0;\n\n    return;\n}\n\n"
 
 export const GLSL_FRAGMENT = "#version 300 es\nprecision highp float;\nprecision highp int;\nlayout(location = 0)\nout vec4 entryPointParam_fs_main_0;\n\n\nin vec4 v_color;\n\n\nvoid main()\n{\n\n    entryPointParam_fs_main_0 = v_color;\n\n    return;\n}\n\n"

package/src/shaders/coverageInterbase.slang

@@ -1,9 +1,9 @@
 // Interbase histogram: slim 1px-wide bars indicating insertion / softclip /
-// hardclip density across the coverage band, hanging from just below the
-// indicator triangles. See coverageBand.slang.
+// hardclip density across the coverage band, hanging from the indicator
+// triangles' tip row. See coverageBand.slang.
 //! targets: wgsl, glsl
 //! layout-out: packages/alignments-core/src/interbaseHistogramLayout.generated.ts
-//! js-export: interbaseEdgePx
+//! js-export: interbaseEdgePx, interbaseBarLeftPx
 //! js-export-out: packages/alignments-core/src/interbaseEdge.generated.ts
 
 import coverageBand;
@@ -20,20 +20,24 @@ struct Instance {
 };
 
 // Screen-pixel Y of one edge of a bar, `offset` bars' worth of stacked height
-// below the coverage-area top.
+// below the coverage-area top, snapped to a whole pixel. The backends disagree
+// on a half-covered row: Canvas2D saturates the ~40 bars stacked in one column
+// at 40bp/px to opaque, while this pass resolves their union once and stays at
+// 50%. drawInterbaseSegments `//! js-export`s this so both sides round alike.
 //
-// Rounded to a whole pixel, and that rounding is the whole reason this is a
-// named function rather than two expressions inline. TRIANGLE_H is 4.5, so an
-// unsnapped edge half-covers a row — and the backends do not agree on a
-// half-covered row: Canvas2D composites each of the ~40 bars that stack in one
-// column at 40bp/px separately and saturates to opaque, while this pass resolves
-// the union coverage once and stays at 50%. That read as a full row of 128,0,128
-// against 191,127,191 and was the largest remaining cross-backend divergence.
-// Whole-pixel edges leave no partial coverage to disagree about — which only
-// holds while both sides round the same numbers the same way, so
-// drawInterbaseSegments `//! js-export`s this rather than repeating it.
+// The bar starts on the row holding the triangle's tip, floor(TRIANGLE_H),
+// rather than the row after it: the tip covers only a sliver of that row, so
+// starting below it leaves a visible gap between triangle and bar.
 float interbaseEdgePx(float offset, float interbaseHeight) {
-  return floor(TRIANGLE_H + offset * interbaseHeight + 0.5);
+  return floor(TRIANGLE_H) + floor(offset * interbaseHeight + 0.5);
+}
+
+// Left edge of the 1px bar centred on `centerPx`, snapped to the pixel column
+// whose centre is nearest. Unsnapped, the bar straddles two columns: Canvas2D
+// paints both at partial alpha while the GPU fills one solid, so the bar read
+// paler than its triangle and the backends disagreed.
+float interbaseBarLeftPx(float centerPx) {
+  return floor(centerPx);
 }
 
 [shader("vertex")]
@@ -41,18 +45,17 @@ CoverageVsOut vs_main(Instance inst, uint vid : SV_VertexID) {
   CoverageVsOut o;
   float2 local = quadLocal(vid);
 
-  float cx = covBpToClipX(inst.position, u);
-  float barW = pxToClipLen(1.0, u.canvasW);
+  float cx = covFlipX(covBpToClipX(inst.position, u), u);
+  float x0 = pxToClipX(interbaseBarLeftPx(clipXToPx(cx, u.canvasW)), u.canvasW);
+  float x1 = x0 + pxToClipLen(1.0, u.canvasW);
 
-  // Hang from the coverage-area top, starting just below the indicator
-  // triangle row (TRIANGLE_H is the same const coverageIndicator.slang draws
-  // with). Both edges snap to whole pixels BEFORE the clip-space conversion.
+  // Both edges snap to whole pixels BEFORE the clip-space conversion.
   float topPx = interbaseEdgePx(inst.yOffset, u.interbaseHeight);
   float botPx = interbaseEdgePx(inst.yOffset + inst.segHeight, u.interbaseHeight);
   float sTop = covAreaTop(u) - pxToClipLen(topPx, u.canvasH);
   float sBot = covAreaTop(u) - pxToClipLen(botPx, u.canvasH);
 
-  o.position = covFlippedQuad(cx - barW * 0.5, cx + barW * 0.5, sBot, sTop, local, u);
+  o.position = float4(lerp(x0, x1, local.x), lerp(sBot, sTop, local.y), 0.0, 1.0);
   o.color = float4(covClipKindColor(inst.colorType, u), 1.0);
   return o;
 }

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

@@ -1,4 +1,4 @@
 // AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageInterbase.slang.
 // Do not edit. Run `pnpm gen:shaders` to regenerate.
 
-export const WGSL_SOURCE = "struct CoverageBandUniforms_std140_0\n{\n    @align(16) bpHi_0 : f32,\n    @align(4) bpLo_0 : f32,\n    @align(8) bpLen_0 : f32,\n    @align(4) hpZero_0 : f32,\n    @align(16) canvasW_0 : f32,\n    @align(4) canvasH_0 : f32,\n    @align(8) covHeight_0 : f32,\n    @align(4) covYOffset_0 : f32,\n    @align(16) covTop_0 : f32,\n    @align(4) regionMaxDepth_0 : f32,\n    @align(8) depthDomainMax_0 : f32,\n    @align(4) depthDomainMin_0 : f32,\n    @align(16) coverageSymlogConstant_0 : f32,\n    @align(4) binSize_0 : f32,\n    @align(8) interbaseHeight_0 : f32,\n    @align(4) snpMinFreq_0 : f32,\n    @align(16) reversed_0 : f32,\n    @align(4) coverageScaleType_0 : i32,\n    @align(8) colorCoverage_0 : u32,\n    @align(4) colorBaseA_0 : u32,\n    @align(16) colorBaseC_0 : u32,\n    @align(4) colorBaseG_0 : u32,\n    @align(8) colorBaseT_0 : u32,\n    @align(4) colorBaseN_0 : u32,\n    @align(16) colorInsertionIndicator_0 : u32,\n    @align(4) colorSoftclipIndicator_0 : u32,\n    @align(8) colorHardclipIndicator_0 : u32,\n};\n\n@binding(1) @group(0) var<uniform> u_0 : CoverageBandUniforms_std140_0;\nfn quadLocal_0( vid_0 : u32) -> vec2<f32>\n{\n    var v_0 : u32 = vid_0 % u32(6);\n    var _S1 : bool = v_0 == u32(0);\n    var _S2 : bool;\n    if(_S1)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(2);\n    }\n    if(_S2)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(3);\n    }\n    var _S3 : f32;\n    if(_S2)\n    {\n        _S3 = 0.0f;\n    }\n    else\n    {\n        _S3 = 1.0f;\n    }\n    if(_S1)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(1);\n    }\n    if(_S2)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(4);\n    }\n    var _S4 : f32;\n    if(_S2)\n    {\n        _S4 = 0.0f;\n    }\n    else\n    {\n        _S4 = 1.0f;\n    }\n    return vec2<f32>(_S3, _S4);\n}\n\nfn hpSplitUint_0( value_0 : u32) -> vec2<f32>\n{\n    var lo_0 : u32 = (value_0 & (u32(4095)));\n    return vec2<f32>(f32(value_0 - lo_0), f32(lo_0));\n}\n\nfn hpToClipX_0( splitPos_0 : vec2<f32>,  bpRange_0 : vec3<f32>,  hpZero_1 : f32) -> f32\n{\n    var floorBp_0 : f32 = -3.00000000549775576e+38f + hpZero_1;\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 covBpToClipX_0( bp_0 : u32,  u_1 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return hpToClipX_0(hpSplitUint_0(bp_0), vec3<f32>((*u_1).bpHi_0, (*u_1).bpLo_0, (*u_1).bpLen_0), (*u_1).hpZero_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 interbaseEdgePx_0( offset_0 : f32,  interbaseHeight_1 : f32) -> f32\n{\n    return floor(4.5f + offset_0 * interbaseHeight_1 + 0.5f);\n}\n\nfn yPxToClipY_0( yPx_0 : f32,  canvasHeight_0 : f32) -> f32\n{\n    return 1.0f - yPx_0 / canvasHeight_0 * 2.0f;\n}\n\nfn covAreaTop_0( u_2 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return yPxToClipY_0((*u_2).covTop_0, (*u_2).canvasH_0);\n}\n\nfn covFlipX_0( x_0 : f32,  u_3 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return mix(x_0, - x_0, (*u_3).reversed_0);\n}\n\nfn covFlippedQuad_0( x0_0 : f32,  x1_0 : f32,  yBot_0 : f32,  yTop_0 : f32,  local_0 : vec2<f32>,  u_4 : ptr<function, CoverageBandUniforms_std140_0>) -> vec4<f32>\n{\n    var _S5 : f32 = covFlipX_0(mix(x0_0, x1_0, local_0.x), &((*u_4)));\n    return vec4<f32>(_S5, mix(yBot_0, yTop_0, local_0.y), 0.0f, 1.0f);\n}\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 covClipKindColor_0( ct_0 : f32,  u_5 : ptr<function, CoverageBandUniforms_std140_0>) -> vec3<f32>\n{\n    var ci_0 : i32 = i32(ct_0);\n    if(ci_0 == i32(2))\n    {\n        return unpackRGBA_0((*u_5).colorSoftclipIndicator_0).xyz;\n    }\n    if(ci_0 == i32(3))\n    {\n        return unpackRGBA_0((*u_5).colorHardclipIndicator_0).xyz;\n    }\n    return unpackRGBA_0((*u_5).colorInsertionIndicator_0).xyz;\n}\n\nstruct CoverageVsOut_0\n{\n    @builtin(position) position_0 : vec4<f32>,\n    @location(0) color_0 : vec4<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) position_1 : u32,\n    @location(1) yOffset_0 : f32,\n    @location(2) segHeight_0 : f32,\n    @location(3) colorType_0 : f32,\n};\n\n@vertex\nfn vs_main( _S6 : vertexInput_0, @builtin(vertex_index) vid_1 : u32) -> CoverageVsOut_0\n{\n    var local_1 : vec2<f32> = quadLocal_0(vid_1);\n    var _S7 : CoverageBandUniforms_std140_0 = u_0;\n    var _S8 : f32 = covBpToClipX_0(_S6.position_1, &(_S7));\n    var barW_0 : f32 = pxToClipLen_0(1.0f, u_0.canvasW_0);\n    var topPx_0 : f32 = interbaseEdgePx_0(_S6.yOffset_0, u_0.interbaseHeight_0);\n    var botPx_0 : f32 = interbaseEdgePx_0(_S6.yOffset_0 + _S6.segHeight_0, u_0.interbaseHeight_0);\n    var _S9 : f32 = covAreaTop_0(&(_S7));\n    var o_0 : CoverageVsOut_0;\n    var _S10 : f32 = barW_0 * 0.5f;\n    var _S11 : vec4<f32> = covFlippedQuad_0(_S8 - _S10, _S8 + _S10, _S9 - pxToClipLen_0(botPx_0, u_0.canvasH_0), _S9 - pxToClipLen_0(topPx_0, u_0.canvasH_0), local_1, &(_S7));\n    o_0.position_0 = _S11;\n    var _S12 : vec3<f32> = covClipKindColor_0(_S6.colorType_0, &(_S7));\n    o_0.color_0 = vec4<f32>(_S12, 1.0f);\n    return o_0;\n}\n\nstruct pixelOutput_0\n{\n    @location(0) output_0 : vec4<f32>,\n};\n\nstruct pixelInput_0\n{\n    @location(0) color_1 : vec4<f32>,\n};\n\n@fragment\nfn fs_main( _S13 : pixelInput_0, @builtin(position) position_2 : vec4<f32>) -> pixelOutput_0\n{\n    var _S14 : pixelOutput_0 = pixelOutput_0( _S13.color_1 );\n    return _S14;\n}\n\n"
+export const WGSL_SOURCE = "struct CoverageBandUniforms_std140_0\n{\n    @align(16) bpHi_0 : f32,\n    @align(4) bpLo_0 : f32,\n    @align(8) bpLen_0 : f32,\n    @align(4) hpZero_0 : f32,\n    @align(16) canvasW_0 : f32,\n    @align(4) canvasH_0 : f32,\n    @align(8) covHeight_0 : f32,\n    @align(4) covYOffset_0 : f32,\n    @align(16) covTop_0 : f32,\n    @align(4) regionMaxDepth_0 : f32,\n    @align(8) depthDomainMax_0 : f32,\n    @align(4) depthDomainMin_0 : f32,\n    @align(16) coverageSymlogConstant_0 : f32,\n    @align(4) binSize_0 : f32,\n    @align(8) interbaseHeight_0 : f32,\n    @align(4) snpMinFreq_0 : f32,\n    @align(16) reversed_0 : f32,\n    @align(4) coverageScaleType_0 : i32,\n    @align(8) colorCoverage_0 : u32,\n    @align(4) colorBaseA_0 : u32,\n    @align(16) colorBaseC_0 : u32,\n    @align(4) colorBaseG_0 : u32,\n    @align(8) colorBaseT_0 : u32,\n    @align(4) colorBaseN_0 : u32,\n    @align(16) colorInsertionIndicator_0 : u32,\n    @align(4) colorSoftclipIndicator_0 : u32,\n    @align(8) colorHardclipIndicator_0 : u32,\n};\n\n@binding(1) @group(0) var<uniform> u_0 : CoverageBandUniforms_std140_0;\nfn quadLocal_0( vid_0 : u32) -> vec2<f32>\n{\n    var v_0 : u32 = vid_0 % u32(6);\n    var _S1 : bool = v_0 == u32(0);\n    var _S2 : bool;\n    if(_S1)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(2);\n    }\n    if(_S2)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(3);\n    }\n    var _S3 : f32;\n    if(_S2)\n    {\n        _S3 = 0.0f;\n    }\n    else\n    {\n        _S3 = 1.0f;\n    }\n    if(_S1)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(1);\n    }\n    if(_S2)\n    {\n        _S2 = true;\n    }\n    else\n    {\n        _S2 = v_0 == u32(4);\n    }\n    var _S4 : f32;\n    if(_S2)\n    {\n        _S4 = 0.0f;\n    }\n    else\n    {\n        _S4 = 1.0f;\n    }\n    return vec2<f32>(_S3, _S4);\n}\n\nfn hpSplitUint_0( value_0 : u32) -> vec2<f32>\n{\n    var lo_0 : u32 = (value_0 & (u32(4095)));\n    return vec2<f32>(f32(value_0 - lo_0), f32(lo_0));\n}\n\nfn hpToClipX_0( splitPos_0 : vec2<f32>,  bpRange_0 : vec3<f32>,  hpZero_1 : f32) -> f32\n{\n    var floorBp_0 : f32 = -3.00000000549775576e+38f + hpZero_1;\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 covBpToClipX_0( bp_0 : u32,  u_1 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return hpToClipX_0(hpSplitUint_0(bp_0), vec3<f32>((*u_1).bpHi_0, (*u_1).bpLo_0, (*u_1).bpLen_0), (*u_1).hpZero_0);\n}\n\nfn covFlipX_0( x_0 : f32,  u_2 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return mix(x_0, - x_0, (*u_2).reversed_0);\n}\n\nfn clipXToPx_0( clipX_0 : f32,  canvasWidth_0 : f32) -> f32\n{\n    return (clipX_0 + 1.0f) * 0.5f * canvasWidth_0;\n}\n\nfn interbaseBarLeftPx_0( centerPx_0 : f32) -> f32\n{\n    return floor(centerPx_0);\n}\n\nfn pxToClipX_0( xPx_0 : f32,  canvasWidth_1 : f32) -> f32\n{\n    return xPx_0 / canvasWidth_1 * 2.0f - 1.0f;\n}\n\nfn pxToClipLen_0( px_0 : f32,  size_0 : f32) -> f32\n{\n    return px_0 * 2.0f / size_0;\n}\n\nfn interbaseEdgePx_0( offset_0 : f32,  interbaseHeight_1 : f32) -> f32\n{\n    return floor(4.5f) + floor(offset_0 * interbaseHeight_1 + 0.5f);\n}\n\nfn yPxToClipY_0( yPx_0 : f32,  canvasHeight_0 : f32) -> f32\n{\n    return 1.0f - yPx_0 / canvasHeight_0 * 2.0f;\n}\n\nfn covAreaTop_0( u_3 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return yPxToClipY_0((*u_3).covTop_0, (*u_3).canvasH_0);\n}\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 covClipKindColor_0( ct_0 : f32,  u_4 : ptr<function, CoverageBandUniforms_std140_0>) -> vec3<f32>\n{\n    var ci_0 : i32 = i32(ct_0);\n    if(ci_0 == i32(2))\n    {\n        return unpackRGBA_0((*u_4).colorSoftclipIndicator_0).xyz;\n    }\n    if(ci_0 == i32(3))\n    {\n        return unpackRGBA_0((*u_4).colorHardclipIndicator_0).xyz;\n    }\n    return unpackRGBA_0((*u_4).colorInsertionIndicator_0).xyz;\n}\n\nstruct CoverageVsOut_0\n{\n    @builtin(position) position_0 : vec4<f32>,\n    @location(0) color_0 : vec4<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) position_1 : u32,\n    @location(1) yOffset_0 : f32,\n    @location(2) segHeight_0 : f32,\n    @location(3) colorType_0 : f32,\n};\n\n@vertex\nfn vs_main( _S5 : vertexInput_0, @builtin(vertex_index) vid_1 : u32) -> CoverageVsOut_0\n{\n    var local_0 : vec2<f32> = quadLocal_0(vid_1);\n    var _S6 : CoverageBandUniforms_std140_0 = u_0;\n    var _S7 : f32 = covBpToClipX_0(_S5.position_1, &(_S6));\n    var _S8 : f32 = covFlipX_0(_S7, &(_S6));\n    var x0_0 : f32 = pxToClipX_0(interbaseBarLeftPx_0(clipXToPx_0(_S8, u_0.canvasW_0)), u_0.canvasW_0);\n    var x1_0 : f32 = x0_0 + pxToClipLen_0(1.0f, u_0.canvasW_0);\n    var topPx_0 : f32 = interbaseEdgePx_0(_S5.yOffset_0, u_0.interbaseHeight_0);\n    var botPx_0 : f32 = interbaseEdgePx_0(_S5.yOffset_0 + _S5.segHeight_0, u_0.interbaseHeight_0);\n    var _S9 : f32 = covAreaTop_0(&(_S6));\n    var o_0 : CoverageVsOut_0;\n    o_0.position_0 = vec4<f32>(mix(x0_0, x1_0, local_0.x), mix(_S9 - pxToClipLen_0(botPx_0, u_0.canvasH_0), _S9 - pxToClipLen_0(topPx_0, u_0.canvasH_0), local_0.y), 0.0f, 1.0f);\n    var _S10 : vec3<f32> = covClipKindColor_0(_S5.colorType_0, &(_S6));\n    o_0.color_0 = vec4<f32>(_S10, 1.0f);\n    return o_0;\n}\n\nstruct pixelOutput_0\n{\n    @location(0) output_0 : vec4<f32>,\n};\n\nstruct pixelInput_0\n{\n    @location(0) color_1 : vec4<f32>,\n};\n\n@fragment\nfn fs_main( _S11 : pixelInput_0, @builtin(position) position_2 : vec4<f32>) -> pixelOutput_0\n{\n    var _S12 : pixelOutput_0 = pixelOutput_0( _S11.color_1 );\n    return _S12;\n}\n\n"

package/src/shaders/coverageSnp.slang

@@ -15,7 +15,7 @@ public static const uint VERTS_PER_INSTANCE = 6u;
 // via covNormalizeDepth, so each segment is a linear slice of a possibly
 // log-scaled bar without any per-position un-mixing at draw time.
 //
-// Which means that on a LOG axis the coloured fraction reads as the allele
+// Which means that on a LOG axis the colored fraction reads as the allele
 // proportion and not as a count off the y-axis: a 50% allele is half the bar
 // whatever the scale, so its top edge does not land on the depth the axis would
 // read there. Deliberate — the alternative bakes the axis into the buffer and