@libraz/libcantus 0.9.2 → 0.9.3

package/dist/generate/index.cjs

@@ -42,6 +42,86 @@ __export(generate_exports, {
 });
 module.exports = __toCommonJS(generate_exports);
 
+// src/core/validation/index.ts
+var DEFAULT_GENERATION_BUDGET = 1e6;
+function assertFiniteNumber(value, name) {
+  if (!Number.isFinite(value)) {
+    throw new RangeError(`${name} must be finite; received ${value}`);
+  }
+  return value;
+}
+function assertInteger(value, name, min = Number.MIN_SAFE_INTEGER, max = Number.MAX_SAFE_INTEGER) {
+  assertFiniteNumber(value, name);
+  if (!Number.isSafeInteger(value) || value < min || value > max) {
+    throw new RangeError(`${name} must be an integer in [${min}, ${max}]; received ${value}`);
+  }
+  return value;
+}
+function assertPositiveInt(value, name, max = DEFAULT_GENERATION_BUDGET) {
+  return assertInteger(value, name, 1, max);
+}
+function assertRange(value, min, max, name) {
+  assertFiniteNumber(value, name);
+  if (value < min || value > max) {
+    throw new RangeError(`${name} must be in [${min}, ${max}]; received ${value}`);
+  }
+  return value;
+}
+function assertGenerationBudget(estimated, name, limit = DEFAULT_GENERATION_BUDGET) {
+  assertFiniteNumber(estimated, name);
+  assertPositiveInt(limit, `${name} limit`, Number.MAX_SAFE_INTEGER);
+  if (estimated < 0 || estimated > limit) {
+    throw new RangeError(`${name} exceeds the generation budget ${limit}; received ${estimated}`);
+  }
+  return estimated;
+}
+function assertTimeSignature(ts, name = "time signature") {
+  assertPositiveInt(ts.numerator, `${name}.numerator`);
+  assertPositiveInt(ts.denominator, `${name}.denominator`);
+  if (ts.grouping !== void 0) {
+    if (ts.grouping.length === 0) {
+      throw new RangeError(`${name}.grouping must not be empty`);
+    }
+    let sum = 0;
+    for (let index = 0; index < ts.grouping.length; index += 1) {
+      sum += assertPositiveInt(ts.grouping[index] ?? Number.NaN, `${name}.grouping[${index}]`);
+    }
+    const compound = ts.numerator % 3 === 0 && ts.numerator > 3;
+    const pulses = compound ? ts.numerator / 3 : ts.numerator;
+    if (sum !== pulses) {
+      throw new RangeError(`${name}.grouping must sum to ${pulses}; received ${sum}`);
+    }
+  }
+  return ts;
+}
+function assertNoteEvent(event, name = "note event", options = {}) {
+  assertFiniteNumber(event.pitch, `${name}.pitch`);
+  assertRange(event.startBeat, 0, Number.MAX_SAFE_INTEGER, `${name}.startBeat`);
+  assertRange(
+    event.durationBeat,
+    options.allowNonPositiveDuration ? -Number.MAX_SAFE_INTEGER : Number.MIN_VALUE,
+    Number.MAX_SAFE_INTEGER,
+    `${name}.durationBeat`
+  );
+  if (!options.allowNonPositiveDuration && event.durationBeat <= 0) {
+    throw new RangeError(`${name}.durationBeat must be positive; received ${event.durationBeat}`);
+  }
+  if (event.velocity !== void 0) {
+    assertRange(event.velocity, 0, 127, `${name}.velocity`);
+  }
+  return event;
+}
+function assertNoteEvents(events, name = "note events", options = {}) {
+  assertGenerationBudget(events.length, `${name} count`, options.budget);
+  for (let index = 0; index < events.length; index += 1) {
+    const event = events[index];
+    if (event !== void 0) {
+      assertNoteEvent(event, `${name}[${index}]`, options);
+    }
+  }
+  return events;
+}
+
 // src/core/meter/index.ts
 var EPS = 1e-9;
 function isMultiple(value, unit) {
@@ -52,15 +132,18 @@ function isMultiple(value, unit) {
   return Math.abs(ratio - Math.round(ratio)) < EPS;
 }
 function formatTimeSignature(ts) {
+  assertTimeSignature(ts);
   return `${ts.numerator}/${ts.denominator}`;
 }
 function isCompound(ts) {
+  assertTimeSignature(ts);
   return ts.numerator % 3 === 0 && ts.numerator > 3;
 }
 function unitBeats(ts) {
   return 4 / ts.denominator;
 }
 function beatsPerBar(ts) {
+  assertTimeSignature(ts);
   return ts.numerator * unitBeats(ts);
 }
 function pulseBeats(ts) {
@@ -93,11 +176,14 @@ function isGroupHead(grouping, pulseIndex) {
   return false;
 }
 function beatToBarPosition(beatInQuarters, ts) {
+  assertFiniteNumber(beatInQuarters, "beat");
   const barLen = beatsPerBar(ts);
   const bar = Math.floor(beatInQuarters / barLen);
   return { bar, beat: beatInQuarters - bar * barLen };
 }
 function metricWeight(beatInQuarters, ts) {
+  assertFiniteNumber(beatInQuarters, "beat");
+  assertTimeSignature(ts);
   const { beat } = beatToBarPosition(beatInQuarters, ts);
   const pulse = pulseBeats(ts);
   if (!isMultiple(beat, pulse)) {
@@ -128,6 +214,7 @@ function isStrongBeat(beatInQuarters, ts) {
 
 // src/core/random/index.ts
 function createRng(seed) {
+  assertFiniteNumber(seed, "seed");
   let state = seed >>> 0;
   const next = () => {
     state = state + 1831565813 >>> 0;
@@ -138,9 +225,29 @@ function createRng(seed) {
   };
   return {
     next,
-    prob: (p) => next() < p,
-    range: (lo, hi) => lo + Math.floor(next() * (hi - lo + 1)),
-    float: (lo, hi) => lo + next() * (hi - lo)
+    prob: (p) => next() < assertRange(p, 0, 1, "probability"),
+    range: (lo, hi) => {
+      assertInteger(lo, "range lower bound");
+      assertInteger(hi, "range upper bound");
+      if (lo > hi) {
+        throw new RangeError(
+          `range lower bound must not exceed upper bound; received ${lo} > ${hi}`
+        );
+      }
+      assertInteger(hi - lo, "range span", 0);
+      return lo + Math.floor(next() * (hi - lo + 1));
+    },
+    float: (lo, hi) => {
+      assertFiniteNumber(lo, "float lower bound");
+      assertFiniteNumber(hi, "float upper bound");
+      if (lo > hi) {
+        throw new RangeError(
+          `float lower bound must not exceed upper bound; received ${lo} > ${hi}`
+        );
+      }
+      assertFiniteNumber(hi - lo, "float span");
+      return lo + next() * (hi - lo);
+    }
   };
 }
 
@@ -340,13 +447,30 @@ function buildWalking(ctx, seg, next) {
   });
 }
 function generateBassLine(opts) {
+  assertGenerationBudget(opts.segments.length, "bass segments");
+  for (let index = 0; index < opts.segments.length; index += 1) {
+    const segment = opts.segments[index];
+    if (!segment) continue;
+    assertRange(segment.startBeat, 0, Number.MAX_SAFE_INTEGER, `segments[${index}].startBeat`);
+    assertRange(segment.endBeat, 0, Number.MAX_SAFE_INTEGER, `segments[${index}].endBeat`);
+    if (segment.endBeat <= segment.startBeat) {
+      throw new RangeError(`segments[${index}] must have a positive duration`);
+    }
+  }
   const segments = [...opts.segments].sort((a, b) => a.startBeat - b.startBeat);
   if (segments.length === 0) {
     return [];
   }
   const ts = opts.ts ?? DEFAULT_TS;
+  assertTimeSignature(ts);
   const style = opts.style ?? DEFAULT_STYLE;
   const octave = opts.octave ?? DEFAULT_OCTAVE;
+  assertInteger(octave, "bass octave", -1, 9);
+  const estimatedNotes = segments.reduce(
+    (count, segment) => count + Math.max(1, Math.ceil(segment.endBeat - segment.startBeat)),
+    0
+  );
+  assertGenerationBudget(estimatedNotes, "bass notes");
   const low = octave * 12 + 12;
   const ctx = {
     ts,
@@ -402,6 +526,73 @@ function generateBassLine(opts) {
   return out;
 }
 
+// src/core/event-index/index.ts
+var EPS3 = 1e-9;
+function upperBound(values, beat) {
+  let low = 0;
+  let high = values.length;
+  while (low < high) {
+    const middle = Math.floor((low + high) / 2);
+    if ((values[middle]?.note.startBeat ?? Number.POSITIVE_INFINITY) <= beat) {
+      low = middle + 1;
+    } else {
+      high = middle;
+    }
+  }
+  return low;
+}
+function createNoteEventIndex(events, options = {}) {
+  assertNoteEvents(events, "note events", options);
+  const notes = events.map((note, originalIndex) => ({
+    note,
+    originalIndex,
+    endBeat: note.startBeat + note.durationBeat
+  })).sort((a, b) => a.note.startBeat - b.note.startBeat || a.originalIndex - b.originalIndex);
+  const prefixMaxEnd = [];
+  let maxEnd = Number.NEGATIVE_INFINITY;
+  for (const indexed of notes) {
+    maxEnd = Math.max(maxEnd, indexed.endBeat);
+    prefixMaxEnd.push(maxEnd);
+  }
+  return {
+    notes,
+    at(beat) {
+      let index = upperBound(notes, beat + EPS3) - 1;
+      while (index >= 0) {
+        if ((prefixMaxEnd[index] ?? Number.NEGATIVE_INFINITY) <= beat + EPS3) {
+          return void 0;
+        }
+        const indexed = notes[index];
+        if (indexed && indexed.note.durationBeat > 0 && indexed.note.startBeat - EPS3 <= beat && beat < indexed.endBeat - EPS3) {
+          return indexed;
+        }
+        index -= 1;
+      }
+      return void 0;
+    },
+    attacksAt(beat) {
+      const firstAfter = upperBound(notes, beat + EPS3);
+      const candidate = notes[firstAfter - 1];
+      return candidate !== void 0 && Math.abs(candidate.note.startBeat - beat) < EPS3;
+    },
+    onsetsBetween(startBeat, endBeat) {
+      const result = [];
+      let index = upperBound(notes, startBeat + EPS3);
+      while (index < notes.length) {
+        const onset = notes[index]?.note.startBeat;
+        if (onset === void 0 || onset >= endBeat - EPS3) {
+          break;
+        }
+        if (result[result.length - 1] !== onset) {
+          result.push(onset);
+        }
+        index += 1;
+      }
+      return result;
+    }
+  };
+}
+
 // src/theory/chord/index.ts
 var QUALITY_INTERVALS = {
   maj: [0, 4, 7],
@@ -538,8 +729,10 @@ function classifySeventh(thirdIc, fifthIc, seventhIc) {
 function stackThirds(degree, key, size) {
   const tones = scaleTonesInDegreeOrder(key);
   const length = tones.length;
-  if (length === 0) {
-    return { rootPc: pitchClass3(key.rootPc), quality: "maj", intervals: [0, 4, 7] };
+  if (length !== 7) {
+    throw new Error(
+      `diatonic chord stacking requires a heptatonic scale (received ${length} tones)`
+    );
   }
   const idx = (step) => ((degree + step) % length + length) % length;
   const rootPc = tones[idx(0)] ?? pitchClass3(key.rootPc);
@@ -617,6 +810,31 @@ function createsHiddenParallelPerfect(aPrev, aCur, bPrev, bCur) {
 // src/theory/safety/index.ts
 var SUGGESTION_WINDOW = 12;
 var MAX_SUGGESTIONS = 3;
+var REASON_SEVERITY = [
+  [2 /* Tension */, 1 /* Warning */, 1 /* Warning */],
+  [4 /* AvoidNote */, 1 /* Warning */, 2 /* Dissonant */],
+  [8 /* ScaleTone */, 1 /* Warning */, 1 /* Warning */],
+  [16 /* NonScale */, 2 /* Dissonant */, 2 /* Dissonant */],
+  [32 /* OutOfRange */, 1 /* Warning */, 1 /* Warning */],
+  [64 /* Tritone */, 2 /* Dissonant */, 2 /* Dissonant */],
+  [128 /* LargeLeap */, 1 /* Warning */, 2 /* Dissonant */],
+  [1024 /* VerticalDissonance */, 2 /* Dissonant */, 2 /* Dissonant */],
+  [2048 /* ParallelPerfect */, 1 /* Warning */, 2 /* Dissonant */],
+  [4096 /* HiddenParallel */, 1 /* Warning */, 2 /* Dissonant */],
+  [8192 /* VoiceCrossing */, 1 /* Warning */, 2 /* Dissonant */],
+  [16384 /* Suspension */, 2 /* Dissonant */, 2 /* Dissonant */],
+  [32768 /* NeedsResolution */, 1 /* Warning */, 2 /* Dissonant */]
+];
+function minimumSafetyForReasons(reasons, profile) {
+  let minimum = 0 /* Safe */;
+  const severityIndex = profile === "strict" ? 2 : 1;
+  for (const entry of REASON_SEVERITY) {
+    if (reasons & entry[0]) {
+      minimum = Math.max(minimum, entry[severityIndex]);
+    }
+  }
+  return minimum;
+}
 function pitchClass4(pitch) {
   return (Math.trunc(pitch) % 12 + 12) % 12;
 }
@@ -640,6 +858,28 @@ function stepResolution(pitch, chord) {
   }
   return void 0;
 }
+function evaluateSafety(q) {
+  assertSafetyContext(q);
+  assertFiniteNumber(q.candidatePitch, "candidatePitch");
+  return evaluateInternal(q, true);
+}
+function assertSafetyContext(q) {
+  if (q.prevPitch !== void 0) assertFiniteNumber(q.prevPitch, "prevPitch");
+  if (q.vocalLow !== void 0) assertFiniteNumber(q.vocalLow, "vocalLow");
+  if (q.vocalHigh !== void 0) assertFiniteNumber(q.vocalHigh, "vocalHigh");
+  if (q.vocalLow !== void 0 && q.vocalHigh !== void 0 && q.vocalLow > q.vocalHigh) {
+    throw new RangeError("vocalLow must not exceed vocalHigh");
+  }
+  assertGenerationBudget(q.otherVoices.length, "other voices");
+  for (let index = 0; index < q.otherVoices.length; index += 1) {
+    const voice = q.otherVoices[index];
+    if (!voice) continue;
+    assertFiniteNumber(voice.pitch, `otherVoices[${index}].pitch`);
+    if (voice.prevPitch !== void 0) {
+      assertFiniteNumber(voice.prevPitch, `otherVoices[${index}].prevPitch`);
+    }
+  }
+}
 function evaluateInternal(q, collectSuggestions) {
   let reasons = 0;
   let level = 0 /* Safe */;
@@ -700,9 +940,11 @@ function evaluateInternal(q, collectSuggestions) {
     const pc = Math.abs(pitch - q.prevPitch) % 12;
     if (pc === 6) {
       reasons |= 64 /* Tritone */;
+      raise(2 /* Dissonant */);
     }
     if (isForbiddenMelodicLeap(q.prevPitch, pitch)) {
       reasons |= 128 /* LargeLeap */;
+      raise(parallelSeverity);
     }
     if (pc === 1) {
       reasons |= 256 /* MinorSecond */;
@@ -747,6 +989,7 @@ function evaluateInternal(q, collectSuggestions) {
       }
     }
   }
+  raise(minimumSafetyForReasons(reasons, q.profile));
   const result = { safety: level, reasons, rationale: describe(reasons, q) };
   if (resolveTo !== void 0) {
     result.resolveTo = resolveTo;
@@ -810,11 +1053,15 @@ function describe(reasons, q) {
   return q.chord ? "Placeable" : "No chord context";
 }
 function enumerateSafePitches(q, pitchLow, pitchHigh) {
+  assertSafetyContext(q);
   const chordTones = [];
   const others = [];
-  if (!Number.isFinite(pitchLow) || !Number.isFinite(pitchHigh)) {
-    return [];
+  assertInteger(pitchLow, "pitchLow");
+  assertInteger(pitchHigh, "pitchHigh");
+  if (pitchLow > pitchHigh) {
+    throw new RangeError(`pitchLow must not exceed pitchHigh; received ${pitchLow} > ${pitchHigh}`);
   }
+  assertGenerationBudget(pitchHigh - pitchLow + 1, "safe pitch candidates");
   for (let pitch = pitchHigh; pitch >= pitchLow; pitch -= 1) {
     const result = evaluateInternal({ ...q, candidatePitch: pitch }, false);
     if (result.safety === 2 /* Dissonant */) {
@@ -831,7 +1078,7 @@ function enumerateSafePitches(q, pitchLow, pitchHigh) {
 
 // src/generate/countermelody/index.ts
 var DEFAULT_TS2 = { numerator: 4, denominator: 4 };
-var EPS3 = 1e-9;
+var EPS4 = 1e-9;
 var GRID_STEP = 0.5;
 var MIN_ONSET_GAP = 1;
 var REGISTER_OFFSET = 12;
@@ -842,30 +1089,35 @@ var REINFORCE_PROB = 0.35;
 var DEFAULT_MELODY_VELOCITY = 96;
 var VELOCITY_DROP = 16;
 function isWholeBeat(beat) {
-  return Math.abs(beat - Math.round(beat)) < EPS3;
+  return Math.abs(beat - Math.round(beat)) < EPS4;
 }
 function melodyNoteAt(melody, beat) {
-  let found;
-  for (const note of melody) {
-    if (note.startBeat - EPS3 <= beat && beat < note.startBeat + note.durationBeat - EPS3) {
-      if (!found || note.startBeat > found.startBeat) {
-        found = note;
-      }
-    }
-  }
-  return found;
+  return melody.at(beat)?.note;
 }
 function melodyAttacksAt(melody, beat) {
-  return melody.some((note) => Math.abs(note.startBeat - beat) < EPS3);
+  return melody.attacksAt(beat);
+}
+function melodySnapshotAt(melody, beat) {
+  const current = melodyNoteAt(melody, beat);
+  if (current === void 0) {
+    return void 0;
+  }
+  const attacked = Math.abs(current.startBeat - beat) < EPS4;
+  const previous = attacked ? melodyNoteAt(melody, beat - EPS4 * 2) : current;
+  const snapshot = { pitch: current.pitch };
+  if (previous !== void 0) {
+    snapshot.prevPitch = previous.pitch;
+  }
+  return snapshot;
 }
 function followOnsets(melody) {
-  return [...new Set(melody.map((note) => note.startBeat))].sort((a, b) => a - b);
+  return [...new Set(melody.notes.map(({ note }) => note.startBeat))];
 }
 function complementOnsets(melody, ts, spanStart, spanEnd, rng) {
   const onsets = [];
   const gridStart = Math.floor(spanStart / GRID_STEP) * GRID_STEP;
   let last = Number.NEGATIVE_INFINITY;
-  for (let beat = gridStart; beat < spanEnd - EPS3; beat += GRID_STEP) {
+  for (let beat = gridStart; beat < spanEnd - EPS4; beat += GRID_STEP) {
     const attacked = melodyAttacksAt(melody, beat);
     const sounding = melodyNoteAt(melody, beat) !== void 0;
     let place = false;
@@ -876,7 +1128,7 @@ function complementOnsets(melody, ts, spanStart, spanEnd, rng) {
     } else if (isStrongBeat(beat, ts)) {
       place = rng.prob(REINFORCE_PROB);
     }
-    if (place && beat - last >= MIN_ONSET_GAP - EPS3) {
+    if (place && beat - last >= MIN_ONSET_GAP - EPS4) {
       onsets.push(beat);
       last = beat;
     }
@@ -922,54 +1174,102 @@ function scoreCandidate(pitch, melPitch, melPrev, prevPitch, chord, center) {
   score -= Math.abs(pitch - center) * 0.02;
   return score;
 }
-function fallbackPitch(melPitch, prevPitch, center, register, low, high, key) {
-  let target = prevPitch ?? Math.round(center);
-  if (melPitch !== void 0) {
-    if (register === "below" && target >= melPitch) {
-      target = melPitch - 5;
-    } else if (register === "above" && target <= melPitch) {
-      target = melPitch + 5;
-    }
-  }
-  target = Math.min(Math.max(target, low), high);
-  let pitch = nearestScaleTone(target, key);
-  if (melPitch !== void 0) {
-    while (register === "below" && pitch >= melPitch && pitch - 12 >= low) {
-      pitch -= 12;
-    }
-    while (register === "above" && pitch <= melPitch && pitch + 12 <= high) {
-      pitch += 12;
+function heldNoteBoundaries(melody, startBeat, endBeat, chordChangeBeats) {
+  const boundaries = /* @__PURE__ */ new Set([startBeat]);
+  for (const onset of melody.onsetsBetween(startBeat, endBeat)) {
+    boundaries.add(onset);
+  }
+  for (const beat of chordChangeBeats) {
+    if (beat > startBeat + EPS4 && beat < endBeat - EPS4) {
+      boundaries.add(beat);
+    }
+  }
+  const firstGrid = Math.floor(startBeat / GRID_STEP + 1) * GRID_STEP;
+  for (let beat = firstGrid; beat < endBeat - EPS4; beat += GRID_STEP) {
+    boundaries.add(beat);
+  }
+  return [...boundaries].sort((a, b) => a - b);
+}
+function heldPitchSafety(pitch, startBeat, endBeat, prevPitch, melody, opts, profile, register, low, high, ts) {
+  let worst = 0 /* Safe */;
+  const boundaries = heldNoteBoundaries(melody, startBeat, endBeat, opts.chordChangeBeats ?? []);
+  for (const boundary of boundaries) {
+    const melodyVoice = melodySnapshotAt(melody, boundary);
+    if (melodyVoice !== void 0 && (register === "below" && pitch >= melodyVoice.pitch || register === "above" && pitch <= melodyVoice.pitch)) {
+      return null;
+    }
+    const atCounterOnset = Math.abs(boundary - startBeat) < EPS4;
+    const result = evaluateSafety({
+      profile,
+      candidatePitch: pitch,
+      prevPitch: atCounterOnset ? prevPitch : pitch,
+      chord: opts.chordAt(boundary),
+      key: opts.key,
+      otherVoices: melodyVoice ? [melodyVoice] : [],
+      strongBeat: isStrongBeat(boundary, ts),
+      vocalLow: low,
+      vocalHigh: high
+    });
+    worst = Math.max(worst, result.safety);
+    if (worst === 2 /* Dissonant */) {
+      return worst;
     }
   }
-  return Math.min(Math.max(pitch, low), high);
+  return worst;
 }
 function generateCounterMelody(opts) {
-  const melody = [...opts.melody].sort((a, b) => a.startBeat - b.startBeat);
-  if (melody.length === 0) {
+  const melody = createNoteEventIndex(opts.melody);
+  if (melody.notes.length === 0) {
     return [];
   }
   const ts = opts.ts ?? DEFAULT_TS2;
+  assertTimeSignature(ts);
   const register = opts.register ?? "below";
   const rhythm = opts.rhythm ?? "complement";
   const profile = opts.profile ?? "pop";
   const rng = createRng(opts.seed ?? 0);
-  const spanStart = melody[0]?.startBeat ?? 0;
-  const spanEnd = melody.reduce((end, n) => Math.max(end, n.startBeat + n.durationBeat), spanStart);
-  const meanPitch = melody.reduce((sum, n) => sum + n.pitch, 0) / melody.length;
+  const spanStart = melody.notes[0]?.note.startBeat ?? 0;
+  const spanEnd = melody.notes.reduce((end, n) => Math.max(end, n.endBeat), spanStart);
+  assertGenerationBudget(Math.ceil((spanEnd - spanStart) / GRID_STEP), "countermelody grid");
+  const meanPitch = melody.notes.reduce((sum, n) => sum + n.note.pitch, 0) / melody.notes.length;
   const defaultCenter = register === "below" ? meanPitch - REGISTER_OFFSET : meanPitch + REGISTER_OFFSET;
   const low = opts.pitchLow ?? Math.round(defaultCenter) - DEFAULT_HALF_RANGE;
   const high = opts.pitchHigh ?? Math.round(defaultCenter) + DEFAULT_HALF_RANGE;
+  assertInteger(low, "countermelody pitchLow");
+  assertInteger(high, "countermelody pitchHigh");
+  if (low > high) {
+    throw new RangeError(
+      `countermelody pitchLow must not exceed pitchHigh; received ${low} > ${high}`
+    );
+  }
+  assertGenerationBudget(high - low + 1, "countermelody pitch candidates");
+  assertGenerationBudget(opts.chordChangeBeats?.length ?? 0, "chord change beats");
+  for (let index = 0; index < (opts.chordChangeBeats?.length ?? 0); index += 1) {
+    assertRange(
+      opts.chordChangeBeats?.[index] ?? Number.NaN,
+      0,
+      Number.MAX_SAFE_INTEGER,
+      `chordChangeBeats[${index}]`
+    );
+  }
   const center = (low + high) / 2;
   const onsets = rhythm === "follow" ? followOnsets(melody) : complementOnsets(melody, ts, spanStart, spanEnd, rng);
+  assertGenerationBudget(onsets.length * (high - low + 1), "countermelody search");
   const out = [];
   let prevPitch;
-  let prevOnset;
-  for (const beat of onsets) {
+  for (let onsetIndex = 0; onsetIndex < onsets.length; onsetIndex += 1) {
+    const beat = onsets[onsetIndex];
+    if (beat === void 0) {
+      continue;
+    }
+    const nextOnset = onsets[onsetIndex + 1];
+    const endBeat = nextOnset ?? Math.max(spanEnd, beat + GRID_STEP);
     const melNote = melodyNoteAt(melody, beat);
     const melPitch = melNote?.pitch;
-    const melPrev = prevOnset !== void 0 ? melodyNoteAt(melody, prevOnset)?.pitch : void 0;
+    const melodyVoice = melodySnapshotAt(melody, beat);
+    const melPrev = melodyVoice?.prevPitch;
     const chord = opts.chordAt(beat);
-    const otherVoices = melPitch !== void 0 ? [{ pitch: melPitch, prevPitch: melPrev }] : [];
+    const otherVoices = melodyVoice !== void 0 ? [melodyVoice] : [];
     const candidates = enumerateSafePitches(
       {
         profile,
@@ -985,15 +1285,26 @@ function generateCounterMelody(opts) {
       high
     );
     let bestPitch;
+    let bestWorstSafety = 2 /* Dissonant */;
     let bestScore = Number.NEGATIVE_INFINITY;
     for (const candidate of candidates) {
+      const worstSafety = heldPitchSafety(
+        candidate,
+        beat,
+        endBeat,
+        prevPitch,
+        melody,
+        opts,
+        profile,
+        register,
+        low,
+        high,
+        ts
+      );
+      if (worstSafety === null || worstSafety === 2 /* Dissonant */) {
+        continue;
+      }
       if (melPitch !== void 0) {
-        if (register === "below" && candidate >= melPitch) {
-          continue;
-        }
-        if (register === "above" && candidate <= melPitch) {
-          continue;
-        }
         if (prevPitch !== void 0 && melPrev !== void 0 && // Parallel octaves are the perfect-class-zero case of this predicate,
         // so a single check covers both parallel fifths and octaves.
         createsParallelPerfect(prevPitch, candidate, melPrev, melPitch)) {
@@ -1001,29 +1312,24 @@ function generateCounterMelody(opts) {
         }
       }
       const score = scoreCandidate(candidate, melPitch, melPrev, prevPitch, chord, center);
-      const tie = Math.abs(score - bestScore) <= EPS3 && candidate < (bestPitch ?? Number.POSITIVE_INFINITY);
-      if (score > bestScore + EPS3 || tie) {
+      const tie = Math.abs(score - bestScore) <= EPS4 && candidate < (bestPitch ?? Number.POSITIVE_INFINITY);
+      if (worstSafety < bestWorstSafety || worstSafety === bestWorstSafety && (score > bestScore + EPS4 || tie)) {
+        bestWorstSafety = worstSafety;
         bestScore = score;
         bestPitch = candidate;
       }
     }
-    const pitch = bestPitch ?? fallbackPitch(melPitch, prevPitch, center, register, low, high, opts.key);
+    if (bestPitch === void 0) {
+      prevPitch = void 0;
+      continue;
+    }
+    const pitch = bestPitch;
     const velocity = Math.max(
       1,
       Math.min(127, Math.round((melNote?.velocity ?? DEFAULT_MELODY_VELOCITY) - VELOCITY_DROP))
     );
-    out.push({ pitch, startBeat: beat, durationBeat: GRID_STEP, velocity });
+    out.push({ pitch, startBeat: beat, durationBeat: endBeat - beat, velocity });
     prevPitch = pitch;
-    prevOnset = beat;
-  }
-  for (let i = 0; i < out.length; i += 1) {
-    const note = out[i];
-    if (!note) {
-      continue;
-    }
-    const next = out[i + 1];
-    const end = next ? next.startBeat : Math.max(spanEnd, note.startBeat + GRID_STEP);
-    note.durationBeat = end - note.startBeat;
   }
   return out;
 }
@@ -1743,9 +2049,10 @@ function generateHiHatForBeat(ctx, sec) {
 
 // src/generate/drums/euclid.ts
 function euclideanRhythm(pulses, steps, rotation = 0) {
-  if (steps < 1) {
-    return [];
-  }
+  assertPositiveInt(steps, "euclidean steps");
+  assertInteger(pulses, "euclidean pulses");
+  assertInteger(rotation, "euclidean rotation");
+  assertGenerationBudget(steps, "euclidean pattern steps");
   const p = Math.max(0, Math.min(pulses, steps));
   const base = bjorklund(p, steps);
   const shift = (rotation % steps + steps) % steps;
@@ -1754,18 +2061,6 @@ function euclideanRhythm(pulses, steps, rotation = 0) {
   }
   return base.map((_, i) => base[(i - shift + steps) % steps] ?? false);
 }
-function hasHit(mask, step) {
-  return (mask >> step & 1) === 1;
-}
-function patternToMask(pattern) {
-  let mask = 0;
-  for (let i = 0; i < pattern.length; i += 1) {
-    if (pattern[i]) {
-      mask |= 1 << i;
-    }
-  }
-  return mask;
-}
 function bjorklund(pulses, steps) {
   if (pulses === 0) {
     return Array.from({ length: steps }, () => false);
@@ -2097,19 +2392,6 @@ function isInPreChorusLift(section, bar, sectionBars, nextIsChorus) {
   }
   return bar >= sectionBars - 2;
 }
-function euclideanToKickPattern(mask) {
-  const slot = (step) => hasHit(mask, step) || hasHit(mask, (step + 15) % 16);
-  return {
-    beat1: slot(0),
-    beat1and: slot(2),
-    beat2: slot(4),
-    beat2and: slot(6),
-    beat3: slot(8),
-    beat3and: slot(10),
-    beat4: slot(12),
-    beat4and: slot(14)
-  };
-}
 function getKickPattern(section, style, bar, rng) {
   const p = { ...EMPTY };
   if (section === "outro") {
@@ -2231,7 +2513,7 @@ function getPercussionConfig(mood, section) {
   };
 }
 function generateAuxPercussionForBar(track, barStart, config, role, densityMult, rng, bpm) {
-  if (role === "minimal" || role === "fxOnly") {
+  if (role === "minimal") {
     return;
   }
   if (config.tambourine) {
@@ -2280,6 +2562,16 @@ function generateAuxPercussionForBar(track, barStart, config, role, densityMult,
 
 // src/generate/drums/index.ts
 function generateDrums(opts) {
+  assertPositiveInt(opts.bars, "drum bars");
+  assertRange(opts.bpm, Number.MIN_VALUE, 1e3, "drum bpm");
+  assertRange(opts.density, 0, 1, "drum density");
+  assertGenerationBudget(opts.bars * 128, "drum hits");
+  if (opts.euclideanKick) {
+    const steps = opts.euclideanKick.steps ?? 16;
+    assertInteger(steps, "euclidean steps", 1, 16);
+    assertInteger(opts.euclideanKick.pulses, "euclidean pulses");
+    assertInteger(opts.euclideanKick.rotation ?? 0, "euclidean rotation");
+  }
   const track = new HitList();
   const rng = createRng(opts.seed ?? 0);
   const mapping = mapStyle(opts.style);
@@ -2316,15 +2608,8 @@ function generateDrums(opts) {
   const percMood = percMoodCategory(style);
   const nextSection = opts.nextSection ? mapSection(opts.nextSection) : section;
   const nextEnergy = sectionEnergy(nextSection);
-  const euclidKick = opts.euclideanKick ? euclideanToKickPattern(
-    patternToMask(
-      euclideanRhythm(
-        opts.euclideanKick.pulses,
-        opts.euclideanKick.steps ?? 16,
-        opts.euclideanKick.rotation ?? 0
-      )
-    )
-  ) : void 0;
+  const euclidSteps = opts.euclideanKick?.steps ?? 16;
+  const euclidKick = opts.euclideanKick ? euclideanRhythm(opts.euclideanKick.pulses, euclidSteps, opts.euclideanKick.rotation ?? 0) : void 0;
   const reuseSectionKick = (section === "b" || section === "chorus") && style !== "sparse";
   let sectionKick;
   for (let bar = 0; bar < opts.bars; bar += 1) {
@@ -2340,12 +2625,10 @@ function generateDrums(opts) {
       barHasOpenHh = !track.hasCrashNear(barStart + openHhBeatIndex);
     }
     let kick;
-    if (euclidKick) {
-      kick = euclidKick;
-    } else if (reuseSectionKick) {
+    if (!euclidKick && reuseSectionKick) {
       sectionKick ??= getKickPattern(section, style, 0, rng);
       kick = sectionKick;
-    } else {
+    } else if (!euclidKick) {
       kick = getKickPattern(section, style, bar, rng);
     }
     const inLift = isInPreChorusLift(section, bar, opts.bars, section === "b");
@@ -2386,7 +2669,23 @@ function generateDrums(opts) {
         if (inLift) {
           generatePreChorusBuildup(ctx, sec, isLastBar);
         }
-        generateKickForBeat(ctx, sec, kick);
+        if (euclidKick && !inLift) {
+          const stepLength = 4 / euclidSteps;
+          for (let step = 0; step < euclidKick.length; step += 1) {
+            const onset = step * stepLength;
+            if (euclidKick[step] && Math.floor(onset) === beat) {
+              const wholeBeat = Math.abs(onset - Math.round(onset)) < 1e-9;
+              track.add(
+                GM.BD,
+                barStart + onset,
+                Math.min(0.5, stepLength),
+                velocity * (wholeBeat ? 1 : 0.85)
+              );
+            }
+          }
+        } else if (kick) {
+          generateKickForBeat(ctx, sec, kick);
+        }
         generateSnareForBeat(ctx, sec, section === "intro" && bar === 0);
         if (sec.useGhostNotes && !inLift) {
           generateGhostNotesForBeat(ctx, sec);
@@ -2403,10 +2702,19 @@ function generateDrums(opts) {
         }
       }
     }
+    let percussionConfig = getPercussionConfig(percMood, section);
+    if (role === "fxOnly" && !percussionConfig.tambourine && !percussionConfig.shaker && !percussionConfig.handclap) {
+      percussionConfig = {
+        tambourine: false,
+        shaker: true,
+        handclap: false,
+        shaker16th: false
+      };
+    }
     generateAuxPercussionForBar(
       track,
       barStart,
-      getPercussionConfig(percMood, section),
+      percussionConfig,
       role,
       densityMult,
       rng,
@@ -2427,10 +2735,16 @@ var MIN_VELOCITY = 1;
 var MAX_VELOCITY = 127;
 function humanize(events, opts = {}) {
   const ts = opts.ts ?? DEFAULT_TS3;
+  assertTimeSignature(ts);
+  assertNoteEvents(events, "humanize events");
   const timing = opts.timing ?? DEFAULT_TIMING;
   const velocityJitter = opts.velocity ?? DEFAULT_VELOCITY_JITTER;
   const accent = opts.accent ?? DEFAULT_ACCENT;
   const baseVelocity = opts.baseVelocity ?? DEFAULT_BASE_VELOCITY;
+  assertRange(timing, 0, Number.MAX_SAFE_INTEGER, "humanize timing");
+  assertRange(velocityJitter, 0, 127, "humanize velocity jitter");
+  assertRange(accent, 0, 127, "humanize accent");
+  assertRange(baseVelocity, 0, 127, "humanize base velocity");
   const rng = createRng(opts.seed ?? 0);
   return events.map((event) => {
     const timingOffset = rng.float(-timing, timing);
@@ -2463,8 +2777,12 @@ function quantizeToGrid(beat, barBeats, subdivision, slotsPerBar) {
   };
 }
 function extractGrooveTemplate(events, ts, subdivision) {
+  assertTimeSignature(ts);
+  assertNoteEvents(events, "groove source events");
+  assertPositiveInt(subdivision, "groove subdivision");
   const barBeats = beatsPerBar(ts);
   const slotsPerBar = Math.round(barBeats * subdivision);
+  assertGenerationBudget(slotsPerBar, "groove template slots");
   const offsetSums = new Array(slotsPerBar).fill(0);
   const velocitySums = new Array(slotsPerBar).fill(0);
   const counts = new Array(slotsPerBar).fill(0);
@@ -2499,6 +2817,14 @@ function extractGrooveTemplate(events, ts, subdivision) {
   return { subdivision, slotsPerBar, slots, ts };
 }
 function applyGrooveTemplate(events, template, ts) {
+  assertTimeSignature(ts);
+  assertNoteEvents(events, "groove target events");
+  assertPositiveInt(template.subdivision, "template subdivision");
+  assertPositiveInt(template.slotsPerBar, "template slotsPerBar");
+  assertGenerationBudget(template.slots.length, "template slots");
+  if (template.slots.length !== template.slotsPerBar) {
+    throw new RangeError("template slots length must equal slotsPerBar");
+  }
   if (template.ts && (template.ts.numerator !== ts.numerator || template.ts.denominator !== ts.denominator)) {
     throw new Error(
       `Groove template meter ${formatTimeSignature(template.ts)} does not match apply-time meter ${formatTimeSignature(ts)}`
@@ -2524,107 +2850,291 @@ function applyGrooveTemplate(events, template, ts) {
   });
 }
 
-// src/theory/harmony/index.ts
-function roleOf(pitch, chord, chordId = 0) {
-  const interval = ((Math.trunc(pitch) - chord.rootPc) % 12 + 12) % 12;
-  const tones = new Set(chord.intervals.map((i) => (i % 12 + 12) % 12));
-  const hasHigherSeventh = tones.has(10) || tones.has(11);
-  const hasThird = tones.has(3) || tones.has(4);
-  const isSuspendedTone = !hasThird && (interval === 5 && tones.has(5) || interval === 2 && tones.has(2));
-  let role;
-  let lock;
-  if (interval === 0) {
-    role = "root";
-    lock = "identity";
-  } else if (isSuspendedTone) {
-    role = "third";
-    lock = "quality";
-  } else if (interval === 3 || interval === 4) {
-    role = "third";
-    lock = "quality";
-  } else if (interval === 6 || interval === 7 || interval === 8) {
-    role = "fifth";
-    lock = "voicing";
-  } else if (interval === 9 && tones.has(3) && tones.has(6) && !hasHigherSeventh) {
-    role = "seventh";
-    lock = "voicing";
-  } else if (interval === 9 && tones.has(9) && !hasHigherSeventh) {
-    role = "sixth";
-    lock = "voicing";
-  } else if (interval === 10 || interval === 11) {
-    role = "seventh";
-    lock = "voicing";
-  } else {
-    role = "tension";
-    lock = "voicing";
-  }
-  return { role, lock, belongsToChordId: chordId };
+// src/analyze/functional/internal.ts
+function mod12(n) {
+  return (n % 12 + 12) % 12;
 }
-
-// src/generate/harmonize/index.ts
-var FALLBACK = { rootPc: 0, quality: "maj", secondaryDominant: false, base: 0 };
-var DEFAULT_METER = { numerator: 4, denominator: 4 };
-var COMFORT_LOW = 55;
-var COMFORT_HIGH = 79;
-var TESSITURA_WEIGHT = 1e-3;
-var TIE_BREAK_JITTER = 1e-6;
-function pitchClass5(pitch) {
-  return (Math.trunc(pitch) % 12 + 12) % 12;
+function isNeapolitan(chord, key) {
+  return mod12(chord.rootPc - key.rootPc) === 1 && chord.quality === "maj";
 }
-function mulberry32(seed) {
-  let state = seed >>> 0;
-  return () => {
-    state = state + 1831565813 >>> 0;
-    let t = state;
-    t = Math.imul(t ^ t >>> 15, t | 1);
-    t ^= t + Math.imul(t ^ t >>> 7, t | 61);
-    return ((t ^ t >>> 14) >>> 0) / 4294967296;
-  };
+
+// src/analyze/functional/roman.ts
+var ROMAN = ["I", "II", "III", "IV", "V", "VI", "VII"];
+var OFFSET_SPELLING = [
+  [1, ""],
+  [2, "b"],
+  [2, ""],
+  [3, "b"],
+  [3, ""],
+  [4, ""],
+  [4, "#"],
+  [5, ""],
+  [6, "b"],
+  [6, ""],
+  [7, "b"],
+  [7, ""]
+];
+var ROMAN_STYLE = {
+  maj: { lower: false, suffix: "" },
+  min: { lower: true, suffix: "" },
+  dim: { lower: true, suffix: "o" },
+  aug: { lower: false, suffix: "+" },
+  maj7: { lower: false, suffix: "maj7" },
+  min7: { lower: true, suffix: "7" },
+  dom7: { lower: false, suffix: "7" },
+  dim7: { lower: true, suffix: "o7" },
+  m7b5: { lower: true, suffix: "\xF87" },
+  minMaj7: { lower: true, suffix: "maj7" },
+  aug7: { lower: false, suffix: "+7" },
+  augMaj7: { lower: false, suffix: "+maj7" },
+  majb5: { lower: false, suffix: "b5" },
+  "6": { lower: false, suffix: "add6" },
+  min6: { lower: true, suffix: "add6" },
+  "6/9": { lower: false, suffix: "69" },
+  sus2: { lower: false, suffix: "sus2" },
+  sus4: { lower: false, suffix: "sus4" },
+  add9: { lower: false, suffix: "add9" },
+  add11: { lower: false, suffix: "add11" },
+  maj9: { lower: false, suffix: "maj9" },
+  min9: { lower: true, suffix: "9" },
+  dom9: { lower: false, suffix: "9" },
+  "7b9": { lower: false, suffix: "7b9" },
+  "7#9": { lower: false, suffix: "7#9" },
+  "7#11": { lower: false, suffix: "7#11" },
+  "7b13": { lower: false, suffix: "7b13" },
+  "11": { lower: false, suffix: "11" },
+  "13": { lower: false, suffix: "13" },
+  "5": { lower: false, suffix: "5" }
+};
+var SUFFIX_QUALITY = new Map(
+  Object.entries(ROMAN_STYLE).map(
+    ([quality, style]) => [`${style.lower ? "l" : "u"}:${style.suffix}`, quality]
+  )
+);
+function romanStyle(quality) {
+  return ROMAN_STYLE[quality];
 }
-function inferKey(melody) {
-  let best = majorKey(0);
-  let bestScore = Number.NEGATIVE_INFINITY;
-  const candidates = [];
-  for (let tonic = 0; tonic < 12; tonic += 1) {
-    candidates.push(majorKey(tonic), minorKey(tonic));
+function romanSpelling(rootPc, key) {
+  const tones = scaleTonesInDegreeOrder(key);
+  const diatonic = tones.indexOf(mod12(rootPc));
+  if (diatonic >= 0) {
+    return { degreeNumber: diatonic + 1, accidental: "" };
   }
-  for (const key of candidates) {
-    const tonic = pitchClass5(key.rootPc);
-    const isMinor = key.modeMask12 === NATURAL_MINOR_MASK;
-    const leadingTone = (tonic + 11) % 12;
-    let score = 0;
-    for (const n of melody) {
-      const w = Math.max(0.25, n.durationBeat);
-      const pc = pitchClass5(n.pitch);
-      if (isScaleTone(n.pitch, key) || isMinor && pc === leadingTone) {
-        score += w;
-      }
-      if (pc === tonic) {
-        score += 0.5 * w;
-      }
+  for (let i = 0; i < tones.length; i += 1) {
+    const degreeNumber2 = i + 1;
+    if (degreeNumber2 === 1 || degreeNumber2 === 5) {
+      continue;
     }
-    if (score > bestScore) {
-      bestScore = score;
-      best = key;
+    if (mod12((tones[i] ?? 0) - 1) === mod12(rootPc)) {
+      return { degreeNumber: degreeNumber2, accidental: "b" };
     }
   }
-  return best;
+  for (let i = 0; i < tones.length; i += 1) {
+    if (mod12((tones[i] ?? 0) + 1) === mod12(rootPc)) {
+      return { degreeNumber: i + 1, accidental: "#" };
+    }
+  }
+  const [degreeNumber, accidental] = OFFSET_SPELLING[mod12(rootPc - key.rootPc)] ?? [1, ""];
+  return { degreeNumber, accidental };
 }
-function buildCandidates(key, reharmonize) {
-  const tones = scaleTonesInDegreeOrder(key);
-  const candidates = tones.map((rootPc, degree) => ({
-    rootPc,
-    quality: diatonicTriad(degree, key).quality,
-    degree,
-    secondaryDominant: false,
-    base: 0
-  }));
-  const tonic = tones[0] ?? 0;
-  if (reharmonize !== "diatonic") {
-    for (const target of [1, 3, 4, 5]) {
-      const targetRoot = tones[target] ?? 0;
-      candidates.push({
-        rootPc: (targetRoot + 7) % 12,
+function baseMarker(quality) {
+  if (quality === "dim" || quality === "dim7") {
+    return "o";
+  }
+  if (quality === "m7b5") {
+    return "\xF8";
+  }
+  if (quality === "aug" || quality === "aug7") {
+    return "+";
+  }
+  if (quality === "augMaj7") {
+    return "+maj7";
+  }
+  if (quality === "maj7" || quality === "minMaj7") {
+    return "maj7";
+  }
+  return "";
+}
+var TRIAD_FIGURES = { 1: "6", 2: "64" };
+var SEVENTH_FIGURES = { 1: "65", 2: "43", 3: "42" };
+var TRIAD_FIGURE_QUALITIES = /* @__PURE__ */ new Set(["maj", "min", "dim", "aug"]);
+var SEVENTH_FIGURE_QUALITIES = /* @__PURE__ */ new Set([
+  "maj7",
+  "min7",
+  "dom7",
+  "dim7",
+  "m7b5",
+  "minMaj7",
+  "aug7",
+  "augMaj7"
+]);
+function chordToRoman(chord, key) {
+  const { degreeNumber, accidental } = romanSpelling(chord.rootPc, key);
+  const { lower, suffix } = romanStyle(chord.quality);
+  const numeral = ROMAN[degreeNumber - 1] ?? "I";
+  const cased = lower ? numeral.toLowerCase() : numeral;
+  let inversion = 0;
+  if (chord.bassPc !== void 0) {
+    const idx = chord.intervals.findIndex((iv) => mod12(chord.rootPc + iv) === chord.bassPc);
+    if (idx > 0) {
+      inversion = idx;
+    }
+  }
+  if (inversion > 0) {
+    const figure = SEVENTH_FIGURE_QUALITIES.has(chord.quality) ? SEVENTH_FIGURES[inversion] : TRIAD_FIGURE_QUALITIES.has(chord.quality) ? TRIAD_FIGURES[inversion] : void 0;
+    if (figure !== void 0) {
+      return `${accidental}${cased}${baseMarker(chord.quality)}${figure}`;
+    }
+  }
+  return `${accidental}${cased}${suffix}`;
+}
+
+// src/analyze/functional/function.ts
+var FUNCTION_BY_OFFSET = [
+  "tonic",
+  // 0  I
+  "subdominant",
+  // 1  bII (Neapolitan)
+  "subdominant",
+  // 2  ii
+  "tonic",
+  // 3  bIII
+  "tonic",
+  // 4  iii
+  "subdominant",
+  // 5  IV
+  "dominant",
+  // 6  #IV / bV
+  "dominant",
+  // 7  V
+  "subdominant",
+  // 8  bVI
+  "tonic",
+  // 9  vi
+  "subdominant",
+  // 10 bVII
+  "dominant"
+  // 11 vii
+];
+function isMinorKey(key) {
+  const hasMinorThird = key.modeMask12 >> 3 & 1;
+  const hasMajorThird = key.modeMask12 >> 4 & 1;
+  return Boolean(hasMinorThird) && !hasMajorThird;
+}
+function functionOf(chord, key) {
+  return FUNCTION_BY_OFFSET[mod12(chord.rootPc - key.rootPc)] ?? "tonic";
+}
+function isDiatonic(chord, key) {
+  return chordPitchClasses(chord).every((pc) => isScaleTone(pc, key));
+}
+function parallelKey(key) {
+  return {
+    rootPc: mod12(key.rootPc),
+    modeMask12: isMinorKey(key) ? MAJOR_MASK : NATURAL_MINOR_MASK
+  };
+}
+
+// src/analyze/functional/borrowed.ts
+function isHarmonicMinorAlteration(chord, key) {
+  return isMinorKey(key) && isDiatonic(chord, { rootPc: mod12(key.rootPc), modeMask12: HARMONIC_MINOR_MASK });
+}
+function isBorrowedChord(chord, key) {
+  if (isDiatonic(chord, key)) {
+    return false;
+  }
+  if (isHarmonicMinorAlteration(chord, key)) {
+    return false;
+  }
+  return isDiatonic(chord, parallelKey(key));
+}
+function borrowedSource(chord, key) {
+  if (isDiatonic(chord, key)) {
+    return null;
+  }
+  if (isNeapolitan(chord, key)) {
+    return "neapolitan";
+  }
+  if (isBorrowedChord(chord, key)) {
+    return isMinorKey(key) ? "parallel-major" : "parallel-minor";
+  }
+  return null;
+}
+
+// src/theory/harmony/index.ts
+function roleOf(pitch, chord, chordId = 0) {
+  const interval = ((Math.trunc(pitch) - chord.rootPc) % 12 + 12) % 12;
+  const tones = new Set(chord.intervals.map((i) => (i % 12 + 12) % 12));
+  const hasThird = tones.has(3) || tones.has(4);
+  const isSuspendedTone = !hasThird && (interval === 5 && tones.has(5) || interval === 2 && tones.has(2));
+  const chordRole = chordToneRole(pitch, chord);
+  let role;
+  let lock;
+  if (isSuspendedTone) {
+    role = "third";
+    lock = "quality";
+  } else if (chordRole !== null) {
+    role = chordRole;
+    lock = chordRole === "root" ? "identity" : chordRole === "third" ? "quality" : "voicing";
+  } else {
+    role = "tension";
+    lock = "voicing";
+  }
+  return { role, lock, belongsToChordId: chordId };
+}
+
+// src/generate/harmonize/index.ts
+var FALLBACK = { rootPc: 0, quality: "maj", secondaryDominant: false, base: 0 };
+var DEFAULT_METER = { numerator: 4, denominator: 4 };
+var COMFORT_LOW = 55;
+var COMFORT_HIGH = 79;
+var TESSITURA_WEIGHT = 1e-3;
+var TIE_BREAK_JITTER = 1e-6;
+function pitchClass5(pitch) {
+  return (Math.trunc(pitch) % 12 + 12) % 12;
+}
+function inferKey(melody) {
+  let best = majorKey(0);
+  let bestScore = Number.NEGATIVE_INFINITY;
+  const candidates = [];
+  for (let tonic = 0; tonic < 12; tonic += 1) {
+    candidates.push(majorKey(tonic), minorKey(tonic));
+  }
+  for (const key of candidates) {
+    const tonic = pitchClass5(key.rootPc);
+    const isMinor = key.modeMask12 === NATURAL_MINOR_MASK;
+    const leadingTone = (tonic + 11) % 12;
+    let score = 0;
+    for (const n of melody) {
+      const w = Math.max(0.25, n.durationBeat);
+      const pc = pitchClass5(n.pitch);
+      if (isScaleTone(n.pitch, key) || isMinor && pc === leadingTone) {
+        score += w;
+      }
+      if (pc === tonic) {
+        score += 0.5 * w;
+      }
+    }
+    if (score > bestScore) {
+      bestScore = score;
+      best = key;
+    }
+  }
+  return best;
+}
+function buildCandidates(key, reharmonize) {
+  const tones = scaleTonesInDegreeOrder(key);
+  const candidates = tones.map((rootPc, degree) => ({
+    rootPc,
+    quality: diatonicTriad(degree, key).quality,
+    degree,
+    secondaryDominant: false,
+    base: 0
+  }));
+  if (reharmonize !== "diatonic") {
+    for (const target of [1, 3, 4, 5]) {
+      const targetRoot = tones[target] ?? 0;
+      candidates.push({
+        rootPc: (targetRoot + 7) % 12,
         quality: "dom7",
         secondaryDominant: true,
         targetDegree: target,
@@ -2633,26 +3143,29 @@ function buildCandidates(key, reharmonize) {
     }
   }
   if (reharmonize === "borrowed") {
-    candidates.push({
-      rootPc: (tonic + 10) % 12,
-      quality: "maj",
-      secondaryDominant: false,
-      base: 2
-    });
-    candidates.push({
-      rootPc: (tonic + 8) % 12,
-      quality: "maj",
-      secondaryDominant: false,
-      base: 2
-    });
-    candidates.push({
-      rootPc: (tonic + 5) % 12,
-      quality: "min",
-      secondaryDominant: false,
-      base: 2
-    });
+    const parallel = parallelKey(key);
+    const parallelTones = scaleTonesInDegreeOrder(parallel);
+    for (let degree = 0; degree < parallelTones.length; degree += 1) {
+      const chord = diatonicTriad(degree, parallel);
+      if (!isDiatonic(chord, key)) {
+        candidates.push({
+          rootPc: chord.rootPc,
+          quality: chord.quality,
+          secondaryDominant: false,
+          base: 1.25
+        });
+      }
+    }
   }
-  return candidates;
+  const seen = /* @__PURE__ */ new Set();
+  return candidates.filter((candidate) => {
+    const id = `${candidate.rootPc}:${candidate.quality}`;
+    if (seen.has(id)) {
+      return false;
+    }
+    seen.add(id);
+    return true;
+  });
 }
 function emissionCost(seg, cand, melody, key) {
   const pcs = chordPitchClasses(makeChord(cand.rootPc, cand.quality));
@@ -2756,20 +3269,37 @@ function harmonizeOnce(melody, key, candidates, segments, jitter) {
   return { cost: bestCost, path };
 }
 function harmonizeMelody(opts) {
+  const noteIndex = createNoteEventIndex(opts.melody);
   const key = opts.key === "infer" ? inferKey(opts.melody) : opts.key;
   const candidates = buildCandidates(key, opts.reharmonize);
   const candAt = (i) => candidates[i] ?? FALLBACK;
-  const rng = mulberry32(opts.seed ?? 0);
-  const jitter = candidates.map(() => rng() * TIE_BREAK_JITTER);
-  const melodyEnd = opts.melody.reduce((m, n) => Math.max(m, n.startBeat + n.durationBeat), 0);
-  const hr = Math.max(0.25, opts.harmonicRhythm);
+  const rng = createRng(opts.seed ?? 0);
+  const jitter = candidates.map(() => rng.next() * TIE_BREAK_JITTER);
+  const melodyEnd = noteIndex.notes.reduce((m, n) => Math.max(m, n.endBeat), 0);
+  const requestedHr = assertRange(
+    opts.harmonicRhythm,
+    Number.MIN_VALUE,
+    Number.MAX_SAFE_INTEGER,
+    "harmonic rhythm"
+  );
+  const hr = Math.max(0.25, requestedHr);
   const segCount = Math.max(1, Math.ceil(melodyEnd / hr));
-  const segments = [];
-  for (let s = 0; s < segCount; s += 1) {
-    const startBeat = s * hr;
-    const endBeat = startBeat + hr;
-    const noteIndices = opts.melody.map((n, i) => ({ n, i })).filter(({ n }) => n.startBeat < endBeat && n.startBeat + n.durationBeat > startBeat).map(({ i }) => i);
-    segments.push({ startBeat, endBeat, noteIndices });
+  assertGenerationBudget(segCount, "harmonic segments");
+  assertGenerationBudget(segCount * candidates.length * candidates.length, "harmonization search");
+  const segments = Array.from({ length: segCount }, (_, s) => ({
+    startBeat: s * hr,
+    endBeat: (s + 1) * hr,
+    noteIndices: []
+  }));
+  let memberships = 0;
+  for (const indexed of noteIndex.notes) {
+    const first = Math.max(0, Math.floor(indexed.note.startBeat / hr));
+    const lastExclusive = Math.min(segCount, Math.ceil(indexed.endBeat / hr - Number.EPSILON));
+    memberships += Math.max(0, lastExclusive - first);
+    assertGenerationBudget(memberships, "note-to-segment memberships");
+    for (let segment = first; segment < lastExclusive; segment += 1) {
+      segments[segment]?.noteIndices.push(indexed.originalIndex);
+    }
   }
   const transposes = [0];
   if (opts.placement.transposeSearch) {
@@ -2816,11 +3346,11 @@ function harmonizeMelody(opts) {
     }
     return chord;
   });
-  const melodyRoles = bestMelody.map((note, noteIndex) => {
+  const melodyRoles = bestMelody.map((note, noteIndex2) => {
     const segIdx = Math.min(segments.length - 1, Math.max(0, Math.floor(note.startBeat / hr)));
     const cand = candAt(bestPath[segIdx] ?? 0);
     const chord = makeChord(cand.rootPc, cand.quality);
-    return { noteIndex, role: roleOf(note.pitch, chord).role };
+    return { noteIndex: noteIndex2, role: roleOf(note.pitch, chord).role };
   });
   return { transposeSemitones: bestTs, key, chords, melodyRoles };
 }
@@ -2867,29 +3397,30 @@ function nearestChordTone(pitch, chord) {
   return pitch;
 }
 function cellOrigin(notes) {
-  return notes.length === 0 ? 0 : Math.min(...notes.map((n) => n.startBeat));
+  let origin = Number.POSITIVE_INFINITY;
+  for (const note of notes) {
+    origin = Math.min(origin, note.startBeat);
+  }
+  return Number.isFinite(origin) ? origin : 0;
+}
+function cellEnd(notes) {
+  let end = Number.NEGATIVE_INFINITY;
+  for (const note of notes) {
+    end = Math.max(end, note.startBeat + note.durationBeat);
+  }
+  return Number.isFinite(end) ? end : 0;
 }
 function cellSpan(cell) {
   if (cell.notes.length === 0) {
     return 0;
   }
   const start = cellOrigin(cell.notes);
-  const end = Math.max(...cell.notes.map((n) => n.startBeat + n.durationBeat));
+  const end = cellEnd(cell.notes);
   return end - start;
 }
 function clone(cell) {
   return { notes: cell.notes.map((n) => ({ ...n })) };
 }
-function mulberry322(seed) {
-  let state = seed >>> 0;
-  return () => {
-    state = state + 1831565813 >>> 0;
-    let t = state;
-    t = Math.imul(t ^ t >>> 15, t | 1);
-    t ^= t + Math.imul(t ^ t >>> 7, t | 61);
-    return ((t ^ t >>> 14) >>> 0) / 4294967296;
-  };
-}
 function contourOffsets(contour, count) {
   const offsets = [];
   for (let i = 0; i < count; i += 1) {
@@ -2912,19 +3443,21 @@ function contourOffsets(contour, count) {
 }
 function generateMotif(opts) {
   const contour = opts.contour ?? "arch";
-  const bars = Math.max(1, Math.trunc(opts.bars));
+  const bars = assertPositiveInt(opts.bars, "motif bars");
   const totalBeats = bars * 4;
   const noteCount = Math.max(3, bars * 2);
+  assertGenerationBudget(noteCount, "motif notes");
   const beatsPerNote = totalBeats / noteCount;
-  const rng = mulberry322(opts.seed ?? 0);
-  const jitterProb = Math.min(1, Math.max(0, opts.jitter ?? 0));
+  const rng = createRng(opts.seed ?? 0);
+  const requestedJitter = assertFiniteNumber(opts.jitter ?? 0, "motif jitter");
+  const jitterProb = Math.min(1, Math.max(0, requestedJitter));
   const tonic = pitchClass6(opts.key.rootPc) + 60;
   const offsets = contourOffsets(contour, noteCount);
   const notes = [];
   for (let i = 0; i < noteCount; i += 1) {
     let jitter = 0;
-    if (jitterProb > 0 && rng() < jitterProb) {
-      jitter = rng() < 0.5 ? 1 : -1;
+    if (jitterProb > 0 && rng.next() < jitterProb) {
+      jitter = rng.next() < 0.5 ? 1 : -1;
     }
     let pitch = stepDiatonic(tonic, (offsets[i] ?? 0) + jitter, opts.key);
     const startBeat = i * beatsPerNote;
@@ -2936,6 +3469,16 @@ function generateMotif(opts) {
   return { notes };
 }
 function transformMotif(cell, t, amount, key) {
+  assertNoteEvents(cell.notes, "motif notes");
+  if (amount !== void 0) {
+    assertFiniteNumber(amount, "motif transform amount");
+  }
+  if ((t === "transposeDiatonic" || t === "sequence") && amount !== void 0) {
+    assertInteger(amount, "diatonic transform amount");
+  }
+  if ((t === "augment" || t === "diminish") && amount !== void 0 && amount <= 0) {
+    throw new RangeError("time transform amount must be positive");
+  }
   const notes = cell.notes;
   switch (t) {
     case "transposeChromatic": {
@@ -2962,7 +3505,7 @@ function transformMotif(cell, t, amount, key) {
         return { notes: [] };
       }
       const start = cellOrigin(notes);
-      const end = Math.max(...notes.map((n) => n.startBeat + n.durationBeat));
+      const end = cellEnd(notes);
       return {
         notes: notes.map((n) => ({
           pitch: n.pitch,
@@ -2995,8 +3538,10 @@ function scaleTime(cell, factor) {
   };
 }
 function developMotif(cell, timeline, key, bars) {
+  assertPositiveInt(bars, "development bars");
+  assertNoteEvents(cell.notes, "motif notes");
   const span = cellSpan(cell);
-  const totalBeats = Math.max(1, Math.trunc(bars)) * 4;
+  const totalBeats = bars * 4;
   const origin = cellOrigin(cell.notes);
   const out = [];
   if (span <= 0) {
@@ -3004,6 +3549,7 @@ function developMotif(cell, timeline, key, bars) {
   }
   const tileSpan = Math.max(span, MIN_TILE_SPAN);
   const tileCount = Math.ceil(totalBeats / tileSpan);
+  assertGenerationBudget(tileCount * cell.notes.length, "developed motif notes");
   for (let k = 0; k < tileCount; k += 1) {
     const offset = k * tileSpan;
     for (const n of cell.notes) {
@@ -3202,16 +3748,6 @@ function progressions() {
 function progressionsByStyle(style) {
   return progressions().filter((p) => p.styles.includes(style));
 }
-function mulberry323(seed) {
-  let state = seed >>> 0;
-  return () => {
-    state = state + 1831565813 >>> 0;
-    let t = state;
-    t = Math.imul(t ^ t >>> 15, t | 1);
-    t ^= t + Math.imul(t ^ t >>> 7, t | 61);
-    return ((t ^ t >>> 14) >>> 0) / 4294967296;
-  };
-}
 function degreeToRootPc(degree, key) {
   if (degree >= 0 && degree <= 6) {
     const tones = scaleTonesInDegreeOrder(key);
@@ -3233,6 +3769,8 @@ function autoQuality(degree, key) {
   return "maj";
 }
 function generateProgression(opts) {
+  assertPositiveInt(opts.bars, "progression bars");
+  assertGenerationBudget(opts.bars, "progression chords");
   const seed = opts.seed ?? 0;
   let preset;
   if (opts.presetId !== void 0) {
@@ -3244,8 +3782,8 @@ function generateProgression(opts) {
   if (preset === void 0) {
     const pool = PRESETS.filter((p) => p.styles.includes(opts.style));
     const candidates = pool.length > 0 ? pool : PRESETS;
-    const rng = mulberry323(seed);
-    const index = Math.floor(rng() * candidates.length) % candidates.length;
+    const rng = createRng(seed);
+    const index = Math.floor(rng.next() * candidates.length) % candidates.length;
     preset = candidates[index] ?? PRESETS[0];
   }
   const degrees = preset?.degrees ?? [0];
@@ -3263,7 +3801,7 @@ function generateProgression(opts) {
   }
   if (opts.reharmonize) {
     const tonicPc = (opts.key.rootPc % 12 + 12) % 12;
-    const rrng = mulberry323((seed ^ 2654435769) >>> 0);
+    const rrng = createRng((seed ^ 2654435769) >>> 0);
     for (let i = 0; i < chords.length - 1; i += 1) {
       const cur = chords[i];
       const next = chords[i + 1];
@@ -3277,7 +3815,7 @@ function generateProgression(opts) {
       if (targetsTonic || alreadySecondary || isResolutionTarget) {
         continue;
       }
-      if (rrng() < 0.5) {
+      if (rrng.next() < 0.5) {
         chords[i] = {
           rootPc: domRoot,
           quality: "dom7",
@@ -3297,12 +3835,21 @@ function pitchClass7(value) {
 function uniquePitchClasses(pitches) {
   return [...new Set(pitches.map(pitchClass7))].sort((a, b) => a - b);
 }
-function detectChord(pitches) {
+function assertPitches(pitches) {
+  assertGenerationBudget(pitches.length, "detection pitches");
+  for (let index = 0; index < pitches.length; index += 1) {
+    assertFiniteNumber(pitches[index] ?? Number.NaN, `pitches[${index}]`);
+  }
+}
+function detectChord(pitches, opts = {}) {
+  assertPitches(pitches);
+  if (opts.bassPc !== void 0) assertFiniteNumber(opts.bassPc, "bassPc");
   const input = uniquePitchClasses(pitches);
   if (input.length === 0) {
     return [];
   }
-  const bassPc = pitchClass7(Math.min(...pitches));
+  const inputKind = opts.input === void 0 || opts.input === "auto" ? pitches.every((pitch) => Number.isInteger(pitch) && pitch >= 0 && pitch <= 11) ? "pitchClass" : "midi" : opts.input;
+  const bassPc = opts.bassPc !== void 0 ? pitchClass7(opts.bassPc) : inputKind === "midi" ? pitchClass7(pitches.reduce((lowest, pitch) => Math.min(lowest, pitch), Infinity)) : void 0;
   const inputSet = new Set(input);
   const matches = [];
   const qualities = chordQualities();
@@ -3317,8 +3864,8 @@ function detectChord(pitches) {
         continue;
       }
       const extraPcs = input.filter((pc) => !toneSet.has(pc));
-      const bassIndex = chord.intervals.findIndex((iv) => pitchClass7(rootPc + iv) === bassPc);
-      const inversion = bassIndex > 0 ? bassIndex : 0;
+      const bassIndex = bassPc === void 0 ? -1 : chord.intervals.findIndex((iv) => pitchClass7(rootPc + iv) === bassPc);
+      const inversion = bassPc === void 0 ? null : bassIndex > 0 ? bassIndex : 0;
       const match = {
         rootPc,
         quality,
@@ -3327,7 +3874,7 @@ function detectChord(pitches) {
         exact: extraPcs.length === 0 && missingPcs.length === 0,
         inversion
       };
-      if (inversion > 0) {
+      if (bassPc !== void 0 && bassPc !== rootPc) {
         match.bassPc = bassPc;
       }
       matches.push(match);
@@ -3349,224 +3896,14 @@ function detectChord(pitches) {
   });
   return matches;
 }
-function detectChordBest(pitches) {
-  const best = detectChord(pitches)[0];
+function detectChordBest(pitches, opts = {}) {
+  const best = detectChord(pitches, opts)[0];
   if (!best) {
     return null;
   }
   return makeChord(best.rootPc, best.quality, best.bassPc);
 }
 
-// src/analyze/functional/internal.ts
-function mod12(n) {
-  return (n % 12 + 12) % 12;
-}
-function isNeapolitan(chord, key) {
-  return mod12(chord.rootPc - key.rootPc) === 1 && chord.quality === "maj";
-}
-
-// src/analyze/functional/roman.ts
-var ROMAN = ["I", "II", "III", "IV", "V", "VI", "VII"];
-var OFFSET_SPELLING = [
-  [1, ""],
-  [2, "b"],
-  [2, ""],
-  [3, "b"],
-  [3, ""],
-  [4, ""],
-  [4, "#"],
-  [5, ""],
-  [6, "b"],
-  [6, ""],
-  [7, "b"],
-  [7, ""]
-];
-var ROMAN_STYLE = {
-  maj: { lower: false, suffix: "" },
-  min: { lower: true, suffix: "" },
-  dim: { lower: true, suffix: "o" },
-  aug: { lower: false, suffix: "+" },
-  maj7: { lower: false, suffix: "maj7" },
-  min7: { lower: true, suffix: "7" },
-  dom7: { lower: false, suffix: "7" },
-  dim7: { lower: true, suffix: "o7" },
-  m7b5: { lower: true, suffix: "\xF87" },
-  minMaj7: { lower: true, suffix: "maj7" },
-  aug7: { lower: false, suffix: "+7" },
-  augMaj7: { lower: false, suffix: "+maj7" },
-  majb5: { lower: false, suffix: "b5" },
-  "6": { lower: false, suffix: "add6" },
-  min6: { lower: true, suffix: "add6" },
-  "6/9": { lower: false, suffix: "69" },
-  sus2: { lower: false, suffix: "sus2" },
-  sus4: { lower: false, suffix: "sus4" },
-  add9: { lower: false, suffix: "add9" },
-  add11: { lower: false, suffix: "add11" },
-  maj9: { lower: false, suffix: "maj9" },
-  min9: { lower: true, suffix: "9" },
-  dom9: { lower: false, suffix: "9" },
-  "7b9": { lower: false, suffix: "7b9" },
-  "7#9": { lower: false, suffix: "7#9" },
-  "7#11": { lower: false, suffix: "7#11" },
-  "7b13": { lower: false, suffix: "7b13" },
-  "11": { lower: false, suffix: "11" },
-  "13": { lower: false, suffix: "13" },
-  "5": { lower: false, suffix: "5" }
-};
-var SUFFIX_QUALITY = new Map(
-  Object.entries(ROMAN_STYLE).map(
-    ([quality, style]) => [`${style.lower ? "l" : "u"}:${style.suffix}`, quality]
-  )
-);
-function romanStyle(quality) {
-  return ROMAN_STYLE[quality];
-}
-function romanSpelling(rootPc, key) {
-  const tones = scaleTonesInDegreeOrder(key);
-  const diatonic = tones.indexOf(mod12(rootPc));
-  if (diatonic >= 0) {
-    return { degreeNumber: diatonic + 1, accidental: "" };
-  }
-  for (let i = 0; i < tones.length; i += 1) {
-    const degreeNumber2 = i + 1;
-    if (degreeNumber2 === 1 || degreeNumber2 === 5) {
-      continue;
-    }
-    if (mod12((tones[i] ?? 0) - 1) === mod12(rootPc)) {
-      return { degreeNumber: degreeNumber2, accidental: "b" };
-    }
-  }
-  for (let i = 0; i < tones.length; i += 1) {
-    if (mod12((tones[i] ?? 0) + 1) === mod12(rootPc)) {
-      return { degreeNumber: i + 1, accidental: "#" };
-    }
-  }
-  const [degreeNumber, accidental] = OFFSET_SPELLING[mod12(rootPc - key.rootPc)] ?? [1, ""];
-  return { degreeNumber, accidental };
-}
-function baseMarker(quality) {
-  if (quality === "dim" || quality === "dim7") {
-    return "o";
-  }
-  if (quality === "m7b5") {
-    return "\xF8";
-  }
-  if (quality === "aug" || quality === "aug7") {
-    return "+";
-  }
-  if (quality === "augMaj7") {
-    return "+maj7";
-  }
-  if (quality === "maj7" || quality === "minMaj7") {
-    return "maj7";
-  }
-  return "";
-}
-var TRIAD_FIGURES = { 1: "6", 2: "64" };
-var SEVENTH_FIGURES = { 1: "65", 2: "43", 3: "42" };
-var TRIAD_FIGURE_QUALITIES = /* @__PURE__ */ new Set(["maj", "min", "dim", "aug"]);
-var SEVENTH_FIGURE_QUALITIES = /* @__PURE__ */ new Set([
-  "maj7",
-  "min7",
-  "dom7",
-  "dim7",
-  "m7b5",
-  "minMaj7",
-  "aug7",
-  "augMaj7"
-]);
-function chordToRoman(chord, key) {
-  const { degreeNumber, accidental } = romanSpelling(chord.rootPc, key);
-  const { lower, suffix } = romanStyle(chord.quality);
-  const numeral = ROMAN[degreeNumber - 1] ?? "I";
-  const cased = lower ? numeral.toLowerCase() : numeral;
-  let inversion = 0;
-  if (chord.bassPc !== void 0) {
-    const idx = chord.intervals.findIndex((iv) => mod12(chord.rootPc + iv) === chord.bassPc);
-    if (idx > 0) {
-      inversion = idx;
-    }
-  }
-  if (inversion > 0) {
-    const figure = SEVENTH_FIGURE_QUALITIES.has(chord.quality) ? SEVENTH_FIGURES[inversion] : TRIAD_FIGURE_QUALITIES.has(chord.quality) ? TRIAD_FIGURES[inversion] : void 0;
-    if (figure !== void 0) {
-      return `${accidental}${cased}${baseMarker(chord.quality)}${figure}`;
-    }
-  }
-  return `${accidental}${cased}${suffix}`;
-}
-
-// src/analyze/functional/function.ts
-var FUNCTION_BY_OFFSET = [
-  "tonic",
-  // 0  I
-  "subdominant",
-  // 1  bII (Neapolitan)
-  "subdominant",
-  // 2  ii
-  "tonic",
-  // 3  bIII
-  "tonic",
-  // 4  iii
-  "subdominant",
-  // 5  IV
-  "dominant",
-  // 6  #IV / bV
-  "dominant",
-  // 7  V
-  "subdominant",
-  // 8  bVI
-  "tonic",
-  // 9  vi
-  "subdominant",
-  // 10 bVII
-  "dominant"
-  // 11 vii
-];
-function isMinorKey(key) {
-  const hasMinorThird = key.modeMask12 >> 3 & 1;
-  const hasMajorThird = key.modeMask12 >> 4 & 1;
-  return Boolean(hasMinorThird) && !hasMajorThird;
-}
-function functionOf(chord, key) {
-  return FUNCTION_BY_OFFSET[mod12(chord.rootPc - key.rootPc)] ?? "tonic";
-}
-function isDiatonic(chord, key) {
-  return chordPitchClasses(chord).every((pc) => isScaleTone(pc, key));
-}
-function parallelKey(key) {
-  return {
-    rootPc: mod12(key.rootPc),
-    modeMask12: isMinorKey(key) ? MAJOR_MASK : NATURAL_MINOR_MASK
-  };
-}
-
-// src/analyze/functional/borrowed.ts
-function isHarmonicMinorAlteration(chord, key) {
-  return isMinorKey(key) && isDiatonic(chord, { rootPc: mod12(key.rootPc), modeMask12: HARMONIC_MINOR_MASK });
-}
-function isBorrowedChord(chord, key) {
-  if (isDiatonic(chord, key)) {
-    return false;
-  }
-  if (isHarmonicMinorAlteration(chord, key)) {
-    return false;
-  }
-  return isDiatonic(chord, parallelKey(key));
-}
-function borrowedSource(chord, key) {
-  if (isDiatonic(chord, key)) {
-    return null;
-  }
-  if (isNeapolitan(chord, key)) {
-    return "neapolitan";
-  }
-  if (isBorrowedChord(chord, key)) {
-    return isMinorKey(key) ? "parallel-major" : "parallel-minor";
-  }
-  return null;
-}
-
 // src/generate/reharmony/index.ts
 function mod122(n) {
   return (n % 12 + 12) % 12;
@@ -3582,11 +3919,6 @@ function isStableTarget(chord) {
   const hasThird = hasInterval(chord, 3) || hasInterval(chord, 4);
   return hasFifth && hasThird && !isDominantType(chord);
 }
-function triadPitchClasses(chord) {
-  const third = hasInterval(chord, 4) ? 4 : hasInterval(chord, 3) ? 3 : 4;
-  const fifth = hasInterval(chord, 8) ? 8 : hasInterval(chord, 6) ? 6 : 7;
-  return [mod122(chord.rootPc), mod122(chord.rootPc + third), mod122(chord.rootPc + fifth)];
-}
 function commonToneCount(a, b) {
   const set = new Set(b);
   return a.filter((pc) => set.has(pc)).length;
@@ -3604,7 +3936,7 @@ function sameChord(a, b) {
 }
 var THIRD_OFFSETS = [3, 4, 8, 9];
 function substituteChord(chord, key, opts) {
-  const original = triadPitchClasses(chord);
+  const original = chordPitchClasses(chord);
   const candidates = [];
   if (isDominantType(chord)) {
     candidates.push({ chord: makeChord(mod122(chord.rootPc + 6), "dom7"), type: "tritone" });
@@ -3725,11 +4057,15 @@ function generateRhythm(ts, opts = {}) {
   const seed = opts.seed ?? 0;
   const bars = opts.bars ?? DEFAULT_BARS;
   const subdivision = opts.subdivision ?? DEFAULT_SUBDIVISION;
-  const density = Math.min(1, Math.max(0, opts.density ?? DEFAULT_DENSITY));
+  const requestedDensity = assertFiniteNumber(opts.density ?? DEFAULT_DENSITY, "rhythm density");
+  const density = Math.min(1, Math.max(0, requestedDensity));
+  assertPositiveInt(bars, "rhythm bars");
+  assertPositiveInt(subdivision, "rhythm subdivision");
   const barBeats = beatsPerBar(ts);
   const spanBeats = barBeats * bars;
   const step = 1 / subdivision;
   const slotCount = Math.round(spanBeats * subdivision);
+  assertGenerationBudget(slotCount, "rhythm grid slots");
   const slotsPerBar = Math.round(barBeats * subdivision);
   const rng = createRng(seed);
   const positions = [];