velo-plot 2.0.0 → 2.20.0

package/dist/plugins/analysis.js

@@ -1,18 +1,18 @@
-import { a as A, s as W, d as _, i as z } from "../math-Ddm2EZjj.js";
-import { c as on } from "../math-Ddm2EZjj.js";
+import { a as F, s as W, d as _, i as z } from "../math-Ddm2EZjj.js";
+import { c as cn } from "../math-Ddm2EZjj.js";
 import { g as E, f as T, d as H, s as X, m as Y, a as N, b as O, v as Z, c as G } from "../utils-BK5UUXxj.js";
-import { h as sn, e as cn, i as an } from "../utils-BK5UUXxj.js";
+import { h as ln, e as fn, i as un } from "../utils-BK5UUXxj.js";
 import { r as K, e as Q, s as U } from "../indicators--eJGmVox.js";
-import { i as fn, j as un, g as hn, f as mn, k as gn, d as dn, m as Mn, c as wn, n as pn, o as yn, p as Fn, b as An, h as bn, a as In, t as qn, v as xn, w as Pn } from "../indicators--eJGmVox.js";
+import { i as mn, j as gn, g as dn, f as Mn, k as wn, d as pn, m as yn, c as bn, n as Fn, o as An, p as xn, b as In, h as qn, a as Pn, t as vn, v as kn, w as $n } from "../indicators--eJGmVox.js";
 function J(e, r, o, t = {}) {
   const n = e.series.get(r);
   if (!n) throw new Error(`Series ${r} not found`);
-  const s = n.getData();
-  if (!s || s.x.length < 2) return "";
-  const c = L(s.x, s.y, o, t), i = n.getBounds(), l = (i == null ? void 0 : i.xMin) ?? 0, a = ((i == null ? void 0 : i.xMax) ?? 1) - l, u = 200, h = new Float32Array(u), m = new Float32Array(u);
+  const c = n.getData();
+  if (!c || c.x.length < 2) return "";
+  const s = C(c.x, c.y, o, t), a = n.getBounds(), l = (a == null ? void 0 : a.xMin) ?? 0, i = ((a == null ? void 0 : a.xMax) ?? 1) - l, u = 200, h = new Float32Array(u), m = new Float32Array(u);
   for (let d = 0; d < u; d++) {
-    const w = l + d / (u - 1) * a;
-    h[d] = w, m[d] = c.predict(w);
+    const w = l + d / (u - 1) * i;
+    h[d] = w, m[d] = s.predict(w);
   }
   const g = `${r}-fit-${Date.now()}`, M = n.getStyle();
   return e.addSeries({
@@ -28,10 +28,10 @@ function J(e, r, o, t = {}) {
     }
   }), e.addAnnotation({
     type: "text",
-    x: l + a * 0.05,
-    y: c.predict(l + a * 0.05),
-    text: `${c.equation}
-(R² = ${c.rSquared.toFixed(4)})`,
+    x: l + i * 0.05,
+    y: s.predict(l + i * 0.05),
+    text: `${s.equation}
+(R² = ${s.rSquared.toFixed(4)})`,
     fontSize: 12,
     backgroundColor: "rgba(0,0,0,0.7)",
     color: M.color || "#ffffff",
@@ -40,50 +40,50 @@ function J(e, r, o, t = {}) {
     interactive: !0
   }), g;
 }
-function x(e) {
-  const r = e.length, o = C(r), t = new Float64Array(o);
-  for (let s = 0; s < r; s++)
-    t[s] = e[s];
+function q(e) {
+  const r = e.length, o = R(r), t = new Float64Array(o);
+  for (let c = 0; c < r; c++)
+    t[c] = e[c];
   const n = new Array(o);
-  for (let s = 0; s < o; s++)
-    n[s] = { re: t[s], im: 0 };
-  b(n);
-  for (let s = 2; s <= o; s *= 2) {
-    const c = s / 2, i = -2 * Math.PI / s;
-    for (let l = 0; l < o; l += s)
-      for (let f = 0; f < c; f++) {
-        const a = i * f, u = Math.cos(a), h = Math.sin(a), m = l + f, g = l + f + c, M = u * n[g].re - h * n[g].im, d = u * n[g].im + h * n[g].re;
+  for (let c = 0; c < o; c++)
+    n[c] = { re: t[c], im: 0 };
+  A(n);
+  for (let c = 2; c <= o; c *= 2) {
+    const s = c / 2, a = -2 * Math.PI / c;
+    for (let l = 0; l < o; l += c)
+      for (let f = 0; f < s; f++) {
+        const i = a * f, u = Math.cos(i), h = Math.sin(i), m = l + f, g = l + f + s, M = u * n[g].re - h * n[g].im, d = u * n[g].im + h * n[g].re;
         n[g].re = n[m].re - M, n[g].im = n[m].im - d, n[m].re = n[m].re + M, n[m].im = n[m].im + d;
       }
   }
   return n;
 }
 function V(e) {
-  const r = e.length, o = e.map((s) => ({ re: s.re, im: -s.im })), t = tt(o), n = new Float32Array(r);
-  for (let s = 0; s < r; s++)
-    n[s] = t[s].re / r;
+  const r = e.length, o = e.map((c) => ({ re: c.re, im: -c.im })), t = tt(o), n = new Float32Array(r);
+  for (let c = 0; c < r; c++)
+    n[c] = t[c].re / r;
   return n;
 }
 function tt(e) {
   const r = e.length, o = e.map((t) => ({ ...t }));
-  b(o);
+  A(o);
   for (let t = 2; t <= r; t *= 2) {
-    const n = t / 2, s = -2 * Math.PI / t;
-    for (let c = 0; c < r; c += t)
-      for (let i = 0; i < n; i++) {
-        const l = s * i, f = Math.cos(l), a = Math.sin(l), u = c + i, h = c + i + n, m = f * o[h].re - a * o[h].im, g = f * o[h].im + a * o[h].re;
+    const n = t / 2, c = -2 * Math.PI / t;
+    for (let s = 0; s < r; s += t)
+      for (let a = 0; a < n; a++) {
+        const l = c * a, f = Math.cos(l), i = Math.sin(l), u = s + a, h = s + a + n, m = f * o[h].re - i * o[h].im, g = f * o[h].im + i * o[h].re;
         o[h].re = o[u].re - m, o[h].im = o[u].im - g, o[u].re = o[u].re + m, o[u].im = o[u].im + g;
       }
   }
   return o;
 }
-function b(e) {
+function A(e) {
   const r = e.length, o = Math.log2(r);
   for (let t = 0; t < r; t++) {
     const n = nt(t, o);
     if (n > t) {
-      const s = e[t];
-      e[t] = e[n], e[n] = s;
+      const c = e[t];
+      e[t] = e[n], e[n] = c;
     }
   }
 }
@@ -93,49 +93,49 @@ function nt(e, r) {
     o = o << 1 | e & 1, e >>= 1;
   return o;
 }
-function C(e) {
+function R(e) {
   let r = 1;
   for (; r < e; ) r *= 2;
   return r;
 }
-function j(e, r = 1) {
-  const o = x(e), t = o.length, n = t / 2, s = new Float32Array(n), c = new Float32Array(n), i = new Float32Array(n), l = new Float32Array(t), f = new Float32Array(t);
-  for (let a = 0; a < t; a++)
-    l[a] = o[a].re, f[a] = o[a].im;
-  for (let a = 0; a < n; a++)
-    s[a] = a * r / t, c[a] = Math.sqrt(o[a].re * o[a].re + o[a].im * o[a].im) / t, i[a] = Math.atan2(o[a].im, o[a].re);
-  for (let a = 1; a < n; a++)
-    c[a] *= 2;
+function L(e, r = 1) {
+  const o = q(e), t = o.length, n = t / 2, c = new Float32Array(n), s = new Float32Array(n), a = new Float32Array(n), l = new Float32Array(t), f = new Float32Array(t);
+  for (let i = 0; i < t; i++)
+    l[i] = o[i].re, f[i] = o[i].im;
+  for (let i = 0; i < n; i++)
+    c[i] = i * r / t, s[i] = Math.sqrt(o[i].re * o[i].re + o[i].im * o[i].im) / t, a[i] = Math.atan2(o[i].im, o[i].re);
+  for (let i = 1; i < n; i++)
+    s[i] *= 2;
   return {
-    frequency: s,
-    magnitude: c,
-    phase: i,
+    frequency: c,
+    magnitude: s,
+    phase: a,
     complex: o,
     real: l,
     imag: f
   };
 }
-function kt(e, r = 1) {
-  const o = j(e, r), t = o.frequency.length, n = new Float32Array(t), s = new Float32Array(t);
-  for (let c = 0; c < t; c++)
-    n[c] = o.magnitude[c] * o.magnitude[c], s[c] = 10 * Math.log10(Math.max(n[c], 1e-20));
+function Lt(e, r = 1) {
+  const o = L(e, r), t = o.frequency.length, n = new Float32Array(t), c = new Float32Array(t);
+  for (let s = 0; s < t; s++)
+    n[s] = o.magnitude[s] * o.magnitude[s], c[s] = 10 * Math.log10(Math.max(n[s], 1e-20));
   return {
     frequency: o.frequency,
     power: n,
-    powerDb: s
+    powerDb: c
   };
 }
-function Rt(e, r = 1, o = 0) {
-  const t = j(e, r);
-  let n = 0, s = 0;
-  for (let c = 0; c < t.magnitude.length; c++)
-    t.frequency[c] >= o && t.magnitude[c] > n && (n = t.magnitude[c], s = c);
+function Ct(e, r = 1, o = 0) {
+  const t = L(e, r);
+  let n = 0, c = 0;
+  for (let s = 0; s < t.magnitude.length; s++)
+    t.frequency[s] >= o && t.magnitude[s] > n && (n = t.magnitude[s], c = s);
   return {
-    frequency: t.frequency[s],
+    frequency: t.frequency[c],
     magnitude: n
   };
 }
-function Ct(e) {
+function St(e) {
   const r = e.length, o = new Float32Array(r);
   for (let t = 0; t < r; t++) {
     const n = 0.5 * (1 - Math.cos(2 * Math.PI * t / (r - 1)));
@@ -143,7 +143,7 @@ function Ct(e) {
   }
   return o;
 }
-function jt(e) {
+function Bt(e) {
   const r = e.length, o = new Float32Array(r);
   for (let t = 0; t < r; t++) {
     const n = 0.54 - 0.46 * Math.cos(2 * Math.PI * t / (r - 1));
@@ -151,61 +151,61 @@ function jt(e) {
   }
   return o;
 }
-function Lt(e) {
-  const r = e.length, o = new Float32Array(r), t = 0.42, n = 0.5, s = 0.08;
-  for (let c = 0; c < r; c++) {
-    const i = t - n * Math.cos(2 * Math.PI * c / (r - 1)) + s * Math.cos(4 * Math.PI * c / (r - 1));
-    o[c] = e[c] * i;
+function Dt(e) {
+  const r = e.length, o = new Float32Array(r), t = 0.42, n = 0.5, c = 0.08;
+  for (let s = 0; s < r; s++) {
+    const a = t - n * Math.cos(2 * Math.PI * s / (r - 1)) + c * Math.cos(4 * Math.PI * s / (r - 1));
+    o[s] = e[s] * a;
   }
   return o;
 }
-function St(e, r = 1) {
-  const o = x(e), t = o.length, n = t / 2, s = new Float32Array(t), c = new Float32Array(t), i = new Float32Array(t), l = new Float32Array(t), f = new Float32Array(t);
-  for (let a = 0; a < t; a++)
-    s[a] = o[a].re, c[a] = o[a].im, i[a] = a * r / t, l[a] = Math.sqrt(o[a].re * o[a].re + o[a].im * o[a].im), f[a] = Math.atan2(o[a].im, o[a].re);
+function Wt(e, r = 1) {
+  const o = q(e), t = o.length, n = t / 2, c = new Float32Array(t), s = new Float32Array(t), a = new Float32Array(t), l = new Float32Array(t), f = new Float32Array(t);
+  for (let i = 0; i < t; i++)
+    c[i] = o[i].re, s[i] = o[i].im, a[i] = i * r / t, l[i] = Math.sqrt(o[i].re * o[i].re + o[i].im * o[i].im), f[i] = Math.atan2(o[i].im, o[i].re);
   return {
-    real: s,
-    imag: c,
-    frequency: i,
+    real: c,
+    imag: s,
+    frequency: a,
     magnitude: l,
     phase: f,
     length: t,
     nyquist: n
   };
 }
-function Bt(e, r) {
+function _t(e, r) {
   const o = e.length;
   if (r.length !== o)
     throw new Error("Real and imaginary arrays must have the same length");
-  const t = C(o), n = new Array(t);
-  for (let a = 0; a < t; a++)
-    n[a] = {
-      re: a < o ? e[a] : 0,
-      im: a < o ? r[a] : 0
+  const t = R(o), n = new Array(t);
+  for (let i = 0; i < t; i++)
+    n[i] = {
+      re: i < o ? e[i] : 0,
+      im: i < o ? r[i] : 0
     };
-  b(n);
-  for (let a = 2; a <= t; a *= 2) {
-    const u = a / 2, h = -2 * Math.PI / a;
-    for (let m = 0; m < t; m += a)
+  A(n);
+  for (let i = 2; i <= t; i *= 2) {
+    const u = i / 2, h = -2 * Math.PI / i;
+    for (let m = 0; m < t; m += i)
       for (let g = 0; g < u; g++) {
-        const M = h * g, d = Math.cos(M), w = Math.sin(M), p = m + g, y = m + g + u, F = d * n[y].re - w * n[y].im, P = d * n[y].im + w * n[y].re;
-        n[y].re = n[p].re - F, n[y].im = n[p].im - P, n[p].re = n[p].re + F, n[p].im = n[p].im + P;
+        const M = h * g, d = Math.cos(M), w = Math.sin(M), p = m + g, y = m + g + u, b = d * n[y].re - w * n[y].im, P = d * n[y].im + w * n[y].re;
+        n[y].re = n[p].re - b, n[y].im = n[p].im - P, n[p].re = n[p].re + b, n[p].im = n[p].im + P;
       }
   }
-  const s = new Float32Array(t), c = new Float32Array(t), i = new Float32Array(t), l = new Float32Array(t), f = new Float32Array(t);
-  for (let a = 0; a < t; a++)
-    s[a] = n[a].re, c[a] = n[a].im, i[a] = a / t, l[a] = Math.sqrt(n[a].re * n[a].re + n[a].im * n[a].im), f[a] = Math.atan2(n[a].im, n[a].re);
+  const c = new Float32Array(t), s = new Float32Array(t), a = new Float32Array(t), l = new Float32Array(t), f = new Float32Array(t);
+  for (let i = 0; i < t; i++)
+    c[i] = n[i].re, s[i] = n[i].im, a[i] = i / t, l[i] = Math.sqrt(n[i].re * n[i].re + n[i].im * n[i].im), f[i] = Math.atan2(n[i].im, n[i].re);
   return {
-    real: s,
-    imag: c,
-    frequency: i,
+    real: c,
+    imag: s,
+    frequency: a,
     magnitude: l,
     phase: f,
     length: t,
     nyquist: t / 2
   };
 }
-function Dt(e) {
+function zt(e) {
   const r = e.length, o = new Float32Array(r), t = new Float32Array(r);
   for (let n = 0; n < r; n++)
     o[n] = e[n].re, t[n] = e[n].im;
@@ -220,27 +220,27 @@ function et(e, r) {
     t[n] = { re: e[n], im: r[n] };
   return t;
 }
-function Wt(e, r) {
+function Et(e, r) {
   const o = et(e, r);
   return V(o);
 }
-function _t(e) {
-  const r = e.length, t = e.map((c) => ({ re: c.re, im: -c.im })).map((c) => ({ ...c }));
-  b(t);
-  for (let c = 2; c <= r; c *= 2) {
-    const i = c / 2, l = -2 * Math.PI / c;
-    for (let f = 0; f < r; f += c)
-      for (let a = 0; a < i; a++) {
-        const u = l * a, h = Math.cos(u), m = Math.sin(u), g = f + a, M = f + a + i, d = h * t[M].re - m * t[M].im, w = h * t[M].im + m * t[M].re;
+function Tt(e) {
+  const r = e.length, t = e.map((s) => ({ re: s.re, im: -s.im })).map((s) => ({ ...s }));
+  A(t);
+  for (let s = 2; s <= r; s *= 2) {
+    const a = s / 2, l = -2 * Math.PI / s;
+    for (let f = 0; f < r; f += s)
+      for (let i = 0; i < a; i++) {
+        const u = l * i, h = Math.cos(u), m = Math.sin(u), g = f + i, M = f + i + a, d = h * t[M].re - m * t[M].im, w = h * t[M].im + m * t[M].re;
         t[M].re = t[g].re - d, t[M].im = t[g].im - w, t[g].re = t[g].re + d, t[g].im = t[g].im + w;
       }
   }
-  const n = new Float32Array(r), s = new Float32Array(r);
-  for (let c = 0; c < r; c++)
-    n[c] = t[c].re / r, s[c] = -t[c].im / r;
-  return { real: n, imag: s };
+  const n = new Float32Array(r), c = new Float32Array(r);
+  for (let s = 0; s < r; s++)
+    n[s] = t[s].re / r, c[s] = -t[s].im / r;
+  return { real: n, imag: c };
 }
-function zt(e) {
+function Ht(e) {
   const r = e.nyquist;
   return {
     real: e.real.slice(0, r),
@@ -250,11 +250,11 @@ function zt(e) {
     phase: e.phase.slice(0, r)
   };
 }
-function L(e, r, o, t = {}) {
+function C(e, r, o, t = {}) {
   if (e.length < 2) throw new Error("At least 2 points are required for fitting");
   switch (o) {
     case "linear":
-      return I(e, r, t);
+      return x(e, r, t);
     case "polynomial":
       return ot(e, r, t.degree ?? 2, t);
     case "exponential":
@@ -267,44 +267,44 @@ function L(e, r, o, t = {}) {
       throw new Error(`Unsupported fit type: ${o}`);
   }
 }
-function I(e, r, o) {
+function x(e, r, o) {
   const t = e.length;
-  let n = 0, s = 0, c = 0, i = 0;
+  let n = 0, c = 0, s = 0, a = 0;
   for (let m = 0; m < t; m++)
-    n += e[m], s += r[m], c += e[m] * r[m], i += e[m] * e[m];
-  const l = (t * c - n * s) / (t * i - n * n), f = (s - l * n) / t, a = (m) => l * m + f, u = o.precision ?? 4, h = o.label || "y";
+    n += e[m], c += r[m], s += e[m] * r[m], a += e[m] * e[m];
+  const l = (t * s - n * c) / (t * a - n * n), f = (c - l * n) / t, i = (m) => l * m + f, u = o.precision ?? 4, h = o.label || "y";
   return {
     type: "linear",
     coefficients: [l, f],
     equation: `${h} = ${l.toFixed(u)}x ${f >= 0 ? "+" : "-"} ${Math.abs(f).toFixed(u)}`,
-    rSquared: A(e, r, a),
-    predict: a
+    rSquared: F(e, r, i),
+    predict: i
   };
 }
 function ot(e, r, o, t) {
-  const n = e.length, s = o + 1, c = Array.from({ length: s }, () => new Array(s).fill(0)), i = new Array(s).fill(0);
-  for (let m = 0; m < s; m++) {
-    for (let M = 0; M < s; M++) {
+  const n = e.length, c = o + 1, s = Array.from({ length: c }, () => new Array(c).fill(0)), a = new Array(c).fill(0);
+  for (let m = 0; m < c; m++) {
+    for (let M = 0; M < c; M++) {
       let d = 0;
       for (let w = 0; w < n; w++)
         d += Math.pow(e[w], m + M);
-      c[m][M] = d;
+      s[m][M] = d;
     }
     let g = 0;
     for (let M = 0; M < n; M++)
       g += r[M] * Math.pow(e[M], m);
-    i[m] = g;
+    a[m] = g;
   }
-  const l = W(c, i), f = (m) => {
+  const l = W(s, a), f = (m) => {
     let g = 0;
     for (let M = 0; M < l.length; M++)
       g += l[M] * Math.pow(m, M);
     return g;
-  }, a = t.precision ?? 4;
+  }, i = t.precision ?? 4;
   let u = "y = ";
   for (let m = l.length - 1; m >= 0; m--) {
     const g = l[m];
-    m < l.length - 1 ? u += g >= 0 ? " + " : " - " : g < 0 && (u += "-"), u += `${Math.abs(g).toFixed(a)}${m > 0 ? m > 1 ? `x^${m}` : "x" : ""}`;
+    m < l.length - 1 ? u += g >= 0 ? " + " : " - " : g < 0 && (u += "-"), u += `${Math.abs(g).toFixed(i)}${m > 0 ? m > 1 ? `x^${m}` : "x" : ""}`;
   }
   const h = t.label || "y";
   return {
@@ -312,127 +312,127 @@ function ot(e, r, o, t) {
     coefficients: l,
     equation: `${h} = ${u.substring(4)}`,
     // replace "y = " with custom label
-    rSquared: A(e, r, f),
+    rSquared: F(e, r, f),
     predict: f
   };
 }
 function rt(e, r, o) {
-  const t = e.length, n = [], s = [];
+  const t = e.length, n = [], c = [];
   for (let h = 0; h < t; h++)
-    r[h] > 0 && (n.push(e[h]), s.push(Math.log(r[h])));
-  const c = I(n, s, o), i = c.coefficients[0], l = Math.exp(c.coefficients[1]), f = (h) => l * Math.exp(i * h), a = o.precision ?? 4, u = o.label || "y";
+    r[h] > 0 && (n.push(e[h]), c.push(Math.log(r[h])));
+  const s = x(n, c, o), a = s.coefficients[0], l = Math.exp(s.coefficients[1]), f = (h) => l * Math.exp(a * h), i = o.precision ?? 4, u = o.label || "y";
   return {
     type: "exponential",
-    coefficients: [l, i],
-    equation: `${u} = ${l.toFixed(a)} * e^(${i.toFixed(a)}x)`,
-    rSquared: A(e, r, f),
+    coefficients: [l, a],
+    equation: `${u} = ${l.toFixed(i)} * e^(${a.toFixed(i)}x)`,
+    rSquared: F(e, r, f),
     predict: f
   };
 }
 function st(e, r, o) {
-  const t = e.length, n = [], s = [];
+  const t = e.length, n = [], c = [];
   for (let h = 0; h < t; h++)
-    e[h] > 0 && (n.push(Math.log(e[h])), s.push(r[h]));
-  const c = I(n, s, o), i = c.coefficients[0], l = c.coefficients[1], f = (h) => l + i * Math.log(h), a = o.precision ?? 4, u = o.label || "y";
+    e[h] > 0 && (n.push(Math.log(e[h])), c.push(r[h]));
+  const s = x(n, c, o), a = s.coefficients[0], l = s.coefficients[1], f = (h) => l + a * Math.log(h), i = o.precision ?? 4, u = o.label || "y";
   return {
     type: "logarithmic",
-    coefficients: [l, i],
-    equation: `${u} = ${l.toFixed(a)} ${i >= 0 ? "+" : "-"} ${Math.abs(i).toFixed(a)} * ln(x)`,
-    rSquared: A(e, r, f),
+    coefficients: [l, a],
+    equation: `${u} = ${l.toFixed(i)} ${a >= 0 ? "+" : "-"} ${Math.abs(a).toFixed(i)} * ln(x)`,
+    rSquared: F(e, r, f),
     predict: f
   };
 }
 function ct(e, r, o) {
-  const t = e.length, n = [], s = [];
+  const t = e.length, n = [], c = [];
   for (let h = 0; h < t; h++)
-    e[h] > 0 && r[h] > 0 && (n.push(Math.log(e[h])), s.push(Math.log(r[h])));
-  const c = I(n, s, o), i = c.coefficients[0], l = Math.exp(c.coefficients[1]), f = (h) => l * Math.pow(h, i), a = o.precision ?? 4, u = o.label || "y";
+    e[h] > 0 && r[h] > 0 && (n.push(Math.log(e[h])), c.push(Math.log(r[h])));
+  const s = x(n, c, o), a = s.coefficients[0], l = Math.exp(s.coefficients[1]), f = (h) => l * Math.pow(h, a), i = o.precision ?? 4, u = o.label || "y";
   return {
     type: "power",
-    coefficients: [l, i],
-    equation: `${u} = ${l.toFixed(a)} * x^(${i.toFixed(a)})`,
-    rSquared: A(e, r, f),
+    coefficients: [l, a],
+    equation: `${u} = ${l.toFixed(i)} * x^(${a.toFixed(i)})`,
+    rSquared: F(e, r, f),
     predict: f
   };
 }
-function Et(e, r, o, t = 5) {
-  const n = r / (o / 2), s = S(Math.min(n, 0.99), t);
-  return q(e, s);
+function Xt(e, r, o, t = 5) {
+  const n = r / (o / 2), c = S(Math.min(n, 0.99), t);
+  return I(e, c);
 }
-function Tt(e, r, o, t = 5) {
-  const n = r / (o / 2), s = it(Math.min(n, 0.99), t);
-  return q(e, s);
+function Yt(e, r, o, t = 5) {
+  const n = r / (o / 2), c = at(Math.min(n, 0.99), t);
+  return I(e, c);
 }
-function Ht(e, r, o, t, n = 5) {
-  const s = r / (t / 2), c = o / (t / 2), i = B(
-    Math.max(s, 0.01),
-    Math.min(c, 0.99),
+function Nt(e, r, o, t, n = 5) {
+  const c = r / (t / 2), s = o / (t / 2), a = B(
+    Math.max(c, 0.01),
+    Math.min(s, 0.99),
     n
   );
-  return q(e, i);
+  return I(e, a);
 }
-function Xt(e, r, o, t, n = 5) {
-  const s = r / (t / 2), c = o / (t / 2), i = at(
-    Math.max(s, 0.01),
-    Math.min(c, 0.99),
+function Ot(e, r, o, t, n = 5) {
+  const c = r / (t / 2), s = o / (t / 2), a = it(
+    Math.max(c, 0.01),
+    Math.min(s, 0.99),
     n
   );
-  return q(e, i);
+  return I(e, a);
 }
 function S(e, r) {
   const o = 2 * r + 1, t = new Float32Array(o), n = e;
-  for (let c = 0; c < o; c++) {
-    const i = c - r;
-    i === 0 ? t[c] = 2 * n : t[c] = Math.sin(2 * Math.PI * n * i) / (Math.PI * i), t[c] *= 0.54 - 0.46 * Math.cos(2 * Math.PI * c / (o - 1));
+  for (let s = 0; s < o; s++) {
+    const a = s - r;
+    a === 0 ? t[s] = 2 * n : t[s] = Math.sin(2 * Math.PI * n * a) / (Math.PI * a), t[s] *= 0.54 - 0.46 * Math.cos(2 * Math.PI * s / (o - 1));
   }
-  let s = 0;
-  for (let c = 0; c < o; c++) s += t[c];
-  for (let c = 0; c < o; c++) t[c] /= s;
+  let c = 0;
+  for (let s = 0; s < o; s++) c += t[s];
+  for (let s = 0; s < o; s++) t[s] /= c;
   return t;
 }
-function it(e, r) {
+function at(e, r) {
   const o = S(e, r), t = o.length, n = new Float32Array(t);
-  for (let s = 0; s < t; s++)
-    n[s] = -o[s];
+  for (let c = 0; c < t; c++)
+    n[c] = -o[c];
   return n[r] += 1, n;
 }
 function B(e, r, o) {
-  const t = 2 * o + 1, n = new Float32Array(t), s = e, c = r;
-  for (let i = 0; i < t; i++) {
-    const l = i - o;
-    l === 0 ? n[i] = 2 * (c - s) : n[i] = (Math.sin(2 * Math.PI * c * l) - Math.sin(2 * Math.PI * s * l)) / (Math.PI * l), n[i] *= 0.54 - 0.46 * Math.cos(2 * Math.PI * i / (t - 1));
+  const t = 2 * o + 1, n = new Float32Array(t), c = e, s = r;
+  for (let a = 0; a < t; a++) {
+    const l = a - o;
+    l === 0 ? n[a] = 2 * (s - c) : n[a] = (Math.sin(2 * Math.PI * s * l) - Math.sin(2 * Math.PI * c * l)) / (Math.PI * l), n[a] *= 0.54 - 0.46 * Math.cos(2 * Math.PI * a / (t - 1));
   }
   return n;
 }
-function at(e, r, o) {
-  const t = B(e, r, o), n = t.length, s = new Float32Array(n);
-  for (let c = 0; c < n; c++)
-    s[c] = -t[c];
-  return s[o] += 1, s;
+function it(e, r, o) {
+  const t = B(e, r, o), n = t.length, c = new Float32Array(n);
+  for (let s = 0; s < n; s++)
+    c[s] = -t[s];
+  return c[o] += 1, c;
 }
-function q(e, r) {
-  const o = e.length, t = r.length, n = Math.floor(t / 2), s = new Float32Array(o);
-  for (let c = 0; c < o; c++) {
-    let i = 0;
+function I(e, r) {
+  const o = e.length, t = r.length, n = Math.floor(t / 2), c = new Float32Array(o);
+  for (let s = 0; s < o; s++) {
+    let a = 0;
     for (let l = 0; l < t; l++) {
-      const f = c - n + l;
-      f >= 0 && f < o && (i += e[f] * r[l]);
+      const f = s - n + l;
+      f >= 0 && f < o && (a += e[f] * r[l]);
     }
-    s[c] = i;
+    c[s] = a;
   }
-  return s;
+  return c;
 }
-function Yt(e, r) {
-  const { type: o, cutoff: t, cutoffHigh: n, sampleRate: s = 1, order: c = 2 } = r, i = 2 * t / s, l = n ? 2 * n / s : void 0, { b: f, a } = ut(o, i, l);
-  return D(e, f, a);
+function Zt(e, r) {
+  const { type: o, cutoff: t, cutoffHigh: n, sampleRate: c = 1, order: s = 2 } = r, a = 2 * t / c, l = n ? 2 * n / c : void 0, { b: f, a: i } = ut(o, a, l);
+  return D(e, f, i);
 }
 function lt(e, r) {
-  const { frequency: o, sampleRate: t, bandwidth: n = 1 } = r, { b: s, a: c } = ft(o, t, n);
-  return D(e, s, c);
+  const { frequency: o, sampleRate: t, bandwidth: n = 1 } = r, { b: c, a: s } = ft(o, t, n);
+  return D(e, c, s);
 }
 function ft(e, r, o) {
-  const t = 2 * Math.PI * e / r, s = 1 - 2 * Math.PI * o / r / 2, c = Math.cos(t), i = [1, -2 * c, 1], l = [1, -2 * s * c, s * s];
-  return { b: i, a: l };
+  const t = 2 * Math.PI * e / r, c = 1 - 2 * Math.PI * o / r / 2, s = Math.cos(t), a = [1, -2 * s, 1], l = [1, -2 * c * s, c * c];
+  return { b: a, a: l };
 }
 function ut(e, r, o, t) {
   const n = Math.tan(Math.PI * r / 2);
@@ -441,158 +441,244 @@ function ut(e, r, o, t) {
   if (e === "highpass")
     return ht(n);
   if (e === "bandpass" && o) {
-    const s = Math.tan(Math.PI * o / 2);
-    return mt(n, s);
+    const c = Math.tan(Math.PI * o / 2);
+    return mt(n, c);
   } else if (e === "bandstop" && o) {
-    const s = Math.tan(Math.PI * o / 2);
-    return gt(n, s);
+    const c = Math.tan(Math.PI * o / 2);
+    return gt(n, c);
   }
   return v(n);
 }
 function v(e, r) {
-  const o = e, t = o * o, n = Math.sqrt(2), s = 1 / (1 + n * o + t);
+  const o = e, t = o * o, n = Math.sqrt(2), c = 1 / (1 + n * o + t);
   return {
-    b: [t * s, 2 * t * s, t * s],
-    a: [1, 2 * (t - 1) * s, (1 - n * o + t) * s]
+    b: [t * c, 2 * t * c, t * c],
+    a: [1, 2 * (t - 1) * c, (1 - n * o + t) * c]
   };
 }
 function ht(e, r) {
-  const o = e, t = o * o, n = Math.sqrt(2), s = 1 / (1 + n * o + t);
+  const o = e, t = o * o, n = Math.sqrt(2), c = 1 / (1 + n * o + t);
   return {
-    b: [s, -2 * s, s],
-    a: [1, 2 * (t - 1) * s, (1 - n * o + t) * s]
+    b: [c, -2 * c, c],
+    a: [1, 2 * (t - 1) * c, (1 - n * o + t) * c]
   };
 }
 function mt(e, r, o) {
-  const t = Math.sqrt(e * r), n = r - e, s = t / n, c = 1 / (1 + t / s + t * t);
+  const t = Math.sqrt(e * r), n = r - e, c = t / n, s = 1 / (1 + t / c + t * t);
   return {
-    b: [t / s * c, 0, -t / s * c],
-    a: [1, 2 * (t * t - 1) * c, (1 - t / s + t * t) * c]
+    b: [t / c * s, 0, -t / c * s],
+    a: [1, 2 * (t * t - 1) * s, (1 - t / c + t * t) * s]
   };
 }
 function gt(e, r, o) {
-  const t = Math.sqrt(e * r), n = r - e, s = t / n, c = 1 / (1 + t / s + t * t);
+  const t = Math.sqrt(e * r), n = r - e, c = t / n, s = 1 / (1 + t / c + t * t);
   return {
-    b: [(1 + t * t) * c, 2 * (t * t - 1) * c, (1 + t * t) * c],
-    a: [1, 2 * (t * t - 1) * c, (1 - t / s + t * t) * c]
+    b: [(1 + t * t) * s, 2 * (t * t - 1) * s, (1 + t * t) * s],
+    a: [1, 2 * (t * t - 1) * s, (1 - t / c + t * t) * s]
   };
 }
-function $(e, r, o) {
+function k(e, r, o) {
   const t = e.length, n = new Float32Array(t);
-  for (let s = 0; s < t; s++) {
-    n[s] = r[0] * e[s];
-    for (let c = 1; c < r.length; c++)
-      s - c >= 0 && (n[s] += r[c] * e[s - c], n[s] -= o[c] * n[s - c]);
+  for (let c = 0; c < t; c++) {
+    n[c] = r[0] * e[c];
+    for (let s = 1; s < r.length; s++)
+      c - s >= 0 && (n[c] += r[s] * e[c - s], n[c] -= o[s] * n[c - s]);
   }
   return n;
 }
 function D(e, r, o) {
-  const t = $(e, r, o), n = new Float32Array(t.length);
-  for (let i = 0; i < t.length; i++)
-    n[i] = t[t.length - 1 - i];
-  const s = $(n, r, o), c = new Float32Array(s.length);
-  for (let i = 0; i < s.length; i++)
-    c[i] = s[s.length - 1 - i];
-  return c;
+  const t = k(e, r, o), n = new Float32Array(t.length);
+  for (let a = 0; a < t.length; a++)
+    n[a] = t[t.length - 1 - a];
+  const c = k(n, r, o), s = new Float32Array(c.length);
+  for (let a = 0; a < c.length; a++)
+    s[a] = c[c.length - 1 - a];
+  return s;
 }
-function Nt(e, r = 0.3) {
+function Gt(e, r = 0.3) {
   const o = e.length, t = new Float32Array(o);
   t[0] = e[0];
   for (let n = 1; n < o; n++)
     t[n] = r * e[n] + (1 - r) * t[n - 1];
   return t;
 }
-function Ot(e, r = 2) {
-  const o = e.length, t = new Float32Array(o), n = Math.ceil(r * 3) * 2 + 1, s = Math.floor(n / 2), c = new Float32Array(n);
-  let i = 0;
+function Kt(e, r = 2) {
+  const o = e.length, t = new Float32Array(o), n = Math.ceil(r * 3) * 2 + 1, c = Math.floor(n / 2), s = new Float32Array(n);
+  let a = 0;
   for (let l = 0; l < n; l++) {
-    const f = l - s;
-    c[l] = Math.exp(-(f * f) / (2 * r * r)), i += c[l];
+    const f = l - c;
+    s[l] = Math.exp(-(f * f) / (2 * r * r)), a += s[l];
   }
   for (let l = 0; l < n; l++)
-    c[l] /= i;
+    s[l] /= a;
   for (let l = 0; l < o; l++) {
     let f = 0;
-    for (let a = 0; a < n; a++) {
-      const u = l - s + a;
-      u >= 0 && u < o && (f += e[u] * c[a]);
+    for (let i = 0; i < n; i++) {
+      const u = l - c + i;
+      u >= 0 && u < o && (f += e[u] * s[i]);
     }
     t[l] = f;
   }
   return t;
 }
-function Zt(e, r = 5, o = 2) {
-  const t = e.length, n = new Float32Array(t), s = Math.floor(r / 2), c = dt(r);
-  for (let i = 0; i < t; i++) {
+function Qt(e, r = 5, o = 2) {
+  const t = e.length, n = new Float32Array(t), c = Math.floor(r / 2), s = dt(r);
+  for (let a = 0; a < t; a++) {
     let l = 0;
-    for (let f = -s; f <= s; f++) {
-      const a = i + f;
-      a >= 0 && a < t ? l += e[a] * c[f + s] : a < 0 ? l += e[0] * c[f + s] : l += e[t - 1] * c[f + s];
+    for (let f = -c; f <= c; f++) {
+      const i = a + f;
+      i >= 0 && i < t ? l += e[i] * s[f + c] : i < 0 ? l += e[0] * s[f + c] : l += e[t - 1] * s[f + c];
     }
-    n[i] = l;
+    n[a] = l;
   }
   return n;
 }
 function dt(e, r) {
   const o = Math.floor(e / 2), t = new Float32Array(e), n = 1 / e;
-  for (let c = 0; c < e; c++) {
-    const i = c - o;
-    t[c] = n * (1 - i * i / ((o + 1) * (o + 1)));
+  for (let s = 0; s < e; s++) {
+    const a = s - o;
+    t[s] = n * (1 - a * a / ((o + 1) * (o + 1)));
   }
-  let s = 0;
-  for (let c = 0; c < e; c++) s += t[c];
-  for (let c = 0; c < e; c++) t[c] /= s;
+  let c = 0;
+  for (let s = 0; s < e; s++) c += t[s];
+  for (let s = 0; s < e; s++) t[s] /= c;
   return t;
 }
-function Gt(e, r = 5) {
+function Ut(e, r = 5) {
   const o = e.length, t = new Float32Array(o), n = Math.floor(r / 2);
-  for (let s = 0; s < o; s++) {
-    const c = [];
-    for (let i = -n; i <= n; i++) {
-      const l = s + i;
-      l >= 0 && l < o && c.push(e[l]);
+  for (let c = 0; c < o; c++) {
+    const s = [];
+    for (let a = -n; a <= n; a++) {
+      const l = c + a;
+      l >= 0 && l < o && s.push(e[l]);
     }
-    c.sort((i, l) => i - l), t[s] = c[Math.floor(c.length / 2)];
+    s.sort((a, l) => a - l), t[c] = s[Math.floor(s.length / 2)];
   }
   return t;
 }
-function Kt(e, r, o, t = {}) {
-  const n = r.length, s = o.length, c = Math.min(...e), i = Math.max(...e), l = t.levels || Mt(c, i, t.numLevels || 10), f = [];
-  for (const a of l) {
+function Jt(e, r, o, t = {}) {
+  const n = r.length, c = o.length, s = Math.min(...e), a = Math.max(...e), l = t.levels || yt(s, a, t.numLevels || 10), f = [];
+  for (const i of l) {
     const u = [];
-    for (let h = 0; h < s - 1; h++)
+    for (let h = 0; h < c - 1; h++)
       for (let m = 0; m < n - 1; m++) {
         const g = e[h * n + m], M = e[h * n + (m + 1)], d = e[(h + 1) * n + (m + 1)], w = e[(h + 1) * n + m];
-        wt(m, h, g, M, d, w, a, r, o, u);
+        bt(m, h, g, M, d, w, i, r, o, u);
       }
     if (u.length > 0) {
       const h = [];
-      for (const m of u)
-        h.push({ x: m[0], y: m[1] }), h.push({ x: m[2], y: m[3] });
-      f.push({ level: a, points: h });
+      for (const g of u)
+        h.push({ x: g[0], y: g[1] }), h.push({ x: g[2], y: g[3] });
+      const m = { level: i, points: h };
+      if (t.joinPaths || t.labels) {
+        const g = Mt(u);
+        m.polylines = g, t.labels && (m.labels = pt(
+          g,
+          i,
+          t.minLabelLength ?? 6,
+          t.labelFormatter ?? wt
+        ));
+      }
+      f.push(m);
     }
   }
   return f;
 }
-function Mt(e, r, o) {
+function Mt(e, r = 1e-6) {
+  const o = e.map((a) => ({
+    a: { x: a[0], y: a[1] },
+    b: { x: a[2], y: a[3] },
+    used: !1
+  })), t = (a) => `${Math.round(a.x / r)}:${Math.round(a.y / r)}`, n = (a, l) => Math.abs(a.x - l.x) <= r && Math.abs(a.y - l.y) <= r, c = /* @__PURE__ */ new Map();
+  o.forEach((a, l) => {
+    for (const f of [a.a, a.b]) {
+      const i = t(f);
+      c.has(i) || c.set(i, []), c.get(i).push(l);
+    }
+  });
+  const s = [];
+  for (let a = 0; a < o.length; a++) {
+    if (o[a].used) continue;
+    o[a].used = !0;
+    const l = [o[a].a, o[a].b];
+    let f = !0;
+    for (; f; ) {
+      f = !1;
+      const i = l[l.length - 1], u = c.get(t(i)) ?? [];
+      for (const h of u) {
+        if (o[h].used) continue;
+        const m = o[h];
+        if (n(m.a, i)) {
+          l.push(m.b), m.used = !0, f = !0;
+          break;
+        } else if (n(m.b, i)) {
+          l.push(m.a), m.used = !0, f = !0;
+          break;
+        }
+      }
+    }
+    for (f = !0; f; ) {
+      f = !1;
+      const i = l[0], u = c.get(t(i)) ?? [];
+      for (const h of u) {
+        if (o[h].used) continue;
+        const m = o[h];
+        if (n(m.b, i)) {
+          l.unshift(m.a), m.used = !0, f = !0;
+          break;
+        } else if (n(m.a, i)) {
+          l.unshift(m.b), m.used = !0, f = !0;
+          break;
+        }
+      }
+    }
+    s.push(l);
+  }
+  return s;
+}
+function wt(e) {
+  const r = Math.abs(e);
+  return r !== 0 && (r < 0.01 || r >= 1e5) ? e.toExponential(1) : e.toFixed(2).replace(/\.00$/, "");
+}
+function pt(e, r, o, t) {
+  const n = [], c = t(r);
+  for (const s of e) {
+    if (s.length < o) continue;
+    let a = 0;
+    for (let d = 1; d < s.length; d++)
+      a += Math.hypot(s[d].x - s[d - 1].x, s[d].y - s[d - 1].y);
+    const l = a / 2;
+    let f = 0, i = Math.floor(s.length / 2);
+    for (let d = 1; d < s.length; d++)
+      if (f += Math.hypot(s[d].x - s[d - 1].x, s[d].y - s[d - 1].y), f >= l) {
+        i = d;
+        break;
+      }
+    const u = s[i], h = Math.min(3, Math.floor(s.length / 4)) || 1, m = s[Math.max(0, i - h)], g = s[Math.min(s.length - 1, i + h)];
+    let M = Math.atan2(g.y - m.y, g.x - m.x);
+    M > Math.PI / 2 && (M -= Math.PI), M < -Math.PI / 2 && (M += Math.PI), n.push({ x: u.x, y: u.y, level: r, text: c, angle: M });
+  }
+  return n;
+}
+function yt(e, r, o) {
   const t = [], n = (r - e) / (o + 1);
-  for (let s = 1; s <= o; s++)
-    t.push(e + s * n);
+  for (let c = 1; c <= o; c++)
+    t.push(e + c * n);
   return t;
 }
-function wt(e, r, o, t, n, s, c, i, l, f) {
-  let a = 0;
-  if (o >= c && (a |= 1), t >= c && (a |= 2), n >= c && (a |= 4), s >= c && (a |= 8), a === 0 || a === 15) return;
-  const u = (M, d, w) => M + (d - M) * w, h = (M, d) => (c - M) / (d - M), m = (M) => {
+function bt(e, r, o, t, n, c, s, a, l, f) {
+  let i = 0;
+  if (o >= s && (i |= 1), t >= s && (i |= 2), n >= s && (i |= 4), c >= s && (i |= 8), i === 0 || i === 15) return;
+  const u = (M, d, w) => M + (d - M) * w, h = (M, d) => (s - M) / (d - M), m = (M) => {
     switch (M) {
       case 0:
-        return [u(i[e], i[e + 1], h(o, t)), l[r]];
+        return [u(a[e], a[e + 1], h(o, t)), l[r]];
       case 1:
-        return [i[e + 1], u(l[r], l[r + 1], h(t, n))];
+        return [a[e + 1], u(l[r], l[r + 1], h(t, n))];
       case 2:
-        return [u(i[e], i[e + 1], h(s, n)), l[r + 1]];
+        return [u(a[e], a[e + 1], h(c, n)), l[r + 1]];
       case 3:
-        return [i[e], u(l[r], l[r + 1], h(o, s))];
+        return [a[e], u(l[r], l[r + 1], h(o, c))];
       default:
         return [0, 0];
     }
@@ -600,7 +686,7 @@ function wt(e, r, o, t, n, s, c, i, l, f) {
     const w = m(M), p = m(d);
     f.push([w[0], w[1], p[0], p[1]]);
   };
-  switch (a) {
+  switch (i) {
     case 1:
     case 14:
       g(0, 3);
@@ -633,114 +719,114 @@ function wt(e, r, o, t, n, s, c, i, l, f) {
       break;
   }
 }
-function pt(e, r, o) {
-  const t = e.length, n = r.length, s = o ?? Math.min(t, n) - 1, c = 2 * s + 1, i = new Float32Array(c), l = new Float32Array(c);
-  let f = 0, a = 0;
+function Ft(e, r, o) {
+  const t = e.length, n = r.length, c = o ?? Math.min(t, n) - 1, s = 2 * c + 1, a = new Float32Array(s), l = new Float32Array(s);
+  let f = 0, i = 0;
   for (let d = 0; d < t; d++) f += e[d];
-  for (let d = 0; d < n; d++) a += r[d];
-  f /= t, a /= n;
+  for (let d = 0; d < n; d++) i += r[d];
+  f /= t, i /= n;
   let u = 0, h = 0;
   for (let d = 0; d < t; d++) u += (e[d] - f) ** 2;
-  for (let d = 0; d < n; d++) h += (r[d] - a) ** 2;
+  for (let d = 0; d < n; d++) h += (r[d] - i) ** 2;
   u = Math.sqrt(u / t), h = Math.sqrt(h / n);
   const m = u * h * Math.min(t, n);
   let g = -1 / 0, M = 0;
-  for (let d = 0; d < c; d++) {
-    const w = d - s;
+  for (let d = 0; d < s; d++) {
+    const w = d - c;
     l[d] = w;
     let p = 0;
     for (let y = 0; y < t; y++) {
-      const F = y + w;
-      F >= 0 && F < n && (p += (e[y] - f) * (r[F] - a));
+      const b = y + w;
+      b >= 0 && b < n && (p += (e[y] - f) * (r[b] - i));
     }
-    i[d] = m > 0 ? p / m : 0, i[d] > g && (g = i[d], M = d);
+    a[d] = m > 0 ? p / m : 0, a[d] > g && (g = a[d], M = d);
   }
   return {
-    correlation: i,
+    correlation: a,
     lags: l,
     maxCorrelation: g,
     lagAtMax: l[M]
   };
 }
-function Qt(e, r) {
-  return pt(e, e, r);
+function Vt(e, r) {
+  return Ft(e, e, r);
 }
-function Ut(e, r = {}) {
+function tn(e, r = {}) {
   const { method: o = "zscore", windowSize: t } = r;
   switch (o) {
     case "zscore":
-      return k(e, r.threshold ?? 3, t);
+      return $(e, r.threshold ?? 3, t);
     case "mad":
-      return yt(e, r.threshold ?? 2.5);
+      return At(e, r.threshold ?? 2.5);
     case "iqr":
-      return Ft(e, r.threshold ?? 1.5);
+      return xt(e, r.threshold ?? 1.5);
     case "isolation":
-      return At(e, r.threshold ?? 0.5);
+      return It(e, r.threshold ?? 0.5);
     default:
-      return k(e, 3, t);
+      return $(e, 3, t);
   }
 }
-function k(e, r, o) {
-  const t = e.length, n = new Float32Array(t), s = [];
+function $(e, r, o) {
+  const t = e.length, n = new Float32Array(t), c = [];
   if (o && o < t) {
-    const c = Math.floor(o / 2);
-    for (let i = 0; i < t; i++) {
-      const l = Math.max(0, i - c), f = Math.min(t, i + c + 1);
-      let a = 0;
-      for (let h = l; h < f; h++) a += e[h];
-      a /= f - l;
+    const s = Math.floor(o / 2);
+    for (let a = 0; a < t; a++) {
+      const l = Math.max(0, a - s), f = Math.min(t, a + s + 1);
+      let i = 0;
+      for (let h = l; h < f; h++) i += e[h];
+      i /= f - l;
       let u = 0;
-      for (let h = l; h < f; h++) u += (e[h] - a) ** 2;
-      u = Math.sqrt(u / (f - l)), n[i] = u > 0 ? Math.abs(e[i] - a) / u : 0, n[i] > r && s.push(i);
+      for (let h = l; h < f; h++) u += (e[h] - i) ** 2;
+      u = Math.sqrt(u / (f - l)), n[a] = u > 0 ? Math.abs(e[a] - i) / u : 0, n[a] > r && c.push(a);
     }
   } else {
-    let c = 0;
-    for (let l = 0; l < t; l++) c += e[l];
-    c /= t;
-    let i = 0;
-    for (let l = 0; l < t; l++) i += (e[l] - c) ** 2;
-    i = Math.sqrt(i / t);
+    let s = 0;
+    for (let l = 0; l < t; l++) s += e[l];
+    s /= t;
+    let a = 0;
+    for (let l = 0; l < t; l++) a += (e[l] - s) ** 2;
+    a = Math.sqrt(a / t);
     for (let l = 0; l < t; l++)
-      n[l] = i > 0 ? Math.abs(e[l] - c) / i : 0, n[l] > r && s.push(l);
+      n[l] = a > 0 ? Math.abs(e[l] - s) / a : 0, n[l] > r && c.push(l);
   }
-  return { indices: s, scores: n, threshold: r };
+  return { indices: c, scores: n, threshold: r };
 }
-function yt(e, r, o) {
-  const t = e.length, n = new Float32Array(t), s = [], c = Array.from(e).sort((u, h) => u - h), i = c[Math.floor(t / 2)], a = 1.4826 * c.map((u) => Math.abs(u - i)).sort((u, h) => u - h)[Math.floor(t / 2)];
+function At(e, r, o) {
+  const t = e.length, n = new Float32Array(t), c = [], s = Array.from(e).sort((u, h) => u - h), a = s[Math.floor(t / 2)], i = 1.4826 * s.map((u) => Math.abs(u - a)).sort((u, h) => u - h)[Math.floor(t / 2)];
   for (let u = 0; u < t; u++)
-    n[u] = a > 0 ? Math.abs(e[u] - i) / a : 0, n[u] > r && s.push(u);
-  return { indices: s, scores: n, threshold: r };
+    n[u] = i > 0 ? Math.abs(e[u] - a) / i : 0, n[u] > r && c.push(u);
+  return { indices: c, scores: n, threshold: r };
 }
-function Ft(e, r) {
-  const o = e.length, t = new Float32Array(o), n = [], s = Array.from(e).sort((u, h) => u - h), c = s[Math.floor(o * 0.25)], i = s[Math.floor(o * 0.75)], l = i - c, f = c - r * l, a = i + r * l;
+function xt(e, r) {
+  const o = e.length, t = new Float32Array(o), n = [], c = Array.from(e).sort((u, h) => u - h), s = c[Math.floor(o * 0.25)], a = c[Math.floor(o * 0.75)], l = a - s, f = s - r * l, i = a + r * l;
   for (let u = 0; u < o; u++)
-    e[u] < f ? (t[u] = (f - e[u]) / l, n.push(u)) : e[u] > a ? (t[u] = (e[u] - a) / l, n.push(u)) : t[u] = 0;
+    e[u] < f ? (t[u] = (f - e[u]) / l, n.push(u)) : e[u] > i ? (t[u] = (e[u] - i) / l, n.push(u)) : t[u] = 0;
   return { indices: n, scores: t, threshold: r };
 }
-function At(e, r, o) {
-  const t = e.length, n = new Float32Array(t), s = [], c = bt(e);
-  for (let i = 0; i < t; i++) {
+function It(e, r, o) {
+  const t = e.length, n = new Float32Array(t), c = [], s = qt(e);
+  for (let a = 0; a < t; a++) {
     let l = 0;
     for (let f = 0; f < t; f++)
-      i !== f && Math.abs(e[i] - e[f]) < c && l++;
-    n[i] = 1 - l / t, n[i] > r && s.push(i);
+      a !== f && Math.abs(e[a] - e[f]) < s && l++;
+    n[a] = 1 - l / t, n[a] > r && c.push(a);
   }
-  return { indices: s, scores: n, threshold: r };
+  return { indices: c, scores: n, threshold: r };
 }
-function bt(e) {
+function qt(e) {
   let r = 1 / 0, o = -1 / 0;
   for (let t = 0; t < e.length; t++)
     e[t] < r && (r = e[t]), e[t] > o && (o = e[t]);
   return (o - r) * 0.1;
 }
-function R(e, r) {
+function j(e, r) {
   const o = e.length;
   if (o < 2) return 0;
   let t = 0;
   if (typeof r == "number") {
     const n = r;
-    for (let s = 1; s < o; s++)
-      t += (e[s - 1] + e[s]) * n / 2;
+    for (let c = 1; c < o; c++)
+      t += (e[c - 1] + e[c]) * n / 2;
   } else if (r)
     for (let n = 1; n < o; n++)
       t += (e[n - 1] + e[n]) * (r[n] - r[n - 1]) / 2;
@@ -749,28 +835,28 @@ function R(e, r) {
       t += (e[n - 1] + e[n]) / 2;
   return t;
 }
-function It(e, r) {
+function Pt(e, r) {
   const o = e.length;
-  if (o < 3) return R(e, r);
+  if (o < 3) return j(e, r);
   const t = o % 2 === 0 ? o - 1 : o;
   let n = 0;
   if (typeof r == "number") {
-    const s = r;
+    const c = r;
     n = e[0] + e[t - 1];
-    for (let c = 1; c < t - 1; c++)
-      n += (c % 2 === 0 ? 2 : 4) * e[c];
-    n *= s / 3, o !== t && (n += (e[o - 2] + e[o - 1]) * s / 2);
-  } else return r ? R(e, r) : It(e, 1);
+    for (let s = 1; s < t - 1; s++)
+      n += (s % 2 === 0 ? 2 : 4) * e[s];
+    n *= c / 3, o !== t && (n += (e[o - 2] + e[o - 1]) * c / 2);
+  } else return r ? j(e, r) : Pt(e, 1);
   return n;
 }
-function Jt(e, r) {
+function nn(e, r) {
   const o = e.length, t = new Float32Array(o);
   if (o < 2)
     return t[0] = 0, t;
   if (t[0] = 0, typeof r == "number") {
     const n = r;
-    for (let s = 1; s < o; s++)
-      t[s] = t[s - 1] + (e[s - 1] + e[s]) * n / 2;
+    for (let c = 1; c < o; c++)
+      t[c] = t[c - 1] + (e[c - 1] + e[c]) * n / 2;
   } else if (r)
     for (let n = 1; n < o; n++)
       t[n] = t[n - 1] + (e[n - 1] + e[n]) * (r[n] - r[n - 1]) / 2;
@@ -779,17 +865,17 @@ function Jt(e, r) {
       t[n] = t[n - 1] + (e[n - 1] + e[n]) / 2;
   return t;
 }
-function Vt(e, r) {
+function en(e, r) {
   const o = e.length, t = r.length;
-  let n = 0, s = 0;
+  let n = 0, c = 0;
   for (let g = 0; g < o; g++) n += e[g];
-  for (let g = 0; g < t; g++) s += r[g];
-  n /= o, s /= t;
-  let c = 0, i = 0;
-  for (let g = 0; g < o; g++) c += (e[g] - n) ** 2;
-  for (let g = 0; g < t; g++) i += (r[g] - s) ** 2;
-  c /= o - 1, i /= t - 1;
-  const l = Math.sqrt(c / o + i / t), f = (n - s) / l, a = (c / o + i / t) ** 2, u = (c / o) ** 2 / (o - 1) + (i / t) ** 2 / (t - 1), h = a / u, m = 2 * (1 - qt(Math.abs(f)));
+  for (let g = 0; g < t; g++) c += r[g];
+  n /= o, c /= t;
+  let s = 0, a = 0;
+  for (let g = 0; g < o; g++) s += (e[g] - n) ** 2;
+  for (let g = 0; g < t; g++) a += (r[g] - c) ** 2;
+  s /= o - 1, a /= t - 1;
+  const l = Math.sqrt(s / o + a / t), f = (n - c) / l, i = (s / o + a / t) ** 2, u = (s / o) ** 2 / (o - 1) + (a / t) ** 2 / (t - 1), h = i / u, m = 2 * (1 - vt(Math.abs(f)));
   return {
     tStatistic: f,
     degreesOfFreedom: h,
@@ -797,23 +883,23 @@ function Vt(e, r) {
     significant: m < 0.05
   };
 }
-function qt(e) {
-  const r = 0.254829592, o = -0.284496736, t = 1.421413741, n = -1.453152027, s = 1.061405429, c = 0.3275911, i = e < 0 ? -1 : 1;
+function vt(e) {
+  const r = 0.254829592, o = -0.284496736, t = 1.421413741, n = -1.453152027, c = 1.061405429, s = 0.3275911, a = e < 0 ? -1 : 1;
   e = Math.abs(e) / Math.sqrt(2);
-  const l = 1 / (1 + c * e), f = 1 - ((((s * l + n) * l + t) * l + o) * l + r) * l * Math.exp(-e * e);
-  return 0.5 * (1 + i * f);
+  const l = 1 / (1 + s * e), f = 1 - ((((c * l + n) * l + t) * l + o) * l + r) * l * Math.exp(-e * e);
+  return 0.5 * (1 + a * f);
 }
-const xt = {
+const kt = {
   name: "velo-plot-analysis",
   version: "1.0.0",
   description: "Data analysis and signal processing for velo-plot",
   provides: ["analysis"],
   tags: ["fft", "filters", "statistics", "math"]
 };
-function tn(e = {}) {
+function on(e = {}) {
   let r;
   return {
-    manifest: xt,
+    manifest: kt,
     onInit(o) {
       r = o;
     },
@@ -829,7 +915,7 @@ function tn(e = {}) {
       calculateStats: G,
       validateData: Z,
       // Signal Processing
-      fft: (o) => x(o),
+      fft: (o) => q(o),
       detectPeaks: O,
       detectCycles: N,
       movingAverage: Y,
@@ -837,7 +923,7 @@ function tn(e = {}) {
       downsampleLTTB: H,
       singleFrequencyFilter: lt,
       // Fitting
-      fitData: L,
+      fitData: C,
       // Indicators
       sma: U,
       ema: Q,
@@ -849,79 +935,80 @@ function tn(e = {}) {
   };
 }
 export {
-  tn as PluginAnalysis,
-  fn as adx,
-  St as analyzeComplexSpectrum,
-  j as analyzeSpectrum,
-  un as aroon,
+  on as PluginAnalysis,
+  mn as adx,
+  Wt as analyzeComplexSpectrum,
+  L as analyzeSpectrum,
+  gn as aroon,
   et as arraysToComplex,
-  hn as atr,
-  Qt as autoCorrelation,
-  Ht as bandPassFilter,
-  Xt as bandStopFilter,
-  Lt as blackmanWindow,
-  mn as bollingerBands,
-  Yt as butterworth,
-  A as calculateR2,
+  dn as atr,
+  Vt as autoCorrelation,
+  Nt as bandPassFilter,
+  Ot as bandStopFilter,
+  Dt as blackmanWindow,
+  Mn as bollingerBands,
+  Zt as butterworth,
+  F as calculateR2,
   G as calculateStats,
-  Dt as complexToArrays,
-  pt as crossCorrelation,
-  gn as cumsum,
-  on as cumulativeIntegral,
-  Jt as cumulativeIntegral2,
-  tn as default,
-  dn as dema,
+  zt as complexToArrays,
+  Ft as crossCorrelation,
+  wn as cumsum,
+  cn as cumulativeIntegral,
+  nn as cumulativeIntegral2,
+  on as default,
+  pn as dema,
   _ as derivative,
-  Ut as detectAnomalies,
+  tn as detectAnomalies,
   N as detectCycles,
   O as detectPeaks,
-  Rt as dominantFrequency,
+  Ct as dominantFrequency,
   H as downsampleLTTB,
   Q as ema,
-  Nt as exponentialMovingAverage,
-  x as fft,
-  Bt as fftFromComplexInput,
-  L as fitData,
-  sn as formatScientific,
-  cn as formatValue,
+  Gt as exponentialMovingAverage,
+  q as fft,
+  _t as fftFromComplexInput,
+  C as fitData,
+  ln as formatScientific,
+  fn as formatValue,
   T as formatWithPrefix,
-  Ot as gaussianSmooth,
-  Kt as generateContours,
+  Kt as gaussianSmooth,
+  Jt as generateContours,
   E as generateCycleColors,
-  an as getBestPrefix,
-  zt as getPositiveFrequencies,
-  jt as hammingWindow,
-  Ct as hanningWindow,
-  Tt as highPassFilter,
+  un as getBestPrefix,
+  Ht as getPositiveFrequencies,
+  Bt as hammingWindow,
+  St as hanningWindow,
+  Yt as highPassFilter,
   V as ifft,
-  _t as ifftComplex,
-  Wt as ifftFromArrays,
+  Tt as ifftComplex,
+  Et as ifftFromArrays,
   z as integrate,
-  Et as lowPassFilter,
-  Mn as macd,
-  Gt as medianFilter,
-  wn as momentum,
+  Mt as joinSegments,
+  Xt as lowPassFilter,
+  yn as macd,
+  Ut as medianFilter,
+  bn as momentum,
   Y as movingAverage,
-  C as nextPowerOf2,
-  pn as normalize,
-  yn as obv,
-  Fn as percentChange,
-  kt as powerSpectrum,
-  An as roc,
+  R as nextPowerOf2,
+  Fn as normalize,
+  An as obv,
+  xn as percentChange,
+  Lt as powerSpectrum,
+  In as roc,
   K as rsi,
-  Zt as savitzkyGolay,
-  It as simpsonsIntegration,
+  Qt as savitzkyGolay,
+  Pt as simpsonsIntegration,
   lt as singleFrequencyFilter,
   U as sma,
   W as solveLinearSystem,
-  bn as standardDeviation,
-  In as stochastic,
+  qn as standardDeviation,
+  Pn as stochastic,
   X as subtractBaseline,
-  Vt as tTest,
-  qn as tema,
-  R as trapezoidalIntegration,
+  en as tTest,
+  vn as tema,
+  j as trapezoidalIntegration,
   Z as validateData,
-  xn as vwap,
-  Pn as wma
+  kn as vwap,
+  $n as wma
 };
 //# sourceMappingURL=analysis.js.map