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

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)
+}

package/src/shaders/colorRampLut.slang

@@ -5,7 +5,8 @@
 // remap plus `ldValueComputed` sentinel gate, wiggle's `normalizeScore` +
 // `densityGradientT`) and declares its own `[[vk::binding(2, 0)]]
 // Sampler2D<float4> colorRamp;` — a module can take a sampler as a parameter
-// but not declare the binding for its importer. The TS half of the mechanism is
+// but not declare the binding for its importer, and inherits the `linear`
+// filter the remap below is written against. The TS half of the mechanism is
 // `uploadColorRampLut` (../colorRampLut.ts), which puts the 256-entry
 // `Uint8Array` behind these samplers; `buildColorRampLut`
 // (@jbrowse/core/util/colorRamp) is what builds it on both backends.
@@ -18,9 +19,22 @@
 // default straight-alpha blend, where a premultiplied return would be wrong
 // the day a translucent ramp entry appears. The function name carries the
 // contract, so a consumer picks the form that matches its blend declaration.
+//! texture-filter: colorRamp linear
 //! export-consts: RAMP_LUT_ENTRIES
+//! js-export: rampMidT
 module colorRampLut;
 
+// Where a normalized value `t` reads the ramp when its middle stop sits at
+// `mid`: both sides of the middle on one scale, the farther end reaching its
+// end stop and the nearer stopping short, so equal distances from the middle
+// take equal colours — ggplot2's `rescale_mid`. A `mid` of 0.5 reads `t`
+// itself, exactly, so a ramp with no middle declared reads its table as it
+// always did, and a moving domain moves one uniform and no table.
+public float rampMidT(float t, float mid) {
+  float far = max(mid, 1.0 - mid);
+  return clamp((t + (far - mid)) / (2.0 * far), 0.0, 1.0);
+}
+
 // Entries in the ramp texture. Exported because `buildColorRampLut` sizes the
 // array by it and the mapping below divides by it: a table built to one number
 // and sampled by another is a half-shade of drift at every entry.

package/src/shaders/coverageBand.slang

@@ -21,10 +21,10 @@
 // have to land on the same pixels (adr-051). They live in
 // @jbrowse/alignments-core because that is where those painters are, and it
 // cannot depend on a plugin — nor, being upstream of it, on this package.
-//! js-export: covEffectiveHeightPx, covBottomOffsetPx, normalizeDepthScalar
+//! js-export: covEffectiveHeightPx, covBottomOffsetPx, normalizeDepthScalar, covBarHeightPx, covSegBottomPx, covSegTopPx
 //! js-export-out: packages/alignments-core/src/coverageBandLayout.generated.ts
 //! consts-out: packages/alignments-core/src/indicatorTriangle.generated.ts
-//! js-skip: covExpandMinWidthX — a 1 CSS px floor over hpmath's expandToMinWidthX, which is where the rule and the reason it has no Canvas2D twin are recorded
+//! js-skip: covExpandMinWidthX — the clip-space half of the rule; nothing outside a shader is in clip space, so what a painter wants is the px core, and coverageBar.slang exports that (hpmath's expandToMinWidthLeftPx and RightPx, into alignments-core's spanMinWidth.generated.ts) for fillSpanRect to place its edges with
 module coverageBand;
 
 import colorPack;
@@ -37,8 +37,8 @@ public struct CoverageBandUniforms {
   public float bpHi;
   public float bpLo;
   public float bpLen;
-  // Runtime 0.0 sentinel. HP math uses 1/hpZero to materialize +inf without the
-  // compiler folding the hi/lo combination. MUST always be uploaded as 0.0.
+  // Runtime 0.0 sentinel that keeps the HP clamp floor off the compiler's
+  // constant table; see hpmath.slang. MUST always be uploaded as 0.0.
   public float hpZero;
 
   public float canvasW;
@@ -110,12 +110,23 @@ public struct CoverageVsOut {
   public float4 color    : COLOR0;
 };
 
-// Fold a vertex to the origin (degenerate, invisible), for an instance a gate
-// dropped. Folding rather than shrinking is what leaves the rest of the stack
-// where it was.
+// One depth bar. `relDepth` is depths[bin] / regionMaxDepth, the same packing
+// the SNP and modification segments carry (see coverageSnp.slang);
+// covNormalizeDepth recovers the raw depth from it. The bin spans `u.binSize`
+// bp: 1 for per-bp coverage, wider once the producer downsamples. The bar pass
+// and the clip-strip pass both read this buffer.
+public struct CoverageBarInstance {
+  public uint  position : ATTR0;
+  public float relDepth : ATTR1;
+};
+
+// Send a vertex off clip, for an instance a gate dropped. Dropping the whole
+// instance rather than shrinking its segment is what leaves the rest of the
+// stack where it was. `spanMark` spells its hidden row the same way: a real
+// `w`, so the perspective divide stays a division.
 public CoverageVsOut covDiscardVertex() {
   CoverageVsOut o;
-  o.position = float4(0.0);
+  o.position = float4(2.0, 2.0, 0.0, 1.0);
   o.color = float4(0.0);
   return o;
 }
@@ -186,21 +197,16 @@ public float covAreaTop(CoverageBandUniforms u) {
   return yPxToClipY(u.covTop, u.canvasH);
 }
 
-public float covBottom(CoverageBandUniforms u) {
-  return covAreaTop(u)
-    - (covBottomOffsetPx(u.covHeight, u.covYOffset) / u.canvasH) * 2.0;
+// Screen-px Y of the bars' baseline: covBottomOffsetPx below wherever covTop
+// put the band's top edge.
+public float covBaselinePx(CoverageBandUniforms u) {
+  return u.covTop + covBottomOffsetPx(u.covHeight, u.covYOffset);
 }
 
-// Normalize a raw coverage depth to [0,1] over the domain [min, max], on the
-// current scale type — `scoreScale.slang`'s normalizer under the coverage
-// band's name for it.
-//
-// The name is what the alignments coverage reads it by: it is the scalar twin
-// of makeScoreNormalizer in packages/wiggle-core/src/normalize.ts, which is
-// what the Canvas2D coverage draws and the y-axis ticks go through, and
-// `coverageNormalizeParity.test.ts` pins all of those against each other. The
-// three branches themselves are no longer here — this band and the wiggle
-// renderer had a copy each and they had drifted (see scoreScale.slang).
+// `scoreScale.slang`'s normalizer under the name the alignments coverage reads
+// it by: the scalar twin `coverageNormalizeParity.test.ts` pins against
+// `makeScoreNormalizer` (render-core `scoreScale.ts`), which the Canvas2D
+// coverage draws and the y-axis ticks go through.
 public float normalizeDepthScalar(float rawDepth, float domainMin, float domainMax, int scaleType, float symlogConstant) {
   return normalizeScore(rawDepth, domainMin, domainMax, scaleType, symlogConstant);
 }
@@ -213,7 +219,11 @@ public float normalizeDepthScalar(float rawDepth, float domainMin, float domainM
 // `getNiceDomain` returns verbatim for a typed `Set max score` of 0, made the
 // first factor +inf and the product NaN, reaching `clamp(NaN, 0, 1)`.
 // Split in two so the scalar half can be `//! js-export`ed — the twin emitter
-// cannot carry a uniform struct across.
+// cannot carry a uniform struct across. Hoisting the domain's transformed
+// endpoints into uniforms was declined: it saves two `log2` per vertex on the
+// two scale types that are not the default, and costs a second spelling of
+// `normalizeScore` that `normalizeDepthScalar`, ADR-097's pinned oracle for
+// this path, would no longer be on.
 public float covNormalizeDepth(float relDepth, CoverageBandUniforms u) {
   float rawDepth = relDepth * u.regionMaxDepth;
   return normalizeDepthScalar(
@@ -221,15 +231,32 @@ public float covNormalizeDepth(float relDepth, CoverageBandUniforms u) {
     u.coverageSymlogConstant);
 }
 
-// Bar height scale: maps relDepth to a clip-space height multiplier. Shared by
-// every stacked-segment pass so the formula can't drift between them.
-public float covBarScale(float relDepth, CoverageBandUniforms u) {
-  return covNormalizeDepth(relDepth, u) * covEffHeight(u) / u.canvasH * 2.0;
+// How tall a bar of this normalized depth is drawn, in px — what every stacked
+// fraction in the band is a fraction OF.
+public float covBarHeightPx(float depthFraction, float effectiveHeightPx) {
+  return depthFraction * effectiveHeightPx;
+}
+
+// The two y edges of one stacked segment, in px. `yOffset` and `segHeight` are
+// fractions of THIS position's bar, measured up from the baseline.
+//
+// Two scalars rather than one `float2`, unlike `rectSpanPx` / `linkRadiiPx`: a
+// vec2 twin comes back as a `[number, number]` (adr-051's tuple convention) and
+// these feed loops that run per covered bp under a no-allocation bar.
+public float covSegBottomPx(float baselinePx, float yOffset,
+                            float barHeightPx) {
+  return baselinePx - yOffset * barHeightPx;
+}
+
+public float covSegTopPx(float baselinePx, float yOffset, float segHeight,
+                         float barHeightPx) {
+  return covSegBottomPx(baselinePx, yOffset, barHeightPx)
+    - segHeight * barHeightPx;
 }
 
 // Stacked-segment quad — the geometry of EVERY bar in the band. A `widthBp`-wide
 // rect whose bottom/top are the yOffset / yOffset+segHeight fractions of this
-// position's coverage bar (covBarScale(relDepth)).
+// position's coverage bar.
 //
 // The three passes differ only in what they put in those parameters. The SNP and
 // modification segments are 1bp slices of a bar (widthBp 1, a partial y span);
@@ -238,16 +265,21 @@ public float covBarScale(float relDepth, CoverageBandUniforms u) {
 // so the band could shift under its own segments without anything catching it.
 // Passing `0.0, 1.0` is not a special case here, it is what "the whole bar"
 // spells.
+//
+// The y half is decided in px and converted once, the shape adr-051 asks for —
+// so the three Canvas2D painters of this band call the same two functions.
 public float4 covSegQuad(uint position, uint widthBp, float yOffset,
                          float segHeight, float relDepth, float2 local,
                          CoverageBandUniforms u) {
   float2 xx = covExpandMinWidthX(
     covBpToClipX(position, u), covBpToClipX(position + widthBp, u), u.canvasW);
-  float bottom = covBottom(u);
-  float k = covBarScale(relDepth, u);
-  float segBot = bottom + yOffset * k;
-  float segTop = segBot + segHeight * k;
-  return covFlippedQuad(xx.x, xx.y, segBot, segTop, local, u);
+  float basePx = covBaselinePx(u);
+  float barPx = covBarHeightPx(covNormalizeDepth(relDepth, u), covEffHeight(u));
+  return covFlippedQuad(
+    xx.x, xx.y,
+    yPxToClipY(covSegBottomPx(basePx, yOffset, barPx), u.canvasH),
+    yPxToClipY(covSegTopPx(basePx, yOffset, segHeight, barPx), u.canvasH),
+    local, u);
 }
 
 // SNP segment fill. colorType encodes 1=A 2=C 3=G 4=T 5=N; 5 and any other type

package/src/shaders/coverageBar.generated.ts

@@ -1,10 +1,11 @@
 // AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageBar.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 step_0 : f32 = 2.0f / bpRange_0.z;\n    var _S5 : f32 = - (1.0f / hpZero_1);\n    return dot(vec3<f32>(-1.0f, max(splitPos_0.x - bpRange_0.x, _S5), max(splitPos_0.y - bpRange_0.y, _S5)), 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 expandToMinWidthPx_0( x1_0 : f32,  x2_0 : f32,  minWidth_0 : f32) -> vec2<f32>\n{\n    if((x2_0 - x1_0) < minWidth_0)\n    {\n        var mid_0 : f32 = (x1_0 + x2_0) * 0.5f;\n        var _S6 : f32 = minWidth_0 * 0.5f;\n        return vec2<f32>(mid_0 - _S6, mid_0 + _S6);\n    }\n    return vec2<f32>(x1_0, x2_0);\n}\n\nfn expandToMinWidthX_0( x1_1 : f32,  x2_1 : f32,  minWidthPx_0 : f32,  canvasWidth_0 : f32) -> vec2<f32>\n{\n    return expandToMinWidthPx_0(x1_1, x2_1, minWidthPx_0 * 2.0f / canvasWidth_0);\n}\n\nfn covExpandMinWidthX_0( x1_2 : f32,  x2_2 : f32,  canvasWidth_1 : f32) -> vec2<f32>\n{\n    return expandToMinWidthX_0(x1_2, x2_2, 1.0f, canvasWidth_1);\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 covBottomOffsetPx_0( covHeight_1 : f32,  covYOffset_1 : f32) -> f32\n{\n    return covHeight_1 - covYOffset_1;\n}\n\nfn covBottom_0( u_3 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    var _S7 : f32 = covAreaTop_0(&((*u_3)));\n    return _S7 - covBottomOffsetPx_0((*u_3).covHeight_0, (*u_3).covYOffset_0) / (*u_3).canvasH_0 * 2.0f;\n}\n\nfn symlogTransform_0( x_0 : f32,  c_0 : f32) -> f32\n{\n    return f32((i32(sign((x_0))))) * log(1.0f + abs(x_0 / c_0));\n}\n\nfn normalizeScore_0( score_0 : f32,  domainMin_0 : f32,  domainMax_0 : f32,  scaleType_0 : i32,  symlogConstant_0 : f32) -> f32\n{\n    if(scaleType_0 == i32(2))\n    {\n        var tMin_0 : f32 = symlogTransform_0(domainMin_0, symlogConstant_0);\n        var tRange_0 : f32 = symlogTransform_0(domainMax_0, symlogConstant_0) - tMin_0;\n        if(tRange_0 <= 0.0f)\n        {\n            return 0.0f;\n        }\n        return clamp((symlogTransform_0(score_0, symlogConstant_0) - tMin_0) / tRange_0, 0.0f, 1.0f);\n    }\n    if(scaleType_0 == i32(1))\n    {\n        var floorV_0 : f32;\n        if(domainMin_0 > 0.0f)\n        {\n            floorV_0 = domainMin_0;\n        }\n        else\n        {\n            floorV_0 = 1.0f;\n        }\n        var logMin_0 : f32 = log2(floorV_0);\n        var logRange_0 : f32 = log2(max(domainMax_0, floorV_0)) - logMin_0;\n        if(logRange_0 <= 0.0f)\n        {\n            return 0.0f;\n        }\n        return clamp((log2(max(score_0, floorV_0)) - logMin_0) / logRange_0, 0.0f, 1.0f);\n    }\n    var range_0 : f32 = domainMax_0 - domainMin_0;\n    if(range_0 <= 0.0f)\n    {\n        return 0.0f;\n    }\n    return clamp((score_0 - domainMin_0) / range_0, 0.0f, 1.0f);\n}\n\nfn normalizeDepthScalar_0( rawDepth_0 : f32,  domainMin_1 : f32,  domainMax_1 : f32,  scaleType_1 : i32,  symlogConstant_1 : f32) -> f32\n{\n    return normalizeScore_0(rawDepth_0, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1);\n}\n\nfn covNormalizeDepth_0( relDepth_0 : f32,  u_4 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return normalizeDepthScalar_0(relDepth_0 * (*u_4).regionMaxDepth_0, (*u_4).depthDomainMin_0, (*u_4).depthDomainMax_0, (*u_4).coverageScaleType_0, (*u_4).coverageSymlogConstant_0);\n}\n\nfn covEffectiveHeightPx_0( covHeight_2 : f32,  covYOffset_2 : f32) -> f32\n{\n    return max(covHeight_2 - 2.0f * covYOffset_2, 0.0f);\n}\n\nfn covEffHeight_0( u_5 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return covEffectiveHeightPx_0((*u_5).covHeight_0, (*u_5).covYOffset_0);\n}\n\nfn covBarScale_0( relDepth_1 : f32,  u_6 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    var _S8 : f32 = covNormalizeDepth_0(relDepth_1, &((*u_6)));\n    var _S9 : f32 = covEffHeight_0(&((*u_6)));\n    return _S8 * _S9 / (*u_6).canvasH_0 * 2.0f;\n}\n\nfn covFlipX_0( x_1 : f32,  u_7 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return mix(x_1, - x_1, (*u_7).reversed_0);\n}\n\nfn covFlippedQuad_0( x0_0 : f32,  x1_3 : f32,  yBot_0 : f32,  yTop_0 : f32,  local_0 : vec2<f32>,  u_8 : ptr<function, CoverageBandUniforms_std140_0>) -> vec4<f32>\n{\n    var _S10 : f32 = covFlipX_0(mix(x0_0, x1_3, local_0.x), &((*u_8)));\n    return vec4<f32>(_S10, mix(yBot_0, yTop_0, local_0.y), 0.0f, 1.0f);\n}\n\nfn covSegQuad_0( position_0 : u32,  widthBp_0 : u32,  yOffset_0 : f32,  segHeight_0 : f32,  relDepth_2 : f32,  local_1 : vec2<f32>,  u_9 : ptr<function, CoverageBandUniforms_std140_0>) -> vec4<f32>\n{\n    var _S11 : f32 = covBpToClipX_0(position_0, &((*u_9)));\n    var _S12 : f32 = covBpToClipX_0(position_0 + widthBp_0, &((*u_9)));\n    var xx_0 : vec2<f32> = covExpandMinWidthX_0(_S11, _S12, (*u_9).canvasW_0);\n    var _S13 : f32 = covBottom_0(&((*u_9)));\n    var _S14 : f32 = covBarScale_0(relDepth_2, &((*u_9)));\n    var segBot_0 : f32 = _S13 + yOffset_0 * _S14;\n    var _S15 : vec4<f32> = covFlippedQuad_0(xx_0.x, xx_0.y, segBot_0, segBot_0 + segHeight_0 * _S14, local_1, &((*u_9)));\n    return _S15;\n}\n\nfn unpackRGBA_0( c_1 : u32) -> vec4<f32>\n{\n    return vec4<f32>(f32((((c_1 >> (u32(0)))) & (u32(255)))), f32((((c_1 >> (u32(8)))) & (u32(255)))), f32((((c_1 >> (u32(16)))) & (u32(255)))), f32((((c_1 >> (u32(24)))) & (u32(255))))) / vec4<f32>(255.0f);\n}\n\nstruct CoverageVsOut_0\n{\n    @builtin(position) position_1 : vec4<f32>,\n    @location(0) color_0 : vec4<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) position_2 : u32,\n    @location(1) relDepth_3 : f32,\n};\n\n@vertex\nfn vs_main( _S16 : vertexInput_0, @builtin(vertex_index) vid_1 : u32) -> CoverageVsOut_0\n{\n    var o_0 : CoverageVsOut_0;\n    var _S17 : u32 = u32(u_0.binSize_0);\n    var _S18 : vec2<f32> = quadLocal_0(vid_1);\n    var _S19 : CoverageBandUniforms_std140_0 = u_0;\n    var _S20 : vec4<f32> = covSegQuad_0(_S16.position_2, _S17, 0.0f, 1.0f, _S16.relDepth_3, _S18, &(_S19));\n    o_0.position_1 = _S20;\n    o_0.color_0 = vec4<f32>(unpackRGBA_0(u_0.colorCoverage_0).xyz, 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( _S21 : pixelInput_0, @builtin(position) position_3 : vec4<f32>) -> pixelOutput_0\n{\n    var _S22 : pixelOutput_0 = pixelOutput_0( _S21.color_1 );\n    return _S22;\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    float step_0 = 2.0 / bpRange_0.z;\n    float _S5 = - (1.0 / hpZero_1);\n\n    return dot(vec3(-1.0, max(splitPos_0.x - bpRange_0.x, _S5), max(splitPos_0.y - bpRange_0.y, _S5)), 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\nvec2 expandToMinWidthPx_0(float x1_0, float x2_0, float minWidth_0)\n{\n\n    if((x2_0 - x1_0) < minWidth_0)\n    {\n\n        float mid_0 = (x1_0 + x2_0) * 0.5;\n        float _S6 = minWidth_0 * 0.5;\n\n        return vec2(mid_0 - _S6, mid_0 + _S6);\n    }\n    return vec2(x1_0, x2_0);\n}\n\nvec2 expandToMinWidthX_0(float x1_1, float x2_1, float minWidthPx_0, float canvasWidth_0)\n{\n\n    return expandToMinWidthPx_0(x1_1, x2_1, minWidthPx_0 * 2.0 / canvasWidth_0);\n}\n\n\nvec2 covExpandMinWidthX_0(float x1_2, float x2_2, float canvasWidth_1)\n{\n\n    return expandToMinWidthX_0(x1_2, x2_2, 1.0, canvasWidth_1);\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 covBottomOffsetPx_0(float covHeight_1, float covYOffset_1)\n{\n\n    return covHeight_1 - covYOffset_1;\n}\n\n\nfloat covBottom_0(CoverageBandUniforms_0 u_3)\n{\n\n    return covAreaTop_0(u_3) - covBottomOffsetPx_0(u_3.covHeight_0, u_3.covYOffset_0) / u_3.canvasH_0 * 2.0;\n}\n\n\nfloat symlogTransform_0(float x_0, float c_0)\n{\n\n    return float((int(sign((x_0))))) * log(1.0 + abs(x_0 / c_0));\n}\n\nfloat normalizeScore_0(float score_0, float domainMin_0, float domainMax_0, int scaleType_0, float symlogConstant_0)\n{\n\n    if(scaleType_0 == 2)\n    {\n\n        float tMin_0 = symlogTransform_0(domainMin_0, symlogConstant_0);\n\n        float tRange_0 = symlogTransform_0(domainMax_0, symlogConstant_0) - tMin_0;\n        if(tRange_0 <= 0.0)\n        {\n\n            return 0.0;\n        }\n        return clamp((symlogTransform_0(score_0, symlogConstant_0) - tMin_0) / tRange_0, 0.0, 1.0);\n    }\n    if(scaleType_0 == 1)\n    {\n\n        float floorV_0;\n\n        if(domainMin_0 > 0.0)\n        {\n\n            floorV_0 = domainMin_0;\n\n        }\n        else\n        {\n\n            floorV_0 = 1.0;\n\n        }\n        float logMin_0 = log2(floorV_0);\n\n        float logRange_0 = log2(max(domainMax_0, floorV_0)) - logMin_0;\n        if(logRange_0 <= 0.0)\n        {\n\n            return 0.0;\n        }\n        return clamp((log2(max(score_0, floorV_0)) - logMin_0) / logRange_0, 0.0, 1.0);\n    }\n    float range_0 = domainMax_0 - domainMin_0;\n    if(range_0 <= 0.0)\n    {\n\n        return 0.0;\n    }\n\n    return clamp((score_0 - domainMin_0) / range_0, 0.0, 1.0);\n}\n\n\nfloat normalizeDepthScalar_0(float rawDepth_0, float domainMin_1, float domainMax_1, int scaleType_1, float symlogConstant_1)\n{\n\n    return normalizeScore_0(rawDepth_0, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1);\n}\n\n\nfloat covNormalizeDepth_0(float relDepth_0, CoverageBandUniforms_0 u_4)\n{\n    return normalizeDepthScalar_0(relDepth_0 * u_4.regionMaxDepth_0, u_4.depthDomainMin_0, u_4.depthDomainMax_0, u_4.coverageScaleType_0, u_4.coverageSymlogConstant_0);\n}\n\n\nfloat covEffectiveHeightPx_0(float covHeight_2, float covYOffset_2)\n{\n\n    return max(covHeight_2 - 2.0 * covYOffset_2, 0.0);\n}\n\n\nfloat covEffHeight_0(CoverageBandUniforms_0 u_5)\n{\n\n    return covEffectiveHeightPx_0(u_5.covHeight_0, u_5.covYOffset_0);\n}\n\n\nfloat covBarScale_0(float relDepth_1, CoverageBandUniforms_0 u_6)\n{\n\n    return covNormalizeDepth_0(relDepth_1, u_6) * covEffHeight_0(u_6) / u_6.canvasH_0 * 2.0;\n}\n\n\nfloat covFlipX_0(float x_1, CoverageBandUniforms_0 u_7)\n{\n\n    return mix(x_1, - x_1, u_7.reversed_0);\n}\n\n\nvec4 covFlippedQuad_0(float x0_0, float x1_3, float yBot_0, float yTop_0, vec2 local_0, CoverageBandUniforms_0 u_8)\n{\n    return vec4(covFlipX_0(mix(x0_0, x1_3, local_0.x), u_8), mix(yBot_0, yTop_0, local_0.y), 0.0, 1.0);\n}\n\n\nvec4 covSegQuad_0(uint position_0, uint widthBp_0, float yOffset_0, float segHeight_0, float relDepth_2, vec2 local_1, CoverageBandUniforms_0 u_9)\n{\n\n    vec2 xx_0 = covExpandMinWidthX_0(covBpToClipX_0(position_0, u_9), covBpToClipX_0(position_0 + widthBp_0, u_9), u_9.canvasW_0);\n\n\n    float k_0 = covBarScale_0(relDepth_2, u_9);\n    float segBot_0 = covBottom_0(u_9) + yOffset_0 * k_0;\n\n    return covFlippedQuad_0(xx_0.x, xx_0.y, segBot_0, segBot_0 + segHeight_0 * k_0, local_1, u_9);\n}\n\n\nvec4 unpackRGBA_0(uint c_1)\n{\n\n    return vec4(float((c_1 >> 0U) & 255U), float((c_1 >> 8U) & 255U), float((c_1 >> 16U) & 255U), float((c_1 >> 24U) & 255U)) / 255.0;\n}\n\n\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_position;\n\n\nlayout(location = 1)\nin float a_relDepth;\n\n\nstruct CoverageVsOut_0\n{\n    vec4 position_1;\n    vec4 color_0;\n};\n\n\nvoid main()\n{\n\n    CoverageVsOut_0 o_0;\n\n    CoverageBandUniforms_0 _S7 = 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    o_0.position_1 = covSegQuad_0(a_position, uint(u_0.binSize_0), 0.0, 1.0, a_relDepth, quadLocal_0(uint(gl_VertexID)), _S7);\n\n    o_0.color_0 = vec4(unpackRGBA_0(u_0.colorCoverage_0).xyz, 1.0);\n    CoverageVsOut_0 _S8 = o_0;\n\n    gl_Position = o_0.position_1;\n\n    v_color = _S8.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"
+import type { ShaderSource } from '@jbrowse/render-core/hal'
 
 export * from './coverageBar.iface.generated.ts'
+
+export const SOURCE: ShaderSource = {
+  wgsl: () => import('./coverageBar.wgsl.generated.ts'),
+  glsl: () => import('./coverageBar.glsl.generated.ts'),
+}

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

@@ -0,0 +1,6 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageBar.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 expandToMinWidthLeftPx_0(float x1_0, float x2_0, float minWidth_0)\n{\n\n    float _S5;\n    if((x2_0 - x1_0) < minWidth_0)\n    {\n\n        _S5 = (x1_0 + x2_0) * 0.5 - minWidth_0 * 0.5;\n\n    }\n    else\n    {\n\n        _S5 = x1_0;\n\n    }\n\n    return _S5;\n}\n\nfloat expandToMinWidthRightPx_0(float x1_1, float x2_1, float minWidth_1)\n{\n\n    float _S6;\n    if((x2_1 - x1_1) < minWidth_1)\n    {\n\n        _S6 = (x1_1 + x2_1) * 0.5 + minWidth_1 * 0.5;\n\n    }\n    else\n    {\n\n        _S6 = x2_1;\n\n    }\n\n    return _S6;\n}\n\nvec2 expandToMinWidthPx_0(float x1_2, float x2_2, float minWidth_2)\n{\n\n    return vec2(expandToMinWidthLeftPx_0(x1_2, x2_2, minWidth_2), expandToMinWidthRightPx_0(x1_2, x2_2, minWidth_2));\n}\n\n\n\n\nvec2 expandToMinWidthX_0(float x1_3, float x2_3, float minWidthPx_0, float canvasWidth_0)\n{\n\n    return expandToMinWidthPx_0(x1_3, x2_3, pxToClipLen_0(minWidthPx_0, canvasWidth_0));\n}\n\n\nvec2 covExpandMinWidthX_0(float x1_4, float x2_4, float canvasWidth_1)\n{\n\n    return expandToMinWidthX_0(x1_4, x2_4, 1.0, canvasWidth_1);\n}\n\n\nfloat covBottomOffsetPx_0(float covHeight_1, float covYOffset_1)\n{\n\n    return covHeight_1 - covYOffset_1;\n}\n\n\nfloat covBaselinePx_0(CoverageBandUniforms_0 u_2)\n{\n\n    return u_2.covTop_0 + covBottomOffsetPx_0(u_2.covHeight_0, u_2.covYOffset_0);\n}\n\n\nfloat log1pf_0(float x_0)\n{\n\n    float u_3 = 1.0 + x_0;\n\n    float _S7;\n    if(u_3 == 1.0)\n    {\n\n        _S7 = x_0;\n\n    }\n    else\n    {\n\n        _S7 = log(u_3) * (x_0 / (u_3 - 1.0));\n\n    }\n\n    return _S7;\n}\n\nfloat symlogTransform_0(float x_1, float c_0)\n{\n\n    return float((int(sign((x_1))))) * log1pf_0(abs(x_1 / c_0));\n}\n\n\nfloat normalizeScoreUnclamped_0(float score_0, float domainMin_0, float domainMax_0, int scaleType_0, float symlogConstant_0)\n{\n\n    float tMax_0;\n\n    float tMin_0;\n\n    float t_0;\n\n\n\n    if(scaleType_0 == 2)\n    {\n\n        float _S8 = symlogTransform_0(score_0, symlogConstant_0);\n        float _S9 = symlogTransform_0(domainMin_0, symlogConstant_0);\n\n        tMax_0 = symlogTransform_0(domainMax_0, symlogConstant_0);\n\n        tMin_0 = _S9;\n\n        t_0 = _S8;\n\n    }\n    else\n    {\n\n        if(scaleType_0 == 1)\n        {\n            if(domainMin_0 > 0.0)\n            {\n\n                tMax_0 = domainMin_0;\n\n            }\n            else\n            {\n\n                tMax_0 = 1.0;\n\n            }\n            float _S10 = log2(max(score_0, tMax_0));\n            float _S11 = log2(tMax_0);\n\n            tMax_0 = log2(max(domainMax_0, tMax_0));\n\n            tMin_0 = _S11;\n\n            t_0 = _S10;\n\n        }\n        else\n        {\n\n            tMax_0 = domainMax_0;\n\n            tMin_0 = domainMin_0;\n\n            t_0 = score_0;\n\n        }\n\n    }\n\n    float tRange_0 = tMax_0 - tMin_0;\n    if(tRange_0 <= 0.0)\n    {\n\n        if(t_0 > tMin_0)\n        {\n\n            tMax_0 = 1.0;\n\n        }\n        else\n        {\n\n            tMax_0 = 0.0;\n\n        }\n\n        return tMax_0;\n    }\n    return (t_0 - tMin_0) / tRange_0;\n}\n\nfloat normalizeScore_0(float score_1, float domainMin_1, float domainMax_1, int scaleType_1, float symlogConstant_1)\n{\n\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1), 0.0, 1.0);\n}\n\n\nfloat normalizeDepthScalar_0(float rawDepth_0, float domainMin_2, float domainMax_2, int scaleType_2, float symlogConstant_2)\n{\n\n    return normalizeScore_0(rawDepth_0, domainMin_2, domainMax_2, scaleType_2, symlogConstant_2);\n}\n\n\nfloat covNormalizeDepth_0(float relDepth_0, CoverageBandUniforms_0 u_4)\n{\n    return normalizeDepthScalar_0(relDepth_0 * u_4.regionMaxDepth_0, u_4.depthDomainMin_0, u_4.depthDomainMax_0, u_4.coverageScaleType_0, u_4.coverageSymlogConstant_0);\n}\n\n\nfloat covEffectiveHeightPx_0(float covHeight_2, float covYOffset_2)\n{\n\n    return max(covHeight_2 - 2.0 * covYOffset_2, 0.0);\n}\n\n\nfloat covEffHeight_0(CoverageBandUniforms_0 u_5)\n{\n\n    return covEffectiveHeightPx_0(u_5.covHeight_0, u_5.covYOffset_0);\n}\n\n\nfloat covBarHeightPx_0(float depthFraction_0, float effectiveHeightPx_0)\n{\n\n    return depthFraction_0 * effectiveHeightPx_0;\n}\n\n\nfloat covSegBottomPx_0(float baselinePx_0, float yOffset_0, float barHeightPx_0)\n{\n    return baselinePx_0 - yOffset_0 * barHeightPx_0;\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 covSegTopPx_0(float baselinePx_1, float yOffset_1, float segHeight_0, float barHeightPx_1)\n{\n    return covSegBottomPx_0(baselinePx_1, yOffset_1, barHeightPx_1) - segHeight_0 * barHeightPx_1;\n}\n\n\nfloat covFlipX_0(float x_2, CoverageBandUniforms_0 u_6)\n{\n\n    return mix(x_2, - x_2, u_6.reversed_0);\n}\n\n\nvec4 covFlippedQuad_0(float x0_0, float x1_5, float yBot_0, float yTop_0, vec2 local_0, CoverageBandUniforms_0 u_7)\n{\n    return vec4(covFlipX_0(mix(x0_0, x1_5, local_0.x), u_7), mix(yBot_0, yTop_0, local_0.y), 0.0, 1.0);\n}\n\n\nvec4 covSegQuad_0(uint position_0, uint widthBp_0, float yOffset_2, float segHeight_1, float relDepth_1, vec2 local_1, CoverageBandUniforms_0 u_8)\n{\n\n    vec2 xx_0 = covExpandMinWidthX_0(covBpToClipX_0(position_0, u_8), covBpToClipX_0(position_0 + widthBp_0, u_8), u_8.canvasW_0);\n\n    float basePx_0 = covBaselinePx_0(u_8);\n    float barPx_0 = covBarHeightPx_0(covNormalizeDepth_0(relDepth_1, u_8), covEffHeight_0(u_8));\n    return covFlippedQuad_0(xx_0.x, xx_0.y, yPxToClipY_0(covSegBottomPx_0(basePx_0, yOffset_2, barPx_0), u_8.canvasH_0), yPxToClipY_0(covSegTopPx_0(basePx_0, yOffset_2, segHeight_1, barPx_0), u_8.canvasH_0), local_1, u_8);\n}\n\n\nvec4 unpackRGBA_0(uint c_1)\n{\n\n    return vec4(float((c_1 >> 0U) & 255U), float((c_1 >> 8U) & 255U), float((c_1 >> 16U) & 255U), float((c_1 >> 24U) & 255U)) / 255.0;\n}\n\n\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_position;\n\n\nlayout(location = 1)\nin float a_relDepth;\n\n\nstruct CoverageVsOut_0\n{\n    vec4 position_1;\n    vec4 color_0;\n};\n\n\nvoid main()\n{\n\n    CoverageVsOut_0 o_0;\n\n    CoverageBandUniforms_0 _S12 = 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    o_0.position_1 = covSegQuad_0(a_position, uint(u_0.binSize_0), 0.0, 1.0, a_relDepth, quadLocal_0(uint(gl_VertexID)), _S12);\n\n    o_0.color_0 = vec4(unpackRGBA_0(u_0.colorCoverage_0).xyz, 1.0);\n    CoverageVsOut_0 _S13 = o_0;\n\n    gl_Position = o_0.position_1;\n\n    v_color = _S13.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/coverageBar.iface.generated.ts

@@ -4,7 +4,7 @@
 import type { VertexAttributeLayout, ShaderBinding } from '@jbrowse/render-core/hal'
 
 export const BINDINGS: readonly ShaderBinding[] = [
-  { index: 1, kind: 'uniform', name: 'u' },
+  { index: 1, kind: 'uniform', name: 'u', stages: ['vertex'] },
 ]
 
 export const VERTS_PER_INSTANCE = 6

package/src/shaders/coverageBar.slang

@@ -4,6 +4,15 @@
 // coverageBand.slang.
 //! targets: wgsl, glsl
 //! layout-out: packages/alignments-core/src/coverageLayout.generated.ts
+// hpmath's midpoint min-width rule, lifted from the pass that draws through it
+// (covSegQuad -> covExpandMinWidthX). It is authored in a module whose own
+// js-export-out redirect is spoken for, and its twin is wanted in
+// @jbrowse/alignments-core, which cannot import render-core — the case adr-051
+// keeps "a shader may export a function it only imports" for. Consumer:
+// `fillSpanRect`, the one Canvas2D spelling the band and the pileup's span
+// marks share.
+//! js-export: expandToMinWidthLeftPx, expandToMinWidthRightPx
+//! js-export-out: packages/alignments-core/src/spanMinWidth.generated.ts
 
 import coverageBand;
 import colorPack;
@@ -12,17 +21,8 @@ import hpmath;
 public static const uint VERTS_PER_INSTANCE = 6u;
 [[vk::binding(1, 0)]] ConstantBuffer<CoverageBandUniforms> u;
 
-// `relDepth` is depths[bin] / regionMaxDepth, the same packing the SNP and
-// modification segments carry (see coverageSnp.slang) — covNormalizeDepth
-// recovers the raw depth from it. The bin spans `u.binSize` bp: 1 for per-bp
-// coverage, wider once the producer downsamples.
-struct Instance {
-  uint  position : ATTR0;
-  float relDepth : ATTR1;
-};
-
 [shader("vertex")]
-CoverageVsOut vs_main(Instance inst, uint vid : SV_VertexID) {
+CoverageVsOut vs_main(CoverageBarInstance inst, uint vid : SV_VertexID) {
   CoverageVsOut o;
   // The whole bar: y span 0..1 of itself. Same quad the segments stack inside.
   o.position = covSegQuad(inst.position, uint(u.binSize), 0.0, 1.0,

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

@@ -0,0 +1,4 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageBar.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 expandToMinWidthLeftPx_0( x1_0 : f32,  x2_0 : f32,  minWidth_0 : f32) -> f32\n{\n    var _S5 : f32;\n    if((x2_0 - x1_0) < minWidth_0)\n    {\n        _S5 = (x1_0 + x2_0) * 0.5f - minWidth_0 * 0.5f;\n    }\n    else\n    {\n        _S5 = x1_0;\n    }\n    return _S5;\n}\n\nfn expandToMinWidthRightPx_0( x1_1 : f32,  x2_1 : f32,  minWidth_1 : f32) -> f32\n{\n    var _S6 : f32;\n    if((x2_1 - x1_1) < minWidth_1)\n    {\n        _S6 = (x1_1 + x2_1) * 0.5f + minWidth_1 * 0.5f;\n    }\n    else\n    {\n        _S6 = x2_1;\n    }\n    return _S6;\n}\n\nfn expandToMinWidthPx_0( x1_2 : f32,  x2_2 : f32,  minWidth_2 : f32) -> vec2<f32>\n{\n    return vec2<f32>(expandToMinWidthLeftPx_0(x1_2, x2_2, minWidth_2), expandToMinWidthRightPx_0(x1_2, x2_2, minWidth_2));\n}\n\nfn expandToMinWidthX_0( x1_3 : f32,  x2_3 : f32,  minWidthPx_0 : f32,  canvasWidth_0 : f32) -> vec2<f32>\n{\n    return expandToMinWidthPx_0(x1_3, x2_3, pxToClipLen_0(minWidthPx_0, canvasWidth_0));\n}\n\nfn covExpandMinWidthX_0( x1_4 : f32,  x2_4 : f32,  canvasWidth_1 : f32) -> vec2<f32>\n{\n    return expandToMinWidthX_0(x1_4, x2_4, 1.0f, canvasWidth_1);\n}\n\nfn covBottomOffsetPx_0( covHeight_1 : f32,  covYOffset_1 : f32) -> f32\n{\n    return covHeight_1 - covYOffset_1;\n}\n\nfn covBaselinePx_0( u_2 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return (*u_2).covTop_0 + covBottomOffsetPx_0((*u_2).covHeight_0, (*u_2).covYOffset_0);\n}\n\nfn log1pf_0( x_0 : f32) -> f32\n{\n    var u_3 : f32 = 1.0f + x_0;\n    var _S7 : f32;\n    if(u_3 == 1.0f)\n    {\n        _S7 = x_0;\n    }\n    else\n    {\n        _S7 = log(u_3) * (x_0 / (u_3 - 1.0f));\n    }\n    return _S7;\n}\n\nfn symlogTransform_0( x_1 : f32,  c_0 : f32) -> f32\n{\n    return f32((i32(sign((x_1))))) * log1pf_0(abs(x_1 / c_0));\n}\n\nfn normalizeScoreUnclamped_0( score_0 : f32,  domainMin_0 : f32,  domainMax_0 : f32,  scaleType_0 : i32,  symlogConstant_0 : f32) -> f32\n{\n    var tMax_0 : f32;\n    var tMin_0 : f32;\n    var t_0 : f32;\n    if(scaleType_0 == i32(2))\n    {\n        var _S8 : f32 = symlogTransform_0(score_0, symlogConstant_0);\n        var _S9 : f32 = symlogTransform_0(domainMin_0, symlogConstant_0);\n        tMax_0 = symlogTransform_0(domainMax_0, symlogConstant_0);\n        tMin_0 = _S9;\n        t_0 = _S8;\n    }\n    else\n    {\n        if(scaleType_0 == i32(1))\n        {\n            if(domainMin_0 > 0.0f)\n            {\n                tMax_0 = domainMin_0;\n            }\n            else\n            {\n                tMax_0 = 1.0f;\n            }\n            var _S10 : f32 = log2(max(score_0, tMax_0));\n            var _S11 : f32 = log2(tMax_0);\n            tMax_0 = log2(max(domainMax_0, tMax_0));\n            tMin_0 = _S11;\n            t_0 = _S10;\n        }\n        else\n        {\n            tMax_0 = domainMax_0;\n            tMin_0 = domainMin_0;\n            t_0 = score_0;\n        }\n    }\n    var tRange_0 : f32 = tMax_0 - tMin_0;\n    if(tRange_0 <= 0.0f)\n    {\n        if(t_0 > 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_0 - tMin_0) / tRange_0;\n}\n\nfn normalizeScore_0( score_1 : f32,  domainMin_1 : f32,  domainMax_1 : f32,  scaleType_1 : i32,  symlogConstant_1 : f32) -> f32\n{\n    return clamp(normalizeScoreUnclamped_0(score_1, domainMin_1, domainMax_1, scaleType_1, symlogConstant_1), 0.0f, 1.0f);\n}\n\nfn normalizeDepthScalar_0( rawDepth_0 : f32,  domainMin_2 : f32,  domainMax_2 : f32,  scaleType_2 : i32,  symlogConstant_2 : f32) -> f32\n{\n    return normalizeScore_0(rawDepth_0, domainMin_2, domainMax_2, scaleType_2, symlogConstant_2);\n}\n\nfn covNormalizeDepth_0( relDepth_0 : f32,  u_4 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return normalizeDepthScalar_0(relDepth_0 * (*u_4).regionMaxDepth_0, (*u_4).depthDomainMin_0, (*u_4).depthDomainMax_0, (*u_4).coverageScaleType_0, (*u_4).coverageSymlogConstant_0);\n}\n\nfn covEffectiveHeightPx_0( covHeight_2 : f32,  covYOffset_2 : f32) -> f32\n{\n    return max(covHeight_2 - 2.0f * covYOffset_2, 0.0f);\n}\n\nfn covEffHeight_0( u_5 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return covEffectiveHeightPx_0((*u_5).covHeight_0, (*u_5).covYOffset_0);\n}\n\nfn covBarHeightPx_0( depthFraction_0 : f32,  effectiveHeightPx_0 : f32) -> f32\n{\n    return depthFraction_0 * effectiveHeightPx_0;\n}\n\nfn covSegBottomPx_0( baselinePx_0 : f32,  yOffset_0 : f32,  barHeightPx_0 : f32) -> f32\n{\n    return baselinePx_0 - yOffset_0 * barHeightPx_0;\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 covSegTopPx_0( baselinePx_1 : f32,  yOffset_1 : f32,  segHeight_0 : f32,  barHeightPx_1 : f32) -> f32\n{\n    return covSegBottomPx_0(baselinePx_1, yOffset_1, barHeightPx_1) - segHeight_0 * barHeightPx_1;\n}\n\nfn covFlipX_0( x_2 : f32,  u_6 : ptr<function, CoverageBandUniforms_std140_0>) -> f32\n{\n    return mix(x_2, - x_2, (*u_6).reversed_0);\n}\n\nfn covFlippedQuad_0( x0_0 : f32,  x1_5 : f32,  yBot_0 : f32,  yTop_0 : f32,  local_0 : vec2<f32>,  u_7 : ptr<function, CoverageBandUniforms_std140_0>) -> vec4<f32>\n{\n    var _S12 : f32 = covFlipX_0(mix(x0_0, x1_5, local_0.x), &((*u_7)));\n    return vec4<f32>(_S12, mix(yBot_0, yTop_0, local_0.y), 0.0f, 1.0f);\n}\n\nfn covSegQuad_0( position_0 : u32,  widthBp_0 : u32,  yOffset_2 : f32,  segHeight_1 : f32,  relDepth_1 : f32,  local_1 : vec2<f32>,  u_8 : ptr<function, CoverageBandUniforms_std140_0>) -> vec4<f32>\n{\n    var _S13 : f32 = covBpToClipX_0(position_0, &((*u_8)));\n    var _S14 : f32 = covBpToClipX_0(position_0 + widthBp_0, &((*u_8)));\n    var xx_0 : vec2<f32> = covExpandMinWidthX_0(_S13, _S14, (*u_8).canvasW_0);\n    var _S15 : f32 = covBaselinePx_0(&((*u_8)));\n    var _S16 : f32 = covNormalizeDepth_0(relDepth_1, &((*u_8)));\n    var _S17 : f32 = covEffHeight_0(&((*u_8)));\n    var barPx_0 : f32 = covBarHeightPx_0(_S16, _S17);\n    var _S18 : vec4<f32> = covFlippedQuad_0(xx_0.x, xx_0.y, yPxToClipY_0(covSegBottomPx_0(_S15, yOffset_2, barPx_0), (*u_8).canvasH_0), yPxToClipY_0(covSegTopPx_0(_S15, yOffset_2, segHeight_1, barPx_0), (*u_8).canvasH_0), local_1, &((*u_8)));\n    return _S18;\n}\n\nfn unpackRGBA_0( c_1 : u32) -> vec4<f32>\n{\n    return vec4<f32>(f32((((c_1 >> (u32(0)))) & (u32(255)))), f32((((c_1 >> (u32(8)))) & (u32(255)))), f32((((c_1 >> (u32(16)))) & (u32(255)))), f32((((c_1 >> (u32(24)))) & (u32(255))))) / vec4<f32>(255.0f);\n}\n\nstruct CoverageVsOut_0\n{\n    @builtin(position) position_1 : vec4<f32>,\n    @location(0) color_0 : vec4<f32>,\n};\n\nstruct vertexInput_0\n{\n    @location(0) position_2 : u32,\n    @location(1) relDepth_2 : f32,\n};\n\n@vertex\nfn vs_main( _S19 : vertexInput_0, @builtin(vertex_index) vid_1 : u32) -> CoverageVsOut_0\n{\n    var o_0 : CoverageVsOut_0;\n    var _S20 : u32 = u32(u_0.binSize_0);\n    var _S21 : vec2<f32> = quadLocal_0(vid_1);\n    var _S22 : CoverageBandUniforms_std140_0 = u_0;\n    var _S23 : vec4<f32> = covSegQuad_0(_S19.position_2, _S20, 0.0f, 1.0f, _S19.relDepth_2, _S21, &(_S22));\n    o_0.position_1 = _S23;\n    o_0.color_0 = vec4<f32>(unpackRGBA_0(u_0.colorCoverage_0).xyz, 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( _S24 : pixelInput_0, @builtin(position) position_3 : vec4<f32>) -> pixelOutput_0\n{\n    var _S25 : pixelOutput_0 = pixelOutput_0( _S24.color_1 );\n    return _S25;\n}\n\n"

package/src/shaders/coverageClip.generated.ts

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

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

@@ -0,0 +1,6 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageClip.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\nint clipSide_0(float value_0, float domainMin_0, float domainMax_0, int scaleType_0)\n{\n\n    float tolerance_0 = (domainMax_0 - domainMin_0) * 0.00000999999974738;\n    if(value_0 > (domainMax_0 + tolerance_0))\n    {\n\n        return 1;\n    }\n\n    bool _S1;\n\n    if(value_0 < (domainMin_0 - tolerance_0))\n    {\n\n        _S1 = scaleType_0 != 1;\n\n    }\n    else\n    {\n\n        _S1 = false;\n\n    }\n\n    if(_S1)\n    {\n\n        return -1;\n    }\n    return 0;\n}\n\n\nstruct CoverageVsOut_0\n{\n    vec4 position_0;\n    vec4 color_0;\n};\n\n\nCoverageVsOut_0 covDiscardVertex_0()\n{\n\n    CoverageVsOut_0 o_0;\n    o_0.position_0 = vec4(2.0, 2.0, 0.0, 1.0);\n    o_0.color_0 = vec4(0.0);\n    return o_0;\n}\n\n\nvec2 hpSplitUint_0(uint value_1)\n{\n\n    uint lo_0 = value_1 & 4095U;\n\n    return vec2(float(value_1 - lo_0), float(lo_0));\n}\n\n\nfloat hpToClipX_0(vec2 splitPos_0, vec3 bpRange_0, float hpZero_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 expandToMinWidthLeftPx_0(float x1_0, float x2_0, float minWidth_0)\n{\n\n    float _S2;\n    if((x2_0 - x1_0) < minWidth_0)\n    {\n\n        _S2 = (x1_0 + x2_0) * 0.5 - minWidth_0 * 0.5;\n\n    }\n    else\n    {\n\n        _S2 = x1_0;\n\n    }\n\n    return _S2;\n}\n\nfloat expandToMinWidthRightPx_0(float x1_1, float x2_1, float minWidth_1)\n{\n\n    float _S3;\n    if((x2_1 - x1_1) < minWidth_1)\n    {\n\n        _S3 = (x1_1 + x2_1) * 0.5 + minWidth_1 * 0.5;\n\n    }\n    else\n    {\n\n        _S3 = x2_1;\n\n    }\n\n    return _S3;\n}\n\nvec2 expandToMinWidthPx_0(float x1_2, float x2_2, float minWidth_2)\n{\n\n    return vec2(expandToMinWidthLeftPx_0(x1_2, x2_2, minWidth_2), expandToMinWidthRightPx_0(x1_2, x2_2, minWidth_2));\n}\n\n\n\n\nvec2 expandToMinWidthX_0(float x1_3, float x2_3, float minWidthPx_0, float canvasWidth_0)\n{\n\n    return expandToMinWidthPx_0(x1_3, x2_3, pxToClipLen_0(minWidthPx_0, canvasWidth_0));\n}\n\n\nvec2 covExpandMinWidthX_0(float x1_4, float x2_4, float canvasWidth_1)\n{\n\n    return expandToMinWidthX_0(x1_4, x2_4, 1.0, canvasWidth_1);\n}\n\n\nfloat covBottomOffsetPx_0(float covHeight_1, float covYOffset_1)\n{\n\n    return covHeight_1 - covYOffset_1;\n}\n\n\nfloat covBaselinePx_0(CoverageBandUniforms_0 u_2)\n{\n\n    return u_2.covTop_0 + covBottomOffsetPx_0(u_2.covHeight_0, u_2.covYOffset_0);\n}\n\n\nfloat covEffectiveHeightPx_0(float covHeight_2, float covYOffset_2)\n{\n\n    return max(covHeight_2 - 2.0 * covYOffset_2, 0.0);\n}\n\n\nfloat covEffHeight_0(CoverageBandUniforms_0 u_3)\n{\n\n    return covEffectiveHeightPx_0(u_3.covHeight_0, u_3.covYOffset_0);\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\nvec2 quadLocal_0(uint vid_0)\n{\n\n    uint v_0 = vid_0 % 6U;\n\n    bool _S4 = v_0 == 0U;\n\n    bool _S5;\n\n    if(_S4)\n    {\n\n        _S5 = true;\n\n    }\n    else\n    {\n\n        _S5 = v_0 == 2U;\n\n    }\n\n    if(_S5)\n    {\n\n        _S5 = true;\n\n    }\n    else\n    {\n\n        _S5 = v_0 == 3U;\n\n    }\n\n    float _S6;\n\n    if(_S5)\n    {\n\n        _S6 = 0.0;\n\n    }\n    else\n    {\n\n        _S6 = 1.0;\n\n    }\n    if(_S4)\n    {\n\n        _S5 = true;\n\n    }\n    else\n    {\n\n        _S5 = v_0 == 1U;\n\n    }\n\n    if(_S5)\n    {\n\n        _S5 = true;\n\n    }\n    else\n    {\n\n        _S5 = v_0 == 4U;\n\n    }\n\n    float _S7;\n\n    if(_S5)\n    {\n\n        _S7 = 0.0;\n\n    }\n    else\n    {\n\n        _S7 = 1.0;\n\n    }\n\n    return vec2(_S6, _S7);\n}\n\n\nfloat covFlipX_0(float x_0, CoverageBandUniforms_0 u_4)\n{\n\n    return mix(x_0, - x_0, u_4.reversed_0);\n}\n\n\nvec4 covFlippedQuad_0(float x0_0, float x1_5, float yBot_0, float yTop_0, vec2 local_0, CoverageBandUniforms_0 u_5)\n{\n    return vec4(covFlipX_0(mix(x0_0, x1_5, local_0.x), u_5), 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\nout vec4 v_color;\n\n\nlayout(location = 0)\nin uint a_position;\n\n\nlayout(location = 1)\nin float a_relDepth;\n\n\nvoid main()\n{\n    if((clipSide_0(a_relDepth * u_0.regionMaxDepth_0, u_0.depthDomainMin_0, u_0.depthDomainMax_0, u_0.coverageScaleType_0)) != 1)\n    {\n\n        CoverageVsOut_0 _S8 = covDiscardVertex_0();\n\n        gl_Position = _S8.position_0;\n\n        v_color = _S8.color_0;\n\n        return;\n    }\n\n    CoverageBandUniforms_0 _S9 = 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    vec2 xx_0 = covExpandMinWidthX_0(covBpToClipX_0(a_position, _S9), covBpToClipX_0(a_position + uint(u_0.binSize_0), _S9), u_0.canvasW_0);\n\n    float topPx_0 = covBaselinePx_0(_S9) - covEffHeight_0(_S9);\n    CoverageVsOut_0 o_1;\n    o_1.position_0 = covFlippedQuad_0(xx_0.x, xx_0.y, yPxToClipY_0(topPx_0 + 2.0, u_0.canvasH_0), yPxToClipY_0(topPx_0, u_0.canvasH_0), quadLocal_0(uint(gl_VertexID)), _S9);\n\n    o_1.color_0 = vec4(unpackRGBA_0(4280821958U).xyz, 1.0);\n    CoverageVsOut_0 _S10 = o_1;\n\n    gl_Position = o_1.position_0;\n\n    v_color = _S10.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/coverageClip.iface.generated.ts

@@ -0,0 +1,172 @@
+// AUTO-GENERATED by packages/shader-tools/src/shader-codegen from coverageClip.slang.
+// Do not edit. Run `pnpm gen:shaders` to regenerate.
+
+import type { VertexAttributeLayout, ShaderBinding } from '@jbrowse/render-core/hal'
+
+export const BINDINGS: readonly ShaderBinding[] = [
+  { index: 1, kind: 'uniform', name: 'u', stages: ['vertex'] },
+]
+
+export const VERTS_PER_INSTANCE = 6
+
+export const UNIFORMS_SIZE_BYTES = 112
+
+// Word indices into a Float32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_F32 = {
+  bpHi: 0,
+  bpLo: 1,
+  bpLen: 2,
+  hpZero: 3,
+  canvasW: 4,
+  canvasH: 5,
+  covHeight: 6,
+  covYOffset: 7,
+  covTop: 8,
+  regionMaxDepth: 9,
+  depthDomainMax: 10,
+  depthDomainMin: 11,
+  coverageSymlogConstant: 12,
+  binSize: 13,
+  interbaseHeight: 14,
+  snpMinFreq: 15,
+  reversed: 16,
+} as const
+
+// Word indices into a Uint32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_U32 = {
+  colorCoverage: 18,
+  colorBaseA: 19,
+  colorBaseC: 20,
+  colorBaseG: 21,
+  colorBaseT: 22,
+  colorBaseN: 23,
+  colorInsertionIndicator: 24,
+  colorSoftclipIndicator: 25,
+  colorHardclipIndicator: 26,
+} as const
+
+// Word indices into a Int32Array view over the uniform buffer.
+export const UNIFORM_OFFSET_I32 = {
+  coverageScaleType: 17,
+} as const
+
+
+export interface Uniforms {
+  bpHi: number
+  bpLo: number
+  bpLen: number
+  hpZero: number
+  canvasW: number
+  canvasH: number
+  covHeight: number
+  covYOffset: number
+  covTop: number
+  regionMaxDepth: number
+  depthDomainMax: number
+  depthDomainMin: number
+  coverageSymlogConstant: number
+  binSize: number
+  interbaseHeight: number
+  snpMinFreq: number
+  reversed: number
+  coverageScaleType: number
+  colorCoverage: number
+  colorBaseA: number
+  colorBaseC: number
+  colorBaseG: number
+  colorBaseT: number
+  colorBaseN: number
+  colorInsertionIndicator: number
+  colorSoftclipIndicator: number
+  colorHardclipIndicator: number
+}
+
+export function writeUniforms(buf: ArrayBuffer, uniforms: Uniforms) {
+  const f32 = new Float32Array(buf)
+  const u32 = new Uint32Array(buf)
+  const i32 = new Int32Array(buf)
+  f32[0] = uniforms.bpHi
+  f32[1] = uniforms.bpLo
+  f32[2] = uniforms.bpLen
+  f32[3] = uniforms.hpZero
+  f32[4] = uniforms.canvasW
+  f32[5] = uniforms.canvasH
+  f32[6] = uniforms.covHeight
+  f32[7] = uniforms.covYOffset
+  f32[8] = uniforms.covTop
+  f32[9] = uniforms.regionMaxDepth
+  f32[10] = uniforms.depthDomainMax
+  f32[11] = uniforms.depthDomainMin
+  f32[12] = uniforms.coverageSymlogConstant
+  f32[13] = uniforms.binSize
+  f32[14] = uniforms.interbaseHeight
+  f32[15] = uniforms.snpMinFreq
+  f32[16] = uniforms.reversed
+  i32[17] = uniforms.coverageScaleType
+  u32[18] = uniforms.colorCoverage
+  u32[19] = uniforms.colorBaseA
+  u32[20] = uniforms.colorBaseC
+  u32[21] = uniforms.colorBaseG
+  u32[22] = uniforms.colorBaseT
+  u32[23] = uniforms.colorBaseN
+  u32[24] = uniforms.colorInsertionIndicator
+  u32[25] = uniforms.colorSoftclipIndicator
+  u32[26] = uniforms.colorHardclipIndicator
+}
+
+export const INSTANCE_STRIDE_BYTES = 8
+export const INSTANCE_STRIDE_WORDS = 2
+
+// Word indices into a Float32Array view over the instance buffer.
+export const INSTANCE_OFFSET_F32 = {
+  relDepth: 1,
+} as const
+
+// Word indices into a Uint32Array view over the instance buffer.
+export const INSTANCE_OFFSET_U32 = {
+  position: 0,
+} as const
+
+export const VERTEX_ATTRIBUTES: readonly VertexAttributeLayout[] = [
+  { name: 'a_position', components: 1, type: 'uint', offsetBytes: 0, integer: true },
+  { name: 'a_relDepth', components: 1, type: 'float', offsetBytes: 4, integer: false },
+]
+
+export interface InstanceArrays {
+  position: ArrayLike<number>
+  relDepth: 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 { position, relDepth } = arrays
+  for (let i = 0; i < numInstances; i++) {
+    const o = i * INSTANCE_STRIDE_WORDS
+    u32[o + 0] = position[i]!
+    f32[o + 1] = relDepth[i]!
+  }
+  return buf
+}
+
+// Instance `i`'s `position`.
+export function getInstancePosition(u32: Uint32Array, i: number) {
+  return u32[i * INSTANCE_STRIDE_WORDS]!
+}
+
+export function setInstancePosition(u32: Uint32Array, i: number, v: number) {
+  u32[i * INSTANCE_STRIDE_WORDS] = v
+}
+
+// Instance `i`'s `relDepth`.
+export function getInstanceRelDepth(f32: Float32Array, i: number) {
+  return f32[i * INSTANCE_STRIDE_WORDS + 1]!
+}
+
+export function setInstanceRelDepth(f32: Float32Array, i: number, v: number) {
+  f32[i * INSTANCE_STRIDE_WORDS + 1] = v
+}

package/src/shaders/coverageClip.slang

@@ -0,0 +1,38 @@
+// The coverage band's clip strips: `CLIP_STRIP_PX` of red across the top of
+// every depth bar whose depth the domain cut (ADR-183). Its own pass off the
+// depth bars' buffer, drawn last, since the SNP, modification and interbase
+// layers stack inside the bar and would paint over a strip drawn with it.
+// `drawCoverageClipStrips` in @jbrowse/alignments-core is the Canvas2D twin.
+//! targets: wgsl, glsl
+
+import coverageBand;
+import colorPack;
+import clipStrip;
+import hpmath;
+
+public static const uint VERTS_PER_INSTANCE = 6u;
+[[vk::binding(1, 0)]] ConstantBuffer<CoverageBandUniforms> u;
+
+[shader("vertex")]
+CoverageVsOut vs_main(CoverageBarInstance inst, uint vid : SV_VertexID) {
+  float rawDepth = inst.relDepth * u.regionMaxDepth;
+  if (clipSide(rawDepth, u.depthDomainMin, u.depthDomainMax, u.coverageScaleType) != 1) {
+    return covDiscardVertex();
+  }
+  float2 xx = covExpandMinWidthX(
+    covBpToClipX(inst.position, u), covBpToClipX(inst.position + uint(u.binSize), u), u.canvasW);
+  float topPx = covBaselinePx(u) - covEffHeight(u);
+  CoverageVsOut o;
+  o.position = covFlippedQuad(
+    xx.x, xx.y,
+    yPxToClipY(topPx + CLIP_STRIP_PX, u.canvasH),
+    yPxToClipY(topPx, u.canvasH),
+    quadLocal(vid), u);
+  o.color = float4(unpackRGBA(CLIP_STRIP_COLOR).xyz, 1.0);
+  return o;
+}
+
+[shader("fragment")]
+float4 fs_main(CoverageVsOut fragIn) : SV_Target {
+  return fragIn.color;
+}