edf2csv 0.5.82 → 0.5.84

package/CHANGELOG.md

@@ -3,6 +3,70 @@
 Notable changes to edf2csv. Versions follow [semantic versioning](https://semver.org); while the
 major version is 0, a minor bump may contain breaking changes.
 
+## 0.5.84
+
+### Fixed: a record duration too small to divide into dropped every sample and blamed the window
+
+A sampling rate is samples per record over record duration. EDF's record-duration field is
+eight characters and accepts `1e-308`, so four samples in one of those records is `Infinity` —
+and `1 / Infinity` is zero, which the resolution check reads as "no step to report" rather than
+"no resolution at all".
+
+So a file holding two complete records of four samples wrote none of them, exited 0, and
+printed one warning:
+
+```
+warning: This recording's 2 data records carry no samples in range, so the signal files hold
+         their headers and no data.
+```
+
+Untrue twice: the records carry their eight samples, and no range was asked for. `--info`
+listed the channel at `Infinity Hz` and predicted zero rows, which at least agreed with the
+conversion.
+
+One power of ten away, at 1e-300, the rate is 4e300 and the same file converts all eight rows
+with `TIME_RESOLUTION`. That code covers this now too, in a branch of its own — the existing
+hint promises "Every sample is written, in order", which would have been the third false
+sentence — and `EMPTY_WINDOW` no longer fires over the top of it, since the rate warning is the
+accurate account of the same zero.
+
+The same guard `decimalsAreClamped` had before 0.5.83, one column over: a step of exactly zero
+means the quantity could not be computed, not that there is nothing to say about it.
+
+## 0.5.83
+
+### Fixed: a physical span too small to represent became a flat channel, silently
+
+```
+time_s,MAG
+0.000,0.000
+0.250,0.000
+0.500,0.000
+```
+
+Eight samples spanning digital -16,000 to +12,000, all written as the same number, no
+diagnostic anywhere, `--strict` exiting 0. The header declares -1e-320 to 1e-320 over the full
+16-bit range: 65,536 distinct physical values, none of them equal to another.
+
+The gain is the span over the digital range — 2e-320/65535, or 3e-325, which is smaller than
+the smallest subnormal double and underflows to +0. `makeScaler` tests `gain === 0` and takes
+its flat-range branch, whose comment is correct about the case it was written for: "A flat
+physical range makes every sample the same value ... That mapping is defined, so its constant
+is written." An underflowed gain is not a flat range, and from inside that test the two look
+identical.
+
+The answer was already in the function, eight lines below. Overflow gets it: "the physical span
+overflowed a double, so there is no mapping at all. Returning physicalMin filled the column with
+one enormous constant — every distinct sample rendered as the same 300-digit number — and raised
+nothing." Underflow is the same fact about the same header and now takes the same route: empty
+cells, and `UNUSABLE_PHYSICAL_RANGE` saying the span is too small rather than too large.
+
+A genuinely flat range still writes its constant. That mapping is defined, every sample really
+is that value, and it has `DEGENERATE_PHYSICAL_RANGE` of its own.
+
+One power of ten away, at 1e-319, the same file has always raised `VALUE_RESOLUTION` — this was
+the one gap in a row of neighbours that all report themselves.
+
 ## 0.5.82
 
 ### Fixed: `--stdout --gzip` onto a full destination announced every row and exited 0

package/dist/convert/plan.js

@@ -114,7 +114,34 @@ export function buildPlan(input, options = {}) {
     */
     for (const group of groups) {
         const step = group.rate > 0 ? 1 / group.rate : 0;
-        if (step > 0 && step < 10 ** -group.timeDecimals) {
+        /*
+          The limit of the same failure, which read as the absence of it.
+    
+          `samplesPerRecord / recordDuration` is a double, and a record duration of 1e-308 with
+          four samples in it is Infinity. `1 / Infinity` is 0, so `step > 0` was false and this
+          said nothing — while every sample was dropped, the run exited 0, and the only warning
+          printed was EMPTY_WINDOW's "This recording's 2 data records carry no samples in range",
+          which is untrue twice over: the records carry eight samples and no range was asked for.
+    
+          One power of ten away, at 1e-300, the rate is 4e300 and the file converts with the
+          warning below. Same guard `decimalsAreClamped` had before 0.5.83, in the column next
+          door: a step of exactly zero means no resolution at all, not nothing to report.
+    
+          Its own branch because the hint below is false here — no rows are written at all, so
+          "Every sample is written, in order" would be the third untrue sentence.
+        */
+        if (!Number.isFinite(group.rate)) {
+            diagnostics.push({
+                code: 'TIME_RESOLUTION',
+                severity: 'warning',
+                message: `Channels in ${group.fileName} work out to a sampling rate of ${formatRate(group.rate)} Hz ` +
+                    `— their samples per record over a record duration too small to divide into — so ` +
+                    `their samples cannot be placed in time and no rows are written for them.`,
+                hint: 'Check the record duration in the header. One power of ten larger and the same ' +
+                    'file converts, with consecutive rows carrying the same time_s.',
+            });
+        }
+        else if (step > 0 && step < 10 ** -group.timeDecimals) {
             diagnostics.push({
                 code: 'TIME_RESOLUTION',
                 severity: 'warning',
@@ -173,7 +200,14 @@ export function buildPlan(input, options = {}) {
       crossed with every option set, so raising it from the plan says the same thing the rows
       would have.
     */
-    if (writeSignals && groups.length > 0 && estimate.rows === 0) {
+    /*
+      Not when a rate above already explained it. EMPTY_WINDOW says the records "carry no
+      samples in range", and on a recording whose rate overflowed to Infinity that is untrue
+      twice: the records carry their samples, and no range was asked for. The rate warning is
+      the accurate account of the same zero.
+    */
+    const untimeable = groups.some((group) => !Number.isFinite(group.rate));
+    if (writeSignals && groups.length > 0 && estimate.rows === 0 && !untimeable) {
         diagnostics.push(emptyWindow(range, input.recordCount));
     }
     if (estimate.exceedsSpreadsheetLimit) {

package/dist/convert/plan.js.map

@@ -1 +1 @@
-{"version":3,"file":"plan.js","sourceRoot":"","sources":["../../src/convert/plan.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAIH,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AACxE,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AACpE,OAAO,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,YAAY,EAAE,MAAM,qBAAqB,CAAC;AAC1D,OAAO,EAAE,gBAAgB,EAAE,kBAAkB,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AACrF,OAAO,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAC7C,OAAO,EAAE,mBAAmB,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAiGpE,MAAM,SAAS,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC;AAE9C,oFAAoF;AACpF,MAAM,CAAC,MAAM,qBAAqB,GAAG,SAAS,CAAC;AAE/C,MAAM,UAAU,SAAS,CAAC,KAAgB,EAAE,UAAuB,EAAE;IACnE,oFAAoF;IACpF,yEAAyE;IACzE,aAAa,CAAC,OAAO,CAAC,CAAC;IAEvB,MAAM,WAAW,GAAiB,EAAE,CAAC;IACrC,MAAM,WAAW,GAAG,gBAAgB,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IAEpD,sFAAsF;IACtF,uFAAuF;IACvF,wFAAwF;IACxF,KAAK,MAAM,MAAM,IAAI,kBAAkB,CAAC,KAAK,CAAC,OAAO,EAAE,WAAW,CAAC,EAAE,CAAC;QACpE,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,iBAAiB;YACvB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,UAAU,MAAM,CAAC,KAAK,iBAAiB,MAAM,CAAC,KAAK,mCAAmC;gBACtF,8DAA8D,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI;YACjG,IAAI,EAAE,oFAAoF;SAC3F,CAAC,CAAC;IACL,CAAC;IAED,MAAM,KAAK,GAAG,YAAY,CAAC;QACzB,KAAK,EAAE,OAAO,CAAC,KAAK;QACpB,SAAS,EAAE,OAAO,CAAC,SAAS;QAC5B,QAAQ,EAAE,OAAO,CAAC,QAAQ;QAC1B,GAAG,EAAE,OAAO,CAAC,GAAG;QAChB,OAAO,EAAE,OAAO,CAAC,OAAO;QACxB,cAAc,EAAE,KAAK,CAAC,cAAc;QACpC,WAAW,EAAE,KAAK,CAAC,WAAW;QAC9B,YAAY,EAAE,KAAK,CAAC,YAAY;KACjC,CAAC,CAAC;IAEH,MAAM,YAAY,GAAG,OAAO,CAAC,eAAe,KAAK,IAAI,CAAC;IAEtD,IAAI,MAAM,GAAgB,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC;IAExE;;;;;;;;;MASE;IACF,IAAI,OAAO,CAAC,QAAQ,IAAI,OAAO,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpD,MAAM,SAAS,GAAG,cAAc,CAAC,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;QAClE,IAAI,YAAY;YAAE,MAAM,GAAG,SAAS,CAAC,OAAO,CAAC;QAC7C,KAAK,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,SAAS,CAAC,SAAS,EAAE,CAAC;YACpD,WAAW,CAAC,IAAI,CAAC;gBACf,IAAI,EAAE,iBAAiB;gBACvB,QAAQ,EAAE,SAAS;gBACnB,OAAO,EACL,IAAI,IAAI,aAAa,OAAO,CAAC,MAAM,uBAAuB;oBAC1D,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,+BAA+B;gBAC7E,IAAI,EAAE,oBAAoB,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,IAAI,CAAC,qBAAqB;aACtE,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IAED,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,MAAM,CAAC;IACxC,MAAM,MAAM,GAAG,YAAY;QACzB,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,WAAW,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,IAAI,KAAK,IAAI,EAAE,MAAM,CAAC;QACnF,CAAC,CAAC,EAAE,CAAC;IACP,MAAM,QAAQ,GAAG,cAAc,CAC7B,MAAM,EACN,KAAK,EACL,KAAK,CAAC,cAAc,EACpB,KAAK,CAAC,YAAY,EAClB,OAAO,CAAC,GAAG,KAAK,IAAI,EACpB,MAAM,CACP,CAAC;IAEF;;;;;;;;;;;;MAYE;IACF,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACtB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gBAAgB,MAAM,CAAC,MAAM,4BAA4B;gBACzD,IAAI,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI;YAC9E,IAAI,EACF,MAAM,KAAK,MAAM;gBACf,CAAC,CAAC,mFAAmF;gBACrF,CAAC,CAAC,mEAAmE;SAC1E,CAAC,CAAC;IACL,CAAC;IAED;;;;;;;;;MASE;IACF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACjD,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC;YACjD,WAAW,CAAC,IAAI,CAAC;gBACf,IAAI,EAAE,iBAAiB;gBACvB,QAAQ,EAAE,SAAS;gBACnB,OAAO,EACL,eAAe,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,6CAA6C;oBAClF,uCAAuC,KAAK,CAAC,QAAQ,+BAA+B;gBACtF,IAAI,EACF,mFAAmF;oBACnF,4DAA4D;aAC/D,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IAED;;;;;;;;;;;;;;;;;;;;MAoBE;IACF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,kBAAkB,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;QACzE,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;YAAE,SAAS;QACjC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,kBAAkB;YACxB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,GAAG,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,WAAW;gBACzF,mFAAmF;gBACnF,qBAAqB,KAAK,CAAC,QAAQ,GAAG;YACxC,IAAI,EACF,kFAAkF;gBAClF,iEAAiE;SACpE,CAAC,CAAC;IACL,CAAC;IAED;;;;;;;;;;;;MAYE;IACF,IAAI,YAAY,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,IAAI,QAAQ,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;QAC7D,WAAW,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC;IAC1D,CAAC;IAED,IAAI,QAAQ,CAAC,uBAAuB,EAAE,CAAC;QACrC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,cAAc;YACpB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gDAAgD,qBAAqB,CAAC,cAAc,CAAC,OAAO,CAAC,GAAG;gBAChG,qDAAqD;YACvD,IAAI,EAAE,qFAAqF;SAC5F,CAAC,CAAC;IACL,CAAC;IAED,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,OAAO,CAAC,IAAI,KAAK,IAAI,EAAE,KAAK,EAAE,WAAW,EAAE,YAAY,EAAE,WAAW,EAAE,QAAQ,EAAE,CAAC;AAClH,CAAC;AAED;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,OAA6B,EAC7B,WAAgC,EAChC,cAAkC,EAClC,IAAa,EACb,MAAuB;IAEvB,MAAM,MAAM,GAAG,IAAI,GAAG,EAAuB,CAAC;IAC9C,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE,CAAC;QAC7B,wEAAwE;QACxE,2EAA2E;QAC3E,IAAI,MAAM,CAAC,gBAAgB,KAAK,CAAC;YAAE,SAAS;QAC5C,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;QAC/C,IAAI,MAAM;YAAE,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;;YAC3B,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC;IACjD,CAAC;IAED,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACvD,oFAAoF;IACpF,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,MAAM,KAAK,MAAM,CAAC;IAEvD;;;;;;;;;;;;;;;;;MAiBE;IACF,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC;IACzC,MAAM,KAAK,GAAG,WAAW,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;IAC9E,MAAM,IAAI,GAAG,IAAI,GAAG,EAAU,CAAC;IAC/B,MAAM,UAAU,GAAG,CAAC,KAAa,EAAU,EAAE;QAC3C,MAAM,IAAI,GAAG,WAAW,KAAK,CAAC,KAAK,CAAC,EAAE,CAAC;QACvC,IAAI,IAAI,GAAG,GAAG,IAAI,GAAG,MAAM,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,GAAG,GAAG,IAAI,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;QACpE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;IAEF,OAAO,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE;QAC/B,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACzB,OAAO;YACL,IAAI;YACJ,gBAAgB,EAAE,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;YACpD,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC,UAAU,MAAM,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC;YACzD;;;;cAIE;YACF,YAAY,EAAE,MAAM,KAAK,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,IAAI,CAAC;YAC3F,QAAQ,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;gBACjC,MAAM;gBACN,MAAM,EAAE,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,UAAU,MAAM,CAAC,KAAK,EAAE;gBACjE,QAAQ,EAAE,cAAc,IAAI,iBAAiB,CAAC,MAAM,CAAC;aACtD,CAAC,CAAC;SACJ,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,sBAAsB,CAAC,WAAkC;IACvE,OAAO,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,sBAAsB,CAAC,CAAC;AACtE,CAAC;AAED,gEAAgE;AAChE,MAAM,UAAU,QAAQ,CAAC,IAAY;IACnC,OAAO,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC;AACnD,CAAC;AAED,mFAAmF;AACnF,SAAS,OAAO,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAAM,GAAG,KAAK;IAClE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5B,MAAM,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAEjD;;;;;;;;;MASE;IACF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,IAAI,GAAG,CAAC,GAAG,QAAQ,CAAC;IACvD,IAAI,IAAI,GAAG,IAAI;QAAE,OAAO,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC;IAC5E,OAAO,IAAI,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;AAC9D,CAAC;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,SAAS,WAAW,CAAC,KAAoB,EAAE,WAAmB;IAC5D,MAAM,KAAK,GAAG,CAAC,KAAK,CAAC,gBAAgB,CAAC;IACtC,OAAO;QACL,IAAI,EAAE,cAAc;QACpB,QAAQ,EAAE,SAAS;QACnB,OAAO,EAAE,KAAK;YACZ,CAAC,CAAC,gDAAgD,KAAK,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,OAAO;gBACnF,GAAG,KAAK,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,yDAAyD;YACxF,CAAC,CAAC,oBAAoB,WAAW,kDAAkD;gBACjF,8CAA8C;QAClD,IAAI,EAAE,KAAK;YACT,CAAC,CAAC;;;;;;;;;;gBAUE;gBACF,KAAK,CAAC,UAAU,IAAI,KAAK,CAAC,qBAAqB;oBAC/C,CAAC,CAAC,4BAA4B,KAAK,CAAC,KAAK,CAAC,qBAAqB,EAAE,CAAC,CAAC,kBAAkB;wBACnF,mFAAmF;wBACnF,sCAAsC;oBACxC,CAAC,CAAC,iFAAiF;wBACjF,qFAAqF;wBACrF,2BAA2B;YAC/B,CAAC,CAAC,kDAAkD;KACvD,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CACrB,MAA4B,EAC5B,KAAoB,EACpB,cAAsB,EACtB,YAA6C,EAC7C,GAAY,EACZ,MAAuB;IAEvB,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,OAAO,GAAG,KAAK,CAAC;IACpB,oFAAoF;IACpF,sFAAsF;IACtF,IAAI,QAAQ,GAAG,CAAC,CAAC;IAEjB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,KAAK,IAAI,MAAM,GAAG,KAAK,CAAC,WAAW,EAAE,MAAM,GAAG,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,CAAC;YACxE,MAAM,WAAW,GAAG,YAAY;gBAC9B,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,cAAc,CAAC;gBACnD,CAAC,CAAC,MAAM,GAAG,cAAc,CAAC;YAC5B,SAAS,IAAI,mBAAmB,CAAC;gBAC/B,WAAW;gBACX,IAAI,EAAE,KAAK,CAAC,IAAI;gBAChB,gBAAgB,EAAE,KAAK,CAAC,gBAAgB;gBACxC,YAAY,EAAE,KAAK,CAAC,YAAY;gBAChC,UAAU,EAAE,KAAK,CAAC,UAAU;aAC7B,CAAC,CAAC;QACL,CAAC;QACD,IAAI,MAAM,KAAK,MAAM,EAAE,CAAC;YACtB,kEAAkE;YAClE,MAAM,UAAU,GAAG,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC;YACrD,IAAI,IAAI,UAAU,CAAC;YACnB,QAAQ,IAAI,UAAU,CAAC;YACvB;;;;cAIE;YACF,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,CAAC,UAAU,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;YAChE,KAAK,MAAM,OAAO,IAAI,KAAK,CAAC,QAAQ,EAAE,CAAC;gBACrC,MAAM,UAAU,GAAG,OAAO,CACxB,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EACpF,OAAO,CAAC,QAAQ,EAChB,OAAO,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,IAAI,OAAO,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,CACjE,CAAC;gBACF,MAAM,SAAS,GAAG,MAAM,CAAC,UAAU,CAAC,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;gBACpE,8BAA8B;gBAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,SAAS,GAAG,SAAS,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC;YAChE,CAAC;YACD,SAAS;QACX,CAAC;QAED,IAAI,IAAI,SAAS,CAAC;QAClB,IAAI,SAAS,GAAG,CAAC,GAAG,qBAAqB;YAAE,OAAO,GAAG,IAAI,CAAC;QAE1D;;;;;;;;;;;;;;;;UAgBE;QACF,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,CAAC,UAAU,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;QAChE,MAAM,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CACrC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CACT,GAAG;YACH,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EACxE,CAAC,CAAC,QAAQ,EACV,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,IAAI,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,CACrD,EACH,CAAC,CACF,CAAC;QACF,2CAA2C;QAC3C,KAAK,IAAI,SAAS,GAAG,CAAC,SAAS,GAAG,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACzE;;;;;;;;;;;UAWE;QACF,KAAK,IAAI,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,QAAQ,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC3F,mFAAmF;QACnF,sFAAsF;QACtF,IAAI,GAAG;YAAE,KAAK,IAAI,SAAS,CAAC;IAC9B,CAAC;IAED,IAAI,MAAM,KAAK,MAAM,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC3C,IAAI,QAAQ,GAAG,CAAC,GAAG,qBAAqB;YAAE,OAAO,GAAG,IAAI,CAAC;QACzD,KAAK,IAAI,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACvE,IAAI,GAAG;YAAE,KAAK,IAAI,SAAS,CAAC;IAC9B,CAAC;IAED,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,uBAAuB,EAAE,OAAO,EAAE,CAAC;AAC3D,CAAC","sourcesContent":["/**\n * Turning a request into a concrete conversion plan.\n *\n * The plan is where the tool's central promise is enforced: channels recorded at\n * different sampling rates are never merged into one table. A single wide CSV can\n * only hold mixed rates by inventing samples for the slow channels — MNE, for\n * instance, expands three genuine 1 Hz temperature readings into 768 interpolated\n * values without warning. Instead each distinct rate gets its own file, so every\n * number in every output file is a number that was actually recorded.\n */\n\nimport type { Diagnostic } from '../edf/errors.js';\nimport type { EdfSignal } from '../edf/header.js';\nimport { formatRate, formatRates } from '../edf/header.js';\nimport { decimalsAreClamped, decimalsForSignal } from '../edf/scale.js';\nimport { UTF8_BOM, csvRow, escapeCsvField } from '../format/csv.js';\nimport { listed } from '../format/list.js';\nimport { fixed, timeDecimals } from '../format/number.js';\nimport { buildColumnNames, renamedByCollision, selectChannels } from './channels.js';\nimport { assertOptions } from './options.js';\nimport { countSamplesInRange, resolveRange } from './time-range.js';\nimport type { ResolvedRange } from './time-range.js';\n\nexport interface PlannedChannel {\n  signal: EdfSignal;\n  column: string;\n  decimals: number;\n}\n\nexport interface RateGroup {\n  /** Sampling rate in Hz shared by every channel in this group. */\n  rate: number;\n  samplesPerRecord: number;\n  fileName: string;\n  timeDecimals: number;\n  channels: PlannedChannel[];\n}\n\nexport interface PlanInput {\n  signals: readonly EdfSignal[];\n  recordDuration: number;\n  recordCount: number;\n  hasAnnotationChannel: boolean;\n  /**\n   * True start time of each data record, supplied for discontinuous files. The\n   * requested time window is resolved against these rather than against\n   * `recordCount * recordDuration`, which for a file with gaps is the amount of\n   * data rather than the span of time it covers.\n   */\n  recordStarts?: Float64Array | null | undefined;\n}\n\nexport interface PlanOptions {\n  channels?: readonly string[] | undefined;\n  start?: number | undefined;\n  /** The `--start` value exactly as typed, for error messages. */\n  startText?: string | undefined;\n  duration?: number | undefined;\n  end?: number | undefined;\n  /** The `--end` value exactly as typed, for error messages. */\n  endText?: string | undefined;\n  annotationsOnly?: boolean | undefined;\n  /** Force a fixed number of decimals instead of deriving it per channel. */\n  decimals?: number | undefined;\n  /** Compress each CSV with gzip, giving every one of them a `.gz` name. */\n  gzip?: boolean | undefined;\n  /** Start each CSV with a UTF-8 byte order mark, so Excel reads it as UTF-8. */\n  bom?: boolean | undefined;\n  /**\n   * How the samples are arranged in the CSV.\n   *\n   * `'wide'`, the default, gives one column per channel and one file per sampling rate.\n   * `'long'` gives one file, three columns — `time_s`, `channel`, `value` — and one row per\n   * sample. See ConversionPlan.layout for why that is the only way to put channels recorded\n   * at different rates in one table without inventing samples.\n   */\n  layout?: 'wide' | 'long' | undefined;\n}\n\nexport interface ConversionPlan {\n  groups: RateGroup[];\n  /**\n   * How the samples are arranged. `'wide'` is a column per channel and a file per rate;\n   * `'long'` is `time_s,channel,value`, one row per sample, all rates in one file.\n   *\n   * The wide layout has to split a mixed-rate recording across files: a 100 Hz channel and\n   * a 1 Hz channel share no rows, and putting them in one wide table means either 99 empty\n   * cells out of every hundred or inventing the samples that would fill them. In the long\n   * layout each sample carries its own time, so nothing has to line up and nothing is\n   * invented — which also makes it the one layout `--stdout` can stream for such a file.\n   */\n  layout: 'wide' | 'long';\n\n  /**\n   * Whether the CSVs will be compressed.\n   *\n   * Recorded rather than inferred from the group file names. Under `--annotations-only`\n   * there are no groups to read it off, and `--info` named `annotations.csv` for a run that\n   * wrote `annotations.csv.gz`.\n   */\n  gzip: boolean;\n  range: ResolvedRange;\n  columnNames: Map<number, string>;\n  writeSignals: boolean;\n  diagnostics: Diagnostic[];\n  estimate: OutputEstimate;\n}\n\nexport interface OutputEstimate {\n  /** Total data rows across every signal file. */\n  rows: number;\n  /** Approximate size of the signal CSVs on disk. */\n  bytes: number;\n  /** True when any single file would exceed Excel's row limit. */\n  exceedsSpreadsheetLimit: boolean;\n}\n\nconst BOM_BYTES = Buffer.byteLength(UTF8_BOM);\n\n/** Excel and most spreadsheet tools stop at 1,048,576 rows including the header. */\nexport const SPREADSHEET_ROW_LIMIT = 1_048_576;\n\nexport function buildPlan(input: PlanInput, options: PlanOptions = {}): ConversionPlan {\n  // First, and before a directory is created or a stream opened, so a rejected option\n  // leaves nothing behind. See assertOptions for what used to get through.\n  assertOptions(options);\n\n  const diagnostics: Diagnostic[] = [];\n  const columnNames = buildColumnNames(input.signals);\n\n  // A channel whose own label was taken by another channel's disambiguating suffix. The\n  // duplicate-label warning is about the labels that collided; this is about the channel\n  // that lost its name to them, which is the one whose column no longer matches the file.\n  for (const signal of renamedByCollision(input.signals, columnNames)) {\n    diagnostics.push({\n      code: 'DUPLICATE_LABEL',\n      severity: 'warning',\n      message:\n        `Signal ${signal.index} is labelled \"${signal.label}\", which is also the column name ` +\n        `another channel's \"_ch\" suffix produces, so its column is \"${columnNames.get(signal.index)}\".`,\n      hint: 'Column names are unique; look this channel up in channels.csv by its signal_index.',\n    });\n  }\n\n  const range = resolveRange({\n    start: options.start,\n    startText: options.startText,\n    duration: options.duration,\n    end: options.end,\n    endText: options.endText,\n    recordDuration: input.recordDuration,\n    recordCount: input.recordCount,\n    recordStarts: input.recordStarts,\n  });\n\n  const writeSignals = options.annotationsOnly !== true;\n\n  let chosen: EdfSignal[] = input.signals.filter((s) => !s.isAnnotations);\n\n  /*\n    Channel names are checked even under --annotations-only, where the selection is not\n    otherwise used.\n\n    Skipping the check meant `--channels TYPO --annotations-only` exited 0 in silence while\n    the same typo without the flag was a usage error, and `--channels \"\"` stayed an error\n    in both — so a mistyped name was the one form of bad input the tool accepted quietly.\n    Everywhere else a term matching nothing is reported rather than ignored; a flag that\n    happens not to apply is a poor reason to make an exception.\n  */\n  if (options.channels && options.channels.length > 0) {\n    const selection = selectChannels(input.signals, options.channels);\n    if (writeSignals) chosen = selection.signals;\n    for (const { term, matched } of selection.ambiguous) {\n      diagnostics.push({\n        code: 'DUPLICATE_LABEL',\n        severity: 'warning',\n        message:\n          `\"${term}\" matches ${matched.length} channels (positions ` +\n          `${listed(matched.map((s) => `#${s.index}`))}); all of them were selected.`,\n        hint: `Use --channels \"#${matched[0]?.index ?? 0}\" to pick just one.`,\n      });\n    }\n  }\n\n  const layout = options.layout ?? 'wide';\n  const groups = writeSignals\n    ? groupByRate(chosen, columnNames, options.decimals, options.gzip === true, layout)\n    : [];\n  const estimate = estimateOutput(\n    groups,\n    range,\n    input.recordDuration,\n    input.recordStarts,\n    options.bom === true,\n    layout,\n  );\n\n  /*\n    The mixed-rate warning describes what this conversion does, not what the file holds.\n\n    The header parser raises its own, which is right for `parseHeader` — but it sees every\n    channel and knows nothing about `--channels`. Converting one channel out of a three-rate\n    recording therefore announced \"3 different sampling rates ... written to one file per\n    rate\" over a run that wrote one file, in the same output where `--info` had already\n    marked the other two \"(not selected)\". Selecting two of the three was wrong the other\n    way: still \"3\".\n\n    Callers combining these with a file's own diagnostics drop that copy in favour of this\n    one; see `withoutFileRateWarning`.\n  */\n  if (groups.length > 1) {\n    diagnostics.push({\n      code: 'MIXED_SAMPLING_RATES',\n      severity: 'warning',\n      message:\n        `Channels use ${groups.length} different sampling rates ` +\n        `(${listed(formatRates(groups.map((g) => g.rate)).map((r) => `${r} Hz`))}).`,\n      hint:\n        layout === 'long'\n          ? 'They share one table, each row carrying its own time, so no channel is resampled.'\n          : 'They are written to one file per rate so no channel is resampled.',\n    });\n  }\n\n  /*\n    A time column that cannot tell two samples apart.\n\n    Sample times are written to at most nine decimal places, which separates everything up to\n    a gigahertz. Below that the column repeats: a recording of 1 ns records holding ten\n    samples each writes twenty rows carrying three distinct times, so joining or plotting on\n    `time_s` silently collapses them. Nothing is lost from the file — every sample is there,\n    in order — but the column stops being an identifier, and that is worth saying rather than\n    leaving to be discovered.\n  */\n  for (const group of groups) {\n    const step = group.rate > 0 ? 1 / group.rate : 0;\n    if (step > 0 && step < 10 ** -group.timeDecimals) {\n      diagnostics.push({\n        code: 'TIME_RESOLUTION',\n        severity: 'warning',\n        message:\n          `Channels at ${formatRate(group.rate)} Hz sample faster than the time column can ` +\n          `distinguish, so consecutive rows in ${group.fileName} carry the same time_s value.`,\n        hint:\n          'Every sample is written, in order. Use the row number rather than time_s to tell ' +\n          'them apart, or convert one rate at a time with --channels.',\n      });\n    }\n  }\n\n  /*\n    The same failure as TIME_RESOLUTION, one column over.\n\n    A channel whose quantization step is below 1e-98 needs more decimals than `toFixed` can\n    print, so consecutive digital codes round to the same text and the arithmetic the FAQ\n    gives for recovering them stops working. That used to happen at 1e-20 and silently — see\n    MAX_DERIVED_DECIMALS. It is rare now, but \"rare\" is the reason to say so rather than the\n    reason not to.\n\n    Asked of the ceiling, not of the precision in use, and so asked whatever `--decimals`\n    says. `--decimals 2` on a channel needing 3 is a trade the caller made knowingly, and\n    reporting it was reporting the flag back at the person who typed it — every channel of an\n    ordinary EEG raised this, and since --strict turns any diagnostic into exit 1,\n    `--decimals 2 --strict` could not succeed on any recording at all.\n\n    0.5.10 fixed that by skipping the check whenever `--decimals` was given, which suppressed\n    the real case along with the false one: at `--decimals 20` a channel stepping by 1e-106\n    printed every code it had as `0.00000000000000000000`, and said nothing. The question is\n    not who chose the precision. It is whether any precision the tool can print would\n    separate consecutive codes.\n  */\n  for (const group of groups) {\n    const short = group.channels.filter((c) => decimalsAreClamped(c.signal));\n    if (short.length === 0) continue;\n    diagnostics.push({\n      code: 'VALUE_RESOLUTION',\n      severity: 'warning',\n      message:\n        `${listed(short.map((c) => c.column))} ${short.length === 1 ? 'steps' : 'step'} by less ` +\n        `than any number of decimals this can print, so some consecutive samples round to ` +\n        `the same value in ${group.fileName}.`,\n      hint:\n        'Every sample is written, in order, and the physical values are computed at full ' +\n        'precision either way. What is lost is only in the printed text.',\n    });\n  }\n\n  /*\n    A window that selects nothing is a fact about the plan, so the plan is where it is raised.\n\n    It was pushed by `convert()` from the rows actually written, which meant `--info` never\n    said it: `--info --start 0.31 --end 0.39` on a 10 Hz recording printed \"Would write 0\n    rows\" with no warning and exited 0 under `--strict`, while converting the same window\n    warned and exited 1. The hint says \"Run with --info to see where the records actually\n    sit\" — advising the reader into the one mode that would not tell them.\n\n    The estimate's row count is exact, which `npm run estimate` checks across every fixture\n    crossed with every option set, so raising it from the plan says the same thing the rows\n    would have.\n  */\n  if (writeSignals && groups.length > 0 && estimate.rows === 0) {\n    diagnostics.push(emptyWindow(range, input.recordCount));\n  }\n\n  if (estimate.exceedsSpreadsheetLimit) {\n    diagnostics.push({\n      code: 'LARGE_OUTPUT',\n      severity: 'warning',\n      message:\n        `At least one output file will have more than ${SPREADSHEET_ROW_LIMIT.toLocaleString('en-US')} ` +\n        `rows, which is more than Excel or Numbers can open.`,\n      hint: 'Use --start and --duration to convert a section, or read the file with pandas or R.',\n    });\n  }\n\n  return { groups, layout, gzip: options.gzip === true, range, columnNames, writeSignals, diagnostics, estimate };\n}\n\n/**\n * Partition channels by sampling rate, largest first.\n *\n * The common case — every channel at one rate — collapses to a single group and a\n * single `signals.csv`, so the honest behaviour costs nothing when there is nothing\n * to be honest about.\n */\nfunction groupByRate(\n  signals: readonly EdfSignal[],\n  columnNames: Map<number, string>,\n  forcedDecimals: number | undefined,\n  gzip: boolean,\n  layout: 'wide' | 'long',\n): RateGroup[] {\n  const byRate = new Map<number, EdfSignal[]>();\n  for (const signal of signals) {\n    // A channel with no samples has no sampling rate to group by, and would\n    // otherwise produce an empty \"0hz\" file. The header parser already warned.\n    if (signal.samplesPerRecord === 0) continue;\n    const bucket = byRate.get(signal.samplingRate);\n    if (bucket) bucket.push(signal);\n    else byRate.set(signal.samplingRate, [signal]);\n  }\n\n  const rates = [...byRate.keys()].sort((a, b) => b - a);\n  // The long layout writes one table whatever the rates are, so every group names it.\n  const single = rates.length === 1 || layout === 'long';\n\n  /*\n    Two distinct rates can produce the same slug, because the slug rounds to six decimal\n    places. Rates come from samplesPerRecord / recordDuration and every channel shares the\n    record duration, so the closest two rates can be is 1 / recordDuration — which drops\n    below 1e-6 once a record is longer than about eleven days. Absurd, but the header\n    permits it, and the failure was silent and destructive: both groups opened a write\n    stream on the same path, so the file ended up holding interleaved rows from both\n    channels under a header naming only one of them.\n\n    Distinct rates therefore get distinct files, always. The suffix is only ever reached by\n    a collision, so ordinary recordings keep the names they have always had.\n\n    Naming from the whole set of rates at once removes most of those collisions before the\n    suffix has to. Rounding each rate on its own gave 1e-6 Hz and 1.25e-6 Hz the same slug,\n    and the numbering below then produced signals_0_000001hz.csv and signals_0_000001hz_2.csv\n    — two files that no longer overwrite each other, but of which only one is named for the\n    rate it holds. The suffix stays as the backstop for anything this still cannot separate.\n  */\n  const suffix = gzip ? '.csv.gz' : '.csv';\n  const slugs = formatRates(rates).map((text) => `${text.replace('.', '_')}hz`);\n  const used = new Set<string>();\n  const uniqueName = (index: number): string => {\n    const base = `signals_${slugs[index]}`;\n    let name = `${base}${suffix}`;\n    for (let n = 2; used.has(name); n++) name = `${base}_${n}${suffix}`;\n    used.add(name);\n    return name;\n  };\n\n  return rates.map((rate, index) => {\n    const members = byRate.get(rate) ?? [];\n    const first = members[0];\n    return {\n      rate,\n      samplesPerRecord: first ? first.samplesPerRecord : 0,\n      fileName: single ? `signals${suffix}` : uniqueName(index),\n      /*\n        In the long layout every rate shares a `time_s` column, so they share its precision:\n        the finest any of them needs. Writing 100 Hz at three places and 256 Hz at eight in\n        the same column would make the column's meaning depend on the row.\n      */\n      timeDecimals: layout === 'long' ? Math.max(...rates.map(timeDecimals)) : timeDecimals(rate),\n      channels: members.map((signal) => ({\n        signal,\n        column: columnNames.get(signal.index) ?? `signal_${signal.index}`,\n        decimals: forcedDecimals ?? decimalsForSignal(signal),\n      })),\n    };\n  });\n}\n\n/**\n * A file's diagnostics with the header's mixed-rate warning removed.\n *\n * `buildPlan` raises that warning for the channels actually being converted, so keeping both\n * would either duplicate it or contradict it. The header parser's copy stays where it is, for\n * callers reading a header without planning a conversion.\n */\nexport function withoutFileRateWarning(diagnostics: readonly Diagnostic[]): Diagnostic[] {\n  return diagnostics.filter((d) => d.code !== 'MIXED_SAMPLING_RATES');\n}\n\n/** `256hz`, `12_5hz` — safe in a filename on every platform. */\nexport function rateSlug(rate: number): string {\n  return `${formatRate(rate).replace('.', '_')}hz`;\n}\n\n/** Characters a fixed-decimal number of this magnitude occupies, sign included. */\nfunction widthOf(magnitude: number, decimals: number, signed = false): number {\n  const size = Math.abs(magnitude);\n  const sign = signed ? 1 : 0;\n  const fraction = decimals > 0 ? 1 + decimals : 0;\n\n  /*\n    Cells are written with toFixed, which rounds. Taking the integer digits from the floor of\n    the bound therefore under-counted whenever rounding carried into a new digit: a channel\n    bounded at 9999.999 and written to zero decimals produces \"10000\", five characters where\n    the floor of 9999.999 suggests four. Every cell on such a channel was a byte short, and\n    `--info` reported 127 KB for a file that came out 131 KB.\n\n    Measuring the bound as rendered removes that. toFixed switches to exponential notation\n    past 1e21, so the arithmetic form still covers magnitudes beyond it.\n  */\n  if (!Number.isFinite(size)) return sign + 1 + fraction;\n  if (size < 1e21) return sign + size.toFixed(Math.min(decimals, 100)).length;\n  return sign + (Math.floor(Math.log10(size)) + 1) + fraction;\n}\n\n/**\n * Raised when the conversion had signal tables to fill and put no data rows in any of them.\n *\n * A window can land where there are no samples without being past the end of the recording:\n * between the last sample and the nominal end of the last record, or — on a discontinuous\n * file — inside a gap. `--start 2 --end 10` on a recording whose records sit at 0s, 1s and\n * 10s asks for eight seconds that contain no data at all.\n *\n * What came out was a signals.csv holding its header and nothing else, exit 0, no warning,\n * and `--strict` passing. The closing summary does say \"signals.csv 0 rows\" and --json\n * carries `rows: 0`, so it was not quite invisible — but a header-only file is exactly what\n * a successful extraction of an empty range looks like, and everywhere else that a request\n * produces nothing this tool says so: a --channels term matching nothing is an error, and\n * --annotations-only on a file with no events raises NO_ANNOTATIONS. A warning rather than\n * an error because a batch of five hundred recordings should not stop for one whose gap\n * happens to line up with the window; --strict turns it into a failure for those who want\n * that.\n */\nfunction emptyWindow(range: ResolvedRange, recordCount: number): Diagnostic {\n  const asked = !range.isWholeRecording;\n  return {\n    code: 'EMPTY_WINDOW',\n    severity: 'warning',\n    message: asked\n      ? `No samples fall inside the requested window (${fixed(range.startSeconds, 3)}s to ` +\n        `${fixed(range.endSeconds, 3)}s), so the signal files hold their headers and no data.`\n      : `This recording's ${recordCount} data records carry no samples in range, so the ` +\n        `signal files hold their headers and no data.`,\n    hint: asked\n      ? /*\n          Which of the two it was, rather than the second one always.\n\n          A recording does not have to start at zero: its first record's timekeeping TAL is\n          what it is timed from, so a file whose records begin at 1000s is asked for with\n          `--start 1000`. `--start 0 --end 1` on that file was told \"The window is inside the\n          recording but lands where there is no data — past the last sample, or inside a gap\n          in a discontinuous file\", when the window sits entirely before the recording and\n          neither offered explanation applies to it. A start at or past the *end* is already\n          an error, so the window being outside can only mean it is before the beginning.\n        */\n        range.endSeconds <= range.recordingStartSeconds\n        ? `This recording starts at ${fixed(range.recordingStartSeconds, 3)}s, so the whole ` +\n          'window sits before it. --start and --end are read on the recording\\'s own clock, ' +\n          'which --info prints as \"Timed from\".'\n        : 'The window is inside the recording but lands where there is no data — past the ' +\n          'last sample, or inside a gap in a discontinuous file. Run with --info to see where ' +\n          'the records actually sit.'\n      : 'Run with --info to see what the header declares.',\n  };\n}\n\nfunction estimateOutput(\n  groups: readonly RateGroup[],\n  range: ResolvedRange,\n  recordDuration: number,\n  recordStarts: Float64Array | null | undefined,\n  bom: boolean,\n  layout: 'wide' | 'long',\n): OutputEstimate {\n  let rows = 0;\n  let bytes = 0;\n  let exceeds = false;\n  // One table in the long layout, so the row limit applies to the sum rather than the\n  // largest group, and the header and mark are counted once rather than once per group.\n  let longRows = 0;\n\n  for (const group of groups) {\n    let groupRows = 0;\n    for (let record = range.startRecord; record < range.endRecord; record++) {\n      const recordStart = recordStarts\n        ? (recordStarts[record] ?? record * recordDuration)\n        : record * recordDuration;\n      groupRows += countSamplesInRange({\n        recordStart,\n        rate: group.rate,\n        samplesPerRecord: group.samplesPerRecord,\n        startSeconds: range.startSeconds,\n        endSeconds: range.endSeconds,\n      });\n    }\n    if (layout === 'long') {\n      // A row per sample per channel rather than a row per sample time.\n      const groupCells = groupRows * group.channels.length;\n      rows += groupCells;\n      longRows += groupCells;\n      /*\n        `time_s,channel,value`: the time, the channel name as it will be escaped into the\n        cell, and the widest the value can print. Same over-counting rule as the wide\n        layout — the declared physical range bounds a cell, and most samples sit under it.\n      */\n      const timeWidth = widthOf(range.endSeconds, group.timeDecimals);\n      for (const channel of group.channels) {\n        const valueWidth = widthOf(\n          Math.max(Math.abs(channel.signal.physicalMin), Math.abs(channel.signal.physicalMax)),\n          channel.decimals,\n          channel.signal.physicalMin < 0 || channel.signal.physicalMax < 0,\n        );\n        const nameWidth = Buffer.byteLength(escapeCsvField(channel.column));\n        // Two commas and the newline.\n        bytes += groupRows * (timeWidth + nameWidth + valueWidth + 3);\n      }\n      continue;\n    }\n\n    rows += groupRows;\n    if (groupRows + 1 > SPREADSHEET_ROW_LIMIT) exceeds = true;\n\n    /*\n      Width per cell, from the channel's own calibration rather than a flat allowance.\n\n      The old `decimals + 6` budgeted six characters for the sign, integer part and decimal\n      point on every channel, whatever it actually held. That over-counted a millivolt\n      channel spanning ±5 by four characters a cell and ran 30-55% high across the fixture\n      set — on a number people use to decide whether a conversion is worth starting.\n\n      The channel's declared physical range is what bounds a cell, so that bound is what is\n      used. Most samples sit below it, so this still reads high, which is the direction a\n      size estimate should err in.\n\n      One case is outside the bound rather than under it: nothing obliges a recording to keep\n      its samples inside the digital range it declares, and one that does not maps outside the\n      physical range too. Such a file can convert larger than the estimate. Clamping the data\n      to make the estimate true is not a trade worth making — the samples are what they are.\n    */\n    const timeWidth = widthOf(range.endSeconds, group.timeDecimals);\n    const cellWidth = group.channels.reduce(\n      (sum, c) =>\n        sum +\n        widthOf(\n          Math.max(Math.abs(c.signal.physicalMin), Math.abs(c.signal.physicalMax)),\n          c.decimals,\n          c.signal.physicalMin < 0 || c.signal.physicalMax < 0,\n        ),\n      0,\n    );\n    // One comma per channel, plus the newline.\n    bytes += groupRows * (timeWidth + cellWidth + group.channels.length + 1);\n    /*\n      The header row, measured as it will be written rather than as the labels are stored.\n\n      A column name is quoted when it contains a comma, a quote, a newline or a leading or\n      trailing space, and every quote inside it is doubled. Counting the raw label under-counted\n      that row: three channels labelled `a,b,c,d,e`, `x\"y` and `plain` write a 32-byte header\n      and were budgeted 27. EDF labels are free text, so commas in them are ordinary — a montage\n      written as `EEG Fpz-Cz, ref` is exactly the kind of thing this is for.\n\n      csvRow is the function that writes it, so it is the function that measures it. Nothing\n      else is in a position to stay correct when the quoting rules change.\n    */\n    bytes += Buffer.byteLength(csvRow(['time_s', ...group.channels.map((c) => c.column)])) + 1;\n    // Three bytes per file under --bom. Small, but the estimate promises never to read\n    // under what gets written, and a one-row conversion is small enough for it to matter.\n    if (bom) bytes += BOM_BYTES;\n  }\n\n  if (layout === 'long' && groups.length > 0) {\n    if (longRows + 1 > SPREADSHEET_ROW_LIMIT) exceeds = true;\n    bytes += Buffer.byteLength(csvRow(['time_s', 'channel', 'value'])) + 1;\n    if (bom) bytes += BOM_BYTES;\n  }\n\n  return { rows, bytes, exceedsSpreadsheetLimit: exceeds };\n}\n"]}
+{"version":3,"file":"plan.js","sourceRoot":"","sources":["../../src/convert/plan.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAIH,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC3D,OAAO,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AACxE,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AACpE,OAAO,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,YAAY,EAAE,MAAM,qBAAqB,CAAC;AAC1D,OAAO,EAAE,gBAAgB,EAAE,kBAAkB,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AACrF,OAAO,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAC7C,OAAO,EAAE,mBAAmB,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AAiGpE,MAAM,SAAS,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC;AAE9C,oFAAoF;AACpF,MAAM,CAAC,MAAM,qBAAqB,GAAG,SAAS,CAAC;AAE/C,MAAM,UAAU,SAAS,CAAC,KAAgB,EAAE,UAAuB,EAAE;IACnE,oFAAoF;IACpF,yEAAyE;IACzE,aAAa,CAAC,OAAO,CAAC,CAAC;IAEvB,MAAM,WAAW,GAAiB,EAAE,CAAC;IACrC,MAAM,WAAW,GAAG,gBAAgB,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IAEpD,sFAAsF;IACtF,uFAAuF;IACvF,wFAAwF;IACxF,KAAK,MAAM,MAAM,IAAI,kBAAkB,CAAC,KAAK,CAAC,OAAO,EAAE,WAAW,CAAC,EAAE,CAAC;QACpE,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,iBAAiB;YACvB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,UAAU,MAAM,CAAC,KAAK,iBAAiB,MAAM,CAAC,KAAK,mCAAmC;gBACtF,8DAA8D,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI;YACjG,IAAI,EAAE,oFAAoF;SAC3F,CAAC,CAAC;IACL,CAAC;IAED,MAAM,KAAK,GAAG,YAAY,CAAC;QACzB,KAAK,EAAE,OAAO,CAAC,KAAK;QACpB,SAAS,EAAE,OAAO,CAAC,SAAS;QAC5B,QAAQ,EAAE,OAAO,CAAC,QAAQ;QAC1B,GAAG,EAAE,OAAO,CAAC,GAAG;QAChB,OAAO,EAAE,OAAO,CAAC,OAAO;QACxB,cAAc,EAAE,KAAK,CAAC,cAAc;QACpC,WAAW,EAAE,KAAK,CAAC,WAAW;QAC9B,YAAY,EAAE,KAAK,CAAC,YAAY;KACjC,CAAC,CAAC;IAEH,MAAM,YAAY,GAAG,OAAO,CAAC,eAAe,KAAK,IAAI,CAAC;IAEtD,IAAI,MAAM,GAAgB,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC;IAExE;;;;;;;;;MASE;IACF,IAAI,OAAO,CAAC,QAAQ,IAAI,OAAO,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACpD,MAAM,SAAS,GAAG,cAAc,CAAC,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;QAClE,IAAI,YAAY;YAAE,MAAM,GAAG,SAAS,CAAC,OAAO,CAAC;QAC7C,KAAK,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,SAAS,CAAC,SAAS,EAAE,CAAC;YACpD,WAAW,CAAC,IAAI,CAAC;gBACf,IAAI,EAAE,iBAAiB;gBACvB,QAAQ,EAAE,SAAS;gBACnB,OAAO,EACL,IAAI,IAAI,aAAa,OAAO,CAAC,MAAM,uBAAuB;oBAC1D,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,+BAA+B;gBAC7E,IAAI,EAAE,oBAAoB,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,IAAI,CAAC,qBAAqB;aACtE,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IAED,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,MAAM,CAAC;IACxC,MAAM,MAAM,GAAG,YAAY;QACzB,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,WAAW,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,IAAI,KAAK,IAAI,EAAE,MAAM,CAAC;QACnF,CAAC,CAAC,EAAE,CAAC;IACP,MAAM,QAAQ,GAAG,cAAc,CAC7B,MAAM,EACN,KAAK,EACL,KAAK,CAAC,cAAc,EACpB,KAAK,CAAC,YAAY,EAClB,OAAO,CAAC,GAAG,KAAK,IAAI,EACpB,MAAM,CACP,CAAC;IAEF;;;;;;;;;;;;MAYE;IACF,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACtB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gBAAgB,MAAM,CAAC,MAAM,4BAA4B;gBACzD,IAAI,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI;YAC9E,IAAI,EACF,MAAM,KAAK,MAAM;gBACf,CAAC,CAAC,mFAAmF;gBACrF,CAAC,CAAC,mEAAmE;SAC1E,CAAC,CAAC;IACL,CAAC;IAED;;;;;;;;;MASE;IACF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACjD;;;;;;;;;;;;;;;UAeE;QACF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;YACjC,WAAW,CAAC,IAAI,CAAC;gBACf,IAAI,EAAE,iBAAiB;gBACvB,QAAQ,EAAE,SAAS;gBACnB,OAAO,EACL,eAAe,KAAK,CAAC,QAAQ,mCAAmC,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM;oBAC5F,kFAAkF;oBAClF,0EAA0E;gBAC5E,IAAI,EACF,gFAAgF;oBAChF,gEAAgE;aACnE,CAAC,CAAC;QACL,CAAC;aAAM,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC;YACxD,WAAW,CAAC,IAAI,CAAC;gBACf,IAAI,EAAE,iBAAiB;gBACvB,QAAQ,EAAE,SAAS;gBACnB,OAAO,EACL,eAAe,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,6CAA6C;oBAClF,uCAAuC,KAAK,CAAC,QAAQ,+BAA+B;gBACtF,IAAI,EACF,mFAAmF;oBACnF,4DAA4D;aAC/D,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IAED;;;;;;;;;;;;;;;;;;;;MAoBE;IACF,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,kBAAkB,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;QACzE,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;YAAE,SAAS;QACjC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,kBAAkB;YACxB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,GAAG,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,WAAW;gBACzF,mFAAmF;gBACnF,qBAAqB,KAAK,CAAC,QAAQ,GAAG;YACxC,IAAI,EACF,kFAAkF;gBAClF,iEAAiE;SACpE,CAAC,CAAC;IACL,CAAC;IAED;;;;;;;;;;;;MAYE;IACF;;;;;MAKE;IACF,MAAM,UAAU,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;IACxE,IAAI,YAAY,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,IAAI,QAAQ,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,UAAU,EAAE,CAAC;QAC5E,WAAW,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC;IAC1D,CAAC;IAED,IAAI,QAAQ,CAAC,uBAAuB,EAAE,CAAC;QACrC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,cAAc;YACpB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gDAAgD,qBAAqB,CAAC,cAAc,CAAC,OAAO,CAAC,GAAG;gBAChG,qDAAqD;YACvD,IAAI,EAAE,qFAAqF;SAC5F,CAAC,CAAC;IACL,CAAC;IAED,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,OAAO,CAAC,IAAI,KAAK,IAAI,EAAE,KAAK,EAAE,WAAW,EAAE,YAAY,EAAE,WAAW,EAAE,QAAQ,EAAE,CAAC;AAClH,CAAC;AAED;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,OAA6B,EAC7B,WAAgC,EAChC,cAAkC,EAClC,IAAa,EACb,MAAuB;IAEvB,MAAM,MAAM,GAAG,IAAI,GAAG,EAAuB,CAAC;IAC9C,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE,CAAC;QAC7B,wEAAwE;QACxE,2EAA2E;QAC3E,IAAI,MAAM,CAAC,gBAAgB,KAAK,CAAC;YAAE,SAAS;QAC5C,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;QAC/C,IAAI,MAAM;YAAE,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;;YAC3B,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC;IACjD,CAAC;IAED,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACvD,oFAAoF;IACpF,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,MAAM,KAAK,MAAM,CAAC;IAEvD;;;;;;;;;;;;;;;;;MAiBE;IACF,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC;IACzC,MAAM,KAAK,GAAG,WAAW,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;IAC9E,MAAM,IAAI,GAAG,IAAI,GAAG,EAAU,CAAC;IAC/B,MAAM,UAAU,GAAG,CAAC,KAAa,EAAU,EAAE;QAC3C,MAAM,IAAI,GAAG,WAAW,KAAK,CAAC,KAAK,CAAC,EAAE,CAAC;QACvC,IAAI,IAAI,GAAG,GAAG,IAAI,GAAG,MAAM,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,IAAI,GAAG,GAAG,IAAI,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;QACpE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC,CAAC;IAEF,OAAO,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE;QAC/B,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACzB,OAAO;YACL,IAAI;YACJ,gBAAgB,EAAE,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;YACpD,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC,UAAU,MAAM,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC;YACzD;;;;cAIE;YACF,YAAY,EAAE,MAAM,KAAK,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,IAAI,CAAC;YAC3F,QAAQ,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;gBACjC,MAAM;gBACN,MAAM,EAAE,WAAW,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,UAAU,MAAM,CAAC,KAAK,EAAE;gBACjE,QAAQ,EAAE,cAAc,IAAI,iBAAiB,CAAC,MAAM,CAAC;aACtD,CAAC,CAAC;SACJ,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,sBAAsB,CAAC,WAAkC;IACvE,OAAO,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,sBAAsB,CAAC,CAAC;AACtE,CAAC;AAED,gEAAgE;AAChE,MAAM,UAAU,QAAQ,CAAC,IAAY;IACnC,OAAO,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC;AACnD,CAAC;AAED,mFAAmF;AACnF,SAAS,OAAO,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAAM,GAAG,KAAK;IAClE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5B,MAAM,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAEjD;;;;;;;;;MASE;IACF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,IAAI,GAAG,CAAC,GAAG,QAAQ,CAAC;IACvD,IAAI,IAAI,GAAG,IAAI;QAAE,OAAO,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC;IAC5E,OAAO,IAAI,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;AAC9D,CAAC;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,SAAS,WAAW,CAAC,KAAoB,EAAE,WAAmB;IAC5D,MAAM,KAAK,GAAG,CAAC,KAAK,CAAC,gBAAgB,CAAC;IACtC,OAAO;QACL,IAAI,EAAE,cAAc;QACpB,QAAQ,EAAE,SAAS;QACnB,OAAO,EAAE,KAAK;YACZ,CAAC,CAAC,gDAAgD,KAAK,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC,OAAO;gBACnF,GAAG,KAAK,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,yDAAyD;YACxF,CAAC,CAAC,oBAAoB,WAAW,kDAAkD;gBACjF,8CAA8C;QAClD,IAAI,EAAE,KAAK;YACT,CAAC,CAAC;;;;;;;;;;gBAUE;gBACF,KAAK,CAAC,UAAU,IAAI,KAAK,CAAC,qBAAqB;oBAC/C,CAAC,CAAC,4BAA4B,KAAK,CAAC,KAAK,CAAC,qBAAqB,EAAE,CAAC,CAAC,kBAAkB;wBACnF,mFAAmF;wBACnF,sCAAsC;oBACxC,CAAC,CAAC,iFAAiF;wBACjF,qFAAqF;wBACrF,2BAA2B;YAC/B,CAAC,CAAC,kDAAkD;KACvD,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CACrB,MAA4B,EAC5B,KAAoB,EACpB,cAAsB,EACtB,YAA6C,EAC7C,GAAY,EACZ,MAAuB;IAEvB,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,OAAO,GAAG,KAAK,CAAC;IACpB,oFAAoF;IACpF,sFAAsF;IACtF,IAAI,QAAQ,GAAG,CAAC,CAAC;IAEjB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,KAAK,IAAI,MAAM,GAAG,KAAK,CAAC,WAAW,EAAE,MAAM,GAAG,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,CAAC;YACxE,MAAM,WAAW,GAAG,YAAY;gBAC9B,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,cAAc,CAAC;gBACnD,CAAC,CAAC,MAAM,GAAG,cAAc,CAAC;YAC5B,SAAS,IAAI,mBAAmB,CAAC;gBAC/B,WAAW;gBACX,IAAI,EAAE,KAAK,CAAC,IAAI;gBAChB,gBAAgB,EAAE,KAAK,CAAC,gBAAgB;gBACxC,YAAY,EAAE,KAAK,CAAC,YAAY;gBAChC,UAAU,EAAE,KAAK,CAAC,UAAU;aAC7B,CAAC,CAAC;QACL,CAAC;QACD,IAAI,MAAM,KAAK,MAAM,EAAE,CAAC;YACtB,kEAAkE;YAClE,MAAM,UAAU,GAAG,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CAAC;YACrD,IAAI,IAAI,UAAU,CAAC;YACnB,QAAQ,IAAI,UAAU,CAAC;YACvB;;;;cAIE;YACF,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,CAAC,UAAU,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;YAChE,KAAK,MAAM,OAAO,IAAI,KAAK,CAAC,QAAQ,EAAE,CAAC;gBACrC,MAAM,UAAU,GAAG,OAAO,CACxB,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EACpF,OAAO,CAAC,QAAQ,EAChB,OAAO,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,IAAI,OAAO,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,CACjE,CAAC;gBACF,MAAM,SAAS,GAAG,MAAM,CAAC,UAAU,CAAC,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;gBACpE,8BAA8B;gBAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,SAAS,GAAG,SAAS,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC;YAChE,CAAC;YACD,SAAS;QACX,CAAC;QAED,IAAI,IAAI,SAAS,CAAC;QAClB,IAAI,SAAS,GAAG,CAAC,GAAG,qBAAqB;YAAE,OAAO,GAAG,IAAI,CAAC;QAE1D;;;;;;;;;;;;;;;;UAgBE;QACF,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,CAAC,UAAU,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;QAChE,MAAM,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,CACrC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CACT,GAAG;YACH,OAAO,CACL,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,WAAW,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,WAAW,CAAC,CAAC,EACxE,CAAC,CAAC,QAAQ,EACV,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,IAAI,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,CAAC,CACrD,EACH,CAAC,CACF,CAAC;QACF,2CAA2C;QAC3C,KAAK,IAAI,SAAS,GAAG,CAAC,SAAS,GAAG,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACzE;;;;;;;;;;;UAWE;QACF,KAAK,IAAI,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,QAAQ,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC3F,mFAAmF;QACnF,sFAAsF;QACtF,IAAI,GAAG;YAAE,KAAK,IAAI,SAAS,CAAC;IAC9B,CAAC;IAED,IAAI,MAAM,KAAK,MAAM,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC3C,IAAI,QAAQ,GAAG,CAAC,GAAG,qBAAqB;YAAE,OAAO,GAAG,IAAI,CAAC;QACzD,KAAK,IAAI,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACvE,IAAI,GAAG;YAAE,KAAK,IAAI,SAAS,CAAC;IAC9B,CAAC;IAED,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,uBAAuB,EAAE,OAAO,EAAE,CAAC;AAC3D,CAAC","sourcesContent":["/**\n * Turning a request into a concrete conversion plan.\n *\n * The plan is where the tool's central promise is enforced: channels recorded at\n * different sampling rates are never merged into one table. A single wide CSV can\n * only hold mixed rates by inventing samples for the slow channels — MNE, for\n * instance, expands three genuine 1 Hz temperature readings into 768 interpolated\n * values without warning. Instead each distinct rate gets its own file, so every\n * number in every output file is a number that was actually recorded.\n */\n\nimport type { Diagnostic } from '../edf/errors.js';\nimport type { EdfSignal } from '../edf/header.js';\nimport { formatRate, formatRates } from '../edf/header.js';\nimport { decimalsAreClamped, decimalsForSignal } from '../edf/scale.js';\nimport { UTF8_BOM, csvRow, escapeCsvField } from '../format/csv.js';\nimport { listed } from '../format/list.js';\nimport { fixed, timeDecimals } from '../format/number.js';\nimport { buildColumnNames, renamedByCollision, selectChannels } from './channels.js';\nimport { assertOptions } from './options.js';\nimport { countSamplesInRange, resolveRange } from './time-range.js';\nimport type { ResolvedRange } from './time-range.js';\n\nexport interface PlannedChannel {\n  signal: EdfSignal;\n  column: string;\n  decimals: number;\n}\n\nexport interface RateGroup {\n  /** Sampling rate in Hz shared by every channel in this group. */\n  rate: number;\n  samplesPerRecord: number;\n  fileName: string;\n  timeDecimals: number;\n  channels: PlannedChannel[];\n}\n\nexport interface PlanInput {\n  signals: readonly EdfSignal[];\n  recordDuration: number;\n  recordCount: number;\n  hasAnnotationChannel: boolean;\n  /**\n   * True start time of each data record, supplied for discontinuous files. The\n   * requested time window is resolved against these rather than against\n   * `recordCount * recordDuration`, which for a file with gaps is the amount of\n   * data rather than the span of time it covers.\n   */\n  recordStarts?: Float64Array | null | undefined;\n}\n\nexport interface PlanOptions {\n  channels?: readonly string[] | undefined;\n  start?: number | undefined;\n  /** The `--start` value exactly as typed, for error messages. */\n  startText?: string | undefined;\n  duration?: number | undefined;\n  end?: number | undefined;\n  /** The `--end` value exactly as typed, for error messages. */\n  endText?: string | undefined;\n  annotationsOnly?: boolean | undefined;\n  /** Force a fixed number of decimals instead of deriving it per channel. */\n  decimals?: number | undefined;\n  /** Compress each CSV with gzip, giving every one of them a `.gz` name. */\n  gzip?: boolean | undefined;\n  /** Start each CSV with a UTF-8 byte order mark, so Excel reads it as UTF-8. */\n  bom?: boolean | undefined;\n  /**\n   * How the samples are arranged in the CSV.\n   *\n   * `'wide'`, the default, gives one column per channel and one file per sampling rate.\n   * `'long'` gives one file, three columns — `time_s`, `channel`, `value` — and one row per\n   * sample. See ConversionPlan.layout for why that is the only way to put channels recorded\n   * at different rates in one table without inventing samples.\n   */\n  layout?: 'wide' | 'long' | undefined;\n}\n\nexport interface ConversionPlan {\n  groups: RateGroup[];\n  /**\n   * How the samples are arranged. `'wide'` is a column per channel and a file per rate;\n   * `'long'` is `time_s,channel,value`, one row per sample, all rates in one file.\n   *\n   * The wide layout has to split a mixed-rate recording across files: a 100 Hz channel and\n   * a 1 Hz channel share no rows, and putting them in one wide table means either 99 empty\n   * cells out of every hundred or inventing the samples that would fill them. In the long\n   * layout each sample carries its own time, so nothing has to line up and nothing is\n   * invented — which also makes it the one layout `--stdout` can stream for such a file.\n   */\n  layout: 'wide' | 'long';\n\n  /**\n   * Whether the CSVs will be compressed.\n   *\n   * Recorded rather than inferred from the group file names. Under `--annotations-only`\n   * there are no groups to read it off, and `--info` named `annotations.csv` for a run that\n   * wrote `annotations.csv.gz`.\n   */\n  gzip: boolean;\n  range: ResolvedRange;\n  columnNames: Map<number, string>;\n  writeSignals: boolean;\n  diagnostics: Diagnostic[];\n  estimate: OutputEstimate;\n}\n\nexport interface OutputEstimate {\n  /** Total data rows across every signal file. */\n  rows: number;\n  /** Approximate size of the signal CSVs on disk. */\n  bytes: number;\n  /** True when any single file would exceed Excel's row limit. */\n  exceedsSpreadsheetLimit: boolean;\n}\n\nconst BOM_BYTES = Buffer.byteLength(UTF8_BOM);\n\n/** Excel and most spreadsheet tools stop at 1,048,576 rows including the header. */\nexport const SPREADSHEET_ROW_LIMIT = 1_048_576;\n\nexport function buildPlan(input: PlanInput, options: PlanOptions = {}): ConversionPlan {\n  // First, and before a directory is created or a stream opened, so a rejected option\n  // leaves nothing behind. See assertOptions for what used to get through.\n  assertOptions(options);\n\n  const diagnostics: Diagnostic[] = [];\n  const columnNames = buildColumnNames(input.signals);\n\n  // A channel whose own label was taken by another channel's disambiguating suffix. The\n  // duplicate-label warning is about the labels that collided; this is about the channel\n  // that lost its name to them, which is the one whose column no longer matches the file.\n  for (const signal of renamedByCollision(input.signals, columnNames)) {\n    diagnostics.push({\n      code: 'DUPLICATE_LABEL',\n      severity: 'warning',\n      message:\n        `Signal ${signal.index} is labelled \"${signal.label}\", which is also the column name ` +\n        `another channel's \"_ch\" suffix produces, so its column is \"${columnNames.get(signal.index)}\".`,\n      hint: 'Column names are unique; look this channel up in channels.csv by its signal_index.',\n    });\n  }\n\n  const range = resolveRange({\n    start: options.start,\n    startText: options.startText,\n    duration: options.duration,\n    end: options.end,\n    endText: options.endText,\n    recordDuration: input.recordDuration,\n    recordCount: input.recordCount,\n    recordStarts: input.recordStarts,\n  });\n\n  const writeSignals = options.annotationsOnly !== true;\n\n  let chosen: EdfSignal[] = input.signals.filter((s) => !s.isAnnotations);\n\n  /*\n    Channel names are checked even under --annotations-only, where the selection is not\n    otherwise used.\n\n    Skipping the check meant `--channels TYPO --annotations-only` exited 0 in silence while\n    the same typo without the flag was a usage error, and `--channels \"\"` stayed an error\n    in both — so a mistyped name was the one form of bad input the tool accepted quietly.\n    Everywhere else a term matching nothing is reported rather than ignored; a flag that\n    happens not to apply is a poor reason to make an exception.\n  */\n  if (options.channels && options.channels.length > 0) {\n    const selection = selectChannels(input.signals, options.channels);\n    if (writeSignals) chosen = selection.signals;\n    for (const { term, matched } of selection.ambiguous) {\n      diagnostics.push({\n        code: 'DUPLICATE_LABEL',\n        severity: 'warning',\n        message:\n          `\"${term}\" matches ${matched.length} channels (positions ` +\n          `${listed(matched.map((s) => `#${s.index}`))}); all of them were selected.`,\n        hint: `Use --channels \"#${matched[0]?.index ?? 0}\" to pick just one.`,\n      });\n    }\n  }\n\n  const layout = options.layout ?? 'wide';\n  const groups = writeSignals\n    ? groupByRate(chosen, columnNames, options.decimals, options.gzip === true, layout)\n    : [];\n  const estimate = estimateOutput(\n    groups,\n    range,\n    input.recordDuration,\n    input.recordStarts,\n    options.bom === true,\n    layout,\n  );\n\n  /*\n    The mixed-rate warning describes what this conversion does, not what the file holds.\n\n    The header parser raises its own, which is right for `parseHeader` — but it sees every\n    channel and knows nothing about `--channels`. Converting one channel out of a three-rate\n    recording therefore announced \"3 different sampling rates ... written to one file per\n    rate\" over a run that wrote one file, in the same output where `--info` had already\n    marked the other two \"(not selected)\". Selecting two of the three was wrong the other\n    way: still \"3\".\n\n    Callers combining these with a file's own diagnostics drop that copy in favour of this\n    one; see `withoutFileRateWarning`.\n  */\n  if (groups.length > 1) {\n    diagnostics.push({\n      code: 'MIXED_SAMPLING_RATES',\n      severity: 'warning',\n      message:\n        `Channels use ${groups.length} different sampling rates ` +\n        `(${listed(formatRates(groups.map((g) => g.rate)).map((r) => `${r} Hz`))}).`,\n      hint:\n        layout === 'long'\n          ? 'They share one table, each row carrying its own time, so no channel is resampled.'\n          : 'They are written to one file per rate so no channel is resampled.',\n    });\n  }\n\n  /*\n    A time column that cannot tell two samples apart.\n\n    Sample times are written to at most nine decimal places, which separates everything up to\n    a gigahertz. Below that the column repeats: a recording of 1 ns records holding ten\n    samples each writes twenty rows carrying three distinct times, so joining or plotting on\n    `time_s` silently collapses them. Nothing is lost from the file — every sample is there,\n    in order — but the column stops being an identifier, and that is worth saying rather than\n    leaving to be discovered.\n  */\n  for (const group of groups) {\n    const step = group.rate > 0 ? 1 / group.rate : 0;\n    /*\n      The limit of the same failure, which read as the absence of it.\n\n      `samplesPerRecord / recordDuration` is a double, and a record duration of 1e-308 with\n      four samples in it is Infinity. `1 / Infinity` is 0, so `step > 0` was false and this\n      said nothing — while every sample was dropped, the run exited 0, and the only warning\n      printed was EMPTY_WINDOW's \"This recording's 2 data records carry no samples in range\",\n      which is untrue twice over: the records carry eight samples and no range was asked for.\n\n      One power of ten away, at 1e-300, the rate is 4e300 and the file converts with the\n      warning below. Same guard `decimalsAreClamped` had before 0.5.83, in the column next\n      door: a step of exactly zero means no resolution at all, not nothing to report.\n\n      Its own branch because the hint below is false here — no rows are written at all, so\n      \"Every sample is written, in order\" would be the third untrue sentence.\n    */\n    if (!Number.isFinite(group.rate)) {\n      diagnostics.push({\n        code: 'TIME_RESOLUTION',\n        severity: 'warning',\n        message:\n          `Channels in ${group.fileName} work out to a sampling rate of ${formatRate(group.rate)} Hz ` +\n          `— their samples per record over a record duration too small to divide into — so ` +\n          `their samples cannot be placed in time and no rows are written for them.`,\n        hint:\n          'Check the record duration in the header. One power of ten larger and the same ' +\n          'file converts, with consecutive rows carrying the same time_s.',\n      });\n    } else if (step > 0 && step < 10 ** -group.timeDecimals) {\n      diagnostics.push({\n        code: 'TIME_RESOLUTION',\n        severity: 'warning',\n        message:\n          `Channels at ${formatRate(group.rate)} Hz sample faster than the time column can ` +\n          `distinguish, so consecutive rows in ${group.fileName} carry the same time_s value.`,\n        hint:\n          'Every sample is written, in order. Use the row number rather than time_s to tell ' +\n          'them apart, or convert one rate at a time with --channels.',\n      });\n    }\n  }\n\n  /*\n    The same failure as TIME_RESOLUTION, one column over.\n\n    A channel whose quantization step is below 1e-98 needs more decimals than `toFixed` can\n    print, so consecutive digital codes round to the same text and the arithmetic the FAQ\n    gives for recovering them stops working. That used to happen at 1e-20 and silently — see\n    MAX_DERIVED_DECIMALS. It is rare now, but \"rare\" is the reason to say so rather than the\n    reason not to.\n\n    Asked of the ceiling, not of the precision in use, and so asked whatever `--decimals`\n    says. `--decimals 2` on a channel needing 3 is a trade the caller made knowingly, and\n    reporting it was reporting the flag back at the person who typed it — every channel of an\n    ordinary EEG raised this, and since --strict turns any diagnostic into exit 1,\n    `--decimals 2 --strict` could not succeed on any recording at all.\n\n    0.5.10 fixed that by skipping the check whenever `--decimals` was given, which suppressed\n    the real case along with the false one: at `--decimals 20` a channel stepping by 1e-106\n    printed every code it had as `0.00000000000000000000`, and said nothing. The question is\n    not who chose the precision. It is whether any precision the tool can print would\n    separate consecutive codes.\n  */\n  for (const group of groups) {\n    const short = group.channels.filter((c) => decimalsAreClamped(c.signal));\n    if (short.length === 0) continue;\n    diagnostics.push({\n      code: 'VALUE_RESOLUTION',\n      severity: 'warning',\n      message:\n        `${listed(short.map((c) => c.column))} ${short.length === 1 ? 'steps' : 'step'} by less ` +\n        `than any number of decimals this can print, so some consecutive samples round to ` +\n        `the same value in ${group.fileName}.`,\n      hint:\n        'Every sample is written, in order, and the physical values are computed at full ' +\n        'precision either way. What is lost is only in the printed text.',\n    });\n  }\n\n  /*\n    A window that selects nothing is a fact about the plan, so the plan is where it is raised.\n\n    It was pushed by `convert()` from the rows actually written, which meant `--info` never\n    said it: `--info --start 0.31 --end 0.39` on a 10 Hz recording printed \"Would write 0\n    rows\" with no warning and exited 0 under `--strict`, while converting the same window\n    warned and exited 1. The hint says \"Run with --info to see where the records actually\n    sit\" — advising the reader into the one mode that would not tell them.\n\n    The estimate's row count is exact, which `npm run estimate` checks across every fixture\n    crossed with every option set, so raising it from the plan says the same thing the rows\n    would have.\n  */\n  /*\n    Not when a rate above already explained it. EMPTY_WINDOW says the records \"carry no\n    samples in range\", and on a recording whose rate overflowed to Infinity that is untrue\n    twice: the records carry their samples, and no range was asked for. The rate warning is\n    the accurate account of the same zero.\n  */\n  const untimeable = groups.some((group) => !Number.isFinite(group.rate));\n  if (writeSignals && groups.length > 0 && estimate.rows === 0 && !untimeable) {\n    diagnostics.push(emptyWindow(range, input.recordCount));\n  }\n\n  if (estimate.exceedsSpreadsheetLimit) {\n    diagnostics.push({\n      code: 'LARGE_OUTPUT',\n      severity: 'warning',\n      message:\n        `At least one output file will have more than ${SPREADSHEET_ROW_LIMIT.toLocaleString('en-US')} ` +\n        `rows, which is more than Excel or Numbers can open.`,\n      hint: 'Use --start and --duration to convert a section, or read the file with pandas or R.',\n    });\n  }\n\n  return { groups, layout, gzip: options.gzip === true, range, columnNames, writeSignals, diagnostics, estimate };\n}\n\n/**\n * Partition channels by sampling rate, largest first.\n *\n * The common case — every channel at one rate — collapses to a single group and a\n * single `signals.csv`, so the honest behaviour costs nothing when there is nothing\n * to be honest about.\n */\nfunction groupByRate(\n  signals: readonly EdfSignal[],\n  columnNames: Map<number, string>,\n  forcedDecimals: number | undefined,\n  gzip: boolean,\n  layout: 'wide' | 'long',\n): RateGroup[] {\n  const byRate = new Map<number, EdfSignal[]>();\n  for (const signal of signals) {\n    // A channel with no samples has no sampling rate to group by, and would\n    // otherwise produce an empty \"0hz\" file. The header parser already warned.\n    if (signal.samplesPerRecord === 0) continue;\n    const bucket = byRate.get(signal.samplingRate);\n    if (bucket) bucket.push(signal);\n    else byRate.set(signal.samplingRate, [signal]);\n  }\n\n  const rates = [...byRate.keys()].sort((a, b) => b - a);\n  // The long layout writes one table whatever the rates are, so every group names it.\n  const single = rates.length === 1 || layout === 'long';\n\n  /*\n    Two distinct rates can produce the same slug, because the slug rounds to six decimal\n    places. Rates come from samplesPerRecord / recordDuration and every channel shares the\n    record duration, so the closest two rates can be is 1 / recordDuration — which drops\n    below 1e-6 once a record is longer than about eleven days. Absurd, but the header\n    permits it, and the failure was silent and destructive: both groups opened a write\n    stream on the same path, so the file ended up holding interleaved rows from both\n    channels under a header naming only one of them.\n\n    Distinct rates therefore get distinct files, always. The suffix is only ever reached by\n    a collision, so ordinary recordings keep the names they have always had.\n\n    Naming from the whole set of rates at once removes most of those collisions before the\n    suffix has to. Rounding each rate on its own gave 1e-6 Hz and 1.25e-6 Hz the same slug,\n    and the numbering below then produced signals_0_000001hz.csv and signals_0_000001hz_2.csv\n    — two files that no longer overwrite each other, but of which only one is named for the\n    rate it holds. The suffix stays as the backstop for anything this still cannot separate.\n  */\n  const suffix = gzip ? '.csv.gz' : '.csv';\n  const slugs = formatRates(rates).map((text) => `${text.replace('.', '_')}hz`);\n  const used = new Set<string>();\n  const uniqueName = (index: number): string => {\n    const base = `signals_${slugs[index]}`;\n    let name = `${base}${suffix}`;\n    for (let n = 2; used.has(name); n++) name = `${base}_${n}${suffix}`;\n    used.add(name);\n    return name;\n  };\n\n  return rates.map((rate, index) => {\n    const members = byRate.get(rate) ?? [];\n    const first = members[0];\n    return {\n      rate,\n      samplesPerRecord: first ? first.samplesPerRecord : 0,\n      fileName: single ? `signals${suffix}` : uniqueName(index),\n      /*\n        In the long layout every rate shares a `time_s` column, so they share its precision:\n        the finest any of them needs. Writing 100 Hz at three places and 256 Hz at eight in\n        the same column would make the column's meaning depend on the row.\n      */\n      timeDecimals: layout === 'long' ? Math.max(...rates.map(timeDecimals)) : timeDecimals(rate),\n      channels: members.map((signal) => ({\n        signal,\n        column: columnNames.get(signal.index) ?? `signal_${signal.index}`,\n        decimals: forcedDecimals ?? decimalsForSignal(signal),\n      })),\n    };\n  });\n}\n\n/**\n * A file's diagnostics with the header's mixed-rate warning removed.\n *\n * `buildPlan` raises that warning for the channels actually being converted, so keeping both\n * would either duplicate it or contradict it. The header parser's copy stays where it is, for\n * callers reading a header without planning a conversion.\n */\nexport function withoutFileRateWarning(diagnostics: readonly Diagnostic[]): Diagnostic[] {\n  return diagnostics.filter((d) => d.code !== 'MIXED_SAMPLING_RATES');\n}\n\n/** `256hz`, `12_5hz` — safe in a filename on every platform. */\nexport function rateSlug(rate: number): string {\n  return `${formatRate(rate).replace('.', '_')}hz`;\n}\n\n/** Characters a fixed-decimal number of this magnitude occupies, sign included. */\nfunction widthOf(magnitude: number, decimals: number, signed = false): number {\n  const size = Math.abs(magnitude);\n  const sign = signed ? 1 : 0;\n  const fraction = decimals > 0 ? 1 + decimals : 0;\n\n  /*\n    Cells are written with toFixed, which rounds. Taking the integer digits from the floor of\n    the bound therefore under-counted whenever rounding carried into a new digit: a channel\n    bounded at 9999.999 and written to zero decimals produces \"10000\", five characters where\n    the floor of 9999.999 suggests four. Every cell on such a channel was a byte short, and\n    `--info` reported 127 KB for a file that came out 131 KB.\n\n    Measuring the bound as rendered removes that. toFixed switches to exponential notation\n    past 1e21, so the arithmetic form still covers magnitudes beyond it.\n  */\n  if (!Number.isFinite(size)) return sign + 1 + fraction;\n  if (size < 1e21) return sign + size.toFixed(Math.min(decimals, 100)).length;\n  return sign + (Math.floor(Math.log10(size)) + 1) + fraction;\n}\n\n/**\n * Raised when the conversion had signal tables to fill and put no data rows in any of them.\n *\n * A window can land where there are no samples without being past the end of the recording:\n * between the last sample and the nominal end of the last record, or — on a discontinuous\n * file — inside a gap. `--start 2 --end 10` on a recording whose records sit at 0s, 1s and\n * 10s asks for eight seconds that contain no data at all.\n *\n * What came out was a signals.csv holding its header and nothing else, exit 0, no warning,\n * and `--strict` passing. The closing summary does say \"signals.csv 0 rows\" and --json\n * carries `rows: 0`, so it was not quite invisible — but a header-only file is exactly what\n * a successful extraction of an empty range looks like, and everywhere else that a request\n * produces nothing this tool says so: a --channels term matching nothing is an error, and\n * --annotations-only on a file with no events raises NO_ANNOTATIONS. A warning rather than\n * an error because a batch of five hundred recordings should not stop for one whose gap\n * happens to line up with the window; --strict turns it into a failure for those who want\n * that.\n */\nfunction emptyWindow(range: ResolvedRange, recordCount: number): Diagnostic {\n  const asked = !range.isWholeRecording;\n  return {\n    code: 'EMPTY_WINDOW',\n    severity: 'warning',\n    message: asked\n      ? `No samples fall inside the requested window (${fixed(range.startSeconds, 3)}s to ` +\n        `${fixed(range.endSeconds, 3)}s), so the signal files hold their headers and no data.`\n      : `This recording's ${recordCount} data records carry no samples in range, so the ` +\n        `signal files hold their headers and no data.`,\n    hint: asked\n      ? /*\n          Which of the two it was, rather than the second one always.\n\n          A recording does not have to start at zero: its first record's timekeeping TAL is\n          what it is timed from, so a file whose records begin at 1000s is asked for with\n          `--start 1000`. `--start 0 --end 1` on that file was told \"The window is inside the\n          recording but lands where there is no data — past the last sample, or inside a gap\n          in a discontinuous file\", when the window sits entirely before the recording and\n          neither offered explanation applies to it. A start at or past the *end* is already\n          an error, so the window being outside can only mean it is before the beginning.\n        */\n        range.endSeconds <= range.recordingStartSeconds\n        ? `This recording starts at ${fixed(range.recordingStartSeconds, 3)}s, so the whole ` +\n          'window sits before it. --start and --end are read on the recording\\'s own clock, ' +\n          'which --info prints as \"Timed from\".'\n        : 'The window is inside the recording but lands where there is no data — past the ' +\n          'last sample, or inside a gap in a discontinuous file. Run with --info to see where ' +\n          'the records actually sit.'\n      : 'Run with --info to see what the header declares.',\n  };\n}\n\nfunction estimateOutput(\n  groups: readonly RateGroup[],\n  range: ResolvedRange,\n  recordDuration: number,\n  recordStarts: Float64Array | null | undefined,\n  bom: boolean,\n  layout: 'wide' | 'long',\n): OutputEstimate {\n  let rows = 0;\n  let bytes = 0;\n  let exceeds = false;\n  // One table in the long layout, so the row limit applies to the sum rather than the\n  // largest group, and the header and mark are counted once rather than once per group.\n  let longRows = 0;\n\n  for (const group of groups) {\n    let groupRows = 0;\n    for (let record = range.startRecord; record < range.endRecord; record++) {\n      const recordStart = recordStarts\n        ? (recordStarts[record] ?? record * recordDuration)\n        : record * recordDuration;\n      groupRows += countSamplesInRange({\n        recordStart,\n        rate: group.rate,\n        samplesPerRecord: group.samplesPerRecord,\n        startSeconds: range.startSeconds,\n        endSeconds: range.endSeconds,\n      });\n    }\n    if (layout === 'long') {\n      // A row per sample per channel rather than a row per sample time.\n      const groupCells = groupRows * group.channels.length;\n      rows += groupCells;\n      longRows += groupCells;\n      /*\n        `time_s,channel,value`: the time, the channel name as it will be escaped into the\n        cell, and the widest the value can print. Same over-counting rule as the wide\n        layout — the declared physical range bounds a cell, and most samples sit under it.\n      */\n      const timeWidth = widthOf(range.endSeconds, group.timeDecimals);\n      for (const channel of group.channels) {\n        const valueWidth = widthOf(\n          Math.max(Math.abs(channel.signal.physicalMin), Math.abs(channel.signal.physicalMax)),\n          channel.decimals,\n          channel.signal.physicalMin < 0 || channel.signal.physicalMax < 0,\n        );\n        const nameWidth = Buffer.byteLength(escapeCsvField(channel.column));\n        // Two commas and the newline.\n        bytes += groupRows * (timeWidth + nameWidth + valueWidth + 3);\n      }\n      continue;\n    }\n\n    rows += groupRows;\n    if (groupRows + 1 > SPREADSHEET_ROW_LIMIT) exceeds = true;\n\n    /*\n      Width per cell, from the channel's own calibration rather than a flat allowance.\n\n      The old `decimals + 6` budgeted six characters for the sign, integer part and decimal\n      point on every channel, whatever it actually held. That over-counted a millivolt\n      channel spanning ±5 by four characters a cell and ran 30-55% high across the fixture\n      set — on a number people use to decide whether a conversion is worth starting.\n\n      The channel's declared physical range is what bounds a cell, so that bound is what is\n      used. Most samples sit below it, so this still reads high, which is the direction a\n      size estimate should err in.\n\n      One case is outside the bound rather than under it: nothing obliges a recording to keep\n      its samples inside the digital range it declares, and one that does not maps outside the\n      physical range too. Such a file can convert larger than the estimate. Clamping the data\n      to make the estimate true is not a trade worth making — the samples are what they are.\n    */\n    const timeWidth = widthOf(range.endSeconds, group.timeDecimals);\n    const cellWidth = group.channels.reduce(\n      (sum, c) =>\n        sum +\n        widthOf(\n          Math.max(Math.abs(c.signal.physicalMin), Math.abs(c.signal.physicalMax)),\n          c.decimals,\n          c.signal.physicalMin < 0 || c.signal.physicalMax < 0,\n        ),\n      0,\n    );\n    // One comma per channel, plus the newline.\n    bytes += groupRows * (timeWidth + cellWidth + group.channels.length + 1);\n    /*\n      The header row, measured as it will be written rather than as the labels are stored.\n\n      A column name is quoted when it contains a comma, a quote, a newline or a leading or\n      trailing space, and every quote inside it is doubled. Counting the raw label under-counted\n      that row: three channels labelled `a,b,c,d,e`, `x\"y` and `plain` write a 32-byte header\n      and were budgeted 27. EDF labels are free text, so commas in them are ordinary — a montage\n      written as `EEG Fpz-Cz, ref` is exactly the kind of thing this is for.\n\n      csvRow is the function that writes it, so it is the function that measures it. Nothing\n      else is in a position to stay correct when the quoting rules change.\n    */\n    bytes += Buffer.byteLength(csvRow(['time_s', ...group.channels.map((c) => c.column)])) + 1;\n    // Three bytes per file under --bom. Small, but the estimate promises never to read\n    // under what gets written, and a one-row conversion is small enough for it to matter.\n    if (bom) bytes += BOM_BYTES;\n  }\n\n  if (layout === 'long' && groups.length > 0) {\n    if (longRows + 1 > SPREADSHEET_ROW_LIMIT) exceeds = true;\n    bytes += Buffer.byteLength(csvRow(['time_s', 'channel', 'value'])) + 1;\n    if (bom) bytes += BOM_BYTES;\n  }\n\n  return { rows, bytes, exceedsSpreadsheetLimit: exceeds };\n}\n"]}

package/dist/edf/header.js

@@ -318,12 +318,28 @@ export function parseHeader(buf, fileSize) {
                     hint: 'It is described in channels.csv but left out of the converted data.',
                 });
             }
-            if (!Number.isFinite(physicalMax - physicalMin)) {
+            /*
+              Too large to represent, and too small — the second was silent.
+      
+              The gain is the span divided by the digital range, and a span of 2e-320 over 65,535
+              codes is 3e-325: below the smallest subnormal double, so it underflows to +0. The
+              scaler's flat-range branch then handed every code the same physical value, and a
+              channel of 65,536 distinct readings became one repeated number with nothing raised at
+              all. One power of ten away, at 1e-319, the same file raises VALUE_RESOLUTION.
+      
+              Both are the same fact about the header — the span cannot be turned into a mapping —
+              so both get this code, and both leave the cells empty rather than filling them with a
+              value the header cannot justify.
+            */
+            const span = physicalMax - physicalMin;
+            const underflowed = span !== 0 && span / (digitalMax - digitalMin) === 0;
+            if (!Number.isFinite(span) || underflowed) {
                 diagnostics.push({
                     code: 'UNUSABLE_PHYSICAL_RANGE',
                     severity: 'warning',
                     message: `Signal ${i} ("${label}") declares a physical range from ${physicalMin} to ` +
-                        `${physicalMax}, whose span is too large to represent, so its values cannot be scaled.`,
+                        `${physicalMax}, whose span is too ${underflowed ? 'small' : 'large'} to ` +
+                        `represent, so its values cannot be scaled.`,
                     hint: 'Its cells are left empty rather than filled with a value the header cannot justify.',
                 });
             }

package/dist/edf/header.js.map

@@ -1 +1 @@
-{"version":3,"file":"header.js","sourceRoot":"","sources":["../../src/edf/header.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6BG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAEvC,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,YAAY,EAAE,MAAM,YAAY,CAAC;AAE1C,gEAAgE;AAChE,MAAM,CAAC,MAAM,iBAAiB,GAAG,iBAAiB,CAAC;AACnD,uDAAuD;AACvD,MAAM,CAAC,MAAM,qBAAqB,GAAG,iBAAiB,CAAC;AAEvD,MAAM,CAAC,MAAM,kBAAkB,GAAG,GAAG,CAAC;AACtC,MAAM,CAAC,MAAM,mBAAmB,GAAG,GAAG,CAAC;AA2EvC;;;;;;GAMG;AACH,SAAS,kBAAkB,CAAC,SAAiB;IAC3C,MAAM,IAAI,GAAG,SAAS,CAAC,WAAW,CAAC,CAAC,CAAW,CAAC;IAChD,OAAO,IAAI,IAAI,IAAI,IAAI,CAAC,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,CAAC,CAAC;AACxD,CAAC;AAED,MAAM,GAAG,GAAG,CAAC,GAAe,EAAE,KAAa,EAAE,GAAW,EAAU,EAAE,CAClE,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;AAExC,8FAA8F;AAC9F,MAAM,SAAS,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC,CAAC;AAEzF;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,eAAe,CAAC,KAAiB;IAC/C,MAAM,KAAK,GAAG,MAAM,CAAC,oBAAoB,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACpE,OAAO,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC;AAC7D,CAAC;AAED,4FAA4F;AAC5F,SAAS,oBAAoB,CAAC,GAAW;IACvC,MAAM,IAAI,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;IAC5B,8EAA8E;IAC9E,IAAI,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC;QAC9C,OAAO,EAAE,IAAI,EAAE,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC;IAC1D,CAAC;IACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC;AACnC,CAAC;AAED,SAAS,gBAAgB,CACvB,GAAW,EACX,KAAa,EACb,EAAE,OAAO,GAAG,KAAK,EAAE,QAAQ,KAA2D,EAAE;IAExF,MAAM,UAAU,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC;IAC7C,MAAM,IAAI,GAAG,UAAU,CAAC,IAAI,CAAC;IAC7B,IAAI,UAAU,CAAC,QAAQ,IAAI,QAAQ;QAAE,QAAQ,CAAC,KAAK,GAAG,IAAI,CAAC;IAC3D,IAAI,IAAI,KAAK,EAAE,EAAE,CAAC;QAChB,MAAM,IAAI,QAAQ,CAAC,kBAAkB,EAAE,iBAAiB,KAAK,aAAa,CAAC,CAAC;IAC9E,CAAC;IACD,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,iBAAiB,KAAK,4BAA4B,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,IAAI,EAC1E,qEAAqE,CACtE,CAAC;IACJ,CAAC;IACD,IAAI,OAAO,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;QACpC,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,iBAAiB,KAAK,mCAAmC,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,IAAI,CAClF,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;GAGG;AACH,SAAS,oBAAoB,CAAC,OAAe,EAAE,OAAe;IAC5D,MAAM,CAAC,GAAG,sCAAsC,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;IAC1E,MAAM,CAAC,GAAG,sCAAsC,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;IAC1E,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;QAAE,OAAO,IAAI,CAAC;IAE1B,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAExB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;QAAE,OAAO,IAAI,CAAC;IAEzF,MAAM,IAAI,GAAG,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;IAC9C,MAAM,IAAI,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5E,6CAA6C;IAC7C,IAAI,IAAI,CAAC,WAAW,EAAE,KAAK,EAAE,GAAG,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE;QAAE,OAAO,IAAI,CAAC;IAC3E,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CAAC,GAAe,EAAE,QAAgB;IAC3D,MAAM,WAAW,GAAiB,EAAE,CAAC;IACrC,MAAM,QAAQ,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAElC,IAAI,GAAG,CAAC,MAAM,GAAG,kBAAkB,EAAE,CAAC;QACpC,MAAM,IAAI,QAAQ,CAChB,gBAAgB,EAChB,WAAW,QAAQ,8CAA8C,kBAAkB,GAAG,CACvF,CAAC;IACJ,CAAC;IAED,6EAA6E;IAC7E,8EAA8E;IAC9E,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,SAAS,CAAC;IAC9D,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAE9D,MAAM,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAC7C,MAAM,WAAW,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAChD,MAAM,YAAY,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,YAAY,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,WAAW,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,wBAAwB,EAAE;QAC/E,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,QAAQ,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IAC9C,MAAM,mBAAmB,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,wBAAwB,EAAE;QACvF,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,cAAc,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,2BAA2B,EAAE;QACrF,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,WAAW,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,mBAAmB,EAAE;QAC1E,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IAEH,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,QAAQ,CAChB,sBAAsB,EACtB,mBAAmB,WAAW,gCAAgC,CAC/D,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,CAAC,cAAc,GAAG,CAAC,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,QAAQ,CAChB,yBAAyB,EACzB,6CAA6C,cAAc,gCAAgC,CAC5F,CAAC;IACJ,CAAC;IAED,MAAM,mBAAmB,GAAG,kBAAkB,GAAG,WAAW,GAAG,mBAAmB,CAAC;IACnF,IAAI,GAAG,CAAC,MAAM,GAAG,mBAAmB,EAAE,CAAC;QACrC,MAAM,IAAI,QAAQ,CAChB,gBAAgB;QAChB;;;;;;;UAOE;QACF,iBAAiB,WAAW,2BAA2B,mBAAmB,gBAAgB;YACxF,CAAC,QAAQ,GAAG,mBAAmB;gBAC7B,CAAC,CAAC,wBAAwB,QAAQ,SAAS;gBAC3C,CAAC,CAAC,YAAY,GAAG,CAAC,MAAM,yCAAyC,CAAC,CACvE,CAAC;IACJ,CAAC;IACD,IAAI,WAAW,KAAK,mBAAmB,EAAE,CAAC;QACxC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,uBAAuB;YAC7B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,qBAAqB,WAAW,eAAe,WAAW,mBAAmB;gBAC7E,GAAG,mBAAmB,yDAAyD;SAClF,CAAC,CAAC;IACL,CAAC;IAED,+EAA+E;IAC/E,MAAM,IAAI,GAAG,kBAAkB,CAAC;IAChC,MAAM,SAAS,GAAG,CAAC,WAAmB,EAAE,KAAa,EAAE,CAAS,EAAU,EAAE,CAC1E,GAAG,CAAC,GAAG,EAAE,IAAI,GAAG,WAAW,GAAG,WAAW,GAAG,CAAC,GAAG,KAAK,EAAE,KAAK,CAAC,CAAC;IAEhE,kFAAkF;IAClF,oEAAoE;IACpE,MAAM,aAAa,GAAG,uBAAuB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7D,MAAM,UAAU,GACd,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC;IAE/E,MAAM,OAAO,GAAgB,EAAE,CAAC;IAChC,IAAI,kBAAkB,GAAG,CAAC,CAAC;IAC3B,MAAM,cAAc,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,UAAU,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC/C,MAAM,WAAW,GAAa,EAAE,CAAC;IAEjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,KAAK,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,UAAU,GAAG,SAAS,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACnD,MAAM,iBAAiB,GAAG,SAAS,CAAC,SAAS,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACzD,MAAM,WAAW,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,4BAA4B,CAAC,GAAG,EAAE;YAC3F,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,WAAW,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,4BAA4B,CAAC,GAAG,EAAE;YAC3F,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,UAAU,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,2BAA2B,CAAC,GAAG,EAAE;YACzF,OAAO,EAAE,IAAI;YACb,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,UAAU,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,2BAA2B,CAAC,GAAG,EAAE;YACzF,OAAO,EAAE,IAAI;YACb,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,YAAY,GAAG,SAAS,CAAC,SAAS,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,gBAAgB,GAAG,gBAAgB,CACvC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EACpB,8BAA8B,CAAC,GAAG,EAClC,EAAE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,CAC5B,CAAC;QACF,MAAM,WAAW,GAAG,SAAS,CAAC,SAAS,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAErD,IAAI,gBAAgB,GAAG,CAAC,EAAE,CAAC;YACzB,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,UAAU,CAAC,MAAM,KAAK,eAAe,gBAAgB,sBAAsB,CAC5E,CAAC;QACJ,CAAC;QAED,MAAM,aAAa,GAAG,KAAK,KAAK,iBAAiB,IAAI,KAAK,KAAK,qBAAqB,CAAC;QAErF,OAAO,CAAC,IAAI,CAAC;YACX,KAAK,EAAE,CAAC;YACR,KAAK;YACL,UAAU;YACV,iBAAiB;YACjB,WAAW;YACX,WAAW;YACX,UAAU;YACV,UAAU;YACV,YAAY;YACZ,gBAAgB;YAChB,QAAQ,EAAE,WAAW;YACrB,aAAa;YACb,YAAY,EAAE,gBAAgB,GAAG,cAAc;YAC/C,kBAAkB;SACnB,CAAC,CAAC;QACH,kBAAkB,IAAI,gBAAgB,GAAG,cAAc,CAAC;QAExD,IAAI,CAAC,aAAa,EAAE,CAAC;YACnB;;;;;;;;;;;;;;cAcE;YACF;;;;;;;;;;cAUE;YACF,MAAM,OAAO,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,MAAM,CAAC,kBAAkB,CAAC,CAAC;YACtD,MAAM,MAAM,GAAG,CAAC,GAAG,iBAAiB,CAAC,CAAC,MAAM,CAAC,kBAAkB,CAAC,CAAC;YACjE,MAAM,OAAO,GAAG,CAAC,GAAG,OAAO,EAAE,GAAG,MAAM,CAAC,CAAC;YACxC,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACvB,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,OAAO,CAAC,CAAC;qBAChC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAO,CAAC,CAAC,WAAW,CAAC,CAAC,CAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC;qBAC9E,IAAI,CAAC,IAAI,CAAC,CAAC;gBACd,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;gBAC/C,MAAM,IAAI,GAAG,OAAO,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;gBACrD,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,wBAAwB,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,eAAe,CAAC;gBACxG,kFAAkF;gBAClF,0EAA0E;gBAC1E,MAAM,KAAK,GAAG,IAAI;oBAChB,CAAC,CAAC,2EAA2E;oBAC7E,CAAC,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC;wBAClB,CAAC,CAAC,0CAA0C;wBAC5C,CAAC,CAAC,gDAAgD,CAAC;gBACvD,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,oBAAoB;oBAC1B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,qBAAqB,MAAM,GAAG;wBACtE,IAAI,KAAK,MAAM,KAAK,kCAAkC;oBACxD,IAAI,EACF,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC;wBACjB,CAAC,CAAC,qDAAqD,CAAC,yBAAyB;4BAC/E,kCAAkC;wBACpC,CAAC,CAAC,iDAAiD,KAAK,sBAAsB,CAAC;wBACjF,2DAA2D;iBAC9D,CAAC,CAAC;YACL,CAAC;YAED,IAAI,KAAK,KAAK,EAAE,EAAE,CAAC;gBACjB,qFAAqF;gBACrF,sFAAsF;gBACtF,2DAA2D;gBAC3D,WAAW,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACtB,CAAC;iBAAM,CAAC;gBACN,0EAA0E;gBAC1E,sEAAsE;gBACtE,MAAM,IAAI,GAAG,UAAU,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;gBACnC,IAAI,IAAI;oBAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;oBAClB,UAAU,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAClC,CAAC;YAED,IAAI,gBAAgB,KAAK,CAAC,EAAE,CAAC;gBAC3B,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,YAAY;oBAClB,QAAQ,EAAE,SAAS;oBACnB,OAAO,EAAE,UAAU,CAAC,MAAM,KAAK,mDAAmD;oBAClF,IAAI,EAAE,qEAAqE;iBAC5E,CAAC,CAAC;YACL,CAAC;YAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,GAAG,WAAW,CAAC,EAAE,CAAC;gBAChD,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,yBAAyB;oBAC/B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,qCAAqC,WAAW,MAAM;wBAC5E,GAAG,WAAW,yEAAyE;oBACzF,IAAI,EAAE,qFAAqF;iBAC5F,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,UAAU,KAAK,UAAU,EAAE,CAAC;gBACrC,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,0BAA0B;oBAChC,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,kDAAkD;wBACxE,IAAI,UAAU,oCAAoC;oBACpD,IAAI,EAAE,qFAAqF;iBAC5F,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,WAAW,KAAK,WAAW,EAAE,CAAC;gBACvC,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,2BAA2B;oBACjC,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,oDAAoD;wBAC1E,IAAI,WAAW,gDAAgD;iBAClE,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC;gBACvE;;;;;;;;;;kBAUE;gBACF,MAAM,QAAQ,GACZ,WAAW,GAAG,WAAW;oBACvB,CAAC,CAAC,oBAAoB,WAAW,2BAA2B,WAAW,EAAE;oBACzE,CAAC,CAAC,mBAAmB,UAAU,0BAA0B,UAAU,EAAE,CAAC;gBAC1E,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,yBAAyB;oBAC/B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EAAE,UAAU,CAAC,MAAM,KAAK,eAAe,QAAQ,+BAA+B;oBACrF,IAAI,EAAE,+EAA+E;iBACtF,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IAED;;;;;;;;;;;;;;;;;;MAkBE;IACF,KAAK,MAAM,KAAK,IAAI,WAAW,EAAE,CAAC;QAChC,MAAM,KAAK,GAAG,UAAU,CAAC,GAAG,CAAC,UAAU,KAAK,EAAE,CAAC,CAAC;QAChD,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,aAAa;YACnB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,KAAK,KAAK,SAAS;gBACjB,CAAC,CAAC,UAAU,KAAK,4CAA4C,KAAK,IAAI;gBACtE,CAAC,CAAC,UAAU,KAAK,+CAA+C,KAAK,YAAY;oBAC/E,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,SAAS,IAAI,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,WAAW;oBACpF,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,mCAAmC;oBAC9E,uCAAuC;SAC9C,CAAC,CAAC;IACL,CAAC;IAED,KAAK,MAAM,CAAC,KAAK,EAAE,OAAO,CAAC,IAAI,UAAU,EAAE,CAAC;QAC1C,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC;YAAE,SAAS;QACjC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,iBAAiB;YACvB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,GAAG,OAAO,CAAC,MAAM,6BAA6B,KAAK,gBAAgB,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI;YAClG,IAAI,EAAE,iFAAiF;SACxF,CAAC,CAAC;IACL,CAAC;IAED,MAAM,WAAW,GAAG,kBAAkB,CAAC;IACvC,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,QAAQ,CAChB,YAAY,EACZ,2FAA2F,CAC5F,CAAC;IACJ,CAAC;IAED,IAAI,QAAQ,CAAC,KAAK,EAAE,CAAC;QACnB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,eAAe;YACrB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,uFAAuF;YAChG,IAAI,EAAE,0EAA0E;SACjF,CAAC,CAAC;IACL,CAAC;IAED,MAAM,SAAS,GAAG,QAAQ,GAAG,mBAAmB,CAAC;IACjD,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;QAClB,MAAM,IAAI,QAAQ,CAAC,gBAAgB,EAAE,sCAAsC,CAAC,CAAC;IAC/E,CAAC;IACD,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,WAAW,CAAC,CAAC;IACxD,MAAM,aAAa,GAAG,SAAS,GAAG,WAAW,GAAG,WAAW,CAAC;IAE5D,IAAI,WAAW,KAAK,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,QAAQ,CAChB,iBAAiB,EACjB,yDAAyD,EACzD,qEAAqE,CACtE,CAAC;IACJ,CAAC;IAED,IAAI,mBAAmB,KAAK,CAAC,CAAC,EAAE,CAAC;QAC/B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,yFAAyF;gBACzF,+CAA+C,OAAO,CAAC,WAAW,EAAE,QAAQ,CAAC,8BAA8B;SAC9G,CAAC,CAAC;IACL,CAAC;SAAM,IAAI,mBAAmB,KAAK,WAAW,EAAE,CAAC;QAC/C,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,uBAAuB;YAC7B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,uBAAuB,mBAAmB,sCAAsC;gBAChF,GAAG,WAAW,oBAAoB,OAAO,CAAC,WAAW,EAAE,QAAQ,CAAC,SAAS,WAAW,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,WAAW;YACtH,IAAI,EACF,mBAAmB,GAAG,WAAW;gBAC/B,CAAC,CAAC,mFAAmF;gBACrF,CAAC,CAAC,4CAA4C;SACnD,CAAC,CAAC;IACL,CAAC;IAED,IAAI,aAAa,GAAG,CAAC,EAAE,CAAC;QACtB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,gBAAgB;YACtB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,GAAG,OAAO,CAAC,aAAa,EAAE,MAAM,CAAC,wCAAwC,aAAa,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,WAAW;SAClI,CAAC,CAAC;IACL,CAAC;IAED,MAAM,SAAS,GAAG,UAAU,KAAK,IAAI,CAAC;IACtC,IAAI,UAAU,KAAK,OAAO,EAAE,CAAC;QAC3B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,eAAe;YACrB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,4BAA4B,KAAK,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,oCAAoC;gBACzF,yBAAyB;YAC3B,IAAI,EAAE,yFAAyF;SAChG,CAAC,CAAC;IACL,CAAC;IAED,MAAM,WAAW,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC;IAC5D,IAAI,WAAW,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC7B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,oBAAoB;YAC1B,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,sEAAsE;SAChF,CAAC,CAAC;IACL,CAAC;IAED,uFAAuF;IACvF,wFAAwF;IACxF,sFAAsF;IACtF,wEAAwE;IACxE,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,gBAAgB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC;IACpG,IAAI,KAAK,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACnB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gBAAgB,KAAK,CAAC,IAAI,4BAA4B;gBACtD,IAAI,MAAM,CAAC,WAAW,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI;YACrF,IAAI,EAAE,mEAAmE;SAC1E,CAAC,CAAC;IACL,CAAC;IAED,OAAO;QACL,MAAM,EAAE;YACN,OAAO;YACP,SAAS;YACT,WAAW;YACX,YAAY;YACZ,YAAY;YACZ,aAAa,EAAE,oBAAoB,CAAC,YAAY,EAAE,YAAY,CAAC;YAC/D,WAAW,EAAE,mBAAmB;YAChC,mBAAmB,EAAE,WAAW;YAChC,QAAQ;YACR,SAAS;YACT,KAAK;YACL,UAAU;YACV,mBAAmB;YACnB,cAAc;YACd,WAAW;YACX,OAAO;YACP,cAAc;YACd,WAAW;SACZ;QACD,WAAW;QACX,aAAa;QACb,WAAW;KACZ,CAAC;AACJ,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,eAAe,CAAC,IAAiB;IAC/C,IAAI,CAAC,IAAI;QAAE,OAAO,IAAI,CAAC;IACvB,OAAO,IAAI,CAAC,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;AACzC,CAAC;AAED,oDAAoD;AACpD,MAAM,UAAU,cAAc,CAAC,MAAiB;IAC9C,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC;IAC1C,IAAI,CAAC,MAAM,CAAC,SAAS;QAAE,OAAO,IAAI,CAAC;IACnC,OAAO,GAAG,IAAI,MAAM,MAAM,CAAC,UAAU,KAAK,OAAO,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC;AACxF,CAAC;AAED,qEAAqE;AACrE,MAAM,UAAU,UAAU,CAAC,EAAU;IACnC,IAAI,MAAM,CAAC,SAAS,CAAC,EAAE,CAAC;QAAE,OAAO,MAAM,CAAC,EAAE,CAAC,CAAC;IAC5C,MAAM,OAAO,GAAG,MAAM,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,qFAAqF;IACrF,sFAAsF;IACtF,sFAAsF;IACtF,8EAA8E;IAC9E,IAAI,OAAO,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IAC9C,OAAO,MAAM,CAAC,OAAO,CAAC,CAAC;AACzB,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,WAAW,CAAC,KAAwB;IAClD,MAAM,OAAO,GAAG,KAAK,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;IACrC,OAAO,IAAI,GAAG,CAAC,OAAO,CAAC,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;AACtF,CAAC","sourcesContent":["/**\n * EDF / EDF+ header parsing.\n *\n * Layout (all fields are ASCII, left-justified, space-padded):\n *\n *   fixed header, 256 bytes\n *     0    8   version ('0', or 255 + 'BIOSEMI' for BDF)\n *     8   80   patient identification\n *    88   80   recording identification\n *   168    8   start date  dd.mm.yy\n *   176    8   start time  hh.mm.ss\n *   184    8   number of bytes in the header record\n *   192   44   reserved  ('EDF+C' / 'EDF+D' live here)\n *   236    8   number of data records (-1 if unknown)\n *   244    8   duration of a data record, in seconds (may be fractional)\n *   252    4   number of signals (ns)\n *\n *   signal header, ns * 256 bytes, stored FIELD-major rather than signal-major:\n *   all ns labels, then all ns transducer types, and so on.\n *     ns * 16   label\n *     ns * 80   transducer type\n *     ns *  8   physical dimension\n *     ns *  8   physical minimum\n *     ns *  8   physical maximum\n *     ns *  8   digital minimum\n *     ns *  8   digital maximum\n *     ns * 80   prefiltering\n *     ns *  8   number of samples in each data record\n *     ns * 32   reserved\n */\n\nimport { EdfError } from './errors.js';\nimport type { Diagnostic } from './errors.js';\nimport { counted, listed } from '../format/list.js';\nimport { decodeLatin1 } from './bytes.js';\n\n/** Label the EDF+ spec reserves for the annotations channel. */\nexport const ANNOTATIONS_LABEL = 'EDF Annotations';\n/** BDF+ uses its own spelling for the same channel. */\nexport const BDF_ANNOTATIONS_LABEL = 'BDF Annotations';\n\nexport const FIXED_HEADER_BYTES = 256;\nexport const SIGNAL_HEADER_BYTES = 256;\n\nexport interface EdfSignal {\n  /** Position in the file, 0-based. Stable even when labels collide. */\n  index: number;\n  label: string;\n  transducer: string;\n  physicalDimension: string;\n  physicalMin: number;\n  physicalMax: number;\n  digitalMin: number;\n  digitalMax: number;\n  prefiltering: string;\n  samplesPerRecord: number;\n  reserved: string;\n  /** True for the EDF+ 'EDF Annotations' channel, which carries text, not signal. */\n  isAnnotations: boolean;\n  /** samplesPerRecord / recordDuration, in Hz. */\n  samplingRate: number;\n  /** Byte offset of this signal's samples within one data record. */\n  byteOffsetInRecord: number;\n}\n\nexport interface EdfHeader {\n  version: string;\n  patientId: string;\n  recordingId: string;\n  /** Raw 'dd.mm.yy' as written in the file. */\n  startDateRaw: string;\n  /** Raw 'hh.mm.ss' as written in the file. */\n  startTimeRaw: string;\n  /** Resolved start instant, or null when the file's date/time fields are unusable. */\n  startDateTime: Date | null;\n  /**\n   * Header size in bytes, computed from the signal count rather than read from the field.\n   *\n   * Every data record offset is derived from this, so it has to be the one the layout\n   * actually uses: 256 for the fixed header plus 256 per signal. A writer that fills the\n   * field in carelessly is common enough to have its own warning, HEADER_BYTES_MISMATCH,\n   * and trusting the field over the arithmetic would put every sample at the wrong offset.\n   */\n  headerBytes: number;\n  /**\n   * What the header's own length field says, which need not be the above.\n   *\n   * Exposed for the same reason `declaredRecordCount` is: the two disagreeing is a fact\n   * about the file, and a caller checking how a recording was written should be able to see\n   * what it claimed rather than only what was believed.\n   */\n  declaredHeaderBytes: number;\n  reserved: string;\n  isEdfPlus: boolean;\n  /** True for BioSemi BDF/BDF+ files, whose samples are 3 bytes rather than 2. */\n  isBdf: boolean;\n  /** 'EDF+C' continuous, 'EDF+D' discontinuous, or null for plain EDF. */\n  continuity: 'EDF+C' | 'EDF+D' | null;\n  /** As declared in the header. -1 means \"unknown\", which the spec permits. */\n  declaredRecordCount: number;\n  recordDuration: number;\n  signalCount: number;\n  signals: EdfSignal[];\n  bytesPerSample: number;\n  recordBytes: number;\n}\n\n/** Everything derived by combining the header with the file's real size. */\nexport interface EdfHeaderInfo {\n  header: EdfHeader;\n  /** Record count implied by the actual file size — the one we trust for reading. */\n  recordCount: number;\n  /** Bytes after the last complete data record. */\n  trailingBytes: number;\n  diagnostics: Diagnostic[];\n}\n\n/**\n * A byte the terminal treats as an instruction rather than as text.\n *\n * C0 and C1, plus DEL. Tab is included deliberately: it is harmless to a terminal but it\n * makes a CSV column name that cannot be typed or matched reliably, which is the other half\n * of what this warning is for.\n */\nfunction isControlCharacter(character: string): boolean {\n  const code = character.codePointAt(0) as number;\n  return code <= 0x1f || (code >= 0x7f && code <= 0x9f);\n}\n\nconst dec = (buf: Uint8Array, start: number, len: number): string =>\n  decodeLatin1(buf, start, start + len);\n\n/** EDF fields are space-padded; trailing NULs also occur in files written by sloppy tools. */\nconst trimField = (s: string): string => s.replace(/[\\0\\s]+$/u, '').replace(/^\\s+/u, '');\n\n/**\n * How many signals the fixed header says there are, read exactly as `parseHeader` will.\n *\n * `EdfFile.open` needs this before it can know how much header to read, and it used to work\n * it out with its own `Number(...)` — which was NUL-tolerant but not comma-tolerant, unlike\n * every other numeric field here. A header written with a comma decimal separator, which\n * COMMA_DECIMAL exists to accept and which the documentation lists this field among, was\n * therefore never given its signal headers at all, and the file died on a message that\n * contradicted itself: \"needs a 768-byte header, but the file is only 848 bytes\".\n *\n * Sharing the parse is what keeps the two from disagreeing again about which files are\n * readable. Null means \"not a usable count\", and the caller reads no further header — the\n * real error then comes from `parseHeader`, which is the one place that decides.\n */\nexport function peekSignalCount(fixed: Uint8Array): number | null {\n  const count = Number(normaliseNumberField(dec(fixed, 252, 4)).text);\n  return Number.isInteger(count) && count > 0 ? count : null;\n}\n\n/** A numeric header field, trimmed and with a comma decimal separator turned into a dot. */\nfunction normaliseNumberField(raw: string): { text: string; sawComma: boolean } {\n  const text = trimField(raw);\n  // Some writers emit a comma decimal separator despite the spec requiring '.'.\n  if (text.includes(',') && !text.includes('.')) {\n    return { text: text.replace(',', '.'), sawComma: true };\n  }\n  return { text, sawComma: false };\n}\n\nfunction parseNumberField(\n  raw: string,\n  field: string,\n  { integer = false, sawComma }: { integer?: boolean; sawComma?: { value: boolean } } = {},\n): number {\n  const normalised = normaliseNumberField(raw);\n  const text = normalised.text;\n  if (normalised.sawComma && sawComma) sawComma.value = true;\n  if (text === '') {\n    throw new EdfError('BAD_HEADER_FIELD', `Header field \"${field}\" is empty.`);\n  }\n  const n = Number(text);\n  if (!Number.isFinite(n)) {\n    throw new EdfError(\n      'BAD_HEADER_FIELD',\n      `Header field \"${field}\" is not a number (found ${JSON.stringify(text)}).`,\n      'The file may be truncated, byte-shifted, or not an EDF file at all.',\n    );\n  }\n  if (integer && !Number.isInteger(n)) {\n    throw new EdfError(\n      'BAD_HEADER_FIELD',\n      `Header field \"${field}\" must be a whole number (found ${JSON.stringify(text)}).`,\n    );\n  }\n  return n;\n}\n\n/**\n * EDF stores a two-digit year. The spec pins the century: 85-99 mean 1985-1999\n * and 00-84 mean 2000-2084. Files outside 1985-2084 cannot express their date.\n */\nfunction resolveStartDateTime(dateRaw: string, timeRaw: string): Date | null {\n  const d = /^(\\d{2})[.\\-/](\\d{2})[.\\-/](\\d{2})$/u.exec(trimField(dateRaw));\n  const t = /^(\\d{2})[.:\\-](\\d{2})[.:\\-](\\d{2})$/u.exec(trimField(timeRaw));\n  if (!d || !t) return null;\n\n  const dd = Number(d[1]);\n  const mm = Number(d[2]);\n  const yy = Number(d[3]);\n  const hh = Number(t[1]);\n  const mi = Number(t[2]);\n  const ss = Number(t[3]);\n\n  if (mm < 1 || mm > 12 || dd < 1 || dd > 31 || hh > 23 || mi > 59 || ss > 60) return null;\n\n  const year = yy >= 85 ? 1900 + yy : 2000 + yy;\n  const date = new Date(Date.UTC(year, mm - 1, dd, hh, mi, Math.min(ss, 59)));\n  // Reject dates that rolled over, e.g. 31.02.\n  if (date.getUTCMonth() !== mm - 1 || date.getUTCDate() !== dd) return null;\n  return date;\n}\n\n/**\n * Parse the fixed 256-byte header plus the per-signal header block.\n *\n * @param buf   At least FIXED_HEADER_BYTES + ns * SIGNAL_HEADER_BYTES bytes.\n * @param fileSize Total size of the file on disk, used to derive the real record count.\n */\nexport function parseHeader(buf: Uint8Array, fileSize: number): EdfHeaderInfo {\n  const diagnostics: Diagnostic[] = [];\n  const sawComma = { value: false };\n\n  if (buf.length < FIXED_HEADER_BYTES) {\n    throw new EdfError(\n      'FILE_TOO_SMALL',\n      `File is ${fileSize} bytes; an EDF header alone needs at least ${FIXED_HEADER_BYTES}.`,\n    );\n  }\n\n  // BDF (BioSemi) marks itself with byte 255 followed by 'BIOSEMI', and stores\n  // 3-byte samples instead of 2. Everything else about the layout is identical.\n  const isBdf = buf[0] === 0xff && dec(buf, 1, 7) === 'BIOSEMI';\n  const version = isBdf ? 'BIOSEMI' : trimField(dec(buf, 0, 8));\n\n  const patientId = trimField(dec(buf, 8, 80));\n  const recordingId = trimField(dec(buf, 88, 80));\n  const startDateRaw = trimField(dec(buf, 168, 8));\n  const startTimeRaw = trimField(dec(buf, 176, 8));\n  const headerBytes = parseNumberField(dec(buf, 184, 8), 'number of header bytes', {\n    integer: true,\n    sawComma,\n  });\n  const reserved = trimField(dec(buf, 192, 44));\n  const declaredRecordCount = parseNumberField(dec(buf, 236, 8), 'number of data records', {\n    integer: true,\n    sawComma,\n  });\n  const recordDuration = parseNumberField(dec(buf, 244, 8), 'duration of a data record', {\n    sawComma,\n  });\n  const signalCount = parseNumberField(dec(buf, 252, 4), 'number of signals', {\n    integer: true,\n    sawComma,\n  });\n\n  if (signalCount <= 0) {\n    throw new EdfError(\n      'INVALID_SIGNAL_COUNT',\n      `Header declares ${signalCount} signals; expected at least 1.`,\n    );\n  }\n  if (!(recordDuration > 0)) {\n    throw new EdfError(\n      'INVALID_RECORD_DURATION',\n      `Header declares a data record duration of ${recordDuration}s; expected a positive number.`,\n    );\n  }\n\n  const expectedHeaderBytes = FIXED_HEADER_BYTES + signalCount * SIGNAL_HEADER_BYTES;\n  if (buf.length < expectedHeaderBytes) {\n    throw new EdfError(\n      'FILE_TOO_SMALL',\n      /*\n        Which of the two is actually short.\n\n        The file size was quoted either way, so a caller that had read too little — the\n        signal count parsed one way here and another way there — produced arithmetic that\n        refuted itself: \"needs a 768-byte header, but the file is only 848 bytes\". A reader\n        following that looks for a truncation that is not there.\n      */\n      `File declares ${signalCount} signals, which needs a ${expectedHeaderBytes}-byte header, ` +\n        (fileSize < expectedHeaderBytes\n          ? `but the file is only ${fileSize} bytes.`\n          : `but only ${buf.length} bytes of it were handed to the parser.`),\n    );\n  }\n  if (headerBytes !== expectedHeaderBytes) {\n    diagnostics.push({\n      code: 'HEADER_BYTES_MISMATCH',\n      severity: 'warning',\n      message:\n        `Header says it is ${headerBytes} bytes, but ${signalCount} signals require ` +\n        `${expectedHeaderBytes} bytes. Using the value computed from the signal count.`,\n    });\n  }\n\n  // Signal headers are field-major: all labels, then all transducers, and so on.\n  const base = FIXED_HEADER_BYTES;\n  const readField = (offsetUnits: number, width: number, i: number): string =>\n    dec(buf, base + offsetUnits * signalCount + i * width, width);\n\n  // EDF+ writes 'EDF+C'/'EDF+D' here; BDF+ writes 'BDF+C'/'BDF+D'. The two mean the\n  // same thing, so both are normalised to a single continuity marker.\n  const continuityTag = /^(?:EDF|BDF)\\+([CD])/u.exec(reserved);\n  const continuity: 'EDF+C' | 'EDF+D' | null =\n    continuityTag === null ? null : continuityTag[1] === 'D' ? 'EDF+D' : 'EDF+C';\n\n  const signals: EdfSignal[] = [];\n  let byteOffsetInRecord = 0;\n  const bytesPerSample = isBdf ? 3 : 2;\n  const seenLabels = new Map<string, number[]>();\n  const emptyLabels: number[] = [];\n\n  for (let i = 0; i < signalCount; i++) {\n    const label = trimField(readField(0, 16, i));\n    const transducer = trimField(readField(16, 80, i));\n    const physicalDimension = trimField(readField(96, 8, i));\n    const physicalMin = parseNumberField(readField(104, 8, i), `physical minimum (signal ${i})`, {\n      sawComma,\n    });\n    const physicalMax = parseNumberField(readField(112, 8, i), `physical maximum (signal ${i})`, {\n      sawComma,\n    });\n    const digitalMin = parseNumberField(readField(120, 8, i), `digital minimum (signal ${i})`, {\n      integer: true,\n      sawComma,\n    });\n    const digitalMax = parseNumberField(readField(128, 8, i), `digital maximum (signal ${i})`, {\n      integer: true,\n      sawComma,\n    });\n    const prefiltering = trimField(readField(136, 80, i));\n    const samplesPerRecord = parseNumberField(\n      readField(216, 8, i),\n      `samples per record (signal ${i})`,\n      { integer: true, sawComma },\n    );\n    const sigReserved = trimField(readField(224, 32, i));\n\n    if (samplesPerRecord < 0) {\n      throw new EdfError(\n        'BAD_HEADER_FIELD',\n        `Signal ${i} (\"${label}\") declares ${samplesPerRecord} samples per record.`,\n      );\n    }\n\n    const isAnnotations = label === ANNOTATIONS_LABEL || label === BDF_ANNOTATIONS_LABEL;\n\n    signals.push({\n      index: i,\n      label,\n      transducer,\n      physicalDimension,\n      physicalMin,\n      physicalMax,\n      digitalMin,\n      digitalMax,\n      prefiltering,\n      samplesPerRecord,\n      reserved: sigReserved,\n      isAnnotations,\n      samplingRate: samplesPerRecord / recordDuration,\n      byteOffsetInRecord,\n    });\n    byteOffsetInRecord += samplesPerRecord * bytesPerSample;\n\n    if (!isAnnotations) {\n      /*\n        A label is free text out of the file, and it becomes a column name in signals.csv.\n\n        `--info` has escaped control bytes since it was written, because an ANSI escape in a\n        header can drive the reader's terminal — `\\x1b[2J` clears the screen. The CSV had no\n        such protection and needed none for correctness: quoting makes any byte safe for a\n        parser, and this still passes the label through exactly as the file gives it, because\n        losing what the header says is not an improvement.\n\n        What was missing is the sentence saying so. A recording whose channel is labelled\n        `\\x1b[2Jgone` converted with no warning at all, and `cat signals.csv` then cleared\n        the terminal — while a script referencing that column by name carried an invisible\n        control character in it. NONPRINTABLE_LABEL has been declared and documented as\n        reserved since 0.1; this is it doing its job.\n      */\n      /*\n        Which of the two fields carries them, because the consequences are not the same.\n\n        The message said \"label or unit\", and then said the bytes \"will appear in the CSV\n        column name\" and that \"the name cannot be typed\" — both of which are about the label.\n        A channel labelled plainly `ECG` in a unit of `u\\x07V` got all of it: its column is\n        `ECG`, `--channels ECG` selects it and exits 0, and the byte is in channels.csv's\n        `unit` cell, which the warning never mentioned. Three sentences, none of them true of\n        the file that raised it, on a warning whose whole purpose is to say where an invisible\n        byte went.\n      */\n      const inLabel = [...label].filter(isControlCharacter);\n      const inUnit = [...physicalDimension].filter(isControlCharacter);\n      const control = [...inLabel, ...inUnit];\n      if (control.length > 0) {\n        const shown = [...new Set(control)]\n          .map((c) => `\\\\x${(c.codePointAt(0) as number).toString(16).padStart(2, '0')}`)\n          .join(', ');\n        const plural = control.length === 1 ? '' : 's';\n        const both = inLabel.length > 0 && inUnit.length > 0;\n        const field = both ? 'label and unit contain' : inLabel.length > 0 ? 'label contains' : 'unit contains';\n        // Where they land, which is the question the reader has. A label becomes a column\n        // name in signals.csv; a unit is a cell of channels.csv and nothing else.\n        const lands = both\n          ? 'which will appear in the CSV column name and in channels.csv\\'s unit cell'\n          : inLabel.length > 0\n            ? 'which will appear in the CSV column name'\n            : 'which will appear in channels.csv\\'s unit cell';\n        diagnostics.push({\n          code: 'NONPRINTABLE_LABEL',\n          severity: 'warning',\n          message:\n            `Signal ${i}'s ${field} ${control.length} control character${plural} ` +\n            `(${shown}), ${lands} exactly as the header has them.`,\n          hint:\n            (inLabel.length > 0\n              ? `Address the channel by position with --channels \"#${i}\" rather than by name, ` +\n                'since the name cannot be typed. '\n              : `The column name is unaffected, so --channels \"${label}\" still selects it. `) +\n            'Printing the CSV to a terminal may do more than print it.',\n        });\n      }\n\n      if (label === '') {\n        // Collected, not reported here: what this channel's column ends up called depends on\n        // whether some later channel is literally labelled `signal_<i>`, and inside this loop\n        // the later channels do not exist yet. See the pass below.\n        emptyLabels.push(i);\n      } else {\n        // Collected rather than reported here: a label repeated five times should\n        // produce one warning naming all five, not four near-identical pairs.\n        const seen = seenLabels.get(label);\n        if (seen) seen.push(i);\n        else seenLabels.set(label, [i]);\n      }\n\n      if (samplesPerRecord === 0) {\n        diagnostics.push({\n          code: 'NO_SAMPLES',\n          severity: 'warning',\n          message: `Signal ${i} (\"${label}\") carries no samples at all (0 per data record).`,\n          hint: 'It is described in channels.csv but left out of the converted data.',\n        });\n      }\n\n      if (!Number.isFinite(physicalMax - physicalMin)) {\n        diagnostics.push({\n          code: 'UNUSABLE_PHYSICAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") declares a physical range from ${physicalMin} to ` +\n            `${physicalMax}, whose span is too large to represent, so its values cannot be scaled.`,\n          hint: 'Its cells are left empty rather than filled with a value the header cannot justify.',\n        });\n      } else if (digitalMax === digitalMin) {\n        diagnostics.push({\n          code: 'DEGENERATE_DIGITAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") has digital minimum equal to digital maximum ` +\n            `(${digitalMin}), so its values cannot be scaled.`,\n          hint: 'Its cells are left empty rather than filled with a value the header cannot justify.',\n        });\n      } else if (physicalMax === physicalMin) {\n        diagnostics.push({\n          code: 'DEGENERATE_PHYSICAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") has physical minimum equal to physical maximum ` +\n            `(${physicalMin}), so every sample converts to the same value.`,\n        });\n      } else if ((physicalMax - physicalMin) * (digitalMax - digitalMin) < 0) {\n        /*\n          Polarity is inverted when the gain is negative, and the gain is\n          (physicalMax - physicalMin) / (digitalMax - digitalMin) — so it is the sign of the\n          two spans together that matters, not the physical pair alone.\n\n          Testing only `physicalMax < physicalMin` was wrong in both directions. A file with\n          its DIGITAL bounds reversed is just as inverted and drew no warning at all, handing\n          back sign-flipped EEG with nothing to indicate it. A file with BOTH pairs reversed\n          has a positive gain and is not inverted, yet was warned about — a message that was\n          simply untrue of that recording.\n        */\n        const reversed =\n          physicalMax < physicalMin\n            ? `physical minimum ${physicalMin} above physical maximum ${physicalMax}`\n            : `digital minimum ${digitalMin} above digital maximum ${digitalMax}`;\n        diagnostics.push({\n          code: 'INVERTED_PHYSICAL_RANGE',\n          severity: 'warning',\n          message: `Signal ${i} (\"${label}\") declares ${reversed}, which inverts its polarity.`,\n          hint: 'The values are converted exactly as the header specifies, inversion included.',\n        });\n      }\n    }\n  }\n\n  /*\n    What an unlabelled channel is actually called, which the message used to guess.\n\n    A channel with no label takes `signal_<index>` — unless another channel is literally\n    labelled that, which EDF permits, since labels are free text and nothing enforces anything\n    about them. Then both collide and both are suffixed. The warning said \"It will appear as\n    \"signal_0\"\" while the file's header read `time_s,signal_0_ch0,signal_0_ch1`: the one\n    sentence the run printed named a column that exists in neither signals.csv nor\n    channels.csv.\n\n    The other half was silent. The channel that genuinely carries the label `signal_0` lost\n    its own column name to a collision with a synthesised one, and nothing said so —\n    DUPLICATE_LABEL did not fire, because the two labels are not the same label. Both halves\n    are one sentence here, because they are one event.\n\n    No specific suffixed name is quoted. The suffix rule has a second pass for names that are\n    still shared afterwards, and a message that hard-coded `_ch<index>` would be guessing again\n    in exactly the way this is fixing.\n  */\n  for (const index of emptyLabels) {\n    const taken = seenLabels.get(`signal_${index}`);\n    diagnostics.push({\n      code: 'EMPTY_LABEL',\n      severity: 'warning',\n      message:\n        taken === undefined\n          ? `Signal ${index} has no label. It will appear as \"signal_${index}\".`\n          : `Signal ${index} has no label, so it takes the name \"signal_${index}\" — which ` +\n            `${taken.length === 1 ? 'signal' : 'signals'} ${listed(taken.map(String))} already ` +\n            `${taken.length === 1 ? 'carries' : 'carry'} as a label, so both columns are ` +\n            `suffixed with their position instead.`,\n    });\n  }\n\n  for (const [label, indices] of seenLabels) {\n    if (indices.length < 2) continue;\n    diagnostics.push({\n      code: 'DUPLICATE_LABEL',\n      severity: 'warning',\n      message: `${indices.length} signals share the label \"${label}\" (positions ${indices.join(', ')}).`,\n      hint: 'Their columns are suffixed with the signal number so they stay distinguishable.',\n    });\n  }\n\n  const recordBytes = byteOffsetInRecord;\n  if (recordBytes <= 0) {\n    throw new EdfError(\n      'NO_SAMPLES',\n      'No signal in this file carries any samples (every channel declares 0 samples per record).',\n    );\n  }\n\n  if (sawComma.value) {\n    diagnostics.push({\n      code: 'COMMA_DECIMAL',\n      severity: 'warning',\n      message: 'Some header numbers use a comma decimal separator, which the EDF spec does not allow.',\n      hint: 'They were read as decimal points. Check the values in the channel table.',\n    });\n  }\n\n  const dataBytes = fileSize - expectedHeaderBytes;\n  if (dataBytes < 0) {\n    throw new EdfError('FILE_TOO_SMALL', `File is smaller than its own header.`);\n  }\n  const recordCount = Math.floor(dataBytes / recordBytes);\n  const trailingBytes = dataBytes - recordCount * recordBytes;\n\n  if (recordCount === 0) {\n    throw new EdfError(\n      'NO_DATA_RECORDS',\n      'The file contains a header but no complete data record.',\n      'The recording was probably interrupted before any data was written.',\n    );\n  }\n\n  if (declaredRecordCount === -1) {\n    diagnostics.push({\n      code: 'RECORD_COUNT_UNKNOWN',\n      severity: 'warning',\n      message:\n        `The header does not say how many data records the file has (-1), which the spec allows ` +\n        `for recordings still in progress. Using the ${counted(recordCount, 'record')} the file actually contains.`,\n    });\n  } else if (declaredRecordCount !== recordCount) {\n    diagnostics.push({\n      code: 'RECORD_COUNT_MISMATCH',\n      severity: 'warning',\n      message:\n        `The header declares ${declaredRecordCount} data records but the file contains ` +\n        `${recordCount}. Converting the ${counted(recordCount, 'record')} that ${recordCount === 1 ? 'is' : 'are'} present.`,\n      hint:\n        declaredRecordCount > recordCount\n          ? 'The recording looks truncated. It may have been cut short or copied incompletely.'\n          : 'The file is longer than its header claims.',\n    });\n  }\n\n  if (trailingBytes > 0) {\n    diagnostics.push({\n      code: 'TRAILING_BYTES',\n      severity: 'warning',\n      message: `${counted(trailingBytes, 'byte')} after the last complete data record ${trailingBytes === 1 ? 'was' : 'were'} ignored.`,\n    });\n  }\n\n  const isEdfPlus = continuity !== null;\n  if (continuity === 'EDF+D') {\n    diagnostics.push({\n      code: 'DISCONTINUOUS',\n      severity: 'warning',\n      message:\n        `This is a discontinuous (${isBdf ? 'BDF+D' : 'EDF+D'}) recording: its data records are ` +\n        `not contiguous in time.`,\n      hint: 'Each row carries its true recording time, so gaps stay visible instead of being closed.',\n    });\n  }\n\n  const dataSignals = signals.filter((s) => !s.isAnnotations);\n  if (dataSignals.length === 0) {\n    diagnostics.push({\n      code: 'NO_SIGNAL_CHANNELS',\n      severity: 'warning',\n      message: 'This file has no signal channels; it contains only EDF+ annotations.',\n    });\n  }\n\n  // A channel declaring zero samples per record has no sampling rate to speak of — it is\n  // reported separately as NO_SAMPLES and no file is written for it. Counting its nominal\n  // 0 Hz as a rate made a single-rate recording warn that it used \"2 different sampling\n  // rates (4 Hz, 0 Hz)\" and claim it was splitting output it never split.\n  const rates = new Set(dataSignals.filter((s) => s.samplesPerRecord > 0).map((s) => s.samplingRate));\n  if (rates.size > 1) {\n    diagnostics.push({\n      code: 'MIXED_SAMPLING_RATES',\n      severity: 'warning',\n      message:\n        `Channels use ${rates.size} different sampling rates ` +\n        `(${listed(formatRates([...rates].sort((a, b) => b - a)).map((r) => `${r} Hz`))}).`,\n      hint: 'They are written to one file per rate so no channel is resampled.',\n    });\n  }\n\n  return {\n    header: {\n      version,\n      patientId,\n      recordingId,\n      startDateRaw,\n      startTimeRaw,\n      startDateTime: resolveStartDateTime(startDateRaw, startTimeRaw),\n      headerBytes: expectedHeaderBytes,\n      declaredHeaderBytes: headerBytes,\n      reserved,\n      isEdfPlus,\n      isBdf,\n      continuity,\n      declaredRecordCount,\n      recordDuration,\n      signalCount,\n      signals,\n      bytesPerSample,\n      recordBytes,\n    },\n    recordCount,\n    trailingBytes,\n    diagnostics,\n  };\n}\n\n/**\n * The recording start as a zone-less wall clock, \"YYYY-MM-DDTHH:MM:SS\".\n *\n * EDF stores the start time as local wall-clock digits with no timezone anywhere in\n * the format. `startDateTime` is built with Date.UTC purely so those digits survive a\n * round trip unshifted, which makes it a carrier for the wall clock rather than a\n * real instant. Serialising it with `toISOString()` would append a Z and assert UTC,\n * and any reader converting to local time would then shift the recording by their own\n * offset: 13:43:04 in the file becomes 08:43:04 in New York. The Z is omitted because\n * the file genuinely does not say which zone it meant.\n */\nexport function formatWallClock(date: Date | null): string | null {\n  if (!date) return null;\n  return date.toISOString().slice(0, 19);\n}\n\n/** \"EDF\", \"EDF+ (EDF+D)\", \"BDF\", \"BDF+ (EDF+C)\". */\nexport function describeFormat(header: EdfHeader): string {\n  const base = header.isBdf ? 'BDF' : 'EDF';\n  if (!header.isEdfPlus) return base;\n  return `${base}+ (${header.continuity === 'EDF+D' ? 'discontinuous' : 'continuous'})`;\n}\n\n/** Render a sampling rate without trailing noise: 256, 0.5, 12.5. */\nexport function formatRate(hz: number): string {\n  if (Number.isInteger(hz)) return String(hz);\n  const rounded = Number(hz.toFixed(6));\n  // A rate below 5e-7 rounds away to \"0\", which reads as \"this channel has no sampling\n  // rate\" and made the mixed-rate warning contradict itself: it announced two different\n  // rates and then printed both as \"0 Hz\". Exponent form keeps a real rate legible, and\n  // keeps distinct rates distinct in the channel table and in output filenames.\n  if (rounded === 0) return hz.toExponential(3);\n  return String(rounded);\n}\n\n/**\n * Renders a group of rates so that rates which differ read as differing.\n *\n * `formatRate` rounds to six decimals, which is what keeps an ordinary rate free of\n * float noise — 30 samples in a 0.1-second record is 299.99999999999994 as a double,\n * and belongs on screen as 300. Two rates separated by less than that round to one\n * string, so a file carrying 1e-6 Hz and 1.25e-6 Hz warned that it used \"2 different\n * sampling rates (0.000001 Hz, 0.000001 Hz)\" and named both files the same thing.\n *\n * That is the contradiction the exponent fallback above already removes for rates that\n * round away to zero; this is the same one a step further out. On a collision every rate\n * in the group switches to its shortest exact form, which is unique for distinct values,\n * rather than only the pair that collided — one column in one notation reads better than\n * two.\n */\nexport function formatRates(rates: readonly number[]): string[] {\n  const rounded = rates.map(formatRate);\n  const distinct = new Set(rates).size;\n  return new Set(rounded).size === distinct ? rounded : rates.map((hz) => String(hz));\n}\n"]}
+{"version":3,"file":"header.js","sourceRoot":"","sources":["../../src/edf/header.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6BG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAEvC,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,YAAY,EAAE,MAAM,YAAY,CAAC;AAE1C,gEAAgE;AAChE,MAAM,CAAC,MAAM,iBAAiB,GAAG,iBAAiB,CAAC;AACnD,uDAAuD;AACvD,MAAM,CAAC,MAAM,qBAAqB,GAAG,iBAAiB,CAAC;AAEvD,MAAM,CAAC,MAAM,kBAAkB,GAAG,GAAG,CAAC;AACtC,MAAM,CAAC,MAAM,mBAAmB,GAAG,GAAG,CAAC;AA2EvC;;;;;;GAMG;AACH,SAAS,kBAAkB,CAAC,SAAiB;IAC3C,MAAM,IAAI,GAAG,SAAS,CAAC,WAAW,CAAC,CAAC,CAAW,CAAC;IAChD,OAAO,IAAI,IAAI,IAAI,IAAI,CAAC,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,CAAC,CAAC;AACxD,CAAC;AAED,MAAM,GAAG,GAAG,CAAC,GAAe,EAAE,KAAa,EAAE,GAAW,EAAU,EAAE,CAClE,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;AAExC,8FAA8F;AAC9F,MAAM,SAAS,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC,CAAC;AAEzF;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,eAAe,CAAC,KAAiB;IAC/C,MAAM,KAAK,GAAG,MAAM,CAAC,oBAAoB,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACpE,OAAO,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC;AAC7D,CAAC;AAED,4FAA4F;AAC5F,SAAS,oBAAoB,CAAC,GAAW;IACvC,MAAM,IAAI,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;IAC5B,8EAA8E;IAC9E,IAAI,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC;QAC9C,OAAO,EAAE,IAAI,EAAE,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC;IAC1D,CAAC;IACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC;AACnC,CAAC;AAED,SAAS,gBAAgB,CACvB,GAAW,EACX,KAAa,EACb,EAAE,OAAO,GAAG,KAAK,EAAE,QAAQ,KAA2D,EAAE;IAExF,MAAM,UAAU,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC;IAC7C,MAAM,IAAI,GAAG,UAAU,CAAC,IAAI,CAAC;IAC7B,IAAI,UAAU,CAAC,QAAQ,IAAI,QAAQ;QAAE,QAAQ,CAAC,KAAK,GAAG,IAAI,CAAC;IAC3D,IAAI,IAAI,KAAK,EAAE,EAAE,CAAC;QAChB,MAAM,IAAI,QAAQ,CAAC,kBAAkB,EAAE,iBAAiB,KAAK,aAAa,CAAC,CAAC;IAC9E,CAAC;IACD,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,iBAAiB,KAAK,4BAA4B,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,IAAI,EAC1E,qEAAqE,CACtE,CAAC;IACJ,CAAC;IACD,IAAI,OAAO,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC;QACpC,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,iBAAiB,KAAK,mCAAmC,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,IAAI,CAClF,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;GAGG;AACH,SAAS,oBAAoB,CAAC,OAAe,EAAE,OAAe;IAC5D,MAAM,CAAC,GAAG,sCAAsC,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;IAC1E,MAAM,CAAC,GAAG,sCAAsC,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;IAC1E,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;QAAE,OAAO,IAAI,CAAC;IAE1B,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAExB,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;QAAE,OAAO,IAAI,CAAC;IAEzF,MAAM,IAAI,GAAG,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;IAC9C,MAAM,IAAI,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5E,6CAA6C;IAC7C,IAAI,IAAI,CAAC,WAAW,EAAE,KAAK,EAAE,GAAG,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,KAAK,EAAE;QAAE,OAAO,IAAI,CAAC;IAC3E,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CAAC,GAAe,EAAE,QAAgB;IAC3D,MAAM,WAAW,GAAiB,EAAE,CAAC;IACrC,MAAM,QAAQ,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAElC,IAAI,GAAG,CAAC,MAAM,GAAG,kBAAkB,EAAE,CAAC;QACpC,MAAM,IAAI,QAAQ,CAChB,gBAAgB,EAChB,WAAW,QAAQ,8CAA8C,kBAAkB,GAAG,CACvF,CAAC;IACJ,CAAC;IAED,6EAA6E;IAC7E,8EAA8E;IAC9E,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,SAAS,CAAC;IAC9D,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAE9D,MAAM,SAAS,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAC7C,MAAM,WAAW,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAChD,MAAM,YAAY,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,YAAY,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjD,MAAM,WAAW,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,wBAAwB,EAAE;QAC/E,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,QAAQ,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IAC9C,MAAM,mBAAmB,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,wBAAwB,EAAE;QACvF,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,cAAc,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,2BAA2B,EAAE;QACrF,QAAQ;KACT,CAAC,CAAC;IACH,MAAM,WAAW,GAAG,gBAAgB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,mBAAmB,EAAE;QAC1E,OAAO,EAAE,IAAI;QACb,QAAQ;KACT,CAAC,CAAC;IAEH,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,QAAQ,CAChB,sBAAsB,EACtB,mBAAmB,WAAW,gCAAgC,CAC/D,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,CAAC,cAAc,GAAG,CAAC,CAAC,EAAE,CAAC;QAC1B,MAAM,IAAI,QAAQ,CAChB,yBAAyB,EACzB,6CAA6C,cAAc,gCAAgC,CAC5F,CAAC;IACJ,CAAC;IAED,MAAM,mBAAmB,GAAG,kBAAkB,GAAG,WAAW,GAAG,mBAAmB,CAAC;IACnF,IAAI,GAAG,CAAC,MAAM,GAAG,mBAAmB,EAAE,CAAC;QACrC,MAAM,IAAI,QAAQ,CAChB,gBAAgB;QAChB;;;;;;;UAOE;QACF,iBAAiB,WAAW,2BAA2B,mBAAmB,gBAAgB;YACxF,CAAC,QAAQ,GAAG,mBAAmB;gBAC7B,CAAC,CAAC,wBAAwB,QAAQ,SAAS;gBAC3C,CAAC,CAAC,YAAY,GAAG,CAAC,MAAM,yCAAyC,CAAC,CACvE,CAAC;IACJ,CAAC;IACD,IAAI,WAAW,KAAK,mBAAmB,EAAE,CAAC;QACxC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,uBAAuB;YAC7B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,qBAAqB,WAAW,eAAe,WAAW,mBAAmB;gBAC7E,GAAG,mBAAmB,yDAAyD;SAClF,CAAC,CAAC;IACL,CAAC;IAED,+EAA+E;IAC/E,MAAM,IAAI,GAAG,kBAAkB,CAAC;IAChC,MAAM,SAAS,GAAG,CAAC,WAAmB,EAAE,KAAa,EAAE,CAAS,EAAU,EAAE,CAC1E,GAAG,CAAC,GAAG,EAAE,IAAI,GAAG,WAAW,GAAG,WAAW,GAAG,CAAC,GAAG,KAAK,EAAE,KAAK,CAAC,CAAC;IAEhE,kFAAkF;IAClF,oEAAoE;IACpE,MAAM,aAAa,GAAG,uBAAuB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7D,MAAM,UAAU,GACd,aAAa,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC;IAE/E,MAAM,OAAO,GAAgB,EAAE,CAAC;IAChC,IAAI,kBAAkB,GAAG,CAAC,CAAC;IAC3B,MAAM,cAAc,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,UAAU,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC/C,MAAM,WAAW,GAAa,EAAE,CAAC;IAEjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,KAAK,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,UAAU,GAAG,SAAS,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACnD,MAAM,iBAAiB,GAAG,SAAS,CAAC,SAAS,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACzD,MAAM,WAAW,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,4BAA4B,CAAC,GAAG,EAAE;YAC3F,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,WAAW,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,4BAA4B,CAAC,GAAG,EAAE;YAC3F,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,UAAU,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,2BAA2B,CAAC,GAAG,EAAE;YACzF,OAAO,EAAE,IAAI;YACb,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,UAAU,GAAG,gBAAgB,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,2BAA2B,CAAC,GAAG,EAAE;YACzF,OAAO,EAAE,IAAI;YACb,QAAQ;SACT,CAAC,CAAC;QACH,MAAM,YAAY,GAAG,SAAS,CAAC,SAAS,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,gBAAgB,GAAG,gBAAgB,CACvC,SAAS,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EACpB,8BAA8B,CAAC,GAAG,EAClC,EAAE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,CAC5B,CAAC;QACF,MAAM,WAAW,GAAG,SAAS,CAAC,SAAS,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAErD,IAAI,gBAAgB,GAAG,CAAC,EAAE,CAAC;YACzB,MAAM,IAAI,QAAQ,CAChB,kBAAkB,EAClB,UAAU,CAAC,MAAM,KAAK,eAAe,gBAAgB,sBAAsB,CAC5E,CAAC;QACJ,CAAC;QAED,MAAM,aAAa,GAAG,KAAK,KAAK,iBAAiB,IAAI,KAAK,KAAK,qBAAqB,CAAC;QAErF,OAAO,CAAC,IAAI,CAAC;YACX,KAAK,EAAE,CAAC;YACR,KAAK;YACL,UAAU;YACV,iBAAiB;YACjB,WAAW;YACX,WAAW;YACX,UAAU;YACV,UAAU;YACV,YAAY;YACZ,gBAAgB;YAChB,QAAQ,EAAE,WAAW;YACrB,aAAa;YACb,YAAY,EAAE,gBAAgB,GAAG,cAAc;YAC/C,kBAAkB;SACnB,CAAC,CAAC;QACH,kBAAkB,IAAI,gBAAgB,GAAG,cAAc,CAAC;QAExD,IAAI,CAAC,aAAa,EAAE,CAAC;YACnB;;;;;;;;;;;;;;cAcE;YACF;;;;;;;;;;cAUE;YACF,MAAM,OAAO,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,MAAM,CAAC,kBAAkB,CAAC,CAAC;YACtD,MAAM,MAAM,GAAG,CAAC,GAAG,iBAAiB,CAAC,CAAC,MAAM,CAAC,kBAAkB,CAAC,CAAC;YACjE,MAAM,OAAO,GAAG,CAAC,GAAG,OAAO,EAAE,GAAG,MAAM,CAAC,CAAC;YACxC,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACvB,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,OAAO,CAAC,CAAC;qBAChC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAO,CAAC,CAAC,WAAW,CAAC,CAAC,CAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC;qBAC9E,IAAI,CAAC,IAAI,CAAC,CAAC;gBACd,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;gBAC/C,MAAM,IAAI,GAAG,OAAO,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;gBACrD,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,wBAAwB,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,eAAe,CAAC;gBACxG,kFAAkF;gBAClF,0EAA0E;gBAC1E,MAAM,KAAK,GAAG,IAAI;oBAChB,CAAC,CAAC,2EAA2E;oBAC7E,CAAC,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC;wBAClB,CAAC,CAAC,0CAA0C;wBAC5C,CAAC,CAAC,gDAAgD,CAAC;gBACvD,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,oBAAoB;oBAC1B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,IAAI,OAAO,CAAC,MAAM,qBAAqB,MAAM,GAAG;wBACtE,IAAI,KAAK,MAAM,KAAK,kCAAkC;oBACxD,IAAI,EACF,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC;wBACjB,CAAC,CAAC,qDAAqD,CAAC,yBAAyB;4BAC/E,kCAAkC;wBACpC,CAAC,CAAC,iDAAiD,KAAK,sBAAsB,CAAC;wBACjF,2DAA2D;iBAC9D,CAAC,CAAC;YACL,CAAC;YAED,IAAI,KAAK,KAAK,EAAE,EAAE,CAAC;gBACjB,qFAAqF;gBACrF,sFAAsF;gBACtF,2DAA2D;gBAC3D,WAAW,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACtB,CAAC;iBAAM,CAAC;gBACN,0EAA0E;gBAC1E,sEAAsE;gBACtE,MAAM,IAAI,GAAG,UAAU,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;gBACnC,IAAI,IAAI;oBAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;oBAClB,UAAU,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAClC,CAAC;YAED,IAAI,gBAAgB,KAAK,CAAC,EAAE,CAAC;gBAC3B,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,YAAY;oBAClB,QAAQ,EAAE,SAAS;oBACnB,OAAO,EAAE,UAAU,CAAC,MAAM,KAAK,mDAAmD;oBAClF,IAAI,EAAE,qEAAqE;iBAC5E,CAAC,CAAC;YACL,CAAC;YAED;;;;;;;;;;;;cAYE;YACF,MAAM,IAAI,GAAG,WAAW,GAAG,WAAW,CAAC;YACvC,MAAM,WAAW,GAAG,IAAI,KAAK,CAAC,IAAI,IAAI,GAAG,CAAC,UAAU,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC;YACzE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,WAAW,EAAE,CAAC;gBAC1C,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,yBAAyB;oBAC/B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,qCAAqC,WAAW,MAAM;wBAC5E,GAAG,WAAW,uBAAuB,WAAW,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,MAAM;wBAC1E,4CAA4C;oBAC9C,IAAI,EAAE,qFAAqF;iBAC5F,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,UAAU,KAAK,UAAU,EAAE,CAAC;gBACrC,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,0BAA0B;oBAChC,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,kDAAkD;wBACxE,IAAI,UAAU,oCAAoC;oBACpD,IAAI,EAAE,qFAAqF;iBAC5F,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,WAAW,KAAK,WAAW,EAAE,CAAC;gBACvC,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,2BAA2B;oBACjC,QAAQ,EAAE,SAAS;oBACnB,OAAO,EACL,UAAU,CAAC,MAAM,KAAK,oDAAoD;wBAC1E,IAAI,WAAW,gDAAgD;iBAClE,CAAC,CAAC;YACL,CAAC;iBAAM,IAAI,CAAC,WAAW,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC;gBACvE;;;;;;;;;;kBAUE;gBACF,MAAM,QAAQ,GACZ,WAAW,GAAG,WAAW;oBACvB,CAAC,CAAC,oBAAoB,WAAW,2BAA2B,WAAW,EAAE;oBACzE,CAAC,CAAC,mBAAmB,UAAU,0BAA0B,UAAU,EAAE,CAAC;gBAC1E,WAAW,CAAC,IAAI,CAAC;oBACf,IAAI,EAAE,yBAAyB;oBAC/B,QAAQ,EAAE,SAAS;oBACnB,OAAO,EAAE,UAAU,CAAC,MAAM,KAAK,eAAe,QAAQ,+BAA+B;oBACrF,IAAI,EAAE,+EAA+E;iBACtF,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IAED;;;;;;;;;;;;;;;;;;MAkBE;IACF,KAAK,MAAM,KAAK,IAAI,WAAW,EAAE,CAAC;QAChC,MAAM,KAAK,GAAG,UAAU,CAAC,GAAG,CAAC,UAAU,KAAK,EAAE,CAAC,CAAC;QAChD,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,aAAa;YACnB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,KAAK,KAAK,SAAS;gBACjB,CAAC,CAAC,UAAU,KAAK,4CAA4C,KAAK,IAAI;gBACtE,CAAC,CAAC,UAAU,KAAK,+CAA+C,KAAK,YAAY;oBAC/E,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,SAAS,IAAI,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,WAAW;oBACpF,GAAG,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,mCAAmC;oBAC9E,uCAAuC;SAC9C,CAAC,CAAC;IACL,CAAC;IAED,KAAK,MAAM,CAAC,KAAK,EAAE,OAAO,CAAC,IAAI,UAAU,EAAE,CAAC;QAC1C,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC;YAAE,SAAS;QACjC,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,iBAAiB;YACvB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,GAAG,OAAO,CAAC,MAAM,6BAA6B,KAAK,gBAAgB,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI;YAClG,IAAI,EAAE,iFAAiF;SACxF,CAAC,CAAC;IACL,CAAC;IAED,MAAM,WAAW,GAAG,kBAAkB,CAAC;IACvC,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACrB,MAAM,IAAI,QAAQ,CAChB,YAAY,EACZ,2FAA2F,CAC5F,CAAC;IACJ,CAAC;IAED,IAAI,QAAQ,CAAC,KAAK,EAAE,CAAC;QACnB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,eAAe;YACrB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,uFAAuF;YAChG,IAAI,EAAE,0EAA0E;SACjF,CAAC,CAAC;IACL,CAAC;IAED,MAAM,SAAS,GAAG,QAAQ,GAAG,mBAAmB,CAAC;IACjD,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;QAClB,MAAM,IAAI,QAAQ,CAAC,gBAAgB,EAAE,sCAAsC,CAAC,CAAC;IAC/E,CAAC;IACD,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,WAAW,CAAC,CAAC;IACxD,MAAM,aAAa,GAAG,SAAS,GAAG,WAAW,GAAG,WAAW,CAAC;IAE5D,IAAI,WAAW,KAAK,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,QAAQ,CAChB,iBAAiB,EACjB,yDAAyD,EACzD,qEAAqE,CACtE,CAAC;IACJ,CAAC;IAED,IAAI,mBAAmB,KAAK,CAAC,CAAC,EAAE,CAAC;QAC/B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,yFAAyF;gBACzF,+CAA+C,OAAO,CAAC,WAAW,EAAE,QAAQ,CAAC,8BAA8B;SAC9G,CAAC,CAAC;IACL,CAAC;SAAM,IAAI,mBAAmB,KAAK,WAAW,EAAE,CAAC;QAC/C,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,uBAAuB;YAC7B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,uBAAuB,mBAAmB,sCAAsC;gBAChF,GAAG,WAAW,oBAAoB,OAAO,CAAC,WAAW,EAAE,QAAQ,CAAC,SAAS,WAAW,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,WAAW;YACtH,IAAI,EACF,mBAAmB,GAAG,WAAW;gBAC/B,CAAC,CAAC,mFAAmF;gBACrF,CAAC,CAAC,4CAA4C;SACnD,CAAC,CAAC;IACL,CAAC;IAED,IAAI,aAAa,GAAG,CAAC,EAAE,CAAC;QACtB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,gBAAgB;YACtB,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,GAAG,OAAO,CAAC,aAAa,EAAE,MAAM,CAAC,wCAAwC,aAAa,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,WAAW;SAClI,CAAC,CAAC;IACL,CAAC;IAED,MAAM,SAAS,GAAG,UAAU,KAAK,IAAI,CAAC;IACtC,IAAI,UAAU,KAAK,OAAO,EAAE,CAAC;QAC3B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,eAAe;YACrB,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,4BAA4B,KAAK,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,oCAAoC;gBACzF,yBAAyB;YAC3B,IAAI,EAAE,yFAAyF;SAChG,CAAC,CAAC;IACL,CAAC;IAED,MAAM,WAAW,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC;IAC5D,IAAI,WAAW,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC7B,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,oBAAoB;YAC1B,QAAQ,EAAE,SAAS;YACnB,OAAO,EAAE,sEAAsE;SAChF,CAAC,CAAC;IACL,CAAC;IAED,uFAAuF;IACvF,wFAAwF;IACxF,sFAAsF;IACtF,wEAAwE;IACxE,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,gBAAgB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC;IACpG,IAAI,KAAK,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACnB,WAAW,CAAC,IAAI,CAAC;YACf,IAAI,EAAE,sBAAsB;YAC5B,QAAQ,EAAE,SAAS;YACnB,OAAO,EACL,gBAAgB,KAAK,CAAC,IAAI,4BAA4B;gBACtD,IAAI,MAAM,CAAC,WAAW,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI;YACrF,IAAI,EAAE,mEAAmE;SAC1E,CAAC,CAAC;IACL,CAAC;IAED,OAAO;QACL,MAAM,EAAE;YACN,OAAO;YACP,SAAS;YACT,WAAW;YACX,YAAY;YACZ,YAAY;YACZ,aAAa,EAAE,oBAAoB,CAAC,YAAY,EAAE,YAAY,CAAC;YAC/D,WAAW,EAAE,mBAAmB;YAChC,mBAAmB,EAAE,WAAW;YAChC,QAAQ;YACR,SAAS;YACT,KAAK;YACL,UAAU;YACV,mBAAmB;YACnB,cAAc;YACd,WAAW;YACX,OAAO;YACP,cAAc;YACd,WAAW;SACZ;QACD,WAAW;QACX,aAAa;QACb,WAAW;KACZ,CAAC;AACJ,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,eAAe,CAAC,IAAiB;IAC/C,IAAI,CAAC,IAAI;QAAE,OAAO,IAAI,CAAC;IACvB,OAAO,IAAI,CAAC,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;AACzC,CAAC;AAED,oDAAoD;AACpD,MAAM,UAAU,cAAc,CAAC,MAAiB;IAC9C,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC;IAC1C,IAAI,CAAC,MAAM,CAAC,SAAS;QAAE,OAAO,IAAI,CAAC;IACnC,OAAO,GAAG,IAAI,MAAM,MAAM,CAAC,UAAU,KAAK,OAAO,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC;AACxF,CAAC;AAED,qEAAqE;AACrE,MAAM,UAAU,UAAU,CAAC,EAAU;IACnC,IAAI,MAAM,CAAC,SAAS,CAAC,EAAE,CAAC;QAAE,OAAO,MAAM,CAAC,EAAE,CAAC,CAAC;IAC5C,MAAM,OAAO,GAAG,MAAM,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,qFAAqF;IACrF,sFAAsF;IACtF,sFAAsF;IACtF,8EAA8E;IAC9E,IAAI,OAAO,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IAC9C,OAAO,MAAM,CAAC,OAAO,CAAC,CAAC;AACzB,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,WAAW,CAAC,KAAwB;IAClD,MAAM,OAAO,GAAG,KAAK,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;IACrC,OAAO,IAAI,GAAG,CAAC,OAAO,CAAC,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;AACtF,CAAC","sourcesContent":["/**\n * EDF / EDF+ header parsing.\n *\n * Layout (all fields are ASCII, left-justified, space-padded):\n *\n *   fixed header, 256 bytes\n *     0    8   version ('0', or 255 + 'BIOSEMI' for BDF)\n *     8   80   patient identification\n *    88   80   recording identification\n *   168    8   start date  dd.mm.yy\n *   176    8   start time  hh.mm.ss\n *   184    8   number of bytes in the header record\n *   192   44   reserved  ('EDF+C' / 'EDF+D' live here)\n *   236    8   number of data records (-1 if unknown)\n *   244    8   duration of a data record, in seconds (may be fractional)\n *   252    4   number of signals (ns)\n *\n *   signal header, ns * 256 bytes, stored FIELD-major rather than signal-major:\n *   all ns labels, then all ns transducer types, and so on.\n *     ns * 16   label\n *     ns * 80   transducer type\n *     ns *  8   physical dimension\n *     ns *  8   physical minimum\n *     ns *  8   physical maximum\n *     ns *  8   digital minimum\n *     ns *  8   digital maximum\n *     ns * 80   prefiltering\n *     ns *  8   number of samples in each data record\n *     ns * 32   reserved\n */\n\nimport { EdfError } from './errors.js';\nimport type { Diagnostic } from './errors.js';\nimport { counted, listed } from '../format/list.js';\nimport { decodeLatin1 } from './bytes.js';\n\n/** Label the EDF+ spec reserves for the annotations channel. */\nexport const ANNOTATIONS_LABEL = 'EDF Annotations';\n/** BDF+ uses its own spelling for the same channel. */\nexport const BDF_ANNOTATIONS_LABEL = 'BDF Annotations';\n\nexport const FIXED_HEADER_BYTES = 256;\nexport const SIGNAL_HEADER_BYTES = 256;\n\nexport interface EdfSignal {\n  /** Position in the file, 0-based. Stable even when labels collide. */\n  index: number;\n  label: string;\n  transducer: string;\n  physicalDimension: string;\n  physicalMin: number;\n  physicalMax: number;\n  digitalMin: number;\n  digitalMax: number;\n  prefiltering: string;\n  samplesPerRecord: number;\n  reserved: string;\n  /** True for the EDF+ 'EDF Annotations' channel, which carries text, not signal. */\n  isAnnotations: boolean;\n  /** samplesPerRecord / recordDuration, in Hz. */\n  samplingRate: number;\n  /** Byte offset of this signal's samples within one data record. */\n  byteOffsetInRecord: number;\n}\n\nexport interface EdfHeader {\n  version: string;\n  patientId: string;\n  recordingId: string;\n  /** Raw 'dd.mm.yy' as written in the file. */\n  startDateRaw: string;\n  /** Raw 'hh.mm.ss' as written in the file. */\n  startTimeRaw: string;\n  /** Resolved start instant, or null when the file's date/time fields are unusable. */\n  startDateTime: Date | null;\n  /**\n   * Header size in bytes, computed from the signal count rather than read from the field.\n   *\n   * Every data record offset is derived from this, so it has to be the one the layout\n   * actually uses: 256 for the fixed header plus 256 per signal. A writer that fills the\n   * field in carelessly is common enough to have its own warning, HEADER_BYTES_MISMATCH,\n   * and trusting the field over the arithmetic would put every sample at the wrong offset.\n   */\n  headerBytes: number;\n  /**\n   * What the header's own length field says, which need not be the above.\n   *\n   * Exposed for the same reason `declaredRecordCount` is: the two disagreeing is a fact\n   * about the file, and a caller checking how a recording was written should be able to see\n   * what it claimed rather than only what was believed.\n   */\n  declaredHeaderBytes: number;\n  reserved: string;\n  isEdfPlus: boolean;\n  /** True for BioSemi BDF/BDF+ files, whose samples are 3 bytes rather than 2. */\n  isBdf: boolean;\n  /** 'EDF+C' continuous, 'EDF+D' discontinuous, or null for plain EDF. */\n  continuity: 'EDF+C' | 'EDF+D' | null;\n  /** As declared in the header. -1 means \"unknown\", which the spec permits. */\n  declaredRecordCount: number;\n  recordDuration: number;\n  signalCount: number;\n  signals: EdfSignal[];\n  bytesPerSample: number;\n  recordBytes: number;\n}\n\n/** Everything derived by combining the header with the file's real size. */\nexport interface EdfHeaderInfo {\n  header: EdfHeader;\n  /** Record count implied by the actual file size — the one we trust for reading. */\n  recordCount: number;\n  /** Bytes after the last complete data record. */\n  trailingBytes: number;\n  diagnostics: Diagnostic[];\n}\n\n/**\n * A byte the terminal treats as an instruction rather than as text.\n *\n * C0 and C1, plus DEL. Tab is included deliberately: it is harmless to a terminal but it\n * makes a CSV column name that cannot be typed or matched reliably, which is the other half\n * of what this warning is for.\n */\nfunction isControlCharacter(character: string): boolean {\n  const code = character.codePointAt(0) as number;\n  return code <= 0x1f || (code >= 0x7f && code <= 0x9f);\n}\n\nconst dec = (buf: Uint8Array, start: number, len: number): string =>\n  decodeLatin1(buf, start, start + len);\n\n/** EDF fields are space-padded; trailing NULs also occur in files written by sloppy tools. */\nconst trimField = (s: string): string => s.replace(/[\\0\\s]+$/u, '').replace(/^\\s+/u, '');\n\n/**\n * How many signals the fixed header says there are, read exactly as `parseHeader` will.\n *\n * `EdfFile.open` needs this before it can know how much header to read, and it used to work\n * it out with its own `Number(...)` — which was NUL-tolerant but not comma-tolerant, unlike\n * every other numeric field here. A header written with a comma decimal separator, which\n * COMMA_DECIMAL exists to accept and which the documentation lists this field among, was\n * therefore never given its signal headers at all, and the file died on a message that\n * contradicted itself: \"needs a 768-byte header, but the file is only 848 bytes\".\n *\n * Sharing the parse is what keeps the two from disagreeing again about which files are\n * readable. Null means \"not a usable count\", and the caller reads no further header — the\n * real error then comes from `parseHeader`, which is the one place that decides.\n */\nexport function peekSignalCount(fixed: Uint8Array): number | null {\n  const count = Number(normaliseNumberField(dec(fixed, 252, 4)).text);\n  return Number.isInteger(count) && count > 0 ? count : null;\n}\n\n/** A numeric header field, trimmed and with a comma decimal separator turned into a dot. */\nfunction normaliseNumberField(raw: string): { text: string; sawComma: boolean } {\n  const text = trimField(raw);\n  // Some writers emit a comma decimal separator despite the spec requiring '.'.\n  if (text.includes(',') && !text.includes('.')) {\n    return { text: text.replace(',', '.'), sawComma: true };\n  }\n  return { text, sawComma: false };\n}\n\nfunction parseNumberField(\n  raw: string,\n  field: string,\n  { integer = false, sawComma }: { integer?: boolean; sawComma?: { value: boolean } } = {},\n): number {\n  const normalised = normaliseNumberField(raw);\n  const text = normalised.text;\n  if (normalised.sawComma && sawComma) sawComma.value = true;\n  if (text === '') {\n    throw new EdfError('BAD_HEADER_FIELD', `Header field \"${field}\" is empty.`);\n  }\n  const n = Number(text);\n  if (!Number.isFinite(n)) {\n    throw new EdfError(\n      'BAD_HEADER_FIELD',\n      `Header field \"${field}\" is not a number (found ${JSON.stringify(text)}).`,\n      'The file may be truncated, byte-shifted, or not an EDF file at all.',\n    );\n  }\n  if (integer && !Number.isInteger(n)) {\n    throw new EdfError(\n      'BAD_HEADER_FIELD',\n      `Header field \"${field}\" must be a whole number (found ${JSON.stringify(text)}).`,\n    );\n  }\n  return n;\n}\n\n/**\n * EDF stores a two-digit year. The spec pins the century: 85-99 mean 1985-1999\n * and 00-84 mean 2000-2084. Files outside 1985-2084 cannot express their date.\n */\nfunction resolveStartDateTime(dateRaw: string, timeRaw: string): Date | null {\n  const d = /^(\\d{2})[.\\-/](\\d{2})[.\\-/](\\d{2})$/u.exec(trimField(dateRaw));\n  const t = /^(\\d{2})[.:\\-](\\d{2})[.:\\-](\\d{2})$/u.exec(trimField(timeRaw));\n  if (!d || !t) return null;\n\n  const dd = Number(d[1]);\n  const mm = Number(d[2]);\n  const yy = Number(d[3]);\n  const hh = Number(t[1]);\n  const mi = Number(t[2]);\n  const ss = Number(t[3]);\n\n  if (mm < 1 || mm > 12 || dd < 1 || dd > 31 || hh > 23 || mi > 59 || ss > 60) return null;\n\n  const year = yy >= 85 ? 1900 + yy : 2000 + yy;\n  const date = new Date(Date.UTC(year, mm - 1, dd, hh, mi, Math.min(ss, 59)));\n  // Reject dates that rolled over, e.g. 31.02.\n  if (date.getUTCMonth() !== mm - 1 || date.getUTCDate() !== dd) return null;\n  return date;\n}\n\n/**\n * Parse the fixed 256-byte header plus the per-signal header block.\n *\n * @param buf   At least FIXED_HEADER_BYTES + ns * SIGNAL_HEADER_BYTES bytes.\n * @param fileSize Total size of the file on disk, used to derive the real record count.\n */\nexport function parseHeader(buf: Uint8Array, fileSize: number): EdfHeaderInfo {\n  const diagnostics: Diagnostic[] = [];\n  const sawComma = { value: false };\n\n  if (buf.length < FIXED_HEADER_BYTES) {\n    throw new EdfError(\n      'FILE_TOO_SMALL',\n      `File is ${fileSize} bytes; an EDF header alone needs at least ${FIXED_HEADER_BYTES}.`,\n    );\n  }\n\n  // BDF (BioSemi) marks itself with byte 255 followed by 'BIOSEMI', and stores\n  // 3-byte samples instead of 2. Everything else about the layout is identical.\n  const isBdf = buf[0] === 0xff && dec(buf, 1, 7) === 'BIOSEMI';\n  const version = isBdf ? 'BIOSEMI' : trimField(dec(buf, 0, 8));\n\n  const patientId = trimField(dec(buf, 8, 80));\n  const recordingId = trimField(dec(buf, 88, 80));\n  const startDateRaw = trimField(dec(buf, 168, 8));\n  const startTimeRaw = trimField(dec(buf, 176, 8));\n  const headerBytes = parseNumberField(dec(buf, 184, 8), 'number of header bytes', {\n    integer: true,\n    sawComma,\n  });\n  const reserved = trimField(dec(buf, 192, 44));\n  const declaredRecordCount = parseNumberField(dec(buf, 236, 8), 'number of data records', {\n    integer: true,\n    sawComma,\n  });\n  const recordDuration = parseNumberField(dec(buf, 244, 8), 'duration of a data record', {\n    sawComma,\n  });\n  const signalCount = parseNumberField(dec(buf, 252, 4), 'number of signals', {\n    integer: true,\n    sawComma,\n  });\n\n  if (signalCount <= 0) {\n    throw new EdfError(\n      'INVALID_SIGNAL_COUNT',\n      `Header declares ${signalCount} signals; expected at least 1.`,\n    );\n  }\n  if (!(recordDuration > 0)) {\n    throw new EdfError(\n      'INVALID_RECORD_DURATION',\n      `Header declares a data record duration of ${recordDuration}s; expected a positive number.`,\n    );\n  }\n\n  const expectedHeaderBytes = FIXED_HEADER_BYTES + signalCount * SIGNAL_HEADER_BYTES;\n  if (buf.length < expectedHeaderBytes) {\n    throw new EdfError(\n      'FILE_TOO_SMALL',\n      /*\n        Which of the two is actually short.\n\n        The file size was quoted either way, so a caller that had read too little — the\n        signal count parsed one way here and another way there — produced arithmetic that\n        refuted itself: \"needs a 768-byte header, but the file is only 848 bytes\". A reader\n        following that looks for a truncation that is not there.\n      */\n      `File declares ${signalCount} signals, which needs a ${expectedHeaderBytes}-byte header, ` +\n        (fileSize < expectedHeaderBytes\n          ? `but the file is only ${fileSize} bytes.`\n          : `but only ${buf.length} bytes of it were handed to the parser.`),\n    );\n  }\n  if (headerBytes !== expectedHeaderBytes) {\n    diagnostics.push({\n      code: 'HEADER_BYTES_MISMATCH',\n      severity: 'warning',\n      message:\n        `Header says it is ${headerBytes} bytes, but ${signalCount} signals require ` +\n        `${expectedHeaderBytes} bytes. Using the value computed from the signal count.`,\n    });\n  }\n\n  // Signal headers are field-major: all labels, then all transducers, and so on.\n  const base = FIXED_HEADER_BYTES;\n  const readField = (offsetUnits: number, width: number, i: number): string =>\n    dec(buf, base + offsetUnits * signalCount + i * width, width);\n\n  // EDF+ writes 'EDF+C'/'EDF+D' here; BDF+ writes 'BDF+C'/'BDF+D'. The two mean the\n  // same thing, so both are normalised to a single continuity marker.\n  const continuityTag = /^(?:EDF|BDF)\\+([CD])/u.exec(reserved);\n  const continuity: 'EDF+C' | 'EDF+D' | null =\n    continuityTag === null ? null : continuityTag[1] === 'D' ? 'EDF+D' : 'EDF+C';\n\n  const signals: EdfSignal[] = [];\n  let byteOffsetInRecord = 0;\n  const bytesPerSample = isBdf ? 3 : 2;\n  const seenLabels = new Map<string, number[]>();\n  const emptyLabels: number[] = [];\n\n  for (let i = 0; i < signalCount; i++) {\n    const label = trimField(readField(0, 16, i));\n    const transducer = trimField(readField(16, 80, i));\n    const physicalDimension = trimField(readField(96, 8, i));\n    const physicalMin = parseNumberField(readField(104, 8, i), `physical minimum (signal ${i})`, {\n      sawComma,\n    });\n    const physicalMax = parseNumberField(readField(112, 8, i), `physical maximum (signal ${i})`, {\n      sawComma,\n    });\n    const digitalMin = parseNumberField(readField(120, 8, i), `digital minimum (signal ${i})`, {\n      integer: true,\n      sawComma,\n    });\n    const digitalMax = parseNumberField(readField(128, 8, i), `digital maximum (signal ${i})`, {\n      integer: true,\n      sawComma,\n    });\n    const prefiltering = trimField(readField(136, 80, i));\n    const samplesPerRecord = parseNumberField(\n      readField(216, 8, i),\n      `samples per record (signal ${i})`,\n      { integer: true, sawComma },\n    );\n    const sigReserved = trimField(readField(224, 32, i));\n\n    if (samplesPerRecord < 0) {\n      throw new EdfError(\n        'BAD_HEADER_FIELD',\n        `Signal ${i} (\"${label}\") declares ${samplesPerRecord} samples per record.`,\n      );\n    }\n\n    const isAnnotations = label === ANNOTATIONS_LABEL || label === BDF_ANNOTATIONS_LABEL;\n\n    signals.push({\n      index: i,\n      label,\n      transducer,\n      physicalDimension,\n      physicalMin,\n      physicalMax,\n      digitalMin,\n      digitalMax,\n      prefiltering,\n      samplesPerRecord,\n      reserved: sigReserved,\n      isAnnotations,\n      samplingRate: samplesPerRecord / recordDuration,\n      byteOffsetInRecord,\n    });\n    byteOffsetInRecord += samplesPerRecord * bytesPerSample;\n\n    if (!isAnnotations) {\n      /*\n        A label is free text out of the file, and it becomes a column name in signals.csv.\n\n        `--info` has escaped control bytes since it was written, because an ANSI escape in a\n        header can drive the reader's terminal — `\\x1b[2J` clears the screen. The CSV had no\n        such protection and needed none for correctness: quoting makes any byte safe for a\n        parser, and this still passes the label through exactly as the file gives it, because\n        losing what the header says is not an improvement.\n\n        What was missing is the sentence saying so. A recording whose channel is labelled\n        `\\x1b[2Jgone` converted with no warning at all, and `cat signals.csv` then cleared\n        the terminal — while a script referencing that column by name carried an invisible\n        control character in it. NONPRINTABLE_LABEL has been declared and documented as\n        reserved since 0.1; this is it doing its job.\n      */\n      /*\n        Which of the two fields carries them, because the consequences are not the same.\n\n        The message said \"label or unit\", and then said the bytes \"will appear in the CSV\n        column name\" and that \"the name cannot be typed\" — both of which are about the label.\n        A channel labelled plainly `ECG` in a unit of `u\\x07V` got all of it: its column is\n        `ECG`, `--channels ECG` selects it and exits 0, and the byte is in channels.csv's\n        `unit` cell, which the warning never mentioned. Three sentences, none of them true of\n        the file that raised it, on a warning whose whole purpose is to say where an invisible\n        byte went.\n      */\n      const inLabel = [...label].filter(isControlCharacter);\n      const inUnit = [...physicalDimension].filter(isControlCharacter);\n      const control = [...inLabel, ...inUnit];\n      if (control.length > 0) {\n        const shown = [...new Set(control)]\n          .map((c) => `\\\\x${(c.codePointAt(0) as number).toString(16).padStart(2, '0')}`)\n          .join(', ');\n        const plural = control.length === 1 ? '' : 's';\n        const both = inLabel.length > 0 && inUnit.length > 0;\n        const field = both ? 'label and unit contain' : inLabel.length > 0 ? 'label contains' : 'unit contains';\n        // Where they land, which is the question the reader has. A label becomes a column\n        // name in signals.csv; a unit is a cell of channels.csv and nothing else.\n        const lands = both\n          ? 'which will appear in the CSV column name and in channels.csv\\'s unit cell'\n          : inLabel.length > 0\n            ? 'which will appear in the CSV column name'\n            : 'which will appear in channels.csv\\'s unit cell';\n        diagnostics.push({\n          code: 'NONPRINTABLE_LABEL',\n          severity: 'warning',\n          message:\n            `Signal ${i}'s ${field} ${control.length} control character${plural} ` +\n            `(${shown}), ${lands} exactly as the header has them.`,\n          hint:\n            (inLabel.length > 0\n              ? `Address the channel by position with --channels \"#${i}\" rather than by name, ` +\n                'since the name cannot be typed. '\n              : `The column name is unaffected, so --channels \"${label}\" still selects it. `) +\n            'Printing the CSV to a terminal may do more than print it.',\n        });\n      }\n\n      if (label === '') {\n        // Collected, not reported here: what this channel's column ends up called depends on\n        // whether some later channel is literally labelled `signal_<i>`, and inside this loop\n        // the later channels do not exist yet. See the pass below.\n        emptyLabels.push(i);\n      } else {\n        // Collected rather than reported here: a label repeated five times should\n        // produce one warning naming all five, not four near-identical pairs.\n        const seen = seenLabels.get(label);\n        if (seen) seen.push(i);\n        else seenLabels.set(label, [i]);\n      }\n\n      if (samplesPerRecord === 0) {\n        diagnostics.push({\n          code: 'NO_SAMPLES',\n          severity: 'warning',\n          message: `Signal ${i} (\"${label}\") carries no samples at all (0 per data record).`,\n          hint: 'It is described in channels.csv but left out of the converted data.',\n        });\n      }\n\n      /*\n        Too large to represent, and too small — the second was silent.\n\n        The gain is the span divided by the digital range, and a span of 2e-320 over 65,535\n        codes is 3e-325: below the smallest subnormal double, so it underflows to +0. The\n        scaler's flat-range branch then handed every code the same physical value, and a\n        channel of 65,536 distinct readings became one repeated number with nothing raised at\n        all. One power of ten away, at 1e-319, the same file raises VALUE_RESOLUTION.\n\n        Both are the same fact about the header — the span cannot be turned into a mapping —\n        so both get this code, and both leave the cells empty rather than filling them with a\n        value the header cannot justify.\n      */\n      const span = physicalMax - physicalMin;\n      const underflowed = span !== 0 && span / (digitalMax - digitalMin) === 0;\n      if (!Number.isFinite(span) || underflowed) {\n        diagnostics.push({\n          code: 'UNUSABLE_PHYSICAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") declares a physical range from ${physicalMin} to ` +\n            `${physicalMax}, whose span is too ${underflowed ? 'small' : 'large'} to ` +\n            `represent, so its values cannot be scaled.`,\n          hint: 'Its cells are left empty rather than filled with a value the header cannot justify.',\n        });\n      } else if (digitalMax === digitalMin) {\n        diagnostics.push({\n          code: 'DEGENERATE_DIGITAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") has digital minimum equal to digital maximum ` +\n            `(${digitalMin}), so its values cannot be scaled.`,\n          hint: 'Its cells are left empty rather than filled with a value the header cannot justify.',\n        });\n      } else if (physicalMax === physicalMin) {\n        diagnostics.push({\n          code: 'DEGENERATE_PHYSICAL_RANGE',\n          severity: 'warning',\n          message:\n            `Signal ${i} (\"${label}\") has physical minimum equal to physical maximum ` +\n            `(${physicalMin}), so every sample converts to the same value.`,\n        });\n      } else if ((physicalMax - physicalMin) * (digitalMax - digitalMin) < 0) {\n        /*\n          Polarity is inverted when the gain is negative, and the gain is\n          (physicalMax - physicalMin) / (digitalMax - digitalMin) — so it is the sign of the\n          two spans together that matters, not the physical pair alone.\n\n          Testing only `physicalMax < physicalMin` was wrong in both directions. A file with\n          its DIGITAL bounds reversed is just as inverted and drew no warning at all, handing\n          back sign-flipped EEG with nothing to indicate it. A file with BOTH pairs reversed\n          has a positive gain and is not inverted, yet was warned about — a message that was\n          simply untrue of that recording.\n        */\n        const reversed =\n          physicalMax < physicalMin\n            ? `physical minimum ${physicalMin} above physical maximum ${physicalMax}`\n            : `digital minimum ${digitalMin} above digital maximum ${digitalMax}`;\n        diagnostics.push({\n          code: 'INVERTED_PHYSICAL_RANGE',\n          severity: 'warning',\n          message: `Signal ${i} (\"${label}\") declares ${reversed}, which inverts its polarity.`,\n          hint: 'The values are converted exactly as the header specifies, inversion included.',\n        });\n      }\n    }\n  }\n\n  /*\n    What an unlabelled channel is actually called, which the message used to guess.\n\n    A channel with no label takes `signal_<index>` — unless another channel is literally\n    labelled that, which EDF permits, since labels are free text and nothing enforces anything\n    about them. Then both collide and both are suffixed. The warning said \"It will appear as\n    \"signal_0\"\" while the file's header read `time_s,signal_0_ch0,signal_0_ch1`: the one\n    sentence the run printed named a column that exists in neither signals.csv nor\n    channels.csv.\n\n    The other half was silent. The channel that genuinely carries the label `signal_0` lost\n    its own column name to a collision with a synthesised one, and nothing said so —\n    DUPLICATE_LABEL did not fire, because the two labels are not the same label. Both halves\n    are one sentence here, because they are one event.\n\n    No specific suffixed name is quoted. The suffix rule has a second pass for names that are\n    still shared afterwards, and a message that hard-coded `_ch<index>` would be guessing again\n    in exactly the way this is fixing.\n  */\n  for (const index of emptyLabels) {\n    const taken = seenLabels.get(`signal_${index}`);\n    diagnostics.push({\n      code: 'EMPTY_LABEL',\n      severity: 'warning',\n      message:\n        taken === undefined\n          ? `Signal ${index} has no label. It will appear as \"signal_${index}\".`\n          : `Signal ${index} has no label, so it takes the name \"signal_${index}\" — which ` +\n            `${taken.length === 1 ? 'signal' : 'signals'} ${listed(taken.map(String))} already ` +\n            `${taken.length === 1 ? 'carries' : 'carry'} as a label, so both columns are ` +\n            `suffixed with their position instead.`,\n    });\n  }\n\n  for (const [label, indices] of seenLabels) {\n    if (indices.length < 2) continue;\n    diagnostics.push({\n      code: 'DUPLICATE_LABEL',\n      severity: 'warning',\n      message: `${indices.length} signals share the label \"${label}\" (positions ${indices.join(', ')}).`,\n      hint: 'Their columns are suffixed with the signal number so they stay distinguishable.',\n    });\n  }\n\n  const recordBytes = byteOffsetInRecord;\n  if (recordBytes <= 0) {\n    throw new EdfError(\n      'NO_SAMPLES',\n      'No signal in this file carries any samples (every channel declares 0 samples per record).',\n    );\n  }\n\n  if (sawComma.value) {\n    diagnostics.push({\n      code: 'COMMA_DECIMAL',\n      severity: 'warning',\n      message: 'Some header numbers use a comma decimal separator, which the EDF spec does not allow.',\n      hint: 'They were read as decimal points. Check the values in the channel table.',\n    });\n  }\n\n  const dataBytes = fileSize - expectedHeaderBytes;\n  if (dataBytes < 0) {\n    throw new EdfError('FILE_TOO_SMALL', `File is smaller than its own header.`);\n  }\n  const recordCount = Math.floor(dataBytes / recordBytes);\n  const trailingBytes = dataBytes - recordCount * recordBytes;\n\n  if (recordCount === 0) {\n    throw new EdfError(\n      'NO_DATA_RECORDS',\n      'The file contains a header but no complete data record.',\n      'The recording was probably interrupted before any data was written.',\n    );\n  }\n\n  if (declaredRecordCount === -1) {\n    diagnostics.push({\n      code: 'RECORD_COUNT_UNKNOWN',\n      severity: 'warning',\n      message:\n        `The header does not say how many data records the file has (-1), which the spec allows ` +\n        `for recordings still in progress. Using the ${counted(recordCount, 'record')} the file actually contains.`,\n    });\n  } else if (declaredRecordCount !== recordCount) {\n    diagnostics.push({\n      code: 'RECORD_COUNT_MISMATCH',\n      severity: 'warning',\n      message:\n        `The header declares ${declaredRecordCount} data records but the file contains ` +\n        `${recordCount}. Converting the ${counted(recordCount, 'record')} that ${recordCount === 1 ? 'is' : 'are'} present.`,\n      hint:\n        declaredRecordCount > recordCount\n          ? 'The recording looks truncated. It may have been cut short or copied incompletely.'\n          : 'The file is longer than its header claims.',\n    });\n  }\n\n  if (trailingBytes > 0) {\n    diagnostics.push({\n      code: 'TRAILING_BYTES',\n      severity: 'warning',\n      message: `${counted(trailingBytes, 'byte')} after the last complete data record ${trailingBytes === 1 ? 'was' : 'were'} ignored.`,\n    });\n  }\n\n  const isEdfPlus = continuity !== null;\n  if (continuity === 'EDF+D') {\n    diagnostics.push({\n      code: 'DISCONTINUOUS',\n      severity: 'warning',\n      message:\n        `This is a discontinuous (${isBdf ? 'BDF+D' : 'EDF+D'}) recording: its data records are ` +\n        `not contiguous in time.`,\n      hint: 'Each row carries its true recording time, so gaps stay visible instead of being closed.',\n    });\n  }\n\n  const dataSignals = signals.filter((s) => !s.isAnnotations);\n  if (dataSignals.length === 0) {\n    diagnostics.push({\n      code: 'NO_SIGNAL_CHANNELS',\n      severity: 'warning',\n      message: 'This file has no signal channels; it contains only EDF+ annotations.',\n    });\n  }\n\n  // A channel declaring zero samples per record has no sampling rate to speak of — it is\n  // reported separately as NO_SAMPLES and no file is written for it. Counting its nominal\n  // 0 Hz as a rate made a single-rate recording warn that it used \"2 different sampling\n  // rates (4 Hz, 0 Hz)\" and claim it was splitting output it never split.\n  const rates = new Set(dataSignals.filter((s) => s.samplesPerRecord > 0).map((s) => s.samplingRate));\n  if (rates.size > 1) {\n    diagnostics.push({\n      code: 'MIXED_SAMPLING_RATES',\n      severity: 'warning',\n      message:\n        `Channels use ${rates.size} different sampling rates ` +\n        `(${listed(formatRates([...rates].sort((a, b) => b - a)).map((r) => `${r} Hz`))}).`,\n      hint: 'They are written to one file per rate so no channel is resampled.',\n    });\n  }\n\n  return {\n    header: {\n      version,\n      patientId,\n      recordingId,\n      startDateRaw,\n      startTimeRaw,\n      startDateTime: resolveStartDateTime(startDateRaw, startTimeRaw),\n      headerBytes: expectedHeaderBytes,\n      declaredHeaderBytes: headerBytes,\n      reserved,\n      isEdfPlus,\n      isBdf,\n      continuity,\n      declaredRecordCount,\n      recordDuration,\n      signalCount,\n      signals,\n      bytesPerSample,\n      recordBytes,\n    },\n    recordCount,\n    trailingBytes,\n    diagnostics,\n  };\n}\n\n/**\n * The recording start as a zone-less wall clock, \"YYYY-MM-DDTHH:MM:SS\".\n *\n * EDF stores the start time as local wall-clock digits with no timezone anywhere in\n * the format. `startDateTime` is built with Date.UTC purely so those digits survive a\n * round trip unshifted, which makes it a carrier for the wall clock rather than a\n * real instant. Serialising it with `toISOString()` would append a Z and assert UTC,\n * and any reader converting to local time would then shift the recording by their own\n * offset: 13:43:04 in the file becomes 08:43:04 in New York. The Z is omitted because\n * the file genuinely does not say which zone it meant.\n */\nexport function formatWallClock(date: Date | null): string | null {\n  if (!date) return null;\n  return date.toISOString().slice(0, 19);\n}\n\n/** \"EDF\", \"EDF+ (EDF+D)\", \"BDF\", \"BDF+ (EDF+C)\". */\nexport function describeFormat(header: EdfHeader): string {\n  const base = header.isBdf ? 'BDF' : 'EDF';\n  if (!header.isEdfPlus) return base;\n  return `${base}+ (${header.continuity === 'EDF+D' ? 'discontinuous' : 'continuous'})`;\n}\n\n/** Render a sampling rate without trailing noise: 256, 0.5, 12.5. */\nexport function formatRate(hz: number): string {\n  if (Number.isInteger(hz)) return String(hz);\n  const rounded = Number(hz.toFixed(6));\n  // A rate below 5e-7 rounds away to \"0\", which reads as \"this channel has no sampling\n  // rate\" and made the mixed-rate warning contradict itself: it announced two different\n  // rates and then printed both as \"0 Hz\". Exponent form keeps a real rate legible, and\n  // keeps distinct rates distinct in the channel table and in output filenames.\n  if (rounded === 0) return hz.toExponential(3);\n  return String(rounded);\n}\n\n/**\n * Renders a group of rates so that rates which differ read as differing.\n *\n * `formatRate` rounds to six decimals, which is what keeps an ordinary rate free of\n * float noise — 30 samples in a 0.1-second record is 299.99999999999994 as a double,\n * and belongs on screen as 300. Two rates separated by less than that round to one\n * string, so a file carrying 1e-6 Hz and 1.25e-6 Hz warned that it used \"2 different\n * sampling rates (0.000001 Hz, 0.000001 Hz)\" and named both files the same thing.\n *\n * That is the contradiction the exponent fallback above already removes for rates that\n * round away to zero; this is the same one a step further out. On a collision every rate\n * in the group switches to its shortest exact form, which is unique for distinct values,\n * rather than only the pair that collided — one column in one notation reads better than\n * two.\n */\nexport function formatRates(rates: readonly number[]): string[] {\n  const rounded = rates.map(formatRate);\n  const distinct = new Set(rates).size;\n  return new Set(rounded).size === distinct ? rounded : rates.map((hz) => String(hz));\n}\n"]}

package/dist/edf/scale.js

@@ -37,10 +37,24 @@ export function makeScaler(signal) {
     if (digitalMax === digitalMin)
         return () => NaN;
     const gain = (physicalMax - physicalMin) / (digitalMax - digitalMin);
-    // A flat physical range makes every sample the same value, and would divide by
-    // zero in the offset below. That mapping is defined, so its constant is written.
+    /*
+      A flat physical range makes every sample the same value, and would divide by zero in the
+      offset below. That mapping is defined, so its constant is written.
+  
+      A gain of zero does not always mean flat, and this could not tell the difference. A range
+      of -1e-320 to 1e-320 is not flat — it is 65,536 distinct physical values — but the gain
+      is 2e-320/65535, which is smaller than the smallest subnormal double and underflows to
+      +0. Every distinct sample then took `physicalMin`, so eight codes spanning -16,000 to
+      +12,000 came out as one repeated number, with no diagnostic anywhere and `--strict`
+      exiting 0. At 1e-319, one power of ten away, the same file raises VALUE_RESOLUTION.
+  
+      That is the same situation as the overflow below it, which this codebase already reasoned
+      about and answered: the span cannot be represented, so there is no mapping, so the cells
+      are left empty rather than filled with a value the header cannot justify. Underflow only
+      got the flat-range treatment because `gain === 0` is what both look like from here.
+    */
     if (gain === 0)
-        return () => physicalMin;
+        return physicalMax === physicalMin ? () => physicalMin : () => NaN;
     // A non-finite gain is a different thing: the physical span overflowed a double, so
     // there is no mapping at all. Returning physicalMin filled the column with one enormous
     // constant — every distinct sample rendered as the same 300-digit number — and raised

package/dist/edf/scale.js.map

@@ -1 +1 @@
-{"version":3,"file":"scale.js","sourceRoot":"","sources":["../../src/edf/scale.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAMH,MAAM,UAAU,UAAU,CAAC,MAAiB;IAC1C,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,WAAW,EAAE,GAAG,MAAM,CAAC;IAEpE,oFAAoF;IACpF,gEAAgE;IAChE,EAAE;IACF,oFAAoF;IACpF,oFAAoF;IACpF,qFAAqF;IACrF,oFAAoF;IACpF,iFAAiF;IACjF,oEAAoE;IACpE,IAAI,UAAU,KAAK,UAAU;QAAE,OAAO,GAAG,EAAE,CAAC,GAAG,CAAC;IAEhD,MAAM,IAAI,GAAG,CAAC,WAAW,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,UAAU,CAAC,CAAC;IAErE,+EAA+E;IAC/E,iFAAiF;IACjF,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,GAAG,EAAE,CAAC,WAAW,CAAC;IAEzC,oFAAoF;IACpF,wFAAwF;IACxF,sFAAsF;IACtF,qFAAqF;IACrF,kDAAkD;IAClD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,GAAG,EAAE,CAAC,GAAG,CAAC;IAE7C,mFAAmF;IACnF,kFAAkF;IAClF,6EAA6E;IAC7E,6CAA6C;IAC7C,MAAM,MAAM,GAAG,WAAW,GAAG,IAAI,GAAG,UAAU,CAAC;IAC/C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC7B,OAAO,CAAC,OAAe,EAAU,EAAE,CAAC,CAAC,OAAO,GAAG,UAAU,CAAC,GAAG,IAAI,GAAG,WAAW,CAAC;IAClF,CAAC;IAED,OAAO,CAAC,OAAe,EAAU,EAAE,CAAC,IAAI,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,CAAC;AAChE,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAiB;IAChD,MAAM,WAAW,GAAG,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;IAC1D,IAAI,WAAW,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAChC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,MAAM,CAAC,WAAW,CAAC,GAAG,WAAW,CAAC,CAAC;AAC3E,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,oBAAoB,GAAG,GAAG,CAAC;AAEjC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,iBAAiB,CAAC,MAAiB,EAAE,GAAG,GAAG,oBAAoB;IAC7E,MAAM,IAAI,GAAG,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IACpD,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;IAChD,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;AAC5C,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAiB;IAClD,MAAM,IAAI,GAAG,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,KAAK,CAAC;IACxD,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC;AACjE,CAAC","sourcesContent":["/**\n * Digital-to-physical conversion.\n *\n * EDF defines the mapping by two calibration points, (digitalMin -> physicalMin)\n * and (digitalMax -> physicalMax), which the specification writes as:\n *\n *   gain     = (physicalMax - physicalMin) / (digitalMax - digitalMin)\n *   physical = (digital - digitalMin) * gain + physicalMin\n *\n * That form is evaluated here in EDFlib's algebraically equivalent arrangement:\n *\n *   offset   = physicalMax / gain - digitalMax\n *   physical = gain * (offset + digital)\n *\n * The rearrangement is not cosmetic. Written the first way, a channel spanning\n * +/-800 uV computes a value near 800 and then subtracts 800, and the cancellation\n * throws away low-order bits: digital 0 yields 0.19536019536019467 when the exact\n * value is 0.19536019536019536. EDFlib's form keeps the intermediate small\n * (offset + digital = 0.5 here) and returns the correctly rounded result.\n *\n * Both properties matter. The values are as accurate as a double can express, and\n * they are bit-identical to pyEDFlib and EDFbrowser, which share EDFlib's arithmetic,\n * so the test suite can assert exact equality against a reference implementation\n * rather than settling for a tolerance.\n */\n\nimport type { EdfSignal } from './header.js';\n\nexport type Scaler = (digital: number) => number;\n\nexport function makeScaler(signal: EdfSignal): Scaler {\n  const { digitalMin, digitalMax, physicalMin, physicalMax } = signal;\n\n  // A zero digital span leaves the mapping undefined — the header contradicts itself,\n  // so there is no physical value for any sample on this channel.\n  //\n  // NaN rather than a stand-in number. Writing the physical minimum produces a column\n  // of plausible readings (\"-100.000\" repeated) that is indistinguishable from a real\n  // flat recording once the CSV is opened somewhere else, which is exactly the kind of\n  // invented data this tool exists to avoid. NaN carries through to an empty CSV cell\n  // and reads back as NaN in pandas, matching how a missing annotation duration is\n  // already written. DEGENERATE_DIGITAL_RANGE is raised alongside it.\n  if (digitalMax === digitalMin) return () => NaN;\n\n  const gain = (physicalMax - physicalMin) / (digitalMax - digitalMin);\n\n  // A flat physical range makes every sample the same value, and would divide by\n  // zero in the offset below. That mapping is defined, so its constant is written.\n  if (gain === 0) return () => physicalMin;\n\n  // A non-finite gain is a different thing: the physical span overflowed a double, so\n  // there is no mapping at all. Returning physicalMin filled the column with one enormous\n  // constant — every distinct sample rendered as the same 300-digit number — and raised\n  // nothing. NaN takes the same route as a degenerate digital range: empty cells, plus\n  // UNUSABLE_PHYSICAL_RANGE from the header parser.\n  if (!Number.isFinite(gain)) return () => NaN;\n\n  // Deriving the offset divides by the gain. For every realistic calibration that is\n  // both safe and more accurate, but an absurd header (a huge physical range over a\n  // near-zero gain) could overflow it, so fall back to the specification's own\n  // arrangement rather than emitting Infinity.\n  const offset = physicalMax / gain - digitalMax;\n  if (!Number.isFinite(offset)) {\n    return (digital: number): number => (digital - digitalMin) * gain + physicalMin;\n  }\n\n  return (digital: number): number => gain * (offset + digital);\n}\n\n/**\n * Smallest physical step this channel can express — one digital unit.\n * Used to choose a decimal precision that preserves every distinct sample value.\n */\nexport function quantizationStep(signal: EdfSignal): number {\n  const digitalSpan = signal.digitalMax - signal.digitalMin;\n  if (digitalSpan === 0) return 0;\n  return Math.abs((signal.physicalMax - signal.physicalMin) / digitalSpan);\n}\n\n/**\n * The most `toFixed` accepts. 101 is a RangeError, so this is the ceiling, not a taste.\n *\n * It used to be 20, on the stated grounds that 20 was what `toFixed` allowed. It is not,\n * and the gap was not academic: a magnetometer channel spanning ±1e-16 T over a 16-bit\n * converter has a step of 3.05e-21 and needs 23 places. Clamped to 20, every value landed\n * on a 1e-20 grid — about three digital codes to a printed value — so 69% of the samples\n * could not be recovered, the conversion exited 0, and nothing said a word. The channel\n * type the old comment named as the reason for the ceiling was the one it broke.\n */\nconst MAX_DERIVED_DECIMALS = 100;\n\n/**\n * Decimal places needed so that two adjacent digital codes never round to the same\n * string. Two places past the quantization step keep rounding error far below the\n * resolution the hardware actually recorded, without padding the file with digits\n * that carry no information.\n *\n * Ordinary channels land at three or four: a ±800 µV channel over 12 bits steps by\n * 0.39 µV and needs three. The ceiling is only reached by calibrations whose step is\n * below 1e-98, which an 8-character physical bound can still express — `1e-99` is five\n * characters. Those get VALUE_RESOLUTION rather than silence.\n */\nexport function decimalsForSignal(signal: EdfSignal, max = MAX_DERIVED_DECIMALS): number {\n  const step = quantizationStep(signal);\n  if (!(step > 0) || !Number.isFinite(step)) return 3;\n  const needed = Math.ceil(-Math.log10(step)) + 2;\n  return Math.min(max, Math.max(0, needed));\n}\n\n/**\n * Whether this channel's step is finer than any precision the tool can print.\n *\n * Asked of the ceiling, not of the precision in use. `--decimals 2` on a channel needing 3\n * is a trade the caller made knowingly and is not this warning's business — 0.5.10 fixed a\n * version of this that fired on every ordinary EEG at `--decimals 2` and made\n * `--decimals 2 --strict` impossible. But it fixed it by asking \"did the caller choose the\n * precision\", which suppressed the real case too: at `--decimals 20` a channel stepping by\n * 1e-106 printed every one of its codes as `0.00000000000000000000`, in silence.\n *\n * The question is whether anything the tool can print would separate consecutive codes. When\n * the answer is no, that is a ceiling nobody chose, and it holds whatever `--decimals` says.\n */\nexport function decimalsAreClamped(signal: EdfSignal): boolean {\n  const step = quantizationStep(signal);\n  if (!(step > 0) || !Number.isFinite(step)) return false;\n  return Math.ceil(-Math.log10(step)) + 2 > MAX_DERIVED_DECIMALS;\n}\n"]}
+{"version":3,"file":"scale.js","sourceRoot":"","sources":["../../src/edf/scale.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAMH,MAAM,UAAU,UAAU,CAAC,MAAiB;IAC1C,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,WAAW,EAAE,GAAG,MAAM,CAAC;IAEpE,oFAAoF;IACpF,gEAAgE;IAChE,EAAE;IACF,oFAAoF;IACpF,oFAAoF;IACpF,qFAAqF;IACrF,oFAAoF;IACpF,iFAAiF;IACjF,oEAAoE;IACpE,IAAI,UAAU,KAAK,UAAU;QAAE,OAAO,GAAG,EAAE,CAAC,GAAG,CAAC;IAEhD,MAAM,IAAI,GAAG,CAAC,WAAW,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,GAAG,UAAU,CAAC,CAAC;IAErE;;;;;;;;;;;;;;;MAeE;IACF,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,WAAW,KAAK,WAAW,CAAC,CAAC,CAAC,GAAW,EAAE,CAAC,WAAW,CAAC,CAAC,CAAC,GAAW,EAAE,CAAC,GAAG,CAAC;IAEnG,oFAAoF;IACpF,wFAAwF;IACxF,sFAAsF;IACtF,qFAAqF;IACrF,kDAAkD;IAClD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,GAAG,EAAE,CAAC,GAAG,CAAC;IAE7C,mFAAmF;IACnF,kFAAkF;IAClF,6EAA6E;IAC7E,6CAA6C;IAC7C,MAAM,MAAM,GAAG,WAAW,GAAG,IAAI,GAAG,UAAU,CAAC;IAC/C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC7B,OAAO,CAAC,OAAe,EAAU,EAAE,CAAC,CAAC,OAAO,GAAG,UAAU,CAAC,GAAG,IAAI,GAAG,WAAW,CAAC;IAClF,CAAC;IAED,OAAO,CAAC,OAAe,EAAU,EAAE,CAAC,IAAI,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,CAAC;AAChE,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAiB;IAChD,MAAM,WAAW,GAAG,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC;IAC1D,IAAI,WAAW,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAChC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,WAAW,GAAG,MAAM,CAAC,WAAW,CAAC,GAAG,WAAW,CAAC,CAAC;AAC3E,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,oBAAoB,GAAG,GAAG,CAAC;AAEjC;;;;;;;;;;GAUG;AACH,MAAM,UAAU,iBAAiB,CAAC,MAAiB,EAAE,GAAG,GAAG,oBAAoB;IAC7E,MAAM,IAAI,GAAG,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,CAAC,CAAC;IACpD,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;IAChD,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;AAC5C,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAiB;IAClD,MAAM,IAAI,GAAG,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,OAAO,KAAK,CAAC;IACxD,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC;AACjE,CAAC","sourcesContent":["/**\n * Digital-to-physical conversion.\n *\n * EDF defines the mapping by two calibration points, (digitalMin -> physicalMin)\n * and (digitalMax -> physicalMax), which the specification writes as:\n *\n *   gain     = (physicalMax - physicalMin) / (digitalMax - digitalMin)\n *   physical = (digital - digitalMin) * gain + physicalMin\n *\n * That form is evaluated here in EDFlib's algebraically equivalent arrangement:\n *\n *   offset   = physicalMax / gain - digitalMax\n *   physical = gain * (offset + digital)\n *\n * The rearrangement is not cosmetic. Written the first way, a channel spanning\n * +/-800 uV computes a value near 800 and then subtracts 800, and the cancellation\n * throws away low-order bits: digital 0 yields 0.19536019536019467 when the exact\n * value is 0.19536019536019536. EDFlib's form keeps the intermediate small\n * (offset + digital = 0.5 here) and returns the correctly rounded result.\n *\n * Both properties matter. The values are as accurate as a double can express, and\n * they are bit-identical to pyEDFlib and EDFbrowser, which share EDFlib's arithmetic,\n * so the test suite can assert exact equality against a reference implementation\n * rather than settling for a tolerance.\n */\n\nimport type { EdfSignal } from './header.js';\n\nexport type Scaler = (digital: number) => number;\n\nexport function makeScaler(signal: EdfSignal): Scaler {\n  const { digitalMin, digitalMax, physicalMin, physicalMax } = signal;\n\n  // A zero digital span leaves the mapping undefined — the header contradicts itself,\n  // so there is no physical value for any sample on this channel.\n  //\n  // NaN rather than a stand-in number. Writing the physical minimum produces a column\n  // of plausible readings (\"-100.000\" repeated) that is indistinguishable from a real\n  // flat recording once the CSV is opened somewhere else, which is exactly the kind of\n  // invented data this tool exists to avoid. NaN carries through to an empty CSV cell\n  // and reads back as NaN in pandas, matching how a missing annotation duration is\n  // already written. DEGENERATE_DIGITAL_RANGE is raised alongside it.\n  if (digitalMax === digitalMin) return () => NaN;\n\n  const gain = (physicalMax - physicalMin) / (digitalMax - digitalMin);\n\n  /*\n    A flat physical range makes every sample the same value, and would divide by zero in the\n    offset below. That mapping is defined, so its constant is written.\n\n    A gain of zero does not always mean flat, and this could not tell the difference. A range\n    of -1e-320 to 1e-320 is not flat — it is 65,536 distinct physical values — but the gain\n    is 2e-320/65535, which is smaller than the smallest subnormal double and underflows to\n    +0. Every distinct sample then took `physicalMin`, so eight codes spanning -16,000 to\n    +12,000 came out as one repeated number, with no diagnostic anywhere and `--strict`\n    exiting 0. At 1e-319, one power of ten away, the same file raises VALUE_RESOLUTION.\n\n    That is the same situation as the overflow below it, which this codebase already reasoned\n    about and answered: the span cannot be represented, so there is no mapping, so the cells\n    are left empty rather than filled with a value the header cannot justify. Underflow only\n    got the flat-range treatment because `gain === 0` is what both look like from here.\n  */\n  if (gain === 0) return physicalMax === physicalMin ? (): number => physicalMin : (): number => NaN;\n\n  // A non-finite gain is a different thing: the physical span overflowed a double, so\n  // there is no mapping at all. Returning physicalMin filled the column with one enormous\n  // constant — every distinct sample rendered as the same 300-digit number — and raised\n  // nothing. NaN takes the same route as a degenerate digital range: empty cells, plus\n  // UNUSABLE_PHYSICAL_RANGE from the header parser.\n  if (!Number.isFinite(gain)) return () => NaN;\n\n  // Deriving the offset divides by the gain. For every realistic calibration that is\n  // both safe and more accurate, but an absurd header (a huge physical range over a\n  // near-zero gain) could overflow it, so fall back to the specification's own\n  // arrangement rather than emitting Infinity.\n  const offset = physicalMax / gain - digitalMax;\n  if (!Number.isFinite(offset)) {\n    return (digital: number): number => (digital - digitalMin) * gain + physicalMin;\n  }\n\n  return (digital: number): number => gain * (offset + digital);\n}\n\n/**\n * Smallest physical step this channel can express — one digital unit.\n * Used to choose a decimal precision that preserves every distinct sample value.\n */\nexport function quantizationStep(signal: EdfSignal): number {\n  const digitalSpan = signal.digitalMax - signal.digitalMin;\n  if (digitalSpan === 0) return 0;\n  return Math.abs((signal.physicalMax - signal.physicalMin) / digitalSpan);\n}\n\n/**\n * The most `toFixed` accepts. 101 is a RangeError, so this is the ceiling, not a taste.\n *\n * It used to be 20, on the stated grounds that 20 was what `toFixed` allowed. It is not,\n * and the gap was not academic: a magnetometer channel spanning ±1e-16 T over a 16-bit\n * converter has a step of 3.05e-21 and needs 23 places. Clamped to 20, every value landed\n * on a 1e-20 grid — about three digital codes to a printed value — so 69% of the samples\n * could not be recovered, the conversion exited 0, and nothing said a word. The channel\n * type the old comment named as the reason for the ceiling was the one it broke.\n */\nconst MAX_DERIVED_DECIMALS = 100;\n\n/**\n * Decimal places needed so that two adjacent digital codes never round to the same\n * string. Two places past the quantization step keep rounding error far below the\n * resolution the hardware actually recorded, without padding the file with digits\n * that carry no information.\n *\n * Ordinary channels land at three or four: a ±800 µV channel over 12 bits steps by\n * 0.39 µV and needs three. The ceiling is only reached by calibrations whose step is\n * below 1e-98, which an 8-character physical bound can still express — `1e-99` is five\n * characters. Those get VALUE_RESOLUTION rather than silence.\n */\nexport function decimalsForSignal(signal: EdfSignal, max = MAX_DERIVED_DECIMALS): number {\n  const step = quantizationStep(signal);\n  if (!(step > 0) || !Number.isFinite(step)) return 3;\n  const needed = Math.ceil(-Math.log10(step)) + 2;\n  return Math.min(max, Math.max(0, needed));\n}\n\n/**\n * Whether this channel's step is finer than any precision the tool can print.\n *\n * Asked of the ceiling, not of the precision in use. `--decimals 2` on a channel needing 3\n * is a trade the caller made knowingly and is not this warning's business — 0.5.10 fixed a\n * version of this that fired on every ordinary EEG at `--decimals 2` and made\n * `--decimals 2 --strict` impossible. But it fixed it by asking \"did the caller choose the\n * precision\", which suppressed the real case too: at `--decimals 20` a channel stepping by\n * 1e-106 printed every one of its codes as `0.00000000000000000000`, in silence.\n *\n * The question is whether anything the tool can print would separate consecutive codes. When\n * the answer is no, that is a ceiling nobody chose, and it holds whatever `--decimals` says.\n */\nexport function decimalsAreClamped(signal: EdfSignal): boolean {\n  const step = quantizationStep(signal);\n  if (!(step > 0) || !Number.isFinite(step)) return false;\n  return Math.ceil(-Math.log10(step)) + 2 > MAX_DERIVED_DECIMALS;\n}\n"]}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "edf2csv",
-  "version": "0.5.82",
+  "version": "0.5.84",
   "description": "Convert EDF, EDF+ and BDF biosignal recordings (European Data Format) to CSV from the command line. Local, streaming, and never resamples or alters units.",
   "keywords": [
     "edf",