@libraz/libcantus 1.2.0 → 1.3.0

package/dist/{chunk-HPVSH7SF.cjs → chunk-YFA7X2F6.cjs}

@@ -30,7 +30,6 @@
 
 
 
-var _chunkSGLEQVJHcjs = require('./chunk-SGLEQVJH.cjs');
 
 
 
@@ -41,9 +40,27 @@ var _chunkSGLEQVJHcjs = require('./chunk-SGLEQVJH.cjs');
 
 
 
-var _chunk3R3QVKY5cjs = require('./chunk-3R3QVKY5.cjs');
 
 
+
+
+
+
+
+
+
+
+
+
+var _chunkIFZDWIKMcjs = require('./chunk-IFZDWIKM.cjs');
+
+
+
+
+
+
+
+var _chunk3R3QVKY5cjs = require('./chunk-3R3QVKY5.cjs');
 
 
 
@@ -108,7 +125,7 @@ function assertTrackNotes(tracks, options) {
 }
 function assertGivenTimeline(timeline, name) {
   if (timeline !== void 0) {
-    _chunkSGLEQVJHcjs.assertChordTimeline.call(void 0, timeline, name);
+    _chunkIFZDWIKMcjs.assertChordTimeline.call(void 0, timeline, name);
   }
 }
 function assertGivenKeys(keys, name) {
@@ -195,7 +212,7 @@ function splitIntoSubVoices(ordered) {
         continue;
       }
       const end = last.startBeat + last.durationBeat;
-      if (!_chunkSGLEQVJHcjs.hasEnded.call(void 0, end, note.startBeat)) {
+      if (!_chunkIFZDWIKMcjs.hasEnded.call(void 0, end, note.startBeat)) {
         continue;
       }
       const distance = Math.abs(last.pitch - note.pitch);
@@ -216,13 +233,13 @@ function splitIntoSubVoices(ordered) {
   for (let start = 0; start < ordered.length; ) {
     const onset = _nullishCoalesce(_optionalChain([ordered, 'access', _13 => _13[start], 'optionalAccess', _14 => _14.startBeat]), () => ( 0));
     let end = start + 1;
-    while (end < ordered.length && _chunkSGLEQVJHcjs.sameInstant.call(void 0, _nullishCoalesce(_optionalChain([ordered, 'access', _15 => _15[end], 'optionalAccess', _16 => _16.startBeat]), () => ( 0)), onset)) {
+    while (end < ordered.length && _chunkIFZDWIKMcjs.sameInstant.call(void 0, _nullishCoalesce(_optionalChain([ordered, 'access', _15 => _15[end], 'optionalAccess', _16 => _16.startBeat]), () => ( 0)), onset)) {
       end += 1;
     }
     const block = ordered.slice(start, end);
     if (block.length > 1) {
       const freeLanes = lanes.map((lane, index) => ({ index, last: lane[lane.length - 1] })).filter(
-        (entry) => entry.last !== void 0 && _chunkSGLEQVJHcjs.hasEnded.call(void 0, entry.last.startBeat + entry.last.durationBeat, onset)
+        (entry) => entry.last !== void 0 && _chunkIFZDWIKMcjs.hasEnded.call(void 0, entry.last.startBeat + entry.last.durationBeat, onset)
       ).sort((a, b) => a.last.pitch - b.last.pitch);
       const assignment = pairWithFreeLanes(block, freeLanes);
       const taken = new Set(
@@ -326,7 +343,7 @@ function prepareTracks(tracks) {
       const sounding = voice.map((note, i) => {
         const prev = i > 0 ? voice[i - 1] : void 0;
         const prevEnd = prev === void 0 ? void 0 : prev.startBeat + prev.durationBeat;
-        const contiguous = prevEnd !== void 0 && _chunkSGLEQVJHcjs.adjacent.call(void 0, prevEnd, note.startBeat);
+        const contiguous = prevEnd !== void 0 && _chunkIFZDWIKMcjs.adjacent.call(void 0, prevEnd, note.startBeat);
         return {
           id: note.id,
           originalIndex: note.originalIndex,
@@ -359,7 +376,7 @@ function withGivenKeys(result, given) {
   if (given === void 0) {
     return result;
   }
-  return { ...result, keys: given, prevailingKey: _nullishCoalesce(_chunkSGLEQVJHcjs.prevailingKeyOf.call(void 0, given), () => ( result.prevailingKey)) };
+  return { ...result, keys: given, prevailingKey: _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, given), () => ( result.prevailingKey)) };
 }
 function callerTimeline(timeline, key, given, pooled, meters, budget) {
   const chartEnd = timeline.segments.reduce((end, segment) => Math.max(end, segment.endBeat), 0);
@@ -368,15 +385,15 @@ function callerTimeline(timeline, key, given, pooled, meters, budget) {
     0
   );
   const totalBeats = Math.max(chartEnd, soundingEnd);
-  const { startBeat: musicStart } = _chunkSGLEQVJHcjs.gridForNotes.call(void 0, pooled, _chunkS5BKZZGMcjs.beatsPerBarAt.call(void 0, 0, meters));
-  const keys = _nullishCoalesce(given, () => ( _chunkSGLEQVJHcjs.attachPivots.call(void 0, 
-    key !== void 0 ? [{ startBeat: musicStart, endBeat: totalBeats, key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, key), confidence: 1 }] : _chunkSGLEQVJHcjs.keyTimelineFromNotes.call(void 0, pooled, { meters, totalBeats, budget }),
+  const { startBeat: musicStart } = _chunkIFZDWIKMcjs.gridForNotes.call(void 0, pooled, _chunkS5BKZZGMcjs.beatsPerBarAt.call(void 0, 0, meters));
+  const keys = _nullishCoalesce(given, () => ( _chunkIFZDWIKMcjs.attachPivots.call(void 0, 
+    key !== void 0 ? [{ startBeat: musicStart, endBeat: totalBeats, key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, key), confidence: 1 }] : _chunkIFZDWIKMcjs.keyTimelineFromNotes.call(void 0, pooled, { meters, totalBeats, budget }),
     timeline.segments
   )));
   return {
     timeline,
     keys,
-    prevailingKey: _nullishCoalesce(_chunkSGLEQVJHcjs.prevailingKeyOf.call(void 0, keys), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0)))),
+    prevailingKey: _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, keys), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0)))),
     segmentConfidence: timeline.segments.map(() => 0)
   };
 }
@@ -393,7 +410,7 @@ function soundingVoicesAt(prepared, beat) {
         continue;
       }
       const snap = { pitch: note.pitch };
-      const previous = _chunkSGLEQVJHcjs.sameInstant.call(void 0, note.startBeat, beat) ? note.prevPitch : note.pitch;
+      const previous = _chunkIFZDWIKMcjs.sameInstant.call(void 0, note.startBeat, beat) ? note.prevPitch : note.pitch;
       if (previous !== void 0) {
         snap.prevPitch = previous;
       }
@@ -437,7 +454,7 @@ function analyzePolyphony(notes, chordAtBeat, key) {
     if (subVoice === void 0) {
       continue;
     }
-    const analyzed = _chunkSGLEQVJHcjs.analyzeVoice.call(void 0, 
+    const analyzed = _chunkIFZDWIKMcjs.analyzeVoice.call(void 0, 
       subVoice.voice,
       chordAtBeat,
       key,
@@ -455,7 +472,7 @@ function analyzePolyphony(notes, chordAtBeat, key) {
     if (analyzed !== void 0) {
       return analyzed;
     }
-    const [only] = _chunkSGLEQVJHcjs.analyzeVoice.call(void 0, [{ ...note, id: index }], chordAtBeat, key);
+    const [only] = _chunkIFZDWIKMcjs.analyzeVoice.call(void 0, [{ ...note, id: index }], chordAtBeat, key);
     return _nullishCoalesce(only, () => ( {
       noteId: index,
       pitch: note.pitch,
@@ -522,7 +539,7 @@ function analyzeArrangementWith(tracks, opts = {}, previous, dirty) {
   let evidence;
   let inferred;
   if (opts.timeline === void 0) {
-    const run = _chunkSGLEQVJHcjs.analyzeTimeline.call(void 0, 
+    const run = _chunkIFZDWIKMcjs.analyzeTimeline.call(void 0, 
       pooled,
       {
         key: opts.key,
@@ -540,9 +557,9 @@ function analyzeArrangementWith(tracks, opts = {}, previous, dirty) {
     inferred = callerTimeline(opts.timeline, opts.key, opts.keys, pooled, meters, budget);
   }
   const { timeline, keys, prevailingKey, segmentConfidence } = inferred;
-  const keyAt = _chunkSGLEQVJHcjs.keyLookup.call(void 0, keys, prevailingKey);
+  const keyAt = _chunkIFZDWIKMcjs.keyLookup.call(void 0, keys, prevailingKey);
   const scaleAt = (beat) => _chunk3R3QVKY5cjs.scaleOf.call(void 0, keyAt(beat));
-  const cadences = _chunkSGLEQVJHcjs.detectCadences.call(void 0, timeline, scaleAt);
+  const cadences = _chunkIFZDWIKMcjs.detectCadences.call(void 0, timeline, scaleAt);
   const soundingCache = /* @__PURE__ */ new Map();
   const trackAnalyses = [];
   const conflicts = [];
@@ -561,7 +578,7 @@ function analyzeArrangementWith(tracks, opts = {}, previous, dirty) {
       if (!subVoice) {
         continue;
       }
-      const analyzed = _chunkSGLEQVJHcjs.analyzeVoice.call(void 0, 
+      const analyzed = _chunkIFZDWIKMcjs.analyzeVoice.call(void 0, 
         subVoice.voice,
         timeline.at,
         scaleAt,
@@ -578,7 +595,7 @@ function analyzeArrangementWith(tracks, opts = {}, previous, dirty) {
           continue;
         }
         for (const beat of evaluationBeats(note, timeline)) {
-          const atOnset = _chunkSGLEQVJHcjs.sameInstant.call(void 0, beat, note.startBeat);
+          const atOnset = _chunkIFZDWIKMcjs.sameInstant.call(void 0, beat, note.startBeat);
           const others = otherVoicesSounding(prepared, t, v, beat, soundingCache);
           const safetyQuery = {
             profile,
@@ -735,7 +752,7 @@ function sessionOf(tracks, opts, analysis, evidence) {
       if (edited === null && !reordered) {
         return sessionOf(next, opts, analysis, evidence);
       }
-      const dirty = evidence === void 0 ? void 0 : edited === null ? { from: 0, to: 0 } : _chunkSGLEQVJHcjs.dirtySlotsFor.call(void 0, evidence, edited.startBeat, edited.endBeat);
+      const dirty = evidence === void 0 ? void 0 : edited === null ? { from: 0, to: 0 } : _chunkIFZDWIKMcjs.dirtySlotsFor.call(void 0, evidence, edited.startBeat, edited.endBeat);
       const run = analyzeArrangementWith(next, opts, evidence, dirty);
       return sessionOf(next, opts, run.analysis, run.evidence);
     }
@@ -788,7 +805,7 @@ function tensionCurve(tracks, opts = {}) {
   const pooled = poolNotes(tracks, harmonyTracks);
   const all = poolNotes(tracks);
   const totalBeats = all.reduce((end, n) => Math.max(end, n.startBeat + n.durationBeat), 0);
-  const { origin: sampleStart, startBeat: musicStart } = _chunkSGLEQVJHcjs.gridForNotes.call(void 0, all, step);
+  const { origin: sampleStart, startBeat: musicStart } = _chunkIFZDWIKMcjs.gridForNotes.call(void 0, all, step);
   const suppliedKeys = opts.timeline === void 0 ? void 0 : _nullishCoalesce(opts.keys, () => ( (opts.key !== void 0 ? [
     {
       // The one region a named key holds starts where the music does,
@@ -799,8 +816,8 @@ function tensionCurve(tracks, opts = {}) {
       key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, opts.key),
       confidence: 1
     }
-  ] : _chunkSGLEQVJHcjs.keyTimelineFromNotes.call(void 0, pooled, { meters, totalBeats, budget }))));
-  const { timeline, keys, prevailingKey } = opts.timeline === void 0 || suppliedKeys === void 0 ? _chunkSGLEQVJHcjs.chordTimelineFromNotes.call(void 0, pooled, {
+  ] : _chunkIFZDWIKMcjs.keyTimelineFromNotes.call(void 0, pooled, { meters, totalBeats, budget }))));
+  const { timeline, keys, prevailingKey } = opts.timeline === void 0 || suppliedKeys === void 0 ? _chunkIFZDWIKMcjs.chordTimelineFromNotes.call(void 0, pooled, {
     key: opts.key,
     meters,
     pickupBeats: opts.pickupBeats,
@@ -809,10 +826,10 @@ function tensionCurve(tracks, opts = {}) {
   }) : {
     timeline: opts.timeline,
     keys: suppliedKeys,
-    prevailingKey: _nullishCoalesce(_chunkSGLEQVJHcjs.prevailingKeyOf.call(void 0, suppliedKeys), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0))))
+    prevailingKey: _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, suppliedKeys), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0))))
   };
   const regions = _nullishCoalesce(opts.keys, () => ( keys));
-  const keyAt = _chunkSGLEQVJHcjs.keyLookup.call(void 0, regions, _nullishCoalesce(_chunkSGLEQVJHcjs.prevailingKeyOf.call(void 0, regions), () => ( prevailingKey)));
+  const keyAt = _chunkIFZDWIKMcjs.keyLookup.call(void 0, regions, _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, regions), () => ( prevailingKey)));
   const prepared = prepareTracks(tracks);
   const points = [];
   const sampleCount = Math.max(0, Math.ceil((totalBeats - sampleStart) / step - _chunkS5BKZZGMcjs.BEAT_EPS));
@@ -870,7 +887,7 @@ function sampleTension(prepared, timeline, key, meters, profile, beat) {
     return 0;
   }
   const chord = timeline.at(beat);
-  const functionScore = chord ? FUNCTION_TENSION[_chunkSGLEQVJHcjs.functionOf.call(void 0, chord, key)] : 0;
+  const functionScore = chord ? FUNCTION_TENSION[_chunkIFZDWIKMcjs.functionOf.call(void 0, chord, key)] : 0;
   const chordPcs = chord ? new Set(_chunk3R3QVKY5cjs.chordPitchClasses.call(void 0, chord)) : null;
   const strongBeat = _chunkS5BKZZGMcjs.isStrongBeat.call(void 0, beat, meters);
   let nonChord = 0;
@@ -950,7 +967,7 @@ function accentContrast(novelty, slices, groupBars, phase) {
       continue;
     }
     const value = _nullishCoalesce(novelty[i], () => ( 0));
-    if (_chunkSGLEQVJHcjs.floorMod.call(void 0, slice.index - phase, groupBars) === 0) {
+    if (_chunkIFZDWIKMcjs.floorMod.call(void 0, slice.index - phase, groupBars) === 0) {
       headSum += value;
       headCount += 1;
     } else {
@@ -968,7 +985,7 @@ function cadenceFitOf(cadenceBars, slices, groupBars, phase) {
   let ends = 0;
   let fit = 0;
   for (const slice of slices) {
-    if (slice.index < 0 || _chunkSGLEQVJHcjs.floorMod.call(void 0, slice.index - phase, groupBars) !== groupBars - 1) {
+    if (slice.index < 0 || _chunkIFZDWIKMcjs.floorMod.call(void 0, slice.index - phase, groupBars) !== groupBars - 1) {
       continue;
     }
     ends += 1;
@@ -978,7 +995,7 @@ function cadenceFitOf(cadenceBars, slices, groupBars, phase) {
       fit += ELIDED_FIT;
     }
   }
-  return ends === 0 ? 0 : _chunkSGLEQVJHcjs.clamp01.call(void 0, fit / ends);
+  return ends === 0 ? 0 : _chunkIFZDWIKMcjs.clamp01.call(void 0, fit / ends);
 }
 function undecidedReading(slices) {
   const fullBars = slices.filter((slice) => slice.index >= 0);
@@ -999,12 +1016,12 @@ function hypermeter(notes, meter, opts = {}) {
     budget: asked.budget
   });
   const sounding = notes.filter((note) => note.durationBeat > 0);
-  const firstOnset = _chunkSGLEQVJHcjs.firstSoundingBeat.call(void 0, sounding);
+  const firstOnset = _chunkIFZDWIKMcjs.firstSoundingBeat.call(void 0, sounding);
   const lastEnd = sounding.reduce((end, n) => Math.max(end, n.startBeat + n.durationBeat), 0);
   const spanEnd = _nullishCoalesce(asked.totalBeats, () => ( lastEnd));
   _chunkS5BKZZGMcjs.assertRange.call(void 0, spanEnd, 0, Number.MAX_SAFE_INTEGER, "hypermeter totalBeats");
-  const spanStart = _chunkSGLEQVJHcjs.barGridStart.call(void 0, firstOnset, meters);
-  const slices = _chunkSGLEQVJHcjs.sliceBars.call(void 0, sounding, meters, spanStart, spanEnd, asked.budget);
+  const spanStart = _chunkIFZDWIKMcjs.barGridStart.call(void 0, firstOnset, meters);
+  const slices = _chunkIFZDWIKMcjs.sliceBars.call(void 0, sounding, meters, spanStart, spanEnd, asked.budget);
   const fullBars = slices.reduce((count, slice) => slice.index >= 0 ? count + 1 : count, 0);
   const candidates = GROUP_CANDIDATES.filter((groupBars) => fullBars >= groupBars * MIN_GROUPS);
   if (candidates.length === 0) {
@@ -1015,7 +1032,7 @@ function hypermeter(notes, meter, opts = {}) {
     "hypermeter readings",
     asked.budget
   );
-  const novelty = slices.map((_, index) => _chunkSGLEQVJHcjs.harmonicNovelty.call(void 0, slices, index));
+  const novelty = slices.map((_, index) => _chunkIFZDWIKMcjs.harmonicNovelty.call(void 0, slices, index));
   const cadenceBars = new Set((_nullishCoalesce(asked.cadenceBeats, () => ( []))).map((beat) => _chunkS5BKZZGMcjs.barIndexAt.call(void 0, beat, meters)));
   const hasCadences = cadenceBars.size > 0;
   const readings = [];
@@ -1030,14 +1047,14 @@ function hypermeter(notes, meter, opts = {}) {
       if (hasCadences && cadenceFit < bestFit - _chunkS5BKZZGMcjs.BEAT_EPS) {
         continue;
       }
-      const evidence = _chunkSGLEQVJHcjs.clamp01.call(void 0, 
+      const evidence = _chunkIFZDWIKMcjs.clamp01.call(void 0, 
         hasCadences ? (
           // Floored at zero rather than summed raw: group heads that restate
           // the harmony their group opened with — the commonest way a phrase
           // begins — change less than the bars between them, and a negative
           // contrast subtracted from the cadence term would count that
           // restatement as evidence against the very phase it confirms.
-          CONTRAST_SHARE * _chunkSGLEQVJHcjs.clamp01.call(void 0, contrast) + CADENCE_SHARE * cadenceFit
+          CONTRAST_SHARE * _chunkIFZDWIKMcjs.clamp01.call(void 0, contrast) + CADENCE_SHARE * cadenceFit
         ) : contrast
       );
       readings.push({
@@ -1058,8 +1075,8 @@ function hypermeter(notes, meter, opts = {}) {
     return undecidedReading(slices);
   }
   const runnerUp = readings.find((reading) => reading.groupBars !== best.groupBars);
-  const margin = runnerUp === void 0 || best.score <= _chunkS5BKZZGMcjs.BEAT_EPS ? 1 : _chunkSGLEQVJHcjs.clamp01.call(void 0, (best.score - runnerUp.score) / best.score);
-  const downbeats = slices.filter((slice) => slice.index >= 0 && _chunkSGLEQVJHcjs.floorMod.call(void 0, slice.index - best.phase, best.groupBars) === 0).map((slice) => slice.startBeat);
+  const margin = runnerUp === void 0 || best.score <= _chunkS5BKZZGMcjs.BEAT_EPS ? 1 : _chunkIFZDWIKMcjs.clamp01.call(void 0, (best.score - runnerUp.score) / best.score);
+  const downbeats = slices.filter((slice) => slice.index >= 0 && _chunkIFZDWIKMcjs.floorMod.call(void 0, slice.index - best.phase, best.groupBars) === 0).map((slice) => slice.startBeat);
   const cadencePart = hasCadences ? `${Math.round(best.cadenceFit * 100)}% of group ends closed by a cadence` : "no cadences given";
   return {
     groupBars: best.groupBars,
@@ -1069,554 +1086,11 @@ function hypermeter(notes, meter, opts = {}) {
     // runner-up is nearly always a multiple or a divisor of the winner — four
     // bars nest inside eight — and a reading of the same accents at another
     // level is weaker counter-evidence than an unrelated rival would be.
-    confidence: _chunkSGLEQVJHcjs.clamp01.call(void 0, best.evidence * (MARGIN_FLOOR + (1 - MARGIN_FLOOR) * margin)),
+    confidence: _chunkIFZDWIKMcjs.clamp01.call(void 0, best.evidence * (MARGIN_FLOOR + (1 - MARGIN_FLOOR) * margin)),
     rationale: `${best.groupBars}-bar groups from bar ${best.phase}: harmonic change ${Math.round(best.contrast * 100)}% stronger on group heads, ${cadencePart}`
   };
 }
 
-// src/analyze/melody/internal.ts
-var RHYTHM_GRID = 64;
-function phraseNotes(phrase) {
-  if (Array.isArray(phrase)) {
-    return phrase;
-  }
-  return _chunkS5BKZZGMcjs.assertRecord.call(void 0, phrase, "melodic phrase").notes;
-}
-function orderedNotes(notes, name, budget) {
-  _chunkS5BKZZGMcjs.assertNoteEvents.call(void 0, notes, name, { allowNonPositiveDuration: true, budget });
-  const sounding = notes.filter((note) => note.durationBeat > 0).sort((a, b) => a.startBeat - b.startBeat || a.pitch - b.pitch);
-  const line = [];
-  for (const note of sounding) {
-    const previous = line[line.length - 1];
-    if (previous !== void 0 && Math.abs(note.startBeat - previous.startBeat) <= _chunkS5BKZZGMcjs.BEAT_EPS) {
-      line[line.length - 1] = note;
-      continue;
-    }
-    line.push(note);
-  }
-  return line;
-}
-function intervalsOf(notes) {
-  const out = [];
-  for (let i = 1; i < notes.length; i += 1) {
-    out.push((_nullishCoalesce(_optionalChain([notes, 'access', _43 => _43[i], 'optionalAccess', _44 => _44.pitch]), () => ( 0))) - (_nullishCoalesce(_optionalChain([notes, 'access', _45 => _45[i - 1], 'optionalAccess', _46 => _46.pitch]), () => ( 0))));
-  }
-  return out;
-}
-function onsetGaps(notes) {
-  const out = [];
-  for (let i = 1; i < notes.length; i += 1) {
-    out.push((_nullishCoalesce(_optionalChain([notes, 'access', _47 => _47[i], 'optionalAccess', _48 => _48.startBeat]), () => ( 0))) - (_nullishCoalesce(_optionalChain([notes, 'access', _49 => _49[i - 1], 'optionalAccess', _50 => _50.startBeat]), () => ( 0))));
-  }
-  return out;
-}
-function retrogradeGaps(notes) {
-  const out = [];
-  for (let i = notes.length - 1; i >= 1; i -= 1) {
-    const later = notes[i];
-    const earlier = notes[i - 1];
-    out.push(
-      (_nullishCoalesce(_optionalChain([later, 'optionalAccess', _51 => _51.startBeat]), () => ( 0))) + (_nullishCoalesce(_optionalChain([later, 'optionalAccess', _52 => _52.durationBeat]), () => ( 0))) - ((_nullishCoalesce(_optionalChain([earlier, 'optionalAccess', _53 => _53.startBeat]), () => ( 0))) + (_nullishCoalesce(_optionalChain([earlier, 'optionalAccess', _54 => _54.durationBeat]), () => ( 0))))
-    );
-  }
-  return out;
-}
-function rhythmProfile(gaps) {
-  const unit = _nullishCoalesce(gaps[0], () => ( 0));
-  if (!(unit > _chunkS5BKZZGMcjs.BEAT_EPS)) {
-    return gaps.length === 0 ? [] : null;
-  }
-  return gaps.map((gap) => Math.round(gap / unit * RHYTHM_GRID) / RHYTHM_GRID);
-}
-function sameNumbers(a, b) {
-  if (a.length !== b.length) {
-    return false;
-  }
-  for (let i = 0; i < a.length; i += 1) {
-    if (Math.abs((_nullishCoalesce(a[i], () => ( 0))) - (_nullishCoalesce(b[i], () => ( 0)))) > _chunkS5BKZZGMcjs.BEAT_EPS) {
-      return false;
-    }
-  }
-  return true;
-}
-function roundTo(value, digits) {
-  const scale = 10 ** digits;
-  return Math.round(value * scale) / scale;
-}
-function endBeatOf(notes) {
-  let end = 0;
-  for (const note of notes) {
-    end = Math.max(end, note.startBeat + note.durationBeat);
-  }
-  return end;
-}
-function onsetSpan(notes) {
-  if (notes.length < 2) {
-    return 0;
-  }
-  return (_nullishCoalesce(_optionalChain([notes, 'access', _55 => _55[notes.length - 1], 'optionalAccess', _56 => _56.startBeat]), () => ( 0))) - (_nullishCoalesce(_optionalChain([notes, 'access', _57 => _57[0], 'optionalAccess', _58 => _58.startBeat]), () => ( 0)));
-}
-
-// src/analyze/melody/contour.ts
-var DIRECTIONAL_SHARE = 0.75;
-function extremes(notes) {
-  if (notes.length === 0) {
-    return { peak: -1, trough: -1 };
-  }
-  let peak = 0;
-  let trough = 0;
-  for (let i = 1; i < notes.length; i += 1) {
-    if ((_nullishCoalesce(_optionalChain([notes, 'access', _59 => _59[i], 'optionalAccess', _60 => _60.pitch]), () => ( 0))) > (_nullishCoalesce(_optionalChain([notes, 'access', _61 => _61[peak], 'optionalAccess', _62 => _62.pitch]), () => ( 0)))) {
-      peak = i;
-    }
-    if ((_nullishCoalesce(_optionalChain([notes, 'access', _63 => _63[i], 'optionalAccess', _64 => _64.pitch]), () => ( 0))) < (_nullishCoalesce(_optionalChain([notes, 'access', _65 => _65[trough], 'optionalAccess', _66 => _66.pitch]), () => ( 0)))) {
-      trough = i;
-    }
-  }
-  return { peak, trough };
-}
-function melodicContour(notes) {
-  const line = orderedNotes(phraseNotes(notes), "contour notes");
-  const intervals = intervalsOf(line);
-  const directions = intervals.map(
-    (step) => step > 0 ? "up" : step < 0 ? "down" : "same"
-  );
-  const { peak, trough } = extremes(line);
-  const range = (_nullishCoalesce(_optionalChain([line, 'access', _67 => _67[peak], 'optionalAccess', _68 => _68.pitch]), () => ( 0))) - (_nullishCoalesce(_optionalChain([line, 'access', _69 => _69[trough], 'optionalAccess', _70 => _70.pitch]), () => ( 0)));
-  const motion = intervals.reduce((sum, step) => sum + Math.abs(step), 0);
-  if (motion === 0) {
-    return {
-      directions,
-      shape: "static",
-      peakIndex: peak,
-      troughIndex: trough,
-      range: 0,
-      rationale: line.length < 2 ? "Too short to have a direction, so the line is static" : "Every note at the same pitch, so the line is static"
-    };
-  }
-  const first = _nullishCoalesce(_optionalChain([line, 'access', _71 => _71[0], 'optionalAccess', _72 => _72.pitch]), () => ( 0));
-  const last = _nullishCoalesce(_optionalChain([line, 'access', _73 => _73[line.length - 1], 'optionalAccess', _74 => _74.pitch]), () => ( 0));
-  const net = last - first;
-  const directness = Math.abs(net) / motion;
-  const share = Math.round(directness * 100);
-  if (directness >= DIRECTIONAL_SHARE) {
-    const shape = net > 0 ? "ascending" : "descending";
-    return {
-      directions,
-      shape,
-      peakIndex: peak,
-      troughIndex: trough,
-      range,
-      rationale: `${share}% of the motion runs ${net > 0 ? "up" : "down"}, so the line is ${shape}`
-    };
-  }
-  const peakPitch = _nullishCoalesce(_optionalChain([line, 'access', _75 => _75[peak], 'optionalAccess', _76 => _76.pitch]), () => ( 0));
-  const archMotion = peakPitch - first + (peakPitch - last);
-  if (peak > 0 && peak < line.length - 1 && archMotion >= DIRECTIONAL_SHARE * motion) {
-    return {
-      directions,
-      shape: "arch",
-      peakIndex: peak,
-      troughIndex: trough,
-      range,
-      rationale: `Rises to its peak at note ${peak + 1} of ${line.length} and falls back, so the line is an arch`
-    };
-  }
-  return {
-    directions,
-    shape: "wave",
-    peakIndex: peak,
-    troughIndex: trough,
-    range,
-    rationale: `Changes direction without settling either way (${share}% of the motion nets out), so the line is a wave`
-  };
-}
-
-// src/analyze/melody/motifs.ts
-var MIN_CELL_NOTES = 3;
-var MAX_CELL_NOTES = 8;
-var MIN_OCCURRENCES = 2;
-function windowSignature(notes) {
-  const profile = rhythmProfile(onsetGaps(notes));
-  if (profile === null) {
-    return null;
-  }
-  return `${notes.length}|${intervalsOf(notes).join(",")}|${profile.join(",")}`;
-}
-function nonOverlapping(starts, length) {
-  const picked = [];
-  let lastEnd = Number.NEGATIVE_INFINITY;
-  for (const start of starts) {
-    if (start > lastEnd) {
-      picked.push(start);
-      lastEnd = start + length - 1;
-    }
-  }
-  return picked;
-}
-function subsumes(outer, inner) {
-  let at = 0;
-  for (const start of inner.starts) {
-    while (at + 1 < outer.starts.length && (_nullishCoalesce(outer.starts[at + 1], () => ( 0))) <= start) {
-      at += 1;
-    }
-    const from = outer.starts[at];
-    if (from === void 0 || start < from || start + inner.length > from + outer.length) {
-      return false;
-    }
-  }
-  return true;
-}
-function listBeats(occurrences) {
-  const shown = occurrences.slice(0, 5).map((occurrence) => roundTo(occurrence.startBeat, 3));
-  return occurrences.length > shown.length ? `${shown.join(", ")}, and ${occurrences.length - shown.length} more` : shown.join(", ");
-}
-function motifFromGroup(group, notes) {
-  const first = _nullishCoalesce(group.starts[0], () => ( 0));
-  const prime = notes.slice(first, first + group.length).map((note) => ({ ...note }));
-  const primePitch = _nullishCoalesce(_optionalChain([prime, 'access', _77 => _77[0], 'optionalAccess', _78 => _78.pitch]), () => ( 0));
-  const primeSpan = onsetSpan(prime);
-  const occurrences = group.starts.map((start) => {
-    const statement = notes.slice(start, start + group.length);
-    const span = onsetSpan(statement);
-    return {
-      noteIndex: start,
-      startBeat: _nullishCoalesce(_optionalChain([statement, 'access', _79 => _79[0], 'optionalAccess', _80 => _80.startBeat]), () => ( 0)),
-      endBeat: endBeatOf(statement),
-      transpose: (_nullishCoalesce(_optionalChain([statement, 'access', _81 => _81[0], 'optionalAccess', _82 => _82.pitch]), () => ( 0))) - primePitch,
-      timeRatio: primeSpan > _chunkS5BKZZGMcjs.BEAT_EPS ? span / primeSpan : 1
-    };
-  });
-  const shape = melodicContour(prime).shape;
-  const transposed = occurrences.some((occurrence) => occurrence.transpose !== 0);
-  const stretched = occurrences.some((occurrence) => Math.abs(occurrence.timeRatio - 1) > _chunkS5BKZZGMcjs.BEAT_EPS);
-  const detail = [
-    transposed ? "at more than one pitch level" : "at the same pitch level",
-    stretched ? "and in more than one set of note values" : ""
-  ].filter((part) => part !== "").join(" ");
-  return {
-    notes: prime,
-    intervals: intervalsOf(prime),
-    rhythm: _nullishCoalesce(rhythmProfile(onsetGaps(prime)), () => ( [])),
-    occurrences,
-    rationale: `${group.length}-note ${shape} cell stated ${occurrences.length} ${occurrences.length === 1 ? "time" : "times"} (beats ${listBeats(occurrences)}), ${detail}`
-  };
-}
-function extractMotifs(notes, opts = {}) {
-  const asked = _chunkS5BKZZGMcjs.assertOptions.call(void 0, opts, "opts");
-  const minNotes = _chunkS5BKZZGMcjs.assertInteger.call(void 0, _nullishCoalesce(asked.minNotes, () => ( MIN_CELL_NOTES)), "motif minNotes", 2, 64);
-  const maxNotes = _chunkS5BKZZGMcjs.assertInteger.call(void 0, _nullishCoalesce(asked.maxNotes, () => ( MAX_CELL_NOTES)), "motif maxNotes", 2, 64);
-  if (maxNotes < minNotes) {
-    throw new (0, _chunkS5BKZZGMcjs.InvalidInputError)(
-      `motif maxNotes must be at least minNotes ${minNotes}; received ${maxNotes}`
-    );
-  }
-  const minOccurrences = _chunkS5BKZZGMcjs.assertInteger.call(void 0, 
-    _nullishCoalesce(asked.minOccurrences, () => ( MIN_OCCURRENCES)),
-    "motif minOccurrences",
-    2
-  );
-  const sounding = orderedNotes(notes, "melody notes", asked.budget);
-  if (sounding.length < minNotes) {
-    return [];
-  }
-  _chunkS5BKZZGMcjs.assertGenerationBudget.call(void 0, 
-    sounding.length * (maxNotes - minNotes + 1),
-    "motif windows",
-    asked.budget
-  );
-  const groups = /* @__PURE__ */ new Map();
-  for (let length = minNotes; length <= maxNotes; length += 1) {
-    for (let start = 0; start + length <= sounding.length; start += 1) {
-      const signature = windowSignature(sounding.slice(start, start + length));
-      if (signature === null) {
-        continue;
-      }
-      const group = groups.get(signature);
-      if (group === void 0) {
-        groups.set(signature, { signature, length, starts: [start] });
-      } else {
-        group.starts.push(start);
-      }
-    }
-  }
-  const candidates = [...groups.values()].map((group) => ({ ...group, starts: nonOverlapping(group.starts, group.length) })).filter((group) => group.starts.length >= minOccurrences).sort(
-    (a, b) => b.length - a.length || b.starts.length - a.starts.length || (_nullishCoalesce(a.starts[0], () => ( 0))) - (_nullishCoalesce(b.starts[0], () => ( 0))) || (a.signature < b.signature ? -1 : a.signature > b.signature ? 1 : 0)
-  );
-  const keptByLength = /* @__PURE__ */ new Map();
-  const kept = [];
-  let comparisons = 0;
-  for (const candidate of candidates) {
-    const first = _nullishCoalesce(candidate.starts[0], () => ( 0));
-    let covered = false;
-    for (let length = candidate.length + 1; length <= maxNotes && !covered; length += 1) {
-      const byStart = keptByLength.get(length);
-      if (byStart === void 0) {
-        continue;
-      }
-      for (let from = Math.max(0, first - (length - candidate.length)); from <= first; from += 1) {
-        comparisons += 1;
-        const outer = byStart.get(from);
-        if (outer === void 0 || outer.starts.length < candidate.starts.length) {
-          continue;
-        }
-        comparisons += outer.starts.length + candidate.starts.length;
-        if (subsumes(outer, candidate)) {
-          covered = true;
-          break;
-        }
-      }
-    }
-    _chunkS5BKZZGMcjs.assertGenerationBudget.call(void 0, comparisons, "motif subsumption comparisons", asked.budget);
-    if (!covered) {
-      kept.push(candidate);
-      let byStart = keptByLength.get(candidate.length);
-      if (byStart === void 0) {
-        byStart = /* @__PURE__ */ new Map();
-        keptByLength.set(candidate.length, byStart);
-      }
-      for (const start of candidate.starts) {
-        byStart.set(start, candidate);
-      }
-    }
-  }
-  return kept.map((group) => motifFromGroup(group, sounding));
-}
-function motifFromNotes(notes) {
-  const sounding = orderedNotes(notes, "motif notes").map((note) => ({ ...note }));
-  return motifFromGroup(
-    { signature: _nullishCoalesce(windowSignature(sounding), () => ( "")), length: sounding.length, starts: [0] },
-    sounding
-  );
-}
-
-// src/analyze/melody/relation.ts
-function spellingOf(key) {
-  if (key === void 0) {
-    return "sharp";
-  }
-  return _chunk3R3QVKY5cjs.keySignatureFifths.call(void 0, key.tonic, key.scale) < 0 ? "flat" : "sharp";
-}
-function intervalBetween(from, to, key) {
-  const spelling = spellingOf(key);
-  return _chunkS5BKZZGMcjs.spelledInterval.call(void 0, _chunkS5BKZZGMcjs.midiToNote.call(void 0, from, spelling), _chunkS5BKZZGMcjs.midiToNote.call(void 0, to, spelling));
-}
-function intervalName(interval) {
-  return `${interval.descending && interval.semitones !== 0 ? "-" : ""}${interval.quality}${interval.number}`;
-}
-function ladderIndices(notes, key) {
-  const root = _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, key.rootPc);
-  const offsets = _chunk3R3QVKY5cjs.scaleTonesInDegreeOrder.call(void 0, key).map((pc) => _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, pc - root));
-  if (offsets.length === 0) {
-    return null;
-  }
-  const out = [];
-  for (const note of notes) {
-    const semitones = note.pitch - root;
-    const octave = Math.floor(semitones / 12);
-    const degree = offsets.indexOf(semitones - octave * 12);
-    if (degree < 0) {
-      return null;
-    }
-    out.push(octave * offsets.length + degree);
-  }
-  return out;
-}
-function ladderSteps(indices) {
-  const out = [];
-  for (let i = 1; i < indices.length; i += 1) {
-    out.push((_nullishCoalesce(indices[i], () => ( 0))) - (_nullishCoalesce(indices[i - 1], () => ( 0))));
-  }
-  return out;
-}
-function pitchShapes(a, b) {
-  const shapes = [];
-  if (sameNumbers(a, b)) {
-    shapes.push("exact");
-  }
-  if (sameNumbers(
-    b,
-    a.map((step) => -step)
-  )) {
-    shapes.push("inversion");
-  }
-  const reversed = [...a].reverse();
-  if (sameNumbers(
-    b,
-    reversed.map((step) => -step)
-  )) {
-    shapes.push("retrograde");
-  }
-  if (sameNumbers(b, reversed)) {
-    shapes.push("retrogradeInversion");
-  }
-  return shapes;
-}
-function totalGap(gaps) {
-  return gaps.reduce((sum, gap) => sum + gap, 0);
-}
-function stretchRatio(modelGaps, answerGaps) {
-  const span = totalGap(modelGaps);
-  return span > _chunkS5BKZZGMcjs.BEAT_EPS ? totalGap(answerGaps) / span : 1;
-}
-function stretchPhrase(ratio) {
-  return `note values ${roundTo(ratio, 3)}x the model's`;
-}
-function relateMotifs(a, b, keyLike) {
-  const key = keyLike === void 0 ? void 0 : _chunk3R3QVKY5cjs.resolveKey.call(void 0, keyLike);
-  const model = _chunkS5BKZZGMcjs.assertArray.call(void 0, _chunkS5BKZZGMcjs.assertRecord.call(void 0, a, "a").notes, "a.notes");
-  const answer = _chunkS5BKZZGMcjs.assertArray.call(void 0, _chunkS5BKZZGMcjs.assertRecord.call(void 0, b, "b").notes, "b.notes");
-  if (model.length === 0 || model.length !== answer.length) {
-    return null;
-  }
-  const modelGaps = onsetGaps(model);
-  const answerGaps = onsetGaps(answer);
-  const modelProfile = rhythmProfile(modelGaps);
-  const answerProfile = rhythmProfile(answerGaps);
-  const retroGaps = retrogradeGaps(model);
-  const retroProfile = rhythmProfile(retroGaps);
-  if (modelProfile === null || answerProfile === null) {
-    return null;
-  }
-  const forward = sameNumbers(modelProfile, answerProfile);
-  const backward = retroProfile !== null && sameNumbers(retroProfile, answerProfile);
-  const semitones = (_nullishCoalesce(_optionalChain([answer, 'access', _83 => _83[0], 'optionalAccess', _84 => _84.pitch]), () => ( 0))) - (_nullishCoalesce(_optionalChain([model, 'access', _85 => _85[0], 'optionalAccess', _86 => _86.pitch]), () => ( 0)));
-  const interval = intervalBetween(_nullishCoalesce(_optionalChain([model, 'access', _87 => _87[0], 'optionalAccess', _88 => _88.pitch]), () => ( 0)), _nullishCoalesce(_optionalChain([answer, 'access', _89 => _89[0], 'optionalAccess', _90 => _90.pitch]), () => ( 0)), key);
-  const gap = (_nullishCoalesce(_optionalChain([answer, 'access', _91 => _91[0], 'optionalAccess', _92 => _92.startBeat]), () => ( 0))) - endBeatOf(model);
-  const sequence = gap >= -_chunkSGLEQVJHcjs.HUMANIZE_ADJACENCY && gap <= _chunkSGLEQVJHcjs.HUMANIZE_ADJACENCY;
-  const timeRatio = stretchRatio(modelGaps, answerGaps);
-  const stretched = Math.abs(timeRatio - 1) > _chunkS5BKZZGMcjs.BEAT_EPS;
-  const shapes = pitchShapes(intervalsOf(model), intervalsOf(answer));
-  const base = { sequence, semitones, interval, timeRatio };
-  const backwardBase = { ...base, timeRatio: stretchRatio(retroGaps, answerGaps) };
-  if (forward && shapes.includes("exact")) {
-    if (semitones === 0) {
-      if (!stretched) {
-        return {
-          ...base,
-          kind: "repetition",
-          rationale: "Exact repetition: the same pitches in the same rhythm"
-        };
-      }
-      const augmented = timeRatio > 1;
-      return {
-        ...base,
-        kind: augmented ? "augmentation" : "diminution",
-        rationale: `${augmented ? "Augmentation" : "Diminution"}: the same pitches in ${stretchPhrase(timeRatio)}`
-      };
-    }
-    const name = intervalName(interval);
-    return {
-      ...base,
-      kind: "transposition",
-      rationale: sequence ? `Real sequence: the cell restated ${name} away with every interval kept literally` + (stretched ? `, in ${stretchPhrase(timeRatio)}` : "") : `Transposition by ${name}: every interval kept literally` + (stretched ? `, in ${stretchPhrase(timeRatio)}` : "")
-    };
-  }
-  if (forward && key !== void 0) {
-    const modelLadder = ladderIndices(model, key.scale);
-    const answerLadder = ladderIndices(answer, key.scale);
-    if (modelLadder !== null && answerLadder !== null) {
-      const degrees = (_nullishCoalesce(answerLadder[0], () => ( 0))) - (_nullishCoalesce(modelLadder[0], () => ( 0)));
-      if (degrees !== 0 && sameNumbers(ladderSteps(modelLadder), ladderSteps(answerLadder))) {
-        const direction = degrees > 0 ? "higher" : "lower";
-        const count = Math.abs(degrees);
-        return {
-          ...base,
-          kind: "tonalTransposition",
-          degrees,
-          rationale: `${sequence ? "Tonal sequence" : "Tonal transposition"}: the cell restated ${count} scale ${count === 1 ? "degree" : "degrees"} ${direction} inside the key, so the key resizes its intervals \u2014 a real one would have kept them literal`
-        };
-      }
-    }
-  }
-  if (forward && shapes.includes("inversion")) {
-    return {
-      ...base,
-      kind: "inversion",
-      rationale: "Inversion: every interval turned upside down about the first note" + (stretched ? `, in ${stretchPhrase(timeRatio)}` : "")
-    };
-  }
-  if (backward && shapes.includes("retrograde")) {
-    return {
-      ...backwardBase,
-      kind: "retrograde",
-      rationale: "Retrograde: the cell read back to front, rhythm included"
-    };
-  }
-  if (backward && shapes.includes("retrogradeInversion")) {
-    return {
-      ...backwardBase,
-      kind: "retrogradeInversion",
-      rationale: "Retrograde inversion: the cell read back to front and turned upside down"
-    };
-  }
-  return null;
-}
-
-// src/analyze/melody/similarity.ts
-var INTERVAL_TOLERANCE = 6;
-var PITCH_WEIGHT = 0.75;
-var RHYTHM_WEIGHT = 0.25;
-function gradedEditDistance(a, b, cost) {
-  let previous = Array.from({ length: b.length + 1 }, (_, i) => i);
-  for (let i = 1; i <= a.length; i += 1) {
-    const row = new Array(b.length + 1);
-    row[0] = i;
-    for (let j = 1; j <= b.length; j += 1) {
-      const substitution = (_nullishCoalesce(previous[j - 1], () => ( 0))) + cost(_nullishCoalesce(a[i - 1], () => ( 0)), _nullishCoalesce(b[j - 1], () => ( 0)));
-      const deletion = (_nullishCoalesce(previous[j], () => ( 0))) + 1;
-      const insertion = (_nullishCoalesce(row[j - 1], () => ( 0))) + 1;
-      row[j] = Math.min(substitution, deletion, insertion);
-    }
-    previous = row;
-  }
-  return _nullishCoalesce(previous[b.length], () => ( 0));
-}
-function distanceToSimilarity(distance, a, b) {
-  const longest = Math.max(a, b);
-  if (longest === 0) {
-    return 1;
-  }
-  return Math.min(1, Math.max(0, 1 - distance / longest));
-}
-function intervalCost(x, y) {
-  return Math.min(1, Math.abs(x - y) / INTERVAL_TOLERANCE);
-}
-function gapCost(x, y) {
-  if (!(x > _chunkS5BKZZGMcjs.BEAT_EPS) || !(y > _chunkS5BKZZGMcjs.BEAT_EPS)) {
-    return Math.abs(x - y) > _chunkS5BKZZGMcjs.BEAT_EPS ? 1 : 0;
-  }
-  return Math.min(1, Math.abs(Math.log2(x / y)));
-}
-function compareMelodies(a, b) {
-  const left = orderedNotes(phraseNotes(a), "melody a");
-  const right = orderedNotes(phraseNotes(b), "melody b");
-  const leftIntervals = intervalsOf(left);
-  const rightIntervals = intervalsOf(right);
-  const pitchSimilarity = distanceToSimilarity(
-    gradedEditDistance(leftIntervals, rightIntervals, intervalCost),
-    leftIntervals.length,
-    rightIntervals.length
-  );
-  const leftGaps = _nullishCoalesce(rhythmProfile(onsetGaps(left)), () => ( onsetGaps(left)));
-  const rightGaps = _nullishCoalesce(rhythmProfile(onsetGaps(right)), () => ( onsetGaps(right)));
-  const rhythmSimilarity = distanceToSimilarity(
-    gradedEditDistance(leftGaps, rightGaps, gapCost),
-    leftGaps.length,
-    rightGaps.length
-  );
-  const similarity = PITCH_WEIGHT * pitchSimilarity + RHYTHM_WEIGHT * rhythmSimilarity;
-  return {
-    similarity,
-    pitchSimilarity,
-    rhythmSimilarity,
-    rationale: `Interval shapes ${Math.round(pitchSimilarity * 100)}% alike, rhythms ${Math.round(rhythmSimilarity * 100)}% alike, weighted ${PITCH_WEIGHT}/${RHYTHM_WEIGHT} into ${roundTo(similarity, 3)}`
-  };
-}
-function melodicSimilarity(a, b) {
-  return compareMelodies(a, b).similarity;
-}
-
 // src/analyze/form/phrase.ts
 var CADENCE_BOUNDARY_STRENGTH = {
   authentic: 0.9,
@@ -1668,7 +1142,7 @@ function readBoundary(boundary) {
     (a, b) => b[1] - a[1] || SIGNAL_ORDER.indexOf(a[0]) - SIGNAL_ORDER.indexOf(b[0])
   );
   return {
-    strength: _chunkSGLEQVJHcjs.combineEvidence.call(void 0, signals.map(([, strength]) => strength)),
+    strength: _chunkIFZDWIKMcjs.combineEvidence.call(void 0, signals.map(([, strength]) => strength)),
     signals: signals.map(([signal]) => signal)
   };
 }
@@ -1712,29 +1186,18 @@ function choosePhrasePath(beats, strengths, minPhraseBeats, expected) {
   return path;
 }
 function arrivalEnd(timeline, atBeat) {
-  return _nullishCoalesce(_optionalChain([_chunkSGLEQVJHcjs.segmentStartingAt.call(void 0, timeline.segments, atBeat), 'optionalAccess', _93 => _93.endBeat]), () => ( atBeat));
-}
-function barSpanOf(startBeat, endBeat, meters) {
-  const firstBar = _chunkS5BKZZGMcjs.barIndexAt.call(void 0, startBeat, meters);
-  const lastBar = _chunkS5BKZZGMcjs.barIndexAt.call(void 0, endBeat, meters);
-  const partOf = (beat, bar) => {
-    const barStart = _chunkS5BKZZGMcjs.barPositionToBeat.call(void 0, { bar, beat: 0 }, meters);
-    const barLength = _chunkS5BKZZGMcjs.beatsPerBarAt.call(void 0, barStart, meters);
-    return barLength > 0 ? (beat - barStart) / barLength : 0;
-  };
-  const bars = lastBar - firstBar + partOf(endBeat, lastBar) - partOf(startBeat, firstBar);
-  return Math.round(bars * 100) / 100;
+  return _nullishCoalesce(_optionalChain([_chunkIFZDWIKMcjs.segmentStartingAt.call(void 0, timeline.segments, atBeat), 'optionalAccess', _43 => _43.endBeat]), () => ( atBeat));
 }
 function keyContextFor(opts, notes, meters, totalBeats) {
   if (opts.key !== void 0) {
     return opts.key;
   }
-  const regions = _chunkSGLEQVJHcjs.keyTimelineFromNotes.call(void 0, notes, {
+  const regions = _chunkIFZDWIKMcjs.keyTimelineFromNotes.call(void 0, notes, {
     meters,
     totalBeats,
     budget: opts.budget
   });
-  const keyAt = _chunkSGLEQVJHcjs.keyLookup.call(void 0, regions, _nullishCoalesce(_chunkSGLEQVJHcjs.prevailingKeyOf.call(void 0, regions), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0)))));
+  const keyAt = _chunkIFZDWIKMcjs.keyLookup.call(void 0, regions, _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, regions), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0)))));
   return (beat) => _chunk3R3QVKY5cjs.scaleOf.call(void 0, keyAt(beat));
 }
 function restBoundaries(boundaries, notes, meters, restBeats) {
@@ -1748,7 +1211,7 @@ function restBoundaries(boundaries, notes, meters, restBeats) {
         boundaries,
         note.startBeat,
         "rest",
-        REST_STRENGTH * _chunkSGLEQVJHcjs.clamp01.call(void 0, bar > 0 ? gap / bar : 1)
+        REST_STRENGTH * _chunkIFZDWIKMcjs.clamp01.call(void 0, bar > 0 ? gap / bar : 1)
       );
     }
     covered = Math.max(covered, note.startBeat + note.durationBeat);
@@ -1758,7 +1221,7 @@ function repetitionBoundaries(boundaries, notes, meters, spanStart, spanEnd, win
   if (!(windowBeats > 0)) {
     return;
   }
-  const slices = _chunkSGLEQVJHcjs.sliceBars.call(void 0, notes, meters, spanStart, spanEnd, budget);
+  const slices = _chunkIFZDWIKMcjs.sliceBars.call(void 0, notes, meters, spanStart, spanEnd, budget);
   _chunkS5BKZZGMcjs.assertGenerationBudget.call(void 0, slices.length * notes.length, "form repetition comparisons", budget);
   const onsetsIn = (from, to) => notes.filter((note) => note.startBeat >= from - _chunkS5BKZZGMcjs.BEAT_EPS && note.startBeat < to - _chunkS5BKZZGMcjs.BEAT_EPS);
   let comparisons = 0;
@@ -1774,7 +1237,7 @@ function repetitionBoundaries(boundaries, notes, meters, spanStart, spanEnd, win
     }
     comparisons += before.length * after.length;
     _chunkS5BKZZGMcjs.assertGenerationBudget.call(void 0, comparisons, "form repetition melody comparisons", budget);
-    const similarity = melodicSimilarity(before, after);
+    const similarity = _chunkIFZDWIKMcjs.melodicSimilarity.call(void 0, before, after);
     if (similarity >= REPETITION_THRESHOLD) {
       addSignal(boundaries, beat, "repetition", REPETITION_STRENGTH * similarity);
     }
@@ -1803,7 +1266,7 @@ function structuralWeightOf(cadence, confidence, aligned, isLast) {
   if (aligned) {
     weight += HYPERMETRIC_CLOSE_BONUS;
   }
-  return _chunkSGLEQVJHcjs.clamp01.call(void 0, weight);
+  return _chunkIFZDWIKMcjs.clamp01.call(void 0, weight);
 }
 function describeCadence(cadence) {
   if (cadence === null || cadence.cadence.type === null) {
@@ -1814,7 +1277,7 @@ function describeCadence(cadence) {
   return `a ${graded} cadence at beat ${cadence.atBeat}`;
 }
 function phrasesFromTimeline(timeline, notes, opts = {}) {
-  _chunkSGLEQVJHcjs.assertChordTimeline.call(void 0, timeline);
+  _chunkIFZDWIKMcjs.assertChordTimeline.call(void 0, timeline);
   const meters = _chunkS5BKZZGMcjs.resolveMeters.call(void 0, opts, "phrase meters");
   _chunkS5BKZZGMcjs.assertNoteEvents.call(void 0, notes, "phrase notes", {
     allowNonPositiveDuration: true,
@@ -1830,12 +1293,12 @@ function phrasesFromTimeline(timeline, notes, opts = {}) {
   const lastSegment = segments[segments.length - 1];
   const firstOnset = sounding.reduce(
     (first, n) => Math.min(first, n.startBeat),
-    _nullishCoalesce(_optionalChain([firstSegment, 'optionalAccess', _94 => _94.startBeat]), () => ( 0))
+    _nullishCoalesce(_optionalChain([firstSegment, 'optionalAccess', _44 => _44.startBeat]), () => ( 0))
   );
-  const spanStart = _chunkSGLEQVJHcjs.barGridStart.call(void 0, firstOnset, meters);
+  const spanStart = _chunkIFZDWIKMcjs.barGridStart.call(void 0, firstOnset, meters);
   const notesEnd = sounding.reduce(
     (end, n) => Math.max(end, n.startBeat + n.durationBeat),
-    _nullishCoalesce(_optionalChain([lastSegment, 'optionalAccess', _95 => _95.endBeat]), () => ( 0))
+    _nullishCoalesce(_optionalChain([lastSegment, 'optionalAccess', _45 => _45.endBeat]), () => ( 0))
   );
   const spanEnd = Math.max(spanStart, _nullishCoalesce(opts.totalBeats, () => ( notesEnd)));
   _chunkS5BKZZGMcjs.assertRange.call(void 0, spanEnd, 0, Number.MAX_SAFE_INTEGER, "phrase totalBeats");
@@ -1850,7 +1313,7 @@ function phrasesFromTimeline(timeline, notes, opts = {}) {
   const minPhraseBeats = _nullishCoalesce(opts.minPhraseBeats, () => ( barBeats));
   _chunkS5BKZZGMcjs.assertRange.call(void 0, minPhraseBeats, Number.MIN_VALUE, Number.MAX_SAFE_INTEGER, "phrase minPhraseBeats");
   const key = keyContextFor(opts, sounding, meters, spanEnd);
-  const cadences = _chunkSGLEQVJHcjs.detectCadences.call(void 0, timeline, key);
+  const cadences = _chunkIFZDWIKMcjs.detectCadences.call(void 0, timeline, key);
   const grouping = _nullishCoalesce(opts.hypermeter, () => ( hypermeter(sounding, meters, {
     totalBeats: spanEnd,
     cadenceBeats: cadences.map((hit) => hit.atBeat),
@@ -1916,16 +1379,16 @@ function phrasesFromTimeline(timeline, notes, opts = {}) {
   let from = 0;
   for (let n = 0; n < path.length; n += 1) {
     const to = _nullishCoalesce(path[n], () => ( stops.length - 1));
-    const startBeat = _nullishCoalesce(_optionalChain([stops, 'access', _96 => _96[from], 'optionalAccess', _97 => _97.beat]), () => ( spanStart));
+    const startBeat = _nullishCoalesce(_optionalChain([stops, 'access', _46 => _46[from], 'optionalAccess', _47 => _47.beat]), () => ( spanStart));
     const stop = stops[to];
-    const endBeat = _nullishCoalesce(_optionalChain([stop, 'optionalAccess', _98 => _98.beat]), () => ( spanEnd));
+    const endBeat = _nullishCoalesce(_optionalChain([stop, 'optionalAccess', _48 => _48.beat]), () => ( spanEnd));
     const reading = _nullishCoalesce(readings[to], () => ( { strength: 0, signals: [] }));
     const length = endBeat - startBeat;
-    const confidence = _chunkSGLEQVJHcjs.clamp01.call(void 0, 
+    const confidence = _chunkIFZDWIKMcjs.clamp01.call(void 0, 
       CLOSING_SHARE * reading.strength + (1 - CLOSING_SHARE) * Math.max(0, lengthFit(length, expected))
     );
     const isLast = n === path.length - 1;
-    const closing = _nullishCoalesce(_optionalChain([stop, 'optionalAccess', _99 => _99.cadence]), () => ( null));
+    const closing = _nullishCoalesce(_optionalChain([stop, 'optionalAccess', _49 => _49.cadence]), () => ( null));
     const cadence = closing !== null && closing.atBeat >= startBeat - _chunkS5BKZZGMcjs.BEAT_EPS && closing.atBeat < endBeat - _chunkS5BKZZGMcjs.BEAT_EPS ? closing : null;
     const structuralWeight = structuralWeightOf(
       cadence,
@@ -1934,7 +1397,7 @@ function phrasesFromTimeline(timeline, notes, opts = {}) {
       isLast
     );
     const signals = cadence === null ? reading.signals.filter((signal) => signal !== "cadence") : reading.signals;
-    const bars = barSpanOf(startBeat, endBeat, meters);
+    const bars = Math.round(_chunkIFZDWIKMcjs.barSpanOf.call(void 0, startBeat, endBeat, meters) * 100) / 100;
     phrases.push({
       startBeat,
       endBeat,
@@ -1986,12 +1449,12 @@ function isSilent(unit) {
   return unit.notes.length === 0;
 }
 function unitSimilarity(a, b) {
-  const harmony = _chunkSGLEQVJHcjs.weightSimilarity.call(void 0, a.weights, b.weights);
+  const harmony = _chunkIFZDWIKMcjs.weightSimilarity.call(void 0, a.weights, b.weights);
   if (a.notes.length < 2 || b.notes.length < 2) {
     return a.notes.length === b.notes.length ? harmony : 0;
   }
-  const melody = melodicSimilarity(a.notes, b.notes);
-  return _chunkSGLEQVJHcjs.clamp01.call(void 0, MELODY_SHARE * melody + HARMONY_SHARE * harmony);
+  const melody = _chunkIFZDWIKMcjs.melodicSimilarity.call(void 0, a.notes, b.notes);
+  return _chunkIFZDWIKMcjs.clamp01.call(void 0, MELODY_SHARE * melody + HARMONY_SHARE * harmony);
 }
 function sectionsFromNotes(notes, opts = {}) {
   const meters = _chunkS5BKZZGMcjs.resolveMeters.call(void 0, opts, "section meters");
@@ -2003,7 +1466,7 @@ function sectionsFromNotes(notes, opts = {}) {
   if (sounding.length === 0) {
     return [];
   }
-  const firstOnset = _chunkSGLEQVJHcjs.barGridStart.call(void 0, _chunkSGLEQVJHcjs.firstSoundingBeat.call(void 0, sounding), meters);
+  const firstOnset = _chunkIFZDWIKMcjs.barGridStart.call(void 0, _chunkIFZDWIKMcjs.firstSoundingBeat.call(void 0, sounding), meters);
   const lastEnd = sounding.reduce((end, n) => Math.max(end, n.startBeat + n.durationBeat), 0);
   const spanEnd = Math.max(firstOnset, _nullishCoalesce(opts.totalBeats, () => ( lastEnd)));
   _chunkS5BKZZGMcjs.assertRange.call(void 0, spanEnd, 0, Number.MAX_SAFE_INTEGER, "section totalBeats");
@@ -2012,7 +1475,7 @@ function sectionsFromNotes(notes, opts = {}) {
   const unitBars = _nullishCoalesce(opts.unitBars, () => ( hypermeter(sounding, meters, { totalBeats: spanEnd, budget: opts.budget }).groupBars));
   _chunkS5BKZZGMcjs.assertPositiveInt.call(void 0, unitBars, "section unitBars");
   const firstBar = _chunkS5BKZZGMcjs.barIndexAt.call(void 0, firstOnset, meters);
-  const lastBar = _chunkSGLEQVJHcjs.lastBarOf.call(void 0, meters, firstOnset, spanEnd);
+  const lastBar = _chunkIFZDWIKMcjs.lastBarOf.call(void 0, meters, firstOnset, spanEnd);
   const startBar = Math.max(0, firstBar);
   const unitCount = Math.max(1, Math.ceil((lastBar - startBar + 1) / unitBars));
   _chunkS5BKZZGMcjs.assertGenerationBudget.call(void 0, unitCount * unitCount, "form section comparisons", opts.budget);
@@ -2029,7 +1492,7 @@ function sectionsFromNotes(notes, opts = {}) {
       startBeat,
       endBeat,
       notes: unitNotes,
-      weights: _chunkSGLEQVJHcjs.windowWeights.call(void 0, unitNotes, startBeat, endBeat, meters).weights
+      weights: _chunkIFZDWIKMcjs.windowWeights.call(void 0, unitNotes, startBeat, endBeat, meters).weights
     });
   }
   let comparisons = 0;
@@ -2093,8 +1556,8 @@ function sectionsFromNotes(notes, opts = {}) {
     while (end < units.length && labels[end] === null) {
       end += 1;
     }
-    const startBeat = _nullishCoalesce(_optionalChain([units, 'access', _100 => _100[index], 'optionalAccess', _101 => _101.startBeat]), () => ( 0));
-    const endBeat = _nullishCoalesce(_optionalChain([units, 'access', _102 => _102[end - 1], 'optionalAccess', _103 => _103.endBeat]), () => ( spanEnd));
+    const startBeat = _nullishCoalesce(_optionalChain([units, 'access', _50 => _50[index], 'optionalAccess', _51 => _51.startBeat]), () => ( 0));
+    const endBeat = _nullishCoalesce(_optionalChain([units, 'access', _52 => _52[end - 1], 'optionalAccess', _53 => _53.endBeat]), () => ( spanEnd));
     const sectionFirstBar = startBar + index * unitBars;
     const sectionLastBar = Math.max(
       sectionFirstBar,
@@ -2105,7 +1568,7 @@ function sectionsFromNotes(notes, opts = {}) {
     if (first === void 0) {
       firstByLabel.set(label, sectionIndex);
     }
-    const similarity = first === void 0 ? 1 : _chunkSGLEQVJHcjs.clamp01.call(void 0, scoreSum / Math.max(1, sounded));
+    const similarity = first === void 0 ? 1 : _chunkIFZDWIKMcjs.clamp01.call(void 0, scoreSum / Math.max(1, sounded));
     sections.push({
       label,
       startBeat,
@@ -2113,163 +1576,284 @@ function sectionsFromNotes(notes, opts = {}) {
       bars: sectionLastBar - sectionFirstBar + 1,
       firstOccurrence: _nullishCoalesce(first, () => ( sectionIndex)),
       similarity,
-      rationale: first === void 0 ? `Section ${label}: material first heard at beat ${startBeat}` : `Section ${label}: ${Math.round(similarity * 100)}% alike the ${label} at beat ${_nullishCoalesce(_optionalChain([sections, 'access', _104 => _104[first], 'optionalAccess', _105 => _105.startBeat]), () => ( 0))}`
+      rationale: first === void 0 ? `Section ${label}: material first heard at beat ${startBeat}` : `Section ${label}: ${Math.round(similarity * 100)}% alike the ${label} at beat ${_nullishCoalesce(_optionalChain([sections, 'access', _54 => _54[first], 'optionalAccess', _55 => _55.startBeat]), () => ( 0))}`
     });
     index = end;
   }
   return sections;
 }
 
-// src/analyze/reduction/index.ts
-var TONIC_ARRIVING = /* @__PURE__ */ new Set(["authentic", "plagal", "modal", "deceptive"]);
-function signedStep(fromPc, toPc) {
-  return (_chunkS5BKZZGMcjs.pitchClassOf.call(void 0, toPc) - _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, fromPc) + 18) % 12 - 6;
+// src/analyze/reference/compare.ts
+var RATIONALE_FIELDS = 2;
+function rationaleOf(comparison) {
+  const measured = [];
+  for (const group of [comparison.form, comparison.harmony, comparison.melody, comparison.rhythm]) {
+    for (const [name, value] of Object.entries(group)) {
+      if (value !== null) {
+        measured.push({ name, value, order: measured.length });
+      }
+    }
+  }
+  if (measured.length === 0) {
+    return "Nothing to compare: neither profile has material for any field";
+  }
+  const describe = (entries) => entries.map((entry) => `${entry.name} ${Math.round(entry.value * 100)}%`).join(", ");
+  const highest = [...measured].sort((x, y) => y.value - x.value || x.order - y.order);
+  const lowest = [...measured].sort((x, y) => x.value - y.value || x.order - y.order);
+  return `Most alike: ${describe(highest.slice(0, RATIONALE_FIELDS))}; least alike: ${describe(lowest.slice(0, RATIONALE_FIELDS))}`;
+}
+function compareReferences(a, b, opts) {
+  _chunkIFZDWIKMcjs.assertReferenceProfile.call(void 0, a, "reference profile a");
+  _chunkIFZDWIKMcjs.assertReferenceProfile.call(void 0, b, "reference profile b");
+  const { budget } = _chunkS5BKZZGMcjs.assertOptions.call(void 0, opts, "compare options");
+  const structural = (profile) => profile.harmony.chords.filter((chord) => chord.level === "structural");
+  const cadences = (profile) => profile.form.phrases.map((phrase) => _nullishCoalesce(_optionalChain([phrase, 'access', _56 => _56.cadence, 'optionalAccess', _57 => _57.type]), () => ( null)));
+  const phraseMelodies = (profile) => profile.form.phrases.map((phrase) => phrase.melody);
+  const comparison = {
+    form: {
+      sectionSequenceSimilarity: _chunkIFZDWIKMcjs.sectionSequenceSimilarity.call(void 0, 
+        a.form.sections.map((section) => section.label),
+        b.form.sections.map((section) => section.label)
+      ),
+      phraseLengthSimilarity: _chunkIFZDWIKMcjs.phraseLengthSimilarity.call(void 0, 
+        a.form.phrases.map((phrase) => phrase.bars),
+        b.form.phrases.map((phrase) => phrase.bars)
+      )
+    },
+    harmony: {
+      progressionSimilarity: _chunkIFZDWIKMcjs.progressionSimilarity.call(void 0, structural(a), structural(b), budget),
+      functionSimilarity: _chunkIFZDWIKMcjs.functionSimilarity.call(void 0, a.harmony.chords, b.harmony.chords),
+      cadenceSimilarity: _chunkIFZDWIKMcjs.cadenceSimilarity.call(void 0, cadences(a), cadences(b)),
+      keyPlanSimilarity: _chunkIFZDWIKMcjs.keyPlanSimilarity.call(void 0, 
+        a.harmony.keys.map((region) => region.key),
+        b.harmony.keys.map((region) => region.key)
+      ),
+      harmonicRhythmOnsetSimilarity: _chunkIFZDWIKMcjs.onsetSimilarity.call(void 0, a.harmony.rhythm, b.harmony.rhythm),
+      harmonicRhythmDensitySimilarity: _chunkIFZDWIKMcjs.densitySimilarity.call(void 0, a.harmony.rhythm, b.harmony.rhythm)
+    },
+    melody: {
+      surfaceSimilarity: _chunkIFZDWIKMcjs.surfaceSimilarity.call(void 0, a.melody.motifs, b.melody.motifs),
+      motifStructureSimilarity: _chunkIFZDWIKMcjs.motifStructureSimilarity.call(void 0, 
+        { graph: a.melody.graph, spanBeats: a.span.endBeat - a.span.startBeat },
+        { graph: b.melody.graph, spanBeats: b.span.endBeat - b.span.startBeat }
+      ),
+      contourSimilarity: _chunkIFZDWIKMcjs.contourSimilarity.call(void 0, phraseMelodies(a), phraseMelodies(b)),
+      registerSimilarity: _chunkIFZDWIKMcjs.registerSimilarity.call(void 0, 
+        { phrases: phraseMelodies(a), mean: _nullishCoalesce(_optionalChain([a, 'access', _58 => _58.melody, 'access', _59 => _59.register, 'optionalAccess', _60 => _60.mean]), () => ( null)) },
+        { phrases: phraseMelodies(b), mean: _nullishCoalesce(_optionalChain([b, 'access', _61 => _61.melody, 'access', _62 => _62.register, 'optionalAccess', _63 => _63.mean]), () => ( null)) }
+      )
+    },
+    rhythm: {
+      onsetSimilarity: _chunkIFZDWIKMcjs.onsetSimilarity.call(void 0, a.melody.rhythm, b.melody.rhythm),
+      densitySimilarity: _chunkIFZDWIKMcjs.densitySimilarity.call(void 0, a.melody.rhythm, b.melody.rhythm),
+      durationSimilarity: _chunkIFZDWIKMcjs.durationSimilarity.call(void 0, a.melody.rhythm, b.melody.rhythm),
+      syncopationSimilarity: _chunkIFZDWIKMcjs.syncopationSimilarity.call(void 0, a.melody.rhythm, b.melody.rhythm)
+    }
+  };
+  return { ...comparison, rationale: rationaleOf(comparison) };
+}
+
+// src/analyze/reference/profile.ts
+function spanOf(sounding, timeline, meters, totalBeats) {
+  const segments = timeline.segments;
+  const firstSegment = segments[0];
+  const lastSegment = segments[segments.length - 1];
+  const firstOnset = _chunkIFZDWIKMcjs.firstSoundingBeat.call(void 0, sounding, _nullishCoalesce(_optionalChain([firstSegment, 'optionalAccess', _64 => _64.startBeat]), () => ( 0)));
+  const startBeat = _chunkIFZDWIKMcjs.barGridStart.call(void 0, 
+    Math.min(firstOnset, _nullishCoalesce(_optionalChain([firstSegment, 'optionalAccess', _65 => _65.startBeat]), () => ( firstOnset))),
+    meters
+  );
+  const notesEnd = sounding.reduce(
+    (end, note) => Math.max(end, note.startBeat + note.durationBeat),
+    _nullishCoalesce(_optionalChain([lastSegment, 'optionalAccess', _66 => _66.endBeat]), () => ( 0))
+  );
+  return { startBeat, endBeat: Math.max(startBeat, _nullishCoalesce(totalBeats, () => ( notesEnd))) };
 }
-function isStep(distance) {
-  const size = Math.abs(distance);
-  return size === 1 || size === 2;
+function sameKey(a, b) {
+  return _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, a.scale.rootPc) === _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, b.scale.rootPc) && a.scale.modeMask12 === b.scale.modeMask12;
 }
-function touches(before, after) {
-  return Math.abs(after.startBeat - before.endBeat) <= _chunkS5BKZZGMcjs.BEAT_EPS;
+function keyIndexAt(regions, beat, prevailing) {
+  const holding = regions.findIndex((region) => region.startBeat <= beat && beat < region.endBeat);
+  if (holding >= 0) {
+    return holding;
+  }
+  return Math.max(
+    0,
+    regions.findIndex((region) => sameKey(region.key, prevailing))
+  );
 }
-function salienceFloors(durations) {
-  const order = durations.map((_, index) => index).sort((a, b) => {
-    const left = _nullishCoalesce(durations[a], () => ( 0));
-    const right = _nullishCoalesce(durations[b], () => ( 0));
-    return left - right || a - b;
-  });
-  const rankOf = new Array(durations.length).fill(0);
-  for (let rank = 0; rank < order.length; rank += 1) {
-    rankOf[_nullishCoalesce(order[rank], () => ( 0))] = rank;
-  }
-  const floors = [];
-  const remaining = durations.length - 1;
-  for (let index = 0; index < durations.length; index += 1) {
-    if (remaining <= 0) {
-      floors.push(0);
-      continue;
-    }
-    const skipped = _nullishCoalesce(rankOf[index], () => ( 0));
-    const at = (position) => _nullishCoalesce(durations[_nullishCoalesce(order[position < skipped ? position : position + 1], () => ( 0))], () => ( 0));
-    const middle = remaining >> 1;
-    floors.push(remaining % 2 === 1 ? at(middle) : (at(middle - 1) + at(middle)) / 2);
-  }
-  return floors;
-}
-function isTonicChord(chord, key) {
-  return _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, chord.rootPc - key.rootPc) === 0 && !_chunkSGLEQVJHcjs.isAppliedDominant.call(void 0, chord, key);
-}
-function functionalReason(segment, before, after, keyAt) {
-  const key = keyAt(segment.startBeat);
-  if (after !== void 0 && _chunkSGLEQVJHcjs.isCadentialSixFour.call(void 0, segment.chord, after.chord, key)) {
-    return { kind: "cadentialSixFour" };
-  }
-  if (isTonicChord(segment.chord, key)) {
-    return { kind: "tonic" };
-  }
-  if (_chunk3R3QVKY5cjs.isDominantChordOf.call(void 0, segment.chord, key)) {
-    return { kind: "dominant" };
-  }
-  if (before !== void 0) {
-    const arrival = _chunkSGLEQVJHcjs.detectCadence.call(void 0, before.chord, segment.chord, key).type;
-    if (arrival !== null) {
-      return { kind: "cadenceArrival", cadence: arrival };
-    }
-  }
-  if (after !== void 0) {
-    const agent = _chunkSGLEQVJHcjs.detectCadence.call(void 0, segment.chord, after.chord, keyAt(after.startBeat)).type;
-    if (agent !== null && TONIC_ARRIVING.has(agent)) {
-      return { kind: "cadenceAgent", cadence: agent };
-    }
-  }
-  return null;
-}
-function figureReason(segment, before, after) {
-  if (before === void 0 || after === void 0) {
-    return null;
-  }
-  const incoming = signedStep(before.chord.rootPc, segment.chord.rootPc);
-  const outgoing = signedStep(segment.chord.rootPc, after.chord.rootPc);
-  if (!isStep(incoming) || !isStep(outgoing)) {
-    return null;
-  }
-  if (incoming > 0 === outgoing > 0) {
-    return { kind: "passing" };
-  }
-  return _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, before.chord.rootPc) === _chunkS5BKZZGMcjs.pitchClassOf.call(void 0, after.chord.rootPc) ? { kind: "neighbor" } : null;
-}
-function describe(reason) {
-  switch (reason.kind) {
-    case "tonic":
-      return "Structural: the tonic of the key, which the progression is heard against";
-    case "dominant":
-      return "Structural: the dominant of the key, the chord its tonic is approached from";
-    case "cadentialSixFour":
-      return "Structural: a tonic six-four on the bass of the dominant that follows, which prolongs that dominant rather than framing the tonic";
-    case "cadenceArrival":
-      return `Structural: the arrival of a ${reason.cadence} cadence`;
-    case "cadenceAgent":
-      return `Structural: the chord that drives a ${reason.cadence} cadence`;
-    case "salient":
-      return "Structural: it holds the harmony longer than the chords around it";
-    case "first":
-      return "Structural: the first chord of the progression frames it";
-    case "last":
-      return "Structural: the last chord of the progression frames it";
-    case "unfigured":
-      return "Structural: its root neither steps through to the next chord nor returns to the last, so it embellishes nothing";
-    case "passing":
-      return "Passing: its root steps in from the previous chord and on to the next in the same direction";
-    case "neighbor":
-      return "Neighbor: its root steps away from the surrounding harmony and returns to it";
-  }
-}
-var BASIS_READINGS = {
-  function: ({ segment, before, after, keyAt }) => functionalReason(segment, before, after, keyAt),
-  duration: ({ segment, salienceFloor }) => segment.endBeat - segment.startBeat > salienceFloor ? { kind: "salient" } : null
-};
-var REDUCTION_BASES = Object.keys(BASIS_READINGS);
-function levelOf(reason) {
-  if (reason.kind === "passing") {
-    return "passing";
+function sectionAt(sections, beat) {
+  const index = sections.findIndex(
+    (section) => beat >= section.startBeat - _chunkS5BKZZGMcjs.BEAT_EPS && beat < section.endBeat - _chunkS5BKZZGMcjs.BEAT_EPS
+  );
+  return index < 0 ? null : index;
+}
+function readingsFor(notes, opts = {}) {
+  const asked = _chunkS5BKZZGMcjs.assertOptions.call(void 0, opts, "reference options");
+  const meters = _chunkS5BKZZGMcjs.resolveMeters.call(void 0, asked, "reference meters");
+  const budget = asked.budget;
+  if (asked.timeline === void 0) {
+    return _chunkIFZDWIKMcjs.chordTimelineFromNotes.call(void 0, notes, {
+      meters,
+      key: asked.key,
+      totalBeats: asked.totalBeats,
+      budget
+    });
   }
-  return reason.kind === "neighbor" ? "neighbor" : "structural";
+  const timeline = _chunkIFZDWIKMcjs.assertChordTimeline.call(void 0, asked.timeline, "reference timeline");
+  let keys;
+  if (asked.key === void 0) {
+    keys = _chunkIFZDWIKMcjs.keyTimelineFromNotes.call(void 0, notes, { meters, totalBeats: asked.totalBeats, budget });
+  } else {
+    const sounding = _chunkS5BKZZGMcjs.assertNoteEvents.call(void 0, notes, "reference notes", {
+      allowNonPositiveDuration: true,
+      budget
+    }).filter((note) => note.durationBeat > 0);
+    const span = spanOf(sounding, timeline, meters, asked.totalBeats);
+    keys = [
+      {
+        startBeat: span.startBeat,
+        endBeat: span.endBeat,
+        key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, asked.key),
+        confidence: 1
+      }
+    ];
+  }
+  return { timeline, keys, prevailingKey: _nullishCoalesce(_chunkIFZDWIKMcjs.prevailingKeyOf.call(void 0, keys), () => ( _chunk3R3QVKY5cjs.resolveKey.call(void 0, _chunk3R3QVKY5cjs.majorKey.call(void 0, 0)))) };
 }
-function reduceProgression(timeline, key, opts = {}) {
-  _chunkSGLEQVJHcjs.assertChordTimeline.call(void 0, timeline);
-  const asked = _chunkS5BKZZGMcjs.assertOptions.call(void 0, opts, "opts");
-  const keyAt = _chunkSGLEQVJHcjs.keyScaleAt.call(void 0, key);
-  const basis = _chunkS5BKZZGMcjs.assertOneOf.call(void 0, _nullishCoalesce(asked.basis, () => ( "function")), REDUCTION_BASES, "reduction basis");
-  const reading = BASIS_READINGS[basis];
-  const segments = timeline.segments;
-  const floors = basis === "duration" ? salienceFloors(segments.map((segment) => segment.endBeat - segment.startBeat)) : [];
-  const result = [];
-  for (let index = 0; index < segments.length; index += 1) {
-    const segment = segments[index];
-    if (segment === void 0) {
-      continue;
+function referenceFromReadings(notes, readings, opts = {}) {
+  const asked = _chunkS5BKZZGMcjs.assertOptions.call(void 0, opts, "reference options");
+  const meters = _chunkS5BKZZGMcjs.resolveMeters.call(void 0, asked, "reference meters");
+  const budget = asked.budget;
+  const sounding = _chunkS5BKZZGMcjs.assertNoteEvents.call(void 0, notes, "reference notes", {
+    allowNonPositiveDuration: true,
+    budget
+  }).filter((note) => note.durationBeat > 0);
+  if (asked.totalBeats !== void 0) {
+    _chunkS5BKZZGMcjs.assertRange.call(void 0, asked.totalBeats, 0, Number.MAX_SAFE_INTEGER, "reference totalBeats");
+  }
+  const line = _chunkIFZDWIKMcjs.orderedNotes.call(void 0, _nullishCoalesce(asked.melody, () => ( notes)), "reference melody", budget);
+  const { timeline, prevailingKey } = readings;
+  const span = spanOf(sounding, timeline, meters, asked.totalBeats);
+  const spanEnd = span.endBeat;
+  const keyAt = _chunkIFZDWIKMcjs.keyLookup.call(void 0, readings.keys, prevailingKey);
+  const scaleAt = (beat) => _chunk3R3QVKY5cjs.scaleOf.call(void 0, keyAt(beat));
+  const cadenceBeats = _chunkIFZDWIKMcjs.detectCadences.call(void 0, timeline, scaleAt).map((hit) => hit.atBeat);
+  const grouping = hypermeter(sounding, meters, { cadenceBeats, totalBeats: spanEnd, budget });
+  const phrases = phrasesFromTimeline(timeline, notes, {
+    meters,
+    key: scaleAt,
+    hypermeter: grouping,
+    totalBeats: spanEnd,
+    budget
+  });
+  const cadences = /* @__PURE__ */ new Map();
+  for (const closing of structuralCadences(phrases)) {
+    const type = closing.cadence.type;
+    if (type !== null) {
+      cadences.set(closing.phraseIndex, {
+        atBeat: closing.atBeat,
+        type,
+        weight: closing.weight
+      });
     }
-    const previous = segments[index - 1];
-    const following = segments[index + 1];
-    const before = previous !== void 0 && touches(previous, segment) ? previous : void 0;
-    const after = following !== void 0 && touches(segment, following) ? following : void 0;
-    let reason = reading({ segment, before, after, keyAt, salienceFloor: _nullishCoalesce(floors[index], () => ( 0)) });
-    if (reason === null && index === 0) {
-      reason = { kind: "first" };
-    }
-    if (reason === null && index === segments.length - 1) {
-      reason = { kind: "last" };
-    }
-    reason ??= _nullishCoalesce(figureReason(segment, before, after), () => ( { kind: "unfigured" }));
-    result.push({
-      chord: segment.chord,
-      startBeat: segment.startBeat,
-      endBeat: segment.endBeat,
-      level: levelOf(reason),
-      rationale: describe(reason)
-    });
   }
-  return result;
+  const sections = sectionsFromNotes(notes, {
+    meters,
+    unitBars: grouping.groupBars,
+    totalBeats: spanEnd,
+    budget
+  });
+  const keys = readings.keys.length === 0 && timeline.segments.length > 0 ? [
+    {
+      startBeat: span.startBeat,
+      endBeat: spanEnd,
+      key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, prevailingKey),
+      confidence: 0
+    }
+  ] : readings.keys.map((region) => ({
+    startBeat: region.startBeat,
+    endBeat: region.endBeat,
+    key: _chunk3R3QVKY5cjs.resolveKey.call(void 0, region.key),
+    confidence: region.confidence
+  }));
+  const chords = _chunkIFZDWIKMcjs.reduceProgression.call(void 0, timeline, scaleAt).map((reduced) => {
+    const index = keyIndexAt(keys, (reduced.startBeat + reduced.endBeat) / 2, prevailingKey);
+    const key = keys[index].key;
+    return {
+      startBeat: reduced.startBeat,
+      endBeat: reduced.endBeat,
+      key: index,
+      roman: _chunkIFZDWIKMcjs.chordToRoman.call(void 0, reduced.chord, key),
+      function: _chunkIFZDWIKMcjs.functionOf.call(void 0, reduced.chord, key),
+      level: reduced.level
+    };
+  });
+  const chordOnsets = timeline.segments.map((segment) => ({
+    pitch: 60,
+    startBeat: segment.startBeat,
+    durationBeat: segment.endBeat - segment.startBeat
+  }));
+  const harmonyRhythm = _chunkIFZDWIKMcjs.analyzeRhythm.call(void 0, chordOnsets, { meters, totalBeats: spanEnd, budget });
+  const motifs = _chunkIFZDWIKMcjs.extractMotifs.call(void 0, line, { budget });
+  const graph = _chunkIFZDWIKMcjs.motifGraph.call(void 0, line, motifs, { key: prevailingKey, budget });
+  let register = null;
+  if (line.length > 0) {
+    const pitches = line.map((note) => note.pitch);
+    register = {
+      low: Math.min(...pitches),
+      high: Math.max(...pitches),
+      mean: _chunkIFZDWIKMcjs.weightedMean.call(void 0, line)
+    };
+  }
+  const lineMean = _nullishCoalesce(_optionalChain([register, 'optionalAccess', _67 => _67.mean]), () => ( 0));
+  const referenceMotifs = motifs.map((motif) => {
+    const first = motif.occurrences[0];
+    return {
+      intervals: motif.intervals,
+      rhythm: motif.rhythm,
+      unitBeats: motif.notes.length < 2 ? 0 : motif.notes[1].startBeat - motif.notes[0].startBeat,
+      spanBeats: first === void 0 ? 0 : first.endBeat - first.startBeat,
+      occurrences: motif.occurrences.length
+    };
+  });
+  const referencePhrases = phrases.map((phrase, index) => ({
+    startBeat: phrase.startBeat,
+    endBeat: phrase.endBeat,
+    bars: _chunkIFZDWIKMcjs.barSpanOf.call(void 0, phrase.startBeat, phrase.endBeat, meters),
+    section: sectionAt(sections, phrase.startBeat),
+    confidence: phrase.confidence,
+    cadence: _nullishCoalesce(cadences.get(index), () => ( null)),
+    melody: _chunkIFZDWIKMcjs.phraseMelody.call(void 0, line, phrase.startBeat, phrase.endBeat, lineMean, meters, budget),
+    motifNodes: graph.nodes.flatMap(
+      (node, nodeIndex) => node.startBeat >= phrase.startBeat - _chunkS5BKZZGMcjs.BEAT_EPS && node.startBeat < phrase.endBeat - _chunkS5BKZZGMcjs.BEAT_EPS ? [nodeIndex] : []
+    )
+  }));
+  return {
+    profileVersion: _chunkIFZDWIKMcjs.REFERENCE_PROFILE_VERSION,
+    span: {
+      startBeat: span.startBeat,
+      endBeat: spanEnd,
+      bars: _chunkIFZDWIKMcjs.barSpanOf.call(void 0, span.startBeat, spanEnd, meters)
+    },
+    meters,
+    form: {
+      hypermeter: { groupBars: grouping.groupBars, confidence: grouping.confidence },
+      sections,
+      phrases: referencePhrases
+    },
+    harmony: { keys, chords, rhythm: harmonyRhythm },
+    melody: {
+      register,
+      shape: line.length < 2 ? null : _chunkIFZDWIKMcjs.melodicContour.call(void 0, line).shape,
+      motifs: referenceMotifs,
+      graph,
+      rhythm: _chunkIFZDWIKMcjs.analyzeRhythm.call(void 0, line, { meters, totalBeats: spanEnd, budget })
+    }
+  };
+}
+function analyzeReference(notes, opts) {
+  return referenceFromReadings(notes, readingsFor(notes, opts), opts);
 }
 
 
@@ -2286,9 +1870,5 @@ function reduceProgression(timeline, key, opts = {}) {
 
 
 
-
-
-
-
-exports.trackRoleOf = trackRoleOf; exports.analyzePolyphony = analyzePolyphony; exports.analyzeArrangement = analyzeArrangement; exports.createArrangementSession = createArrangementSession; exports.tensionCurve = tensionCurve; exports.tensionCurveFrom = tensionCurveFrom; exports.hypermeter = hypermeter; exports.melodicContour = melodicContour; exports.extractMotifs = extractMotifs; exports.motifFromNotes = motifFromNotes; exports.relateMotifs = relateMotifs; exports.compareMelodies = compareMelodies; exports.melodicSimilarity = melodicSimilarity; exports.phrasesFromTimeline = phrasesFromTimeline; exports.structuralCadences = structuralCadences; exports.sectionsFromNotes = sectionsFromNotes; exports.reduceProgression = reduceProgression;
-//# sourceMappingURL=chunk-HPVSH7SF.cjs.map
+exports.trackRoleOf = trackRoleOf; exports.analyzePolyphony = analyzePolyphony; exports.analyzeArrangement = analyzeArrangement; exports.createArrangementSession = createArrangementSession; exports.tensionCurve = tensionCurve; exports.tensionCurveFrom = tensionCurveFrom; exports.hypermeter = hypermeter; exports.phrasesFromTimeline = phrasesFromTimeline; exports.structuralCadences = structuralCadences; exports.sectionsFromNotes = sectionsFromNotes; exports.compareReferences = compareReferences; exports.referenceFromReadings = referenceFromReadings; exports.analyzeReference = analyzeReference;
+//# sourceMappingURL=chunk-YFA7X2F6.cjs.map