1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
// Shared utilities for mq_editor
// Evaluate interpolating curve at time t via Lagrange interpolation.
// The curve passes through all 4 control points at their stored time positions.
// Knot positions: u_k = (t_k - t0) / (t3 - t0), k=0..3.
// TODO: support arbitrary number of inner control points
function evalBezier(curve, t) {
const dt = curve.t3 - curve.t0;
if (dt <= 0) return curve.v0;
let u = (t - curve.t0) / dt;
u = Math.max(0, Math.min(1, u));
const u1 = (curve.t1 - curve.t0) / dt;
const u2 = (curve.t2 - curve.t0) / dt;
const d0 = (-u1) * (-u2) * (-1);
const d1 = u1 * (u1 - u2) * (u1 - 1);
const d2 = u2 * (u2 - u1) * (u2 - 1);
const d3 = (1 - u1) * (1 - u2);
if (Math.abs(d0) < 1e-9 || Math.abs(d1) < 1e-9 || Math.abs(d2) < 1e-9 || Math.abs(d3) < 1e-9)
return curve.v0 + (curve.v3 - curve.v0) * u;
const l0 = (u - u1) * (u - u2) * (u - 1) / d0;
const l1 = u * (u - u2) * (u - 1) / d1;
const l2 = u * (u - u1) * (u - 1) / d2;
const l3 = u * (u - u1) * (u - u2) / d3;
return l0 * curve.v0 + l1 * curve.v1 + l2 * curve.v2 + l3 * curve.v3;
}
// Evaluate amplitude component of interpolating curve at time t
function evalBezierAmp(curve, t) {
const dt = curve.t3 - curve.t0;
if (dt <= 0) return curve.a0;
let u = (t - curve.t0) / dt;
u = Math.max(0, Math.min(1, u));
const u1 = (curve.t1 - curve.t0) / dt;
const u2 = (curve.t2 - curve.t0) / dt;
const d0 = (-u1) * (-u2) * (-1);
const d1 = u1 * (u1 - u2) * (u1 - 1);
const d2 = u2 * (u2 - u1) * (u2 - 1);
const d3 = (1 - u1) * (1 - u2);
if (Math.abs(d0) < 1e-9 || Math.abs(d1) < 1e-9 || Math.abs(d2) < 1e-9 || Math.abs(d3) < 1e-9)
return curve.a0 + (curve.a3 - curve.a0) * u;
const l0 = (u - u1) * (u - u2) * (u - 1) / d0;
const l1 = u * (u - u2) * (u - 1) / d1;
const l2 = u * (u - u1) * (u - 1) / d2;
const l3 = u * (u - u1) * (u - u2) / d3;
return l0 * curve.a0 + l1 * curve.a1 + l2 * curve.a2 + l3 * curve.a3;
}
// Get canvas-relative {x, y} from a mouse event
function getCanvasCoords(e, canvas) {
const rect = canvas.getBoundingClientRect();
return { x: e.clientX - rect.left, y: e.clientY - rect.top };
}
// Build upper/lower band point arrays for a frequency curve.
// factorAbove/factorBelow are fractional offsets (e.g. 0.02 = ±2%).
// Returns { upper: [[x,y],...], lower: [[x,y],...] }
function buildBandPoints(viewer, curve, factorAbove, factorBelow) {
const STEPS = 60;
const upper = [], lower = [];
for (let i = 0; i <= STEPS; ++i) {
const t = curve.t0 + (curve.t3 - curve.t0) * i / STEPS;
if (t < viewer.t_view_min - 0.01 || t > viewer.t_view_max + 0.01) continue;
const f = evalBezier(curve, t);
upper.push([viewer.timeToX(t), viewer.freqToY(f * (1 + factorAbove))]);
lower.push([viewer.timeToX(t), viewer.freqToY(f * (1 - factorBelow))]);
}
return { upper, lower };
}
|
