pantograph2d 0.4.0 → 0.6.0

package/dist/Diagram-d848c815.js

@@ -0,0 +1,4252 @@
+function yt(i, t, n, e) {
+  return i <= e && t >= n;
+}
+class it {
+  constructor(t = 1 / 0, n = 1 / 0, e = -1 / 0, s = -1 / 0) {
+    this.xMin = t, this.yMin = n, this.xMax = e, this.yMax = s;
+  }
+  get width() {
+    return this.xMax - this.xMin;
+  }
+  get height() {
+    return this.yMax - this.yMin;
+  }
+  contains(t) {
+    const [n, e] = t;
+    return yt(this.xMin, this.xMax, n, n) && yt(this.yMin, this.yMax, e, e);
+  }
+  overlaps(t) {
+    return yt(this.xMin, this.xMax, t.xMin, t.xMax) && yt(this.yMin, this.yMax, t.yMin, t.yMax);
+  }
+  addPoint(t) {
+    const [n, e] = t;
+    return new it(
+      Math.min(this.xMin, n),
+      Math.min(this.yMin, e),
+      Math.max(this.xMax, n),
+      Math.max(this.yMax, e)
+    );
+  }
+  merge(t) {
+    return new it(
+      Math.min(this.xMin, t.xMin),
+      Math.min(this.yMin, t.yMin),
+      Math.max(this.xMax, t.xMax),
+      Math.max(this.yMax, t.yMax)
+    );
+  }
+}
+function De(i) {
+  return i.reduce((t, n) => t.addPoint(n), new it());
+}
+const Jt = (i, t = 1e-9) => {
+  let n = i;
+  return Math.abs(i) < t && (n = 0), n.toFixed(-Math.log10(t));
+};
+function ct(i, t = 1e-9) {
+  return Array.from(
+    new Map(
+      i.map(([n, e]) => [
+        `[${Jt(n, t)},${Jt(e, t)}]`,
+        [n, e]
+      ])
+    ).values()
+  );
+}
+const Lt = Math.PI / 180, Ue = 180 / Math.PI, T = (i) => `[${i[0]}, ${i[1]}]`, E = ([i, t], [n, e], s = 1e-9) => Math.abs(i - n) <= s && Math.abs(t - e) <= s, O = ([i, t], [n, e]) => [i + n, t + e], A = ([i, t], [n, e]) => [i - n, t - e], jt = ([i, t]) => i * i + t * t, D = ([i, t], n) => [i * n, t * n], mt = ([i, t], [n, e] = [0, 0]) => (i - n) ** 2 + (t - e) ** 2, q = (i, t = [0, 0]) => Math.sqrt(mt(i, t));
+function at([i, t], [n, e]) {
+  return i * e - t * n;
+}
+function $t([i, t], [n, e]) {
+  return i * n + t * e;
+}
+function W([i, t]) {
+  const n = q([i, t]);
+  return [i / n, t / n];
+}
+function ft(i, t) {
+  const n = Math.cos(t) * i, e = Math.sin(t) * i;
+  return [n, e];
+}
+function Ye([i, t]) {
+  return Math.atan2(t, i);
+}
+function Xe(i) {
+  const t = q(i), n = Ye(i);
+  return [t, n];
+}
+function Vt(i, t, n = 1e-9) {
+  const e = at(i, t), s = jt(i), r = jt(t);
+  return e * e < s * r * n * n;
+}
+function J(i) {
+  return [-i[1], i[0]];
+}
+function bt(i) {
+  return [i[1], -i[0]];
+}
+function Ge(i, t, n) {
+  return E(n, i) || Vt(A(n, i), t);
+}
+const ot = (i, t) => {
+  const [n, e, s, r, o, h, u, l, c] = i, [a, g, w, m, d, P, f, p, y] = t;
+  return [
+    n * a + e * m + s * f,
+    n * g + e * d + s * p,
+    n * w + e * P + s * y,
+    r * a + o * m + h * f,
+    r * g + o * d + h * p,
+    r * w + o * P + h * y,
+    u * a + l * m + c * f,
+    u * g + l * d + c * p,
+    u * w + l * P + c * y
+  ];
+}, We = (i) => {
+  const [t, n, e, s, r, o, h, u, l] = i, c = t * (r * l - o * u) - n * (s * l - o * h) + e * (s * u - r * h);
+  return [
+    (r * l - o * u) / c,
+    (e * u - n * l) / c,
+    (n * o - e * r) / c,
+    (o * h - s * l) / c,
+    (t * l - e * h) / c,
+    (e * s - t * o) / c,
+    (s * u - r * h) / c,
+    (n * h - t * u) / c,
+    (t * r - n * s) / c
+  ];
+}, Je = (i) => {
+  const [t, n, e, s, r, o, h, u, l] = i;
+  return [t, s, h, n, r, u, e, o, l];
+};
+class V {
+  constructor(t) {
+    this._matrix = [1, 0, 0, 0, 1, 0, 0, 0, 1], t && (this._matrix = [...t]);
+  }
+  clone() {
+    return new V(this._matrix);
+  }
+  transpose() {
+    return this._matrix = Je(this._matrix), this;
+  }
+  inverse() {
+    return this._matrix = We(this._matrix), this;
+  }
+  translate(t, n) {
+    return this._matrix = ot(this._matrix, [1, 0, t, 0, 1, n, 0, 0, 1]), this;
+  }
+  rotate(t, n) {
+    const e = Math.cos(t), s = Math.sin(t), r = [e, -s, 0, s, e, 0, 0, 0, 1];
+    return n && this.translate(n[0], n[1]), this._matrix = ot(this._matrix, r), n && this.translate(-n[0], -n[1]), this;
+  }
+  mirrorX() {
+    return this._matrix = ot(this._matrix, [1, 0, 0, 0, -1, 0, 0, 0, 1]), this;
+  }
+  mirrorY() {
+    return this._matrix = ot(this._matrix, [-1, 0, 0, 0, 1, 0, 0, 0, 1]), this;
+  }
+  mirrorLine(t, n) {
+    const [e, s] = t, r = Math.atan2(s, e);
+    return n && this.translate(n[0], n[1]), this.rotate(r), this.mirrorX(), this.rotate(-r), n && this.translate(-n[0], -n[1]), this;
+  }
+  mirrorCenter(t) {
+    return t && this.translate(t[0], t[1]), this._matrix = ot(this._matrix, [-1, 0, 0, 0, -1, 0, 0, 0, 1]), t && this.translate(-t[0], -t[1]), this;
+  }
+  scale(t, n) {
+    return n && this.translate(n[0], n[1]), this._matrix = ot(this._matrix, [t, 0, 0, 0, t, 0, 0, 0, 1]), n && this.translate(-n[0], -n[1]), this;
+  }
+  transform(t) {
+    const [n, e] = t, [s, r, o, h, u, l] = this._matrix;
+    return [s * n + r * e + o, h * n + u * e + l];
+  }
+  transformAngle(t) {
+    const [n, e] = this.transform([Math.cos(t), Math.sin(t)]), [s, r] = this.transform([0, 0]);
+    return Math.atan2(e - r, n - s);
+  }
+  keepsOrientation() {
+    const [t, , , , n] = this._matrix;
+    return t * n > 0;
+  }
+  scaleFactor() {
+    const [t, , , n] = this._matrix;
+    return Math.sqrt(t * t + n * n);
+  }
+}
+class _t {
+  translateX(t) {
+    const n = new V().translate(t, 0);
+    return this.transform(n);
+  }
+  translateY(t) {
+    const n = new V().translate(0, t);
+    return this.transform(n);
+  }
+  translate(t, n) {
+    const e = new V().translate(t, n);
+    return this.transform(e);
+  }
+  translateTo([t, n]) {
+    const e = new V().translate(t, n);
+    return this.transform(e);
+  }
+  rotate(t, n) {
+    const e = new V().rotate(
+      t * Lt,
+      n
+    );
+    return this.transform(e);
+  }
+  scale(t, n) {
+    const e = new V().scale(t, n);
+    return this.transform(e);
+  }
+  mirrorCenter(t) {
+    const n = new V().mirrorCenter(t);
+    return this.transform(n);
+  }
+  mirror(t = "x", n) {
+    const e = new V();
+    return t === "x" ? e.mirrorX() : t === "y" ? e.mirrorY() : e.mirrorLine(t, n), this.transform(e);
+  }
+}
+class zt extends _t {
+  constructor(t, n) {
+    super(), this.firstPoint = t, this.lastPoint = n, this.precision = 1e-9, this.firstPoint = t, this.lastPoint = n;
+  }
+  get repr() {
+    return `${this.segmentType} ${T(this.firstPoint)} - ${T(
+      this.lastPoint
+    )}`;
+  }
+  get info() {
+    return this.repr;
+  }
+  [Symbol.for("nodejs.util.inspect.custom")]() {
+    return this.repr;
+  }
+}
+class R extends zt {
+  constructor() {
+    super(...arguments), this.segmentType = "LINE", this._V = null, this._slope = null, this._yIntercept = null, this._boundingBox = null;
+  }
+  isValidParameter(t) {
+    const n = this.length * this.precision;
+    return t >= -n && 1 - t >= -n;
+  }
+  paramPoint(t) {
+    return O(this.firstPoint, D(this.V, t));
+  }
+  get length() {
+    return q(this.firstPoint, this.lastPoint);
+  }
+  get squareLength() {
+    return mt(this.firstPoint, this.lastPoint);
+  }
+  get V() {
+    return this._V === null && (this._V = A(this.lastPoint, this.firstPoint)), this._V;
+  }
+  get slope() {
+    if (this._slope === null) {
+      const [t, n] = this.V;
+      this._slope = n / t;
+    }
+    return this._slope;
+  }
+  get yIntercept() {
+    return this._yIntercept === null && (this._yIntercept = this.firstPoint[1] - this.slope * this.firstPoint[0]), this._yIntercept;
+  }
+  get midPoint() {
+    return O(this.firstPoint, D(this.V, 0.5));
+  }
+  isSame(t) {
+    return t instanceof R ? E(this.firstPoint, t.firstPoint) && E(this.lastPoint, t.lastPoint) || E(this.lastPoint, t.firstPoint) && E(this.firstPoint, t.lastPoint) : !1;
+  }
+  clone() {
+    return new R(this.firstPoint, this.lastPoint);
+  }
+  reverse() {
+    return new R(this.lastPoint, this.firstPoint);
+  }
+  get boundingBox() {
+    return this._boundingBox === null && (this._boundingBox = new it(
+      Math.min(this.firstPoint[0], this.lastPoint[0]) - this.precision,
+      Math.min(this.firstPoint[1], this.lastPoint[1]) - this.precision,
+      Math.max(this.firstPoint[0], this.lastPoint[0]) + this.precision,
+      Math.max(this.firstPoint[1], this.lastPoint[1]) + this.precision
+    )), this._boundingBox;
+  }
+  distanceFrom(t) {
+    const n = A(t, this.firstPoint), e = $t(n, this.V) / this.squareLength;
+    if (e < 0)
+      return q(t, this.firstPoint);
+    if (e > 1)
+      return q(t, this.lastPoint);
+    const s = this.paramPoint(e);
+    return q(t, s);
+  }
+  isOnSegment(t) {
+    if (E(t, this.firstPoint, this.precision))
+      return !0;
+    const n = A(t, this.firstPoint);
+    if (!Vt(this.V, n))
+      return !1;
+    const e = $t(n, this.V) / this.squareLength;
+    return this.isValidParameter(e);
+  }
+  gradientAt(t) {
+    return this.V;
+  }
+  tangentAt(t) {
+    if (!this.isOnSegment(t))
+      throw new Error("Point is not on segment");
+    return W(this.V);
+  }
+  get normalVector() {
+    return J(W(this.V));
+  }
+  get tangentAtFirstPoint() {
+    return W(this.V);
+  }
+  get tangentAtLastPoint() {
+    return W(this.V);
+  }
+  splitAt(t) {
+    let n;
+    if (Array.isArray(t) && t.length === 0)
+      return [this];
+    Array.isArray(t[0]) ? n = t : n = [t], n.forEach((u) => {
+      if (!this.isOnSegment(u))
+        throw new Error(
+          `Point ${T(u)} is not on segment ${this.repr}`
+        );
+    });
+    const e = [this.firstPoint, ...n, this.lastPoint], s = ct(e), r = this.lastPoint[0] - this.firstPoint[0];
+    let o = Math.sign(r), h = 0;
+    return Math.abs(r) < this.precision && (o = Math.sign(this.lastPoint[1] - this.firstPoint[1]), h = 1), s.sort(
+      (u, l) => o * (u[h] - l[h])
+    ), s.flatMap((u, l) => l === s.length - 1 ? [] : new R(u, s[l + 1]));
+  }
+  transform(t) {
+    return new R(
+      t.transform(this.firstPoint),
+      t.transform(this.lastPoint)
+    );
+  }
+}
+function Ke(i) {
+  return Array.from(Array(i).keys());
+}
+function nt(i) {
+  const t = Math.min(...i.map((n) => n.length));
+  return Ke(t).map((n) => i.map((e) => e[n]));
+}
+function st(i, t = 1e-9) {
+  return i < 0 ? i + 2 * Math.PI : i >= 2 * Math.PI ? i % (2 * Math.PI) : i > 2 * Math.PI - t ? 0 : i;
+}
+function vt(i, t, n, e = 1e-9) {
+  let s = t - i;
+  return n && (s = -s), s < 0 && (s += 2 * Math.PI), s > 2 * Math.PI - e ? 0 : s;
+}
+const Ft = (i, t, n) => {
+  const e = at(i.V, t.V), s = jt(i.V), r = jt(t.V), o = n ? n * n : i.precision * t.precision;
+  if (e * e < s * r * o)
+    return "parallel";
+  const h = A(t.firstPoint, i.firstPoint), u = at(h, t.V) / e, l = at(h, i.V) / e;
+  return {
+    intersectionParam1: u,
+    intersectionParam2: l
+  };
+};
+function ge(i, t, n = !1, e) {
+  const s = Ft(i, t, e);
+  if (s === "parallel") {
+    if (!n)
+      return null;
+    if (i.isSame(t))
+      return i;
+    const h = ct(
+      [
+        t.isOnSegment(i.firstPoint) ? i.firstPoint : null,
+        t.isOnSegment(i.lastPoint) ? i.lastPoint : null,
+        i.isOnSegment(t.firstPoint) ? t.firstPoint : null,
+        i.isOnSegment(t.lastPoint) ? t.lastPoint : null
+      ].filter((u) => u !== null)
+    ).sort((u, l) => u[0] - l[0]);
+    if (h.length === 0)
+      return null;
+    if (h.length === 1)
+      return null;
+    if (h.length === 2)
+      return new R(h[0], h[1]);
+    throw console.error(h), new Error(
+      "Unexpected number of points while intersecting parallel lines"
+    );
+  }
+  const { intersectionParam1: r, intersectionParam2: o } = s;
+  return !i.isValidParameter(r) || !t.isValidParameter(o) ? null : i.paramPoint(r);
+}
+const Pt = (i, t) => {
+  const n = A(i, t);
+  return Xe(n);
+};
+class I extends zt {
+  constructor(t, n, e, s = !1, { ignoreChecks: r = !1 } = {}) {
+    if (super(t, n), this.segmentType = "ARC", this._coefficients = null, this._angularLength = null, this._radius = null, this._firstAngle = null, this._lastAngle = null, this._boundingBox = null, this.center = e, this.clockwise = s, !r) {
+      if (E(t, n))
+        throw new Error("Invalid arc, cannot be a full circle");
+      if (Math.abs(this.radius - q(this.lastPoint, this.center)) > this.precision)
+        throw new Error(
+          `Invalid arc, radius does not match between ${T(
+            t
+          )} and ${T(n)}} (center ${T(e)})`
+        );
+    }
+  }
+  get info() {
+    return `ARC(${T(this.firstPoint)}, ${T(
+      this.lastPoint
+    )}, ${T(this.center)}, ${this.clockwise ? "CW" : "CCW"})`;
+  }
+  get coefficients() {
+    if (this._coefficients === null) {
+      const [t, n] = this.center, e = this.radius * this.radius;
+      this._coefficients = {
+        x2: 1 / e,
+        xy: 0,
+        y2: 1 / e,
+        x: -(2 * t) / e,
+        y: -(2 * n) / e,
+        c: (t * t + n * n - e) / e
+      };
+    }
+    return this._coefficients;
+  }
+  isValidParameter(t) {
+    return 1 - t >= -this.precision && t >= -this.precision;
+  }
+  angleToParam(t) {
+    return vt(this.firstAngle, st(t), this.clockwise) / this.angularLength;
+  }
+  get angularLength() {
+    return this._angularLength || (this._angularLength = vt(
+      this.firstAngle,
+      this.lastAngle,
+      this.clockwise
+    )), this._angularLength;
+  }
+  paramPoint(t) {
+    return O(
+      this.center,
+      ft(
+        this.radius,
+        this.firstAngle + t * this.angularLength * (this.clockwise ? -1 : 1)
+      )
+    );
+  }
+  pointToParam(t) {
+    const [n, e] = Pt(t, this.center);
+    if (Math.abs(n - this.radius) > this.precision)
+      throw new Error(
+        `Point ${T(t)} is not on segment ${this.repr}`
+      );
+    const s = this.angleToParam(e);
+    if (!this.isValidParameter(s))
+      throw new Error(
+        `Point ${T(t)} is not on segment ${this.repr}`
+      );
+    return s;
+  }
+  get radius() {
+    return this._radius === null && (this._radius = q(this.firstPoint, this.center)), this._radius;
+  }
+  get firstAngle() {
+    if (this._firstAngle === null) {
+      const [t, n] = A(this.firstPoint, this.center);
+      this._firstAngle = st(Math.atan2(n, t));
+    }
+    return this._firstAngle;
+  }
+  get lastAngle() {
+    if (this._lastAngle === null) {
+      const [t, n] = A(this.lastPoint, this.center);
+      this._lastAngle = st(Math.atan2(n, t));
+    }
+    return this._lastAngle;
+  }
+  get length() {
+    return this.radius * this.angularLength;
+  }
+  get squareLength() {
+    return this.length * this.length;
+  }
+  get midPoint() {
+    return this.paramPoint(0.5);
+  }
+  isSame(t) {
+    return !(t instanceof I) || !E(this.center, t.center) ? !1 : E(this.firstPoint, t.firstPoint) && E(this.lastPoint, t.lastPoint) && this.clockwise === t.clockwise || E(this.lastPoint, t.firstPoint) && E(this.firstPoint, t.lastPoint) && this.clockwise === !t.clockwise;
+  }
+  clone() {
+    return new I(
+      this.firstPoint,
+      this.lastPoint,
+      this.center,
+      this.clockwise
+    );
+  }
+  reverse() {
+    return new I(
+      this.lastPoint,
+      this.firstPoint,
+      this.center,
+      !this.clockwise
+    );
+  }
+  get boundingBox() {
+    if (this._boundingBox === null) {
+      const t = this.radius + this.precision, n = (e) => this.isValidParameter(this.angleToParam(e));
+      this._boundingBox = new it(
+        n(Math.PI) ? this.center[0] - t : Math.min(this.firstPoint[0], this.lastPoint[0]) - this.precision,
+        n(Math.PI * 1.5) ? this.center[1] - t : Math.min(this.firstPoint[1], this.lastPoint[1]) - this.precision,
+        n(0) ? this.center[0] + t : Math.max(this.firstPoint[0], this.lastPoint[0]) + this.precision,
+        n(Math.PI / 2) ? this.center[1] + t : Math.max(this.firstPoint[1], this.lastPoint[1]) + this.precision
+      );
+    }
+    return this._boundingBox;
+  }
+  distanceFrom(t) {
+    const [n, e] = Pt(t, this.center);
+    return this.isValidParameter(this.angleToParam(e)) ? Math.abs(n - this.radius) : Math.sqrt(
+      Math.min(
+        mt(t, this.firstPoint),
+        mt(t, this.lastPoint)
+      )
+    );
+  }
+  isOnSegment(t) {
+    if (E(t, this.firstPoint) || E(t, this.lastPoint))
+      return !0;
+    const [n, e] = Pt(t, this.center);
+    if (Math.abs(n - this.radius) > this.precision)
+      return !1;
+    const s = this.angleToParam(e);
+    return this.isValidParameter(s);
+  }
+  gradientAt(t) {
+    const n = this.firstAngle + t * this.angularLength * (this.clockwise ? -1 : 1), e = this.radius * this.angularLength, s = -e * Math.sin(n), r = e * Math.cos(n);
+    return this.clockwise ? [-s, -r] : [s, r];
+  }
+  tangentAt(t) {
+    const [n, e] = Pt(t, this.center);
+    if (Math.abs(n - this.radius) > this.precision)
+      throw new Error("Point is not on the arc");
+    const s = this.angleToParam(e);
+    if (!this.isValidParameter(s))
+      throw new Error("Point is not on the arc");
+    const r = ft(1, e);
+    return (this.clockwise ? bt : J)(W(r));
+  }
+  get tangentAtFirstPoint() {
+    const t = ft(1, this.firstAngle);
+    return (this.clockwise ? bt : J)(W(t));
+  }
+  get tangentAtLastPoint() {
+    const t = ft(1, this.lastAngle);
+    return (this.clockwise ? bt : J)(W(t));
+  }
+  splitAt(t) {
+    let n;
+    if (Array.isArray(t) && t.length === 0)
+      return [this];
+    Array.isArray(t[0]) ? n = t : n = [t];
+    const s = [0, 1, ...n.map((h) => this.pointToParam(h))], r = new Map(
+      nt([s, [this.firstPoint, this.lastPoint, ...n]])
+    );
+    s.sort((h, u) => h - u);
+    let o = null;
+    return s.flatMap((h, u) => {
+      if (u === s.length - 1)
+        return [];
+      const l = s[u + 1];
+      if (l - h < this.precision)
+        return o === null && (o = h), [];
+      const c = o === null ? h : o, a = new I(
+        r.get(c) || this.paramPoint(c),
+        r.get(l) || this.paramPoint(l),
+        this.center,
+        this.clockwise
+      );
+      return o = null, a;
+    });
+  }
+  transform(t) {
+    return new I(
+      t.transform(this.firstPoint),
+      t.transform(this.lastPoint),
+      t.transform(this.center),
+      t.keepsOrientation() ? this.clockwise : !this.clockwise
+    );
+  }
+}
+function as(i, t, n) {
+  const e = new R(t, i), s = new R(t, n), r = J(e.tangentAtFirstPoint), o = J(s.tangentAtLastPoint), h = Ft(
+    { firstPoint: e.midPoint, V: r, precision: 1e-9 },
+    { firstPoint: s.midPoint, V: o, precision: 1e-9 }
+  );
+  if (h === "parallel")
+    throw new Error("Cannot create an arc from three colinear points");
+  const u = at(
+    A(i, t),
+    A(n, t)
+  ) > 0;
+  return new I(
+    i,
+    n,
+    O(e.midPoint, D(r, h.intersectionParam1)),
+    u,
+    { ignoreChecks: !0 }
+  );
+}
+function cs(i, t, n) {
+  const e = new R(t, i), s = J(e.tangentAtFirstPoint), r = Ft(
+    { firstPoint: e.midPoint, V: s, precision: 1e-9 },
+    {
+      firstPoint: i,
+      V: J(n),
+      precision: 1e-9
+    }
+  );
+  if (r === "parallel")
+    throw new Error("Cannot create an arc from three colinear points");
+  const o = O(
+    e.midPoint,
+    D(s, r.intersectionParam1)
+  ), h = at(
+    A(o, i),
+    A(o, O(i, n))
+  ) < 0;
+  return new I(i, t, o, h, {
+    ignoreChecks: !0
+  });
+}
+function Qe(i, t) {
+  const n = A(t, i.firstPoint), e = $t(n, i.V) / i.squareLength;
+  return i.paramPoint(e);
+}
+function kt(i, t, n) {
+  const e = n || i.precision, s = Qe(i, t.center), r = q(s, t.center);
+  if (r > t.radius + e)
+    return [];
+  if (Math.abs(r - t.radius) < e) {
+    const a = s;
+    return i.isOnSegment(a) && t.isOnSegment(a) ? [a] : [];
+  }
+  const o = [], h = Math.sqrt(
+    t.radius * t.radius - r * r
+  ), u = i.tangentAtFirstPoint, l = O(s, D(u, h));
+  i.isOnSegment(l) && t.isOnSegment(l) && o.push(l);
+  const c = O(s, D(u, -h));
+  return i.isOnSegment(c) && t.isOnSegment(c) && o.push(c), o;
+}
+const Ze = (i) => {
+  const { firstPoint: t, lastPoint: n, center: e, clockwise: s } = i;
+  return new I(n, t, e, s, {
+    ignoreChecks: !0
+  });
+}, He = (i, t) => {
+  if (i.isSame(t))
+    return [i];
+  const n = ct(
+    [
+      t.isOnSegment(i.firstPoint) ? i.firstPoint : null,
+      t.isOnSegment(i.lastPoint) ? i.lastPoint : null,
+      i.isOnSegment(t.firstPoint) ? t.firstPoint : null,
+      i.isOnSegment(t.lastPoint) ? t.lastPoint : null
+    ].filter((e) => e !== null)
+    // We sort by the param value of the first arc. This means that the points
+    // will be sorted with the same orientation than arc1
+  ).sort((e, s) => i.pointToParam(e) - i.pointToParam(s));
+  if (n.length === 0)
+    return [];
+  if (n.length === 1)
+    return [];
+  if (n.length === 2)
+    return i.isSame(Ze(t)) ? [] : [new I(n[0], n[1], i.center, i.clockwise)];
+  if (n.length === 3) {
+    const e = E(n[0], t.lastPoint) || E(n[0], t.firstPoint) ? 1 : 0;
+    return [
+      new I(
+        n[0 + e],
+        n[1 + e],
+        i.center,
+        i.clockwise
+      )
+    ];
+  } else if (n.length === 4)
+    return [
+      new I(n[0], n[1], i.center, i.clockwise),
+      new I(n[2], n[3], i.center, i.clockwise)
+    ];
+  throw new Error("Bug in the arc arc overlap algorithm");
+};
+function me(i, t, n = !1, e) {
+  const s = e || i.precision, r = q(i.center, t.center), o = i.radius + t.radius;
+  if (r > o + s)
+    return [];
+  const h = Math.abs(i.radius - t.radius);
+  if (r < h - s)
+    return [];
+  if (r < s)
+    return h > s ? [] : n ? He(i, t) : [];
+  const u = W(A(t.center, i.center)), l = r > o - s;
+  if (
+    // circles are outside each other
+    l || // circles are inside each other
+    Math.abs(r - h) < s
+  ) {
+    const f = l || i.radius > t.radius ? 1 : -1, p = O(
+      i.center,
+      D(u, f * i.radius)
+    );
+    return i.isOnSegment(p) && t.isOnSegment(p) ? [p] : [];
+  }
+  const c = i.radius * i.radius / (2 * r) - t.radius * t.radius / (2 * r) + r / 2, a = O(
+    i.center,
+    D(u, c)
+  ), g = Math.sqrt(
+    i.radius * i.radius - c * c
+  ), w = J(u), m = O(a, D(w, g)), d = O(a, D(w, -g)), P = [];
+  return i.isOnSegment(m) && t.isOnSegment(m) && P.push(m), i.isOnSegment(d) && t.isOnSegment(d) && P.push(d), P;
+}
+class C extends zt {
+  constructor(t, n, e, s, r, o, h = !1, {
+    ignoreChecks: u = !1,
+    angleUnits: l = "deg"
+  } = {}) {
+    super(t, n), this.segmentType = "ELLIPSE_ARC", this.precision = 1e-6, this._coefficients = null, this.center = e;
+    const c = s >= r;
+    this.majorRadius = c ? s : r, this.minorRadius = c ? r : s;
+    const a = l === "deg" ? o * Lt : o;
+    if (this.tiltAngle = st(
+      c ? a : a + Math.PI / 2
+    ), this.clockwise = h, !u) {
+      if (E(t, n))
+        throw new Error("Invalid arc, cannot be a full circle");
+      if (!this.isPointOnEllipse(t))
+        throw new Error(
+          `First point ${T(
+            t
+          )} not on the ellipse defined by ${this.info}`
+        );
+      if (!this.isPointOnEllipse(n))
+        throw new Error(
+          `Last point ${T(n)} not on the ellipse defined by ${this.info}`
+        );
+      if (Math.abs(this.majorRadius - this.minorRadius) < this.precision)
+        throw new Error(
+          "Both radii should be different, create an arc instead"
+        );
+    }
+  }
+  get coefficients() {
+    if (this._coefficients === null) {
+      const [t, n] = this.center, e = this.majorRadius * this.majorRadius, s = this.minorRadius * this.minorRadius, r = Math.cos(this.tiltAngle), o = Math.sin(this.tiltAngle), h = e * o * o + s * r * r, u = 2 * (s - e) * o * r, l = e * r * r + s * o * o;
+      this._coefficients = {
+        x2: h,
+        xy: u,
+        y2: l,
+        x: -2 * h * t - u * n,
+        y: -2 * l * n - u * t,
+        c: h * t * t + u * t * n + l * n * n - e * s
+      };
+    }
+    return this._coefficients;
+  }
+  get info() {
+    return `ELLIPSE_ARC(${T(this.firstPoint)}, ${T(
+      this.lastPoint
+    )}, ${T(this.center)}, ${this.majorRadius}, ${this.minorRadius}, ${this.tiltAngle * Ue}, ${this.clockwise ? "CW" : "CCW"})`;
+  }
+  reverse() {
+    return new C(
+      this.lastPoint,
+      this.firstPoint,
+      this.center,
+      this.majorRadius,
+      this.minorRadius,
+      this.tiltAngle,
+      !this.clockwise,
+      { ignoreChecks: !0, angleUnits: "rad" }
+    );
+  }
+  clone() {
+    return new C(
+      this.firstPoint,
+      this.lastPoint,
+      this.center,
+      this.majorRadius,
+      this.minorRadius,
+      this.tiltAngle,
+      this.clockwise,
+      { ignoreChecks: !0, angleUnits: "rad" }
+    );
+  }
+  computeEllipseExtrema() {
+    const t = Math.tan(this.tiltAngle), n = this.minorRadius / this.majorRadius, e = -Math.atan(n * t), s = Math.atan(n / t);
+    return [e, Math.PI + e, s, Math.PI + s].map((o) => this.thetaToParam(o)).filter((o) => this.isValidParameter(o)).map((o) => this.paramPoint(o));
+  }
+  get boundingBox() {
+    if (this._boundingBox === void 0) {
+      const t = this.computeEllipseExtrema(), n = [
+        this.firstPoint,
+        this.lastPoint,
+        ...t
+      ];
+      this._boundingBox = De(n);
+    }
+    return this._boundingBox;
+  }
+  get linearExentricity() {
+    return this._linearExentricity === void 0 && (this._linearExentricity = Math.sqrt(
+      this.majorRadius * this.majorRadius - this.minorRadius * this.minorRadius
+    )), this._linearExentricity;
+  }
+  get exentricity() {
+    return this._exentricity === void 0 && (this._exentricity = this.linearExentricity / this.majorRadius), this._exentricity;
+  }
+  get focals() {
+    if (this._focals === void 0) {
+      const t = this.majorAxis;
+      this._focals = [
+        O(this.center, D(t, this.linearExentricity)),
+        O(this.center, D(t, -this.linearExentricity))
+      ];
+    }
+    return this._focals;
+  }
+  get majorAxis() {
+    return ft(1, this.tiltAngle);
+  }
+  paramPoint(t) {
+    const n = this.firstAngle + t * this.deltaAngle * (this.clockwise ? -1 : 1);
+    return this.reverseEllipseReferenceFrameTransform.transform([
+      this.majorRadius * Math.cos(n),
+      this.minorRadius * Math.sin(n)
+    ]);
+  }
+  pointToParam(t) {
+    if (!this.isPointOnEllipse(t))
+      throw new Error(
+        `Point ${T(t)} not on the ellipse defined by ${this.repr}`
+      );
+    const n = this.thetaToParam(this.pointTheta(t));
+    if (!this.isValidParameter(n))
+      throw new Error(
+        `Point ${T(t)} is not on segment ${this.repr}`
+      );
+    return n;
+  }
+  get midPoint() {
+    return this.paramPoint(0.5);
+  }
+  isValidParameter(t) {
+    return 1 - t >= -this.precision && t >= -this.precision;
+  }
+  isSame(t) {
+    return E(this.center, t.center) && Math.abs(this.majorRadius - t.majorRadius) < this.precision && Math.abs(this.minorRadius - t.minorRadius) < this.precision && (Math.abs(this.tiltAngle - t.tiltAngle) < this.precision || Math.abs(Math.abs(this.tiltAngle - t.tiltAngle) - Math.PI) < this.precision) && (E(this.firstPoint, t.firstPoint) && this.clockwise === t.clockwise || E(this.firstPoint, t.lastPoint) && this.clockwise !== t.clockwise);
+  }
+  pointTheta(t) {
+    const n = this.ellipseReferenceFrameTransform.transform(t), e = Math.atan2(
+      n[1] / this.minorRadius,
+      n[0] / this.majorRadius
+    );
+    return st(e);
+  }
+  thetaToParam(t) {
+    return vt(this.firstAngle, st(t), this.clockwise) / this.deltaAngle;
+  }
+  isPointOnEllipse(t) {
+    const [n, e] = this.focals, s = q(t, n), r = q(t, e);
+    return Math.abs(2 * this.majorRadius - s - r) < this.precision;
+  }
+  isOnSegment(t) {
+    return this.isPointOnEllipse(t) ? this.isValidParameter(this.thetaToParam(this.pointTheta(t))) : !1;
+  }
+  distanceFrom(t) {
+    let n;
+    if (E(t, this.center) ? n = O(
+      this.center,
+      D(J(this.majorAxis), this.minorRadius)
+    ) : n = this.reverseEllipseReferenceFrameTransform.transform(
+      en(
+        this.majorRadius,
+        this.minorRadius,
+        this.ellipseReferenceFrameTransform.transform(t)
+      )
+    ), this.isValidParameter(this.thetaToParam(this.pointTheta(t))))
+      return q(t, n);
+    if (Ge(t, this.majorAxis, this.center)) {
+      const e = st(
+        2 * Math.PI - this.pointTheta(t)
+      ), s = this.thetaToParam(e);
+      if (this.isValidParameter(s))
+        return q(t, this.paramPoint(s));
+    }
+    return Math.min(
+      q(t, this.firstPoint),
+      q(t, this.lastPoint)
+    );
+  }
+  get ellipseReferenceFrameTransform() {
+    return this._ellipseReferenceFrameTransform === void 0 && (this._ellipseReferenceFrameTransform = new V().rotate(-this.tiltAngle).translate(-this.center[0], -this.center[1])), this._ellipseReferenceFrameTransform;
+  }
+  get reverseEllipseReferenceFrameTransform() {
+    return this._reverseEllipseReferenceFrameTransform === void 0 && (this._reverseEllipseReferenceFrameTransform = new V().translate(this.center[0], this.center[1]).rotate(this.tiltAngle)), this._reverseEllipseReferenceFrameTransform;
+  }
+  get rotateFromEllipseReferenceFrame() {
+    return this._rotateFromEllipseReferenceFrame === void 0 && (this._rotateFromEllipseReferenceFrame = new V().rotate(
+      this.tiltAngle
+    )), this._rotateFromEllipseReferenceFrame;
+  }
+  get firstAngle() {
+    return this._firstAngle === void 0 && (this._firstAngle = this.pointTheta(this.firstPoint)), this._firstAngle;
+  }
+  get lastAngle() {
+    return this._lastAngle === void 0 && (this._lastAngle = this.pointTheta(this.lastPoint)), this._lastAngle;
+  }
+  get deltaAngle() {
+    return this._deltaAngle === void 0 && (this._deltaAngle = vt(
+      this.firstAngle,
+      this.lastAngle,
+      this.clockwise
+    )), this._deltaAngle;
+  }
+  normalAt(t) {
+    const n = this.tangentAt(t);
+    return this.clockwise ? J(n) : bt(n);
+  }
+  gradientAt(t) {
+    const n = this.firstAngle + t * this.deltaAngle * (this.clockwise ? -1 : 1), e = -this.majorRadius * this.deltaAngle * Math.sin(n), s = this.minorRadius * this.deltaAngle * Math.cos(n), r = this.clockwise ? [-e, -s] : [e, s];
+    return this.rotateFromEllipseReferenceFrame.transform(r);
+  }
+  tangentAt(t) {
+    const n = this.pointTheta(t), e = -this.majorRadius * Math.sin(n), s = this.minorRadius * Math.cos(n), r = this.clockwise ? [-e, -s] : [e, s];
+    return W(this.rotateFromEllipseReferenceFrame.transform(r));
+  }
+  get tangentAtFirstPoint() {
+    const t = -this.majorRadius * Math.sin(this.firstAngle), n = this.minorRadius * Math.cos(this.firstAngle), e = this.clockwise ? [-t, -n] : [t, n];
+    return W(this.rotateFromEllipseReferenceFrame.transform(e));
+  }
+  get tangentAtLastPoint() {
+    const t = -this.majorRadius * Math.sin(this.lastAngle), n = this.minorRadius * Math.cos(this.lastAngle), e = this.clockwise ? [-t, -n] : [t, n];
+    return W(this.rotateFromEllipseReferenceFrame.transform(e));
+  }
+  transform(t) {
+    const n = t.transformAngle(this.tiltAngle), e = t.scaleFactor();
+    return new C(
+      t.transform(this.firstPoint),
+      t.transform(this.lastPoint),
+      t.transform(this.center),
+      this.majorRadius * e,
+      this.minorRadius * e,
+      n,
+      t.keepsOrientation() ? this.clockwise : !this.clockwise,
+      { angleUnits: "rad" }
+    );
+  }
+  splitAt(t) {
+    let n;
+    if (Array.isArray(t) && t.length === 0)
+      return [this];
+    Array.isArray(t[0]) ? n = t : n = [t];
+    const s = [0, 1, ...n.map((h) => this.pointToParam(h))], r = new Map(
+      nt([s, [this.firstPoint, this.lastPoint, ...n]])
+    );
+    s.sort((h, u) => h - u);
+    let o = null;
+    return s.flatMap((h, u) => {
+      if (u === s.length - 1)
+        return [];
+      const l = s[u + 1];
+      if (l - h < this.precision)
+        return o === null && (o = h), [];
+      const c = o === null ? h : o, a = new C(
+        r.get(c) || this.paramPoint(c),
+        r.get(l) || this.paramPoint(l),
+        this.center,
+        this.majorRadius,
+        this.minorRadius,
+        this.tiltAngle,
+        this.clockwise,
+        { angleUnits: "rad" }
+      );
+      return o = null, a;
+    });
+  }
+}
+function fs(i, t, n, e, s, r, o) {
+  const { center: h, rx: u, ry: l } = tn(
+    i,
+    t,
+    n,
+    e,
+    s * Lt,
+    !r,
+    o
+  );
+  return Math.abs(u - l) < 1e-9 ? new I(i, t, h, o) : new C(i, t, h, u, l, s, o);
+}
+function tn([i, t], [n, e], s, r, o, h, u) {
+  if (s < 0 && (s = -s), r < 0 && (r = -r), s == 0 || r == 0)
+    throw Error("rx and ry can not be 0");
+  const l = Math.sin(o), c = Math.cos(o), a = (i - n) / 2, g = (t - e) / 2, w = (i + n) / 2, m = (t + e) / 2, d = c * a + l * g, P = c * g - l * a, f = d * d / (s * s) + P * P / (r * r);
+  f > 1 && (s = s * Math.sqrt(f), r = r * Math.sqrt(f));
+  const p = s * r, y = s * P, M = r * d, b = y * y + M * M;
+  if (!b)
+    throw Error("start point can not be same as end point");
+  let x = Math.sqrt(Math.abs((p * p - b) / b));
+  h == u && (x = -x);
+  const k = x * y / r, S = -x * M / s, _ = c * k - l * S + w, F = l * k + c * S + m;
+  return {
+    center: [_, F],
+    rx: s,
+    ry: r
+  };
+}
+function en(i, t, n) {
+  const e = Math.abs(n[0]), s = Math.abs(n[1]);
+  let r = 0.707, o = 0.707;
+  const h = i, u = t;
+  for (let l = 0; l < 3; l++) {
+    const c = h * r, a = u * o, g = (h * h - u * u) * r ** 3 / h, w = (u * u - h * h) * o ** 3 / u, m = c - g, d = a - w, P = e - g, f = s - w, p = Math.hypot(m, d), y = Math.hypot(P, f);
+    r = Math.min(1, Math.max(0, (P * p / y + g) / h)), o = Math.min(1, Math.max(0, (f * p / y + w) / u));
+    const M = Math.hypot(r, o);
+    r /= M, o /= M;
+  }
+  return [h * r * Math.sign(n[0]), u * o * Math.sign(n[1])];
+}
+function qt(i, t, n = 1e-9) {
+  const e = i.transform(t.ellipseReferenceFrameTransform), s = e.slope, r = e.yIntercept, o = t.majorRadius * t.majorRadius, h = t.minorRadius * t.minorRadius, u = t.majorRadius * t.minorRadius, l = e.slope * e.slope, c = e.yIntercept * e.yIntercept, a = (f) => f.map(
+    (p) => t.reverseEllipseReferenceFrameTransform.transform(p)
+  ).filter((p) => i.isOnSegment(p) && t.isOnSegment(p));
+  if (!Number.isFinite(s)) {
+    const f = e.firstPoint[0];
+    if (Math.abs(f) - t.majorRadius > n)
+      return [];
+    if (Math.abs(Math.abs(f) - t.majorRadius) < n)
+      return a([[f, 0]]);
+    const p = t.minorRadius * Math.sqrt(1 - f * f / o), y = [f, p], M = [f, -p];
+    return a([y, M]);
+  }
+  const g = o * l + h - c;
+  if (g < -n)
+    return [];
+  const w = o * l + h;
+  if (Math.abs(g) < n) {
+    const f = -(o * s * r) / w, p = h * r / w;
+    return a([[f, p]]);
+  }
+  const m = Math.sqrt(g), d = [
+    -(o * s * r + u * m) / w,
+    (h * r - u * s * m) / w
+  ], P = [
+    -(o * s * r - u * m) / w,
+    (h * r + u * s * m) / w
+  ];
+  return a([d, P]);
+}
+const nn = Object.prototype.toString;
+function z(i) {
+  const t = nn.call(i);
+  return t.endsWith("Array]") && !t.includes("Big");
+}
+function sn(i) {
+  var t = arguments.length > 1 && arguments[1] !== void 0 ? arguments[1] : {};
+  if (!z(i))
+    throw new TypeError("input must be an array");
+  if (i.length === 0)
+    throw new TypeError("input must not be empty");
+  var n = t.fromIndex, e = n === void 0 ? 0 : n, s = t.toIndex, r = s === void 0 ? i.length : s;
+  if (e < 0 || e >= i.length || !Number.isInteger(e))
+    throw new Error("fromIndex must be a positive integer smaller than length");
+  if (r <= e || r > i.length || !Number.isInteger(r))
+    throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");
+  for (var o = i[e], h = e + 1; h < r; h++)
+    i[h] > o && (o = i[h]);
+  return o;
+}
+function rn(i) {
+  var t = arguments.length > 1 && arguments[1] !== void 0 ? arguments[1] : {};
+  if (!z(i))
+    throw new TypeError("input must be an array");
+  if (i.length === 0)
+    throw new TypeError("input must not be empty");
+  var n = t.fromIndex, e = n === void 0 ? 0 : n, s = t.toIndex, r = s === void 0 ? i.length : s;
+  if (e < 0 || e >= i.length || !Number.isInteger(e))
+    throw new Error("fromIndex must be a positive integer smaller than length");
+  if (r <= e || r > i.length || !Number.isInteger(r))
+    throw new Error("toIndex must be an integer greater than fromIndex and at most equal to length");
+  for (var o = i[e], h = e + 1; h < r; h++)
+    i[h] < o && (o = i[h]);
+  return o;
+}
+function Kt(i) {
+  var t = arguments.length > 1 && arguments[1] !== void 0 ? arguments[1] : {};
+  if (z(i)) {
+    if (i.length === 0)
+      throw new TypeError("input must not be empty");
+  } else
+    throw new TypeError("input must be an array");
+  var n;
+  if (t.output !== void 0) {
+    if (!z(t.output))
+      throw new TypeError("output option must be an array if specified");
+    n = t.output;
+  } else
+    n = new Array(i.length);
+  var e = rn(i), s = sn(i);
+  if (e === s)
+    throw new RangeError("minimum and maximum input values are equal. Cannot rescale a constant array");
+  var r = t.min, o = r === void 0 ? t.autoMinMax ? e : 0 : r, h = t.max, u = h === void 0 ? t.autoMinMax ? s : 1 : h;
+  if (o >= u)
+    throw new RangeError("min option must be smaller than max option");
+  for (var l = (u - o) / (s - e), c = 0; c < i.length; c++)
+    n[c] = (i[c] - e) * l + o;
+  return n;
+}
+const Mt = " ".repeat(2), pe = " ".repeat(4);
+function on() {
+  return we(this);
+}
+function we(i, t = {}) {
+  const {
+    maxRows: n = 15,
+    maxColumns: e = 10,
+    maxNumSize: s = 8,
+    padMinus: r = "auto"
+  } = t;
+  return `${i.constructor.name} {
+${Mt}[
+${pe}${hn(i, n, e, s, r)}
+${Mt}]
+${Mt}rows: ${i.rows}
+${Mt}columns: ${i.columns}
+}`;
+}
+function hn(i, t, n, e, s) {
+  const { rows: r, columns: o } = i, h = Math.min(r, t), u = Math.min(o, n), l = [];
+  if (s === "auto") {
+    s = !1;
+    t:
+      for (let c = 0; c < h; c++)
+        for (let a = 0; a < u; a++)
+          if (i.get(c, a) < 0) {
+            s = !0;
+            break t;
+          }
+  }
+  for (let c = 0; c < h; c++) {
+    let a = [];
+    for (let g = 0; g < u; g++)
+      a.push(ln(i.get(c, g), e, s));
+    l.push(`${a.join(" ")}`);
+  }
+  return u !== o && (l[l.length - 1] += ` ... ${o - n} more columns`), h !== r && l.push(`... ${r - t} more rows`), l.join(`
+${pe}`);
+}
+function ln(i, t, n) {
+  return (i >= 0 && n ? ` ${Qt(i, t - 1)}` : Qt(i, t)).padEnd(t);
+}
+function Qt(i, t) {
+  let n = i.toString();
+  if (n.length <= t)
+    return n;
+  let e = i.toFixed(t);
+  if (e.length > t && (e = i.toFixed(Math.max(0, t - (e.length - t)))), e.length <= t && !e.startsWith("0.000") && !e.startsWith("-0.000"))
+    return e;
+  let s = i.toExponential(t);
+  return s.length > t && (s = i.toExponential(Math.max(0, t - (s.length - t)))), s.slice(0);
+}
+function un(i, t) {
+  i.prototype.add = function(e) {
+    return typeof e == "number" ? this.addS(e) : this.addM(e);
+  }, i.prototype.addS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) + e);
+    return this;
+  }, i.prototype.addM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) + e.get(s, r));
+    return this;
+  }, i.add = function(e, s) {
+    return new t(e).add(s);
+  }, i.prototype.sub = function(e) {
+    return typeof e == "number" ? this.subS(e) : this.subM(e);
+  }, i.prototype.subS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) - e);
+    return this;
+  }, i.prototype.subM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) - e.get(s, r));
+    return this;
+  }, i.sub = function(e, s) {
+    return new t(e).sub(s);
+  }, i.prototype.subtract = i.prototype.sub, i.prototype.subtractS = i.prototype.subS, i.prototype.subtractM = i.prototype.subM, i.subtract = i.sub, i.prototype.mul = function(e) {
+    return typeof e == "number" ? this.mulS(e) : this.mulM(e);
+  }, i.prototype.mulS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) * e);
+    return this;
+  }, i.prototype.mulM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) * e.get(s, r));
+    return this;
+  }, i.mul = function(e, s) {
+    return new t(e).mul(s);
+  }, i.prototype.multiply = i.prototype.mul, i.prototype.multiplyS = i.prototype.mulS, i.prototype.multiplyM = i.prototype.mulM, i.multiply = i.mul, i.prototype.div = function(e) {
+    return typeof e == "number" ? this.divS(e) : this.divM(e);
+  }, i.prototype.divS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) / e);
+    return this;
+  }, i.prototype.divM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) / e.get(s, r));
+    return this;
+  }, i.div = function(e, s) {
+    return new t(e).div(s);
+  }, i.prototype.divide = i.prototype.div, i.prototype.divideS = i.prototype.divS, i.prototype.divideM = i.prototype.divM, i.divide = i.div, i.prototype.mod = function(e) {
+    return typeof e == "number" ? this.modS(e) : this.modM(e);
+  }, i.prototype.modS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) % e);
+    return this;
+  }, i.prototype.modM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) % e.get(s, r));
+    return this;
+  }, i.mod = function(e, s) {
+    return new t(e).mod(s);
+  }, i.prototype.modulus = i.prototype.mod, i.prototype.modulusS = i.prototype.modS, i.prototype.modulusM = i.prototype.modM, i.modulus = i.mod, i.prototype.and = function(e) {
+    return typeof e == "number" ? this.andS(e) : this.andM(e);
+  }, i.prototype.andS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) & e);
+    return this;
+  }, i.prototype.andM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) & e.get(s, r));
+    return this;
+  }, i.and = function(e, s) {
+    return new t(e).and(s);
+  }, i.prototype.or = function(e) {
+    return typeof e == "number" ? this.orS(e) : this.orM(e);
+  }, i.prototype.orS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) | e);
+    return this;
+  }, i.prototype.orM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) | e.get(s, r));
+    return this;
+  }, i.or = function(e, s) {
+    return new t(e).or(s);
+  }, i.prototype.xor = function(e) {
+    return typeof e == "number" ? this.xorS(e) : this.xorM(e);
+  }, i.prototype.xorS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) ^ e);
+    return this;
+  }, i.prototype.xorM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) ^ e.get(s, r));
+    return this;
+  }, i.xor = function(e, s) {
+    return new t(e).xor(s);
+  }, i.prototype.leftShift = function(e) {
+    return typeof e == "number" ? this.leftShiftS(e) : this.leftShiftM(e);
+  }, i.prototype.leftShiftS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) << e);
+    return this;
+  }, i.prototype.leftShiftM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) << e.get(s, r));
+    return this;
+  }, i.leftShift = function(e, s) {
+    return new t(e).leftShift(s);
+  }, i.prototype.signPropagatingRightShift = function(e) {
+    return typeof e == "number" ? this.signPropagatingRightShiftS(e) : this.signPropagatingRightShiftM(e);
+  }, i.prototype.signPropagatingRightShiftS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) >> e);
+    return this;
+  }, i.prototype.signPropagatingRightShiftM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) >> e.get(s, r));
+    return this;
+  }, i.signPropagatingRightShift = function(e, s) {
+    return new t(e).signPropagatingRightShift(s);
+  }, i.prototype.rightShift = function(e) {
+    return typeof e == "number" ? this.rightShiftS(e) : this.rightShiftM(e);
+  }, i.prototype.rightShiftS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) >>> e);
+    return this;
+  }, i.prototype.rightShiftM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, this.get(s, r) >>> e.get(s, r));
+    return this;
+  }, i.rightShift = function(e, s) {
+    return new t(e).rightShift(s);
+  }, i.prototype.zeroFillRightShift = i.prototype.rightShift, i.prototype.zeroFillRightShiftS = i.prototype.rightShiftS, i.prototype.zeroFillRightShiftM = i.prototype.rightShiftM, i.zeroFillRightShift = i.rightShift, i.prototype.not = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, ~this.get(e, s));
+    return this;
+  }, i.not = function(e) {
+    return new t(e).not();
+  }, i.prototype.abs = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.abs(this.get(e, s)));
+    return this;
+  }, i.abs = function(e) {
+    return new t(e).abs();
+  }, i.prototype.acos = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.acos(this.get(e, s)));
+    return this;
+  }, i.acos = function(e) {
+    return new t(e).acos();
+  }, i.prototype.acosh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.acosh(this.get(e, s)));
+    return this;
+  }, i.acosh = function(e) {
+    return new t(e).acosh();
+  }, i.prototype.asin = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.asin(this.get(e, s)));
+    return this;
+  }, i.asin = function(e) {
+    return new t(e).asin();
+  }, i.prototype.asinh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.asinh(this.get(e, s)));
+    return this;
+  }, i.asinh = function(e) {
+    return new t(e).asinh();
+  }, i.prototype.atan = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.atan(this.get(e, s)));
+    return this;
+  }, i.atan = function(e) {
+    return new t(e).atan();
+  }, i.prototype.atanh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.atanh(this.get(e, s)));
+    return this;
+  }, i.atanh = function(e) {
+    return new t(e).atanh();
+  }, i.prototype.cbrt = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.cbrt(this.get(e, s)));
+    return this;
+  }, i.cbrt = function(e) {
+    return new t(e).cbrt();
+  }, i.prototype.ceil = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.ceil(this.get(e, s)));
+    return this;
+  }, i.ceil = function(e) {
+    return new t(e).ceil();
+  }, i.prototype.clz32 = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.clz32(this.get(e, s)));
+    return this;
+  }, i.clz32 = function(e) {
+    return new t(e).clz32();
+  }, i.prototype.cos = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.cos(this.get(e, s)));
+    return this;
+  }, i.cos = function(e) {
+    return new t(e).cos();
+  }, i.prototype.cosh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.cosh(this.get(e, s)));
+    return this;
+  }, i.cosh = function(e) {
+    return new t(e).cosh();
+  }, i.prototype.exp = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.exp(this.get(e, s)));
+    return this;
+  }, i.exp = function(e) {
+    return new t(e).exp();
+  }, i.prototype.expm1 = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.expm1(this.get(e, s)));
+    return this;
+  }, i.expm1 = function(e) {
+    return new t(e).expm1();
+  }, i.prototype.floor = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.floor(this.get(e, s)));
+    return this;
+  }, i.floor = function(e) {
+    return new t(e).floor();
+  }, i.prototype.fround = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.fround(this.get(e, s)));
+    return this;
+  }, i.fround = function(e) {
+    return new t(e).fround();
+  }, i.prototype.log = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.log(this.get(e, s)));
+    return this;
+  }, i.log = function(e) {
+    return new t(e).log();
+  }, i.prototype.log1p = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.log1p(this.get(e, s)));
+    return this;
+  }, i.log1p = function(e) {
+    return new t(e).log1p();
+  }, i.prototype.log10 = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.log10(this.get(e, s)));
+    return this;
+  }, i.log10 = function(e) {
+    return new t(e).log10();
+  }, i.prototype.log2 = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.log2(this.get(e, s)));
+    return this;
+  }, i.log2 = function(e) {
+    return new t(e).log2();
+  }, i.prototype.round = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.round(this.get(e, s)));
+    return this;
+  }, i.round = function(e) {
+    return new t(e).round();
+  }, i.prototype.sign = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.sign(this.get(e, s)));
+    return this;
+  }, i.sign = function(e) {
+    return new t(e).sign();
+  }, i.prototype.sin = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.sin(this.get(e, s)));
+    return this;
+  }, i.sin = function(e) {
+    return new t(e).sin();
+  }, i.prototype.sinh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.sinh(this.get(e, s)));
+    return this;
+  }, i.sinh = function(e) {
+    return new t(e).sinh();
+  }, i.prototype.sqrt = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.sqrt(this.get(e, s)));
+    return this;
+  }, i.sqrt = function(e) {
+    return new t(e).sqrt();
+  }, i.prototype.tan = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.tan(this.get(e, s)));
+    return this;
+  }, i.tan = function(e) {
+    return new t(e).tan();
+  }, i.prototype.tanh = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.tanh(this.get(e, s)));
+    return this;
+  }, i.tanh = function(e) {
+    return new t(e).tanh();
+  }, i.prototype.trunc = function() {
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.set(e, s, Math.trunc(this.get(e, s)));
+    return this;
+  }, i.trunc = function(e) {
+    return new t(e).trunc();
+  }, i.pow = function(e, s) {
+    return new t(e).pow(s);
+  }, i.prototype.pow = function(e) {
+    return typeof e == "number" ? this.powS(e) : this.powM(e);
+  }, i.prototype.powS = function(e) {
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, Math.pow(this.get(s, r), e));
+    return this;
+  }, i.prototype.powM = function(e) {
+    if (e = t.checkMatrix(e), this.rows !== e.rows || this.columns !== e.columns)
+      throw new RangeError("Matrices dimensions must be equal");
+    for (let s = 0; s < this.rows; s++)
+      for (let r = 0; r < this.columns; r++)
+        this.set(s, r, Math.pow(this.get(s, r), e.get(s, r)));
+    return this;
+  };
+}
+function Y(i, t, n) {
+  let e = n ? i.rows : i.rows - 1;
+  if (t < 0 || t > e)
+    throw new RangeError("Row index out of range");
+}
+function X(i, t, n) {
+  let e = n ? i.columns : i.columns - 1;
+  if (t < 0 || t > e)
+    throw new RangeError("Column index out of range");
+}
+function lt(i, t) {
+  if (t.to1DArray && (t = t.to1DArray()), t.length !== i.columns)
+    throw new RangeError(
+      "vector size must be the same as the number of columns"
+    );
+  return t;
+}
+function ut(i, t) {
+  if (t.to1DArray && (t = t.to1DArray()), t.length !== i.rows)
+    throw new RangeError("vector size must be the same as the number of rows");
+  return t;
+}
+function an(i, t) {
+  if (!z(t))
+    throw new TypeError("row indices must be an array");
+  for (let n = 0; n < t.length; n++)
+    if (t[n] < 0 || t[n] >= i.rows)
+      throw new RangeError("row indices are out of range");
+}
+function cn(i, t) {
+  if (!z(t))
+    throw new TypeError("column indices must be an array");
+  for (let n = 0; n < t.length; n++)
+    if (t[n] < 0 || t[n] >= i.columns)
+      throw new RangeError("column indices are out of range");
+}
+function Zt(i, t, n, e, s) {
+  if (arguments.length !== 5)
+    throw new RangeError("expected 4 arguments");
+  if (Et("startRow", t), Et("endRow", n), Et("startColumn", e), Et("endColumn", s), t > n || e > s || t < 0 || t >= i.rows || n < 0 || n >= i.rows || e < 0 || e >= i.columns || s < 0 || s >= i.columns)
+    throw new RangeError("Submatrix indices are out of range");
+}
+function Tt(i, t = 0) {
+  let n = [];
+  for (let e = 0; e < i; e++)
+    n.push(t);
+  return n;
+}
+function Et(i, t) {
+  if (typeof t != "number")
+    throw new TypeError(`${i} must be a number`);
+}
+function ht(i) {
+  if (i.isEmpty())
+    throw new Error("Empty matrix has no elements to index");
+}
+function fn(i) {
+  let t = Tt(i.rows);
+  for (let n = 0; n < i.rows; ++n)
+    for (let e = 0; e < i.columns; ++e)
+      t[n] += i.get(n, e);
+  return t;
+}
+function gn(i) {
+  let t = Tt(i.columns);
+  for (let n = 0; n < i.rows; ++n)
+    for (let e = 0; e < i.columns; ++e)
+      t[e] += i.get(n, e);
+  return t;
+}
+function mn(i) {
+  let t = 0;
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      t += i.get(n, e);
+  return t;
+}
+function pn(i) {
+  let t = Tt(i.rows, 1);
+  for (let n = 0; n < i.rows; ++n)
+    for (let e = 0; e < i.columns; ++e)
+      t[n] *= i.get(n, e);
+  return t;
+}
+function wn(i) {
+  let t = Tt(i.columns, 1);
+  for (let n = 0; n < i.rows; ++n)
+    for (let e = 0; e < i.columns; ++e)
+      t[e] *= i.get(n, e);
+  return t;
+}
+function dn(i) {
+  let t = 1;
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      t *= i.get(n, e);
+  return t;
+}
+function yn(i, t, n) {
+  const e = i.rows, s = i.columns, r = [];
+  for (let o = 0; o < e; o++) {
+    let h = 0, u = 0, l = 0;
+    for (let c = 0; c < s; c++)
+      l = i.get(o, c) - n[o], h += l, u += l * l;
+    t ? r.push((u - h * h / s) / (s - 1)) : r.push((u - h * h / s) / s);
+  }
+  return r;
+}
+function Pn(i, t, n) {
+  const e = i.rows, s = i.columns, r = [];
+  for (let o = 0; o < s; o++) {
+    let h = 0, u = 0, l = 0;
+    for (let c = 0; c < e; c++)
+      l = i.get(c, o) - n[o], h += l, u += l * l;
+    t ? r.push((u - h * h / e) / (e - 1)) : r.push((u - h * h / e) / e);
+  }
+  return r;
+}
+function Mn(i, t, n) {
+  const e = i.rows, s = i.columns, r = e * s;
+  let o = 0, h = 0, u = 0;
+  for (let l = 0; l < e; l++)
+    for (let c = 0; c < s; c++)
+      u = i.get(l, c) - n, o += u, h += u * u;
+  return t ? (h - o * o / r) / (r - 1) : (h - o * o / r) / r;
+}
+function En(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) - t[n]);
+}
+function Sn(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) - t[e]);
+}
+function bn(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) - t);
+}
+function jn(i) {
+  const t = [];
+  for (let n = 0; n < i.rows; n++) {
+    let e = 0;
+    for (let s = 0; s < i.columns; s++)
+      e += Math.pow(i.get(n, s), 2) / (i.columns - 1);
+    t.push(Math.sqrt(e));
+  }
+  return t;
+}
+function vn(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) / t[n]);
+}
+function kn(i) {
+  const t = [];
+  for (let n = 0; n < i.columns; n++) {
+    let e = 0;
+    for (let s = 0; s < i.rows; s++)
+      e += Math.pow(i.get(s, n), 2) / (i.rows - 1);
+    t.push(Math.sqrt(e));
+  }
+  return t;
+}
+function In(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) / t[e]);
+}
+function xn(i) {
+  const t = i.size - 1;
+  let n = 0;
+  for (let e = 0; e < i.columns; e++)
+    for (let s = 0; s < i.rows; s++)
+      n += Math.pow(i.get(s, e), 2) / t;
+  return Math.sqrt(n);
+}
+function Rn(i, t) {
+  for (let n = 0; n < i.rows; n++)
+    for (let e = 0; e < i.columns; e++)
+      i.set(n, e, i.get(n, e) / t);
+}
+class v {
+  static from1DArray(t, n, e) {
+    if (t * n !== e.length)
+      throw new RangeError("data length does not match given dimensions");
+    let r = new j(t, n);
+    for (let o = 0; o < t; o++)
+      for (let h = 0; h < n; h++)
+        r.set(o, h, e[o * n + h]);
+    return r;
+  }
+  static rowVector(t) {
+    let n = new j(1, t.length);
+    for (let e = 0; e < t.length; e++)
+      n.set(0, e, t[e]);
+    return n;
+  }
+  static columnVector(t) {
+    let n = new j(t.length, 1);
+    for (let e = 0; e < t.length; e++)
+      n.set(e, 0, t[e]);
+    return n;
+  }
+  static zeros(t, n) {
+    return new j(t, n);
+  }
+  static ones(t, n) {
+    return new j(t, n).fill(1);
+  }
+  static rand(t, n, e = {}) {
+    if (typeof e != "object")
+      throw new TypeError("options must be an object");
+    const { random: s = Math.random } = e;
+    let r = new j(t, n);
+    for (let o = 0; o < t; o++)
+      for (let h = 0; h < n; h++)
+        r.set(o, h, s());
+    return r;
+  }
+  static randInt(t, n, e = {}) {
+    if (typeof e != "object")
+      throw new TypeError("options must be an object");
+    const { min: s = 0, max: r = 1e3, random: o = Math.random } = e;
+    if (!Number.isInteger(s))
+      throw new TypeError("min must be an integer");
+    if (!Number.isInteger(r))
+      throw new TypeError("max must be an integer");
+    if (s >= r)
+      throw new RangeError("min must be smaller than max");
+    let h = r - s, u = new j(t, n);
+    for (let l = 0; l < t; l++)
+      for (let c = 0; c < n; c++) {
+        let a = s + Math.round(o() * h);
+        u.set(l, c, a);
+      }
+    return u;
+  }
+  static eye(t, n, e) {
+    n === void 0 && (n = t), e === void 0 && (e = 1);
+    let s = Math.min(t, n), r = this.zeros(t, n);
+    for (let o = 0; o < s; o++)
+      r.set(o, o, e);
+    return r;
+  }
+  static diag(t, n, e) {
+    let s = t.length;
+    n === void 0 && (n = s), e === void 0 && (e = n);
+    let r = Math.min(s, n, e), o = this.zeros(n, e);
+    for (let h = 0; h < r; h++)
+      o.set(h, h, t[h]);
+    return o;
+  }
+  static min(t, n) {
+    t = this.checkMatrix(t), n = this.checkMatrix(n);
+    let e = t.rows, s = t.columns, r = new j(e, s);
+    for (let o = 0; o < e; o++)
+      for (let h = 0; h < s; h++)
+        r.set(o, h, Math.min(t.get(o, h), n.get(o, h)));
+    return r;
+  }
+  static max(t, n) {
+    t = this.checkMatrix(t), n = this.checkMatrix(n);
+    let e = t.rows, s = t.columns, r = new this(e, s);
+    for (let o = 0; o < e; o++)
+      for (let h = 0; h < s; h++)
+        r.set(o, h, Math.max(t.get(o, h), n.get(o, h)));
+    return r;
+  }
+  static checkMatrix(t) {
+    return v.isMatrix(t) ? t : new j(t);
+  }
+  static isMatrix(t) {
+    return t != null && t.klass === "Matrix";
+  }
+  get size() {
+    return this.rows * this.columns;
+  }
+  apply(t) {
+    if (typeof t != "function")
+      throw new TypeError("callback must be a function");
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        t.call(this, n, e);
+    return this;
+  }
+  to1DArray() {
+    let t = [];
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        t.push(this.get(n, e));
+    return t;
+  }
+  to2DArray() {
+    let t = [];
+    for (let n = 0; n < this.rows; n++) {
+      t.push([]);
+      for (let e = 0; e < this.columns; e++)
+        t[n].push(this.get(n, e));
+    }
+    return t;
+  }
+  toJSON() {
+    return this.to2DArray();
+  }
+  isRowVector() {
+    return this.rows === 1;
+  }
+  isColumnVector() {
+    return this.columns === 1;
+  }
+  isVector() {
+    return this.rows === 1 || this.columns === 1;
+  }
+  isSquare() {
+    return this.rows === this.columns;
+  }
+  isEmpty() {
+    return this.rows === 0 || this.columns === 0;
+  }
+  isSymmetric() {
+    if (this.isSquare()) {
+      for (let t = 0; t < this.rows; t++)
+        for (let n = 0; n <= t; n++)
+          if (this.get(t, n) !== this.get(n, t))
+            return !1;
+      return !0;
+    }
+    return !1;
+  }
+  isEchelonForm() {
+    let t = 0, n = 0, e = -1, s = !0, r = !1;
+    for (; t < this.rows && s; ) {
+      for (n = 0, r = !1; n < this.columns && r === !1; )
+        this.get(t, n) === 0 ? n++ : this.get(t, n) === 1 && n > e ? (r = !0, e = n) : (s = !1, r = !0);
+      t++;
+    }
+    return s;
+  }
+  isReducedEchelonForm() {
+    let t = 0, n = 0, e = -1, s = !0, r = !1;
+    for (; t < this.rows && s; ) {
+      for (n = 0, r = !1; n < this.columns && r === !1; )
+        this.get(t, n) === 0 ? n++ : this.get(t, n) === 1 && n > e ? (r = !0, e = n) : (s = !1, r = !0);
+      for (let o = n + 1; o < this.rows; o++)
+        this.get(t, o) !== 0 && (s = !1);
+      t++;
+    }
+    return s;
+  }
+  echelonForm() {
+    let t = this.clone(), n = 0, e = 0;
+    for (; n < t.rows && e < t.columns; ) {
+      let s = n;
+      for (let r = n; r < t.rows; r++)
+        t.get(r, e) > t.get(s, e) && (s = r);
+      if (t.get(s, e) === 0)
+        e++;
+      else {
+        t.swapRows(n, s);
+        let r = t.get(n, e);
+        for (let o = e; o < t.columns; o++)
+          t.set(n, o, t.get(n, o) / r);
+        for (let o = n + 1; o < t.rows; o++) {
+          let h = t.get(o, e) / t.get(n, e);
+          t.set(o, e, 0);
+          for (let u = e + 1; u < t.columns; u++)
+            t.set(o, u, t.get(o, u) - t.get(n, u) * h);
+        }
+        n++, e++;
+      }
+    }
+    return t;
+  }
+  reducedEchelonForm() {
+    let t = this.echelonForm(), n = t.columns, e = t.rows, s = e - 1;
+    for (; s >= 0; )
+      if (t.maxRow(s) === 0)
+        s--;
+      else {
+        let r = 0, o = !1;
+        for (; r < e && o === !1; )
+          t.get(s, r) === 1 ? o = !0 : r++;
+        for (let h = 0; h < s; h++) {
+          let u = t.get(h, r);
+          for (let l = r; l < n; l++) {
+            let c = t.get(h, l) - u * t.get(s, l);
+            t.set(h, l, c);
+          }
+        }
+        s--;
+      }
+    return t;
+  }
+  set() {
+    throw new Error("set method is unimplemented");
+  }
+  get() {
+    throw new Error("get method is unimplemented");
+  }
+  repeat(t = {}) {
+    if (typeof t != "object")
+      throw new TypeError("options must be an object");
+    const { rows: n = 1, columns: e = 1 } = t;
+    if (!Number.isInteger(n) || n <= 0)
+      throw new TypeError("rows must be a positive integer");
+    if (!Number.isInteger(e) || e <= 0)
+      throw new TypeError("columns must be a positive integer");
+    let s = new j(this.rows * n, this.columns * e);
+    for (let r = 0; r < n; r++)
+      for (let o = 0; o < e; o++)
+        s.setSubMatrix(this, this.rows * r, this.columns * o);
+    return s;
+  }
+  fill(t) {
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, t);
+    return this;
+  }
+  neg() {
+    return this.mulS(-1);
+  }
+  getRow(t) {
+    Y(this, t);
+    let n = [];
+    for (let e = 0; e < this.columns; e++)
+      n.push(this.get(t, e));
+    return n;
+  }
+  getRowVector(t) {
+    return j.rowVector(this.getRow(t));
+  }
+  setRow(t, n) {
+    Y(this, t), n = lt(this, n);
+    for (let e = 0; e < this.columns; e++)
+      this.set(t, e, n[e]);
+    return this;
+  }
+  swapRows(t, n) {
+    Y(this, t), Y(this, n);
+    for (let e = 0; e < this.columns; e++) {
+      let s = this.get(t, e);
+      this.set(t, e, this.get(n, e)), this.set(n, e, s);
+    }
+    return this;
+  }
+  getColumn(t) {
+    X(this, t);
+    let n = [];
+    for (let e = 0; e < this.rows; e++)
+      n.push(this.get(e, t));
+    return n;
+  }
+  getColumnVector(t) {
+    return j.columnVector(this.getColumn(t));
+  }
+  setColumn(t, n) {
+    X(this, t), n = ut(this, n);
+    for (let e = 0; e < this.rows; e++)
+      this.set(e, t, n[e]);
+    return this;
+  }
+  swapColumns(t, n) {
+    X(this, t), X(this, n);
+    for (let e = 0; e < this.rows; e++) {
+      let s = this.get(e, t);
+      this.set(e, t, this.get(e, n)), this.set(e, n, s);
+    }
+    return this;
+  }
+  addRowVector(t) {
+    t = lt(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) + t[e]);
+    return this;
+  }
+  subRowVector(t) {
+    t = lt(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) - t[e]);
+    return this;
+  }
+  mulRowVector(t) {
+    t = lt(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) * t[e]);
+    return this;
+  }
+  divRowVector(t) {
+    t = lt(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) / t[e]);
+    return this;
+  }
+  addColumnVector(t) {
+    t = ut(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) + t[n]);
+    return this;
+  }
+  subColumnVector(t) {
+    t = ut(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) - t[n]);
+    return this;
+  }
+  mulColumnVector(t) {
+    t = ut(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) * t[n]);
+    return this;
+  }
+  divColumnVector(t) {
+    t = ut(this, t);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        this.set(n, e, this.get(n, e) / t[n]);
+    return this;
+  }
+  mulRow(t, n) {
+    Y(this, t);
+    for (let e = 0; e < this.columns; e++)
+      this.set(t, e, this.get(t, e) * n);
+    return this;
+  }
+  mulColumn(t, n) {
+    X(this, t);
+    for (let e = 0; e < this.rows; e++)
+      this.set(e, t, this.get(e, t) * n);
+    return this;
+  }
+  max(t) {
+    if (this.isEmpty())
+      return NaN;
+    switch (t) {
+      case "row": {
+        const n = new Array(this.rows).fill(Number.NEGATIVE_INFINITY);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) > n[e] && (n[e] = this.get(e, s));
+        return n;
+      }
+      case "column": {
+        const n = new Array(this.columns).fill(Number.NEGATIVE_INFINITY);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) > n[s] && (n[s] = this.get(e, s));
+        return n;
+      }
+      case void 0: {
+        let n = this.get(0, 0);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) > n && (n = this.get(e, s));
+        return n;
+      }
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  maxIndex() {
+    ht(this);
+    let t = this.get(0, 0), n = [0, 0];
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.get(e, s) > t && (t = this.get(e, s), n[0] = e, n[1] = s);
+    return n;
+  }
+  min(t) {
+    if (this.isEmpty())
+      return NaN;
+    switch (t) {
+      case "row": {
+        const n = new Array(this.rows).fill(Number.POSITIVE_INFINITY);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) < n[e] && (n[e] = this.get(e, s));
+        return n;
+      }
+      case "column": {
+        const n = new Array(this.columns).fill(Number.POSITIVE_INFINITY);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) < n[s] && (n[s] = this.get(e, s));
+        return n;
+      }
+      case void 0: {
+        let n = this.get(0, 0);
+        for (let e = 0; e < this.rows; e++)
+          for (let s = 0; s < this.columns; s++)
+            this.get(e, s) < n && (n = this.get(e, s));
+        return n;
+      }
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  minIndex() {
+    ht(this);
+    let t = this.get(0, 0), n = [0, 0];
+    for (let e = 0; e < this.rows; e++)
+      for (let s = 0; s < this.columns; s++)
+        this.get(e, s) < t && (t = this.get(e, s), n[0] = e, n[1] = s);
+    return n;
+  }
+  maxRow(t) {
+    if (Y(this, t), this.isEmpty())
+      return NaN;
+    let n = this.get(t, 0);
+    for (let e = 1; e < this.columns; e++)
+      this.get(t, e) > n && (n = this.get(t, e));
+    return n;
+  }
+  maxRowIndex(t) {
+    Y(this, t), ht(this);
+    let n = this.get(t, 0), e = [t, 0];
+    for (let s = 1; s < this.columns; s++)
+      this.get(t, s) > n && (n = this.get(t, s), e[1] = s);
+    return e;
+  }
+  minRow(t) {
+    if (Y(this, t), this.isEmpty())
+      return NaN;
+    let n = this.get(t, 0);
+    for (let e = 1; e < this.columns; e++)
+      this.get(t, e) < n && (n = this.get(t, e));
+    return n;
+  }
+  minRowIndex(t) {
+    Y(this, t), ht(this);
+    let n = this.get(t, 0), e = [t, 0];
+    for (let s = 1; s < this.columns; s++)
+      this.get(t, s) < n && (n = this.get(t, s), e[1] = s);
+    return e;
+  }
+  maxColumn(t) {
+    if (X(this, t), this.isEmpty())
+      return NaN;
+    let n = this.get(0, t);
+    for (let e = 1; e < this.rows; e++)
+      this.get(e, t) > n && (n = this.get(e, t));
+    return n;
+  }
+  maxColumnIndex(t) {
+    X(this, t), ht(this);
+    let n = this.get(0, t), e = [0, t];
+    for (let s = 1; s < this.rows; s++)
+      this.get(s, t) > n && (n = this.get(s, t), e[0] = s);
+    return e;
+  }
+  minColumn(t) {
+    if (X(this, t), this.isEmpty())
+      return NaN;
+    let n = this.get(0, t);
+    for (let e = 1; e < this.rows; e++)
+      this.get(e, t) < n && (n = this.get(e, t));
+    return n;
+  }
+  minColumnIndex(t) {
+    X(this, t), ht(this);
+    let n = this.get(0, t), e = [0, t];
+    for (let s = 1; s < this.rows; s++)
+      this.get(s, t) < n && (n = this.get(s, t), e[0] = s);
+    return e;
+  }
+  diag() {
+    let t = Math.min(this.rows, this.columns), n = [];
+    for (let e = 0; e < t; e++)
+      n.push(this.get(e, e));
+    return n;
+  }
+  norm(t = "frobenius") {
+    let n = 0;
+    if (t === "max")
+      return this.max();
+    if (t === "frobenius") {
+      for (let e = 0; e < this.rows; e++)
+        for (let s = 0; s < this.columns; s++)
+          n = n + this.get(e, s) * this.get(e, s);
+      return Math.sqrt(n);
+    } else
+      throw new RangeError(`unknown norm type: ${t}`);
+  }
+  cumulativeSum() {
+    let t = 0;
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        t += this.get(n, e), this.set(n, e, t);
+    return this;
+  }
+  dot(t) {
+    v.isMatrix(t) && (t = t.to1DArray());
+    let n = this.to1DArray();
+    if (n.length !== t.length)
+      throw new RangeError("vectors do not have the same size");
+    let e = 0;
+    for (let s = 0; s < n.length; s++)
+      e += n[s] * t[s];
+    return e;
+  }
+  mmul(t) {
+    t = j.checkMatrix(t);
+    let n = this.rows, e = this.columns, s = t.columns, r = new j(n, s), o = new Float64Array(e);
+    for (let h = 0; h < s; h++) {
+      for (let u = 0; u < e; u++)
+        o[u] = t.get(u, h);
+      for (let u = 0; u < n; u++) {
+        let l = 0;
+        for (let c = 0; c < e; c++)
+          l += this.get(u, c) * o[c];
+        r.set(u, h, l);
+      }
+    }
+    return r;
+  }
+  strassen2x2(t) {
+    t = j.checkMatrix(t);
+    let n = new j(2, 2);
+    const e = this.get(0, 0), s = t.get(0, 0), r = this.get(0, 1), o = t.get(0, 1), h = this.get(1, 0), u = t.get(1, 0), l = this.get(1, 1), c = t.get(1, 1), a = (e + l) * (s + c), g = (h + l) * s, w = e * (o - c), m = l * (u - s), d = (e + r) * c, P = (h - e) * (s + o), f = (r - l) * (u + c), p = a + m - d + f, y = w + d, M = g + m, b = a - g + w + P;
+    return n.set(0, 0, p), n.set(0, 1, y), n.set(1, 0, M), n.set(1, 1, b), n;
+  }
+  strassen3x3(t) {
+    t = j.checkMatrix(t);
+    let n = new j(3, 3);
+    const e = this.get(0, 0), s = this.get(0, 1), r = this.get(0, 2), o = this.get(1, 0), h = this.get(1, 1), u = this.get(1, 2), l = this.get(2, 0), c = this.get(2, 1), a = this.get(2, 2), g = t.get(0, 0), w = t.get(0, 1), m = t.get(0, 2), d = t.get(1, 0), P = t.get(1, 1), f = t.get(1, 2), p = t.get(2, 0), y = t.get(2, 1), M = t.get(2, 2), b = (e + s + r - o - h - c - a) * P, x = (e - o) * (-w + P), k = h * (-g + w + d - P - f - p + M), S = (-e + o + h) * (g - w + P), _ = (o + h) * (-g + w), F = e * g, B = (-e + l + c) * (g - m + f), K = (-e + l) * (m - f), Q = (l + c) * (-g + m), Z = (e + s + r - h - u - l - c) * f, $ = c * (-g + m + d - P - f - p + y), N = (-r + c + a) * (P + p - y), dt = (r - a) * (P - y), H = r * p, rt = (c + a) * (-p + y), L = (-r + h + u) * (f + p - M), Gt = (r - u) * (f - M), Wt = (h + u) * (-p + M), Re = s * d, _e = u * y, Fe = o * m, Te = l * w, Ce = a * M, $e = F + H + Re, qe = b + S + _ + F + N + H + rt, Be = F + B + Q + Z + H + L + Wt, Oe = x + k + S + F + H + L + Gt, Ae = x + S + _ + F + _e, Ne = H + L + Gt + Wt + Fe, Le = F + B + K + $ + N + dt + H, Ve = N + dt + H + rt + Te, ze = F + B + K + Q + Ce;
+    return n.set(0, 0, $e), n.set(0, 1, qe), n.set(0, 2, Be), n.set(1, 0, Oe), n.set(1, 1, Ae), n.set(1, 2, Ne), n.set(2, 0, Le), n.set(2, 1, Ve), n.set(2, 2, ze), n;
+  }
+  mmulStrassen(t) {
+    t = j.checkMatrix(t);
+    let n = this.clone(), e = n.rows, s = n.columns, r = t.rows, o = t.columns;
+    s !== r && console.warn(
+      `Multiplying ${e} x ${s} and ${r} x ${o} matrix: dimensions do not match.`
+    );
+    function h(a, g, w) {
+      let m = a.rows, d = a.columns;
+      if (m === g && d === w)
+        return a;
+      {
+        let P = v.zeros(g, w);
+        return P = P.setSubMatrix(a, 0, 0), P;
+      }
+    }
+    let u = Math.max(e, r), l = Math.max(s, o);
+    n = h(n, u, l), t = h(t, u, l);
+    function c(a, g, w, m) {
+      if (w <= 512 || m <= 512)
+        return a.mmul(g);
+      w % 2 === 1 && m % 2 === 1 ? (a = h(a, w + 1, m + 1), g = h(g, w + 1, m + 1)) : w % 2 === 1 ? (a = h(a, w + 1, m), g = h(g, w + 1, m)) : m % 2 === 1 && (a = h(a, w, m + 1), g = h(g, w, m + 1));
+      let d = parseInt(a.rows / 2, 10), P = parseInt(a.columns / 2, 10), f = a.subMatrix(0, d - 1, 0, P - 1), p = g.subMatrix(0, d - 1, 0, P - 1), y = a.subMatrix(0, d - 1, P, a.columns - 1), M = g.subMatrix(0, d - 1, P, g.columns - 1), b = a.subMatrix(d, a.rows - 1, 0, P - 1), x = g.subMatrix(d, g.rows - 1, 0, P - 1), k = a.subMatrix(d, a.rows - 1, P, a.columns - 1), S = g.subMatrix(d, g.rows - 1, P, g.columns - 1), _ = c(
+        v.add(f, k),
+        v.add(p, S),
+        d,
+        P
+      ), F = c(v.add(b, k), p, d, P), B = c(f, v.sub(M, S), d, P), K = c(k, v.sub(x, p), d, P), Q = c(v.add(f, y), S, d, P), Z = c(
+        v.sub(b, f),
+        v.add(p, M),
+        d,
+        P
+      ), $ = c(
+        v.sub(y, k),
+        v.add(x, S),
+        d,
+        P
+      ), N = v.add(_, K);
+      N.sub(Q), N.add($);
+      let dt = v.add(B, Q), H = v.add(F, K), rt = v.sub(_, F);
+      rt.add(B), rt.add(Z);
+      let L = v.zeros(2 * N.rows, 2 * N.columns);
+      return L = L.setSubMatrix(N, 0, 0), L = L.setSubMatrix(dt, N.rows, 0), L = L.setSubMatrix(H, 0, N.columns), L = L.setSubMatrix(rt, N.rows, N.columns), L.subMatrix(0, w - 1, 0, m - 1);
+    }
+    return c(n, t, u, l);
+  }
+  scaleRows(t = {}) {
+    if (typeof t != "object")
+      throw new TypeError("options must be an object");
+    const { min: n = 0, max: e = 1 } = t;
+    if (!Number.isFinite(n))
+      throw new TypeError("min must be a number");
+    if (!Number.isFinite(e))
+      throw new TypeError("max must be a number");
+    if (n >= e)
+      throw new RangeError("min must be smaller than max");
+    let s = new j(this.rows, this.columns);
+    for (let r = 0; r < this.rows; r++) {
+      const o = this.getRow(r);
+      o.length > 0 && Kt(o, { min: n, max: e, output: o }), s.setRow(r, o);
+    }
+    return s;
+  }
+  scaleColumns(t = {}) {
+    if (typeof t != "object")
+      throw new TypeError("options must be an object");
+    const { min: n = 0, max: e = 1 } = t;
+    if (!Number.isFinite(n))
+      throw new TypeError("min must be a number");
+    if (!Number.isFinite(e))
+      throw new TypeError("max must be a number");
+    if (n >= e)
+      throw new RangeError("min must be smaller than max");
+    let s = new j(this.rows, this.columns);
+    for (let r = 0; r < this.columns; r++) {
+      const o = this.getColumn(r);
+      o.length && Kt(o, {
+        min: n,
+        max: e,
+        output: o
+      }), s.setColumn(r, o);
+    }
+    return s;
+  }
+  flipRows() {
+    const t = Math.ceil(this.columns / 2);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < t; e++) {
+        let s = this.get(n, e), r = this.get(n, this.columns - 1 - e);
+        this.set(n, e, r), this.set(n, this.columns - 1 - e, s);
+      }
+    return this;
+  }
+  flipColumns() {
+    const t = Math.ceil(this.rows / 2);
+    for (let n = 0; n < this.columns; n++)
+      for (let e = 0; e < t; e++) {
+        let s = this.get(e, n), r = this.get(this.rows - 1 - e, n);
+        this.set(e, n, r), this.set(this.rows - 1 - e, n, s);
+      }
+    return this;
+  }
+  kroneckerProduct(t) {
+    t = j.checkMatrix(t);
+    let n = this.rows, e = this.columns, s = t.rows, r = t.columns, o = new j(n * s, e * r);
+    for (let h = 0; h < n; h++)
+      for (let u = 0; u < e; u++)
+        for (let l = 0; l < s; l++)
+          for (let c = 0; c < r; c++)
+            o.set(s * h + l, r * u + c, this.get(h, u) * t.get(l, c));
+    return o;
+  }
+  kroneckerSum(t) {
+    if (t = j.checkMatrix(t), !this.isSquare() || !t.isSquare())
+      throw new Error("Kronecker Sum needs two Square Matrices");
+    let n = this.rows, e = t.rows, s = this.kroneckerProduct(j.eye(e, e)), r = j.eye(n, n).kroneckerProduct(t);
+    return s.add(r);
+  }
+  transpose() {
+    let t = new j(this.columns, this.rows);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        t.set(e, n, this.get(n, e));
+    return t;
+  }
+  sortRows(t = Ht) {
+    for (let n = 0; n < this.rows; n++)
+      this.setRow(n, this.getRow(n).sort(t));
+    return this;
+  }
+  sortColumns(t = Ht) {
+    for (let n = 0; n < this.columns; n++)
+      this.setColumn(n, this.getColumn(n).sort(t));
+    return this;
+  }
+  subMatrix(t, n, e, s) {
+    Zt(this, t, n, e, s);
+    let r = new j(
+      n - t + 1,
+      s - e + 1
+    );
+    for (let o = t; o <= n; o++)
+      for (let h = e; h <= s; h++)
+        r.set(o - t, h - e, this.get(o, h));
+    return r;
+  }
+  subMatrixRow(t, n, e) {
+    if (n === void 0 && (n = 0), e === void 0 && (e = this.columns - 1), n > e || n < 0 || n >= this.columns || e < 0 || e >= this.columns)
+      throw new RangeError("Argument out of range");
+    let s = new j(t.length, e - n + 1);
+    for (let r = 0; r < t.length; r++)
+      for (let o = n; o <= e; o++) {
+        if (t[r] < 0 || t[r] >= this.rows)
+          throw new RangeError(`Row index out of range: ${t[r]}`);
+        s.set(r, o - n, this.get(t[r], o));
+      }
+    return s;
+  }
+  subMatrixColumn(t, n, e) {
+    if (n === void 0 && (n = 0), e === void 0 && (e = this.rows - 1), n > e || n < 0 || n >= this.rows || e < 0 || e >= this.rows)
+      throw new RangeError("Argument out of range");
+    let s = new j(e - n + 1, t.length);
+    for (let r = 0; r < t.length; r++)
+      for (let o = n; o <= e; o++) {
+        if (t[r] < 0 || t[r] >= this.columns)
+          throw new RangeError(`Column index out of range: ${t[r]}`);
+        s.set(o - n, r, this.get(o, t[r]));
+      }
+    return s;
+  }
+  setSubMatrix(t, n, e) {
+    if (t = j.checkMatrix(t), t.isEmpty())
+      return this;
+    let s = n + t.rows - 1, r = e + t.columns - 1;
+    Zt(this, n, s, e, r);
+    for (let o = 0; o < t.rows; o++)
+      for (let h = 0; h < t.columns; h++)
+        this.set(n + o, e + h, t.get(o, h));
+    return this;
+  }
+  selection(t, n) {
+    an(this, t), cn(this, n);
+    let e = new j(t.length, n.length);
+    for (let s = 0; s < t.length; s++) {
+      let r = t[s];
+      for (let o = 0; o < n.length; o++) {
+        let h = n[o];
+        e.set(s, o, this.get(r, h));
+      }
+    }
+    return e;
+  }
+  trace() {
+    let t = Math.min(this.rows, this.columns), n = 0;
+    for (let e = 0; e < t; e++)
+      n += this.get(e, e);
+    return n;
+  }
+  clone() {
+    let t = new j(this.rows, this.columns);
+    for (let n = 0; n < this.rows; n++)
+      for (let e = 0; e < this.columns; e++)
+        t.set(n, e, this.get(n, e));
+    return t;
+  }
+  sum(t) {
+    switch (t) {
+      case "row":
+        return fn(this);
+      case "column":
+        return gn(this);
+      case void 0:
+        return mn(this);
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  product(t) {
+    switch (t) {
+      case "row":
+        return pn(this);
+      case "column":
+        return wn(this);
+      case void 0:
+        return dn(this);
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  mean(t) {
+    const n = this.sum(t);
+    switch (t) {
+      case "row": {
+        for (let e = 0; e < this.rows; e++)
+          n[e] /= this.columns;
+        return n;
+      }
+      case "column": {
+        for (let e = 0; e < this.columns; e++)
+          n[e] /= this.rows;
+        return n;
+      }
+      case void 0:
+        return n / this.size;
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  variance(t, n = {}) {
+    if (typeof t == "object" && (n = t, t = void 0), typeof n != "object")
+      throw new TypeError("options must be an object");
+    const { unbiased: e = !0, mean: s = this.mean(t) } = n;
+    if (typeof e != "boolean")
+      throw new TypeError("unbiased must be a boolean");
+    switch (t) {
+      case "row": {
+        if (!z(s))
+          throw new TypeError("mean must be an array");
+        return yn(this, e, s);
+      }
+      case "column": {
+        if (!z(s))
+          throw new TypeError("mean must be an array");
+        return Pn(this, e, s);
+      }
+      case void 0: {
+        if (typeof s != "number")
+          throw new TypeError("mean must be a number");
+        return Mn(this, e, s);
+      }
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  standardDeviation(t, n) {
+    typeof t == "object" && (n = t, t = void 0);
+    const e = this.variance(t, n);
+    if (t === void 0)
+      return Math.sqrt(e);
+    for (let s = 0; s < e.length; s++)
+      e[s] = Math.sqrt(e[s]);
+    return e;
+  }
+  center(t, n = {}) {
+    if (typeof t == "object" && (n = t, t = void 0), typeof n != "object")
+      throw new TypeError("options must be an object");
+    const { center: e = this.mean(t) } = n;
+    switch (t) {
+      case "row": {
+        if (!z(e))
+          throw new TypeError("center must be an array");
+        return En(this, e), this;
+      }
+      case "column": {
+        if (!z(e))
+          throw new TypeError("center must be an array");
+        return Sn(this, e), this;
+      }
+      case void 0: {
+        if (typeof e != "number")
+          throw new TypeError("center must be a number");
+        return bn(this, e), this;
+      }
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  scale(t, n = {}) {
+    if (typeof t == "object" && (n = t, t = void 0), typeof n != "object")
+      throw new TypeError("options must be an object");
+    let e = n.scale;
+    switch (t) {
+      case "row": {
+        if (e === void 0)
+          e = jn(this);
+        else if (!z(e))
+          throw new TypeError("scale must be an array");
+        return vn(this, e), this;
+      }
+      case "column": {
+        if (e === void 0)
+          e = kn(this);
+        else if (!z(e))
+          throw new TypeError("scale must be an array");
+        return In(this, e), this;
+      }
+      case void 0: {
+        if (e === void 0)
+          e = xn(this);
+        else if (typeof e != "number")
+          throw new TypeError("scale must be a number");
+        return Rn(this, e), this;
+      }
+      default:
+        throw new Error(`invalid option: ${t}`);
+    }
+  }
+  toString(t) {
+    return we(this, t);
+  }
+}
+v.prototype.klass = "Matrix";
+typeof Symbol < "u" && (v.prototype[Symbol.for("nodejs.util.inspect.custom")] = on);
+function Ht(i, t) {
+  return i - t;
+}
+function _n(i) {
+  return i.every((t) => typeof t == "number");
+}
+v.random = v.rand;
+v.randomInt = v.randInt;
+v.diagonal = v.diag;
+v.prototype.diagonal = v.prototype.diag;
+v.identity = v.eye;
+v.prototype.negate = v.prototype.neg;
+v.prototype.tensorProduct = v.prototype.kroneckerProduct;
+class j extends v {
+  constructor(t, n) {
+    if (super(), j.isMatrix(t))
+      return t.clone();
+    if (Number.isInteger(t) && t >= 0)
+      if (this.data = [], Number.isInteger(n) && n >= 0)
+        for (let e = 0; e < t; e++)
+          this.data.push(new Float64Array(n));
+      else
+        throw new TypeError("nColumns must be a positive integer");
+    else if (z(t)) {
+      const e = t;
+      if (t = e.length, n = t ? e[0].length : 0, typeof n != "number")
+        throw new TypeError(
+          "Data must be a 2D array with at least one element"
+        );
+      this.data = [];
+      for (let s = 0; s < t; s++) {
+        if (e[s].length !== n)
+          throw new RangeError("Inconsistent array dimensions");
+        if (!_n(e[s]))
+          throw new TypeError("Input data contains non-numeric values");
+        this.data.push(Float64Array.from(e[s]));
+      }
+    } else
+      throw new TypeError(
+        "First argument must be a positive number or an array"
+      );
+    this.rows = t, this.columns = n;
+  }
+  set(t, n, e) {
+    return this.data[t][n] = e, this;
+  }
+  get(t, n) {
+    return this.data[t][n];
+  }
+  removeRow(t) {
+    return Y(this, t), this.data.splice(t, 1), this.rows -= 1, this;
+  }
+  addRow(t, n) {
+    return n === void 0 && (n = t, t = this.rows), Y(this, t, !0), n = Float64Array.from(lt(this, n)), this.data.splice(t, 0, n), this.rows += 1, this;
+  }
+  removeColumn(t) {
+    X(this, t);
+    for (let n = 0; n < this.rows; n++) {
+      const e = new Float64Array(this.columns - 1);
+      for (let s = 0; s < t; s++)
+        e[s] = this.data[n][s];
+      for (let s = t + 1; s < this.columns; s++)
+        e[s - 1] = this.data[n][s];
+      this.data[n] = e;
+    }
+    return this.columns -= 1, this;
+  }
+  addColumn(t, n) {
+    typeof n > "u" && (n = t, t = this.columns), X(this, t, !0), n = ut(this, n);
+    for (let e = 0; e < this.rows; e++) {
+      const s = new Float64Array(this.columns + 1);
+      let r = 0;
+      for (; r < t; r++)
+        s[r] = this.data[e][r];
+      for (s[r++] = n[e]; r < this.columns + 1; r++)
+        s[r] = this.data[e][r - 1];
+      this.data[e] = s;
+    }
+    return this.columns += 1, this;
+  }
+}
+un(v, j);
+class Fn extends v {
+  constructor(t) {
+    super(), this.data = t, this.rows = t.length, this.columns = t[0].length;
+  }
+  set(t, n, e) {
+    return this.data[t][n] = e, this;
+  }
+  get(t, n) {
+    return this.data[t][n];
+  }
+}
+function te(i, t) {
+  let n = 0;
+  return Math.abs(i) > Math.abs(t) ? (n = t / i, Math.abs(i) * Math.sqrt(1 + n * n)) : t !== 0 ? (n = i / t, Math.abs(t) * Math.sqrt(1 + n * n)) : 0;
+}
+class Tn {
+  constructor(t, n = {}) {
+    const { assumeSymmetric: e = !1 } = n;
+    if (t = Fn.checkMatrix(t), !t.isSquare())
+      throw new Error("Matrix is not a square matrix");
+    if (t.isEmpty())
+      throw new Error("Matrix must be non-empty");
+    let s = t.columns, r = new j(s, s), o = new Float64Array(s), h = new Float64Array(s), u = t, l, c, a = !1;
+    if (e ? a = !0 : a = t.isSymmetric(), a) {
+      for (l = 0; l < s; l++)
+        for (c = 0; c < s; c++)
+          r.set(l, c, u.get(l, c));
+      Cn(s, h, o, r), $n(s, h, o, r);
+    } else {
+      let g = new j(s, s), w = new Float64Array(s);
+      for (c = 0; c < s; c++)
+        for (l = 0; l < s; l++)
+          g.set(l, c, u.get(l, c));
+      qn(s, g, w, r), Bn(s, h, o, r, g);
+    }
+    this.n = s, this.e = h, this.d = o, this.V = r;
+  }
+  get realEigenvalues() {
+    return Array.from(this.d);
+  }
+  get imaginaryEigenvalues() {
+    return Array.from(this.e);
+  }
+  get eigenvectorMatrix() {
+    return this.V;
+  }
+  get diagonalMatrix() {
+    let t = this.n, n = this.e, e = this.d, s = new j(t, t), r, o;
+    for (r = 0; r < t; r++) {
+      for (o = 0; o < t; o++)
+        s.set(r, o, 0);
+      s.set(r, r, e[r]), n[r] > 0 ? s.set(r, r + 1, n[r]) : n[r] < 0 && s.set(r, r - 1, n[r]);
+    }
+    return s;
+  }
+}
+function Cn(i, t, n, e) {
+  let s, r, o, h, u, l, c, a;
+  for (u = 0; u < i; u++)
+    n[u] = e.get(i - 1, u);
+  for (h = i - 1; h > 0; h--) {
+    for (a = 0, o = 0, l = 0; l < h; l++)
+      a = a + Math.abs(n[l]);
+    if (a === 0)
+      for (t[h] = n[h - 1], u = 0; u < h; u++)
+        n[u] = e.get(h - 1, u), e.set(h, u, 0), e.set(u, h, 0);
+    else {
+      for (l = 0; l < h; l++)
+        n[l] /= a, o += n[l] * n[l];
+      for (s = n[h - 1], r = Math.sqrt(o), s > 0 && (r = -r), t[h] = a * r, o = o - s * r, n[h - 1] = s - r, u = 0; u < h; u++)
+        t[u] = 0;
+      for (u = 0; u < h; u++) {
+        for (s = n[u], e.set(u, h, s), r = t[u] + e.get(u, u) * s, l = u + 1; l <= h - 1; l++)
+          r += e.get(l, u) * n[l], t[l] += e.get(l, u) * s;
+        t[u] = r;
+      }
+      for (s = 0, u = 0; u < h; u++)
+        t[u] /= o, s += t[u] * n[u];
+      for (c = s / (o + o), u = 0; u < h; u++)
+        t[u] -= c * n[u];
+      for (u = 0; u < h; u++) {
+        for (s = n[u], r = t[u], l = u; l <= h - 1; l++)
+          e.set(l, u, e.get(l, u) - (s * t[l] + r * n[l]));
+        n[u] = e.get(h - 1, u), e.set(h, u, 0);
+      }
+    }
+    n[h] = o;
+  }
+  for (h = 0; h < i - 1; h++) {
+    if (e.set(i - 1, h, e.get(h, h)), e.set(h, h, 1), o = n[h + 1], o !== 0) {
+      for (l = 0; l <= h; l++)
+        n[l] = e.get(l, h + 1) / o;
+      for (u = 0; u <= h; u++) {
+        for (r = 0, l = 0; l <= h; l++)
+          r += e.get(l, h + 1) * e.get(l, u);
+        for (l = 0; l <= h; l++)
+          e.set(l, u, e.get(l, u) - r * n[l]);
+      }
+    }
+    for (l = 0; l <= h; l++)
+      e.set(l, h + 1, 0);
+  }
+  for (u = 0; u < i; u++)
+    n[u] = e.get(i - 1, u), e.set(i - 1, u, 0);
+  e.set(i - 1, i - 1, 1), t[0] = 0;
+}
+function $n(i, t, n, e) {
+  let s, r, o, h, u, l, c, a, g, w, m, d, P, f, p, y;
+  for (o = 1; o < i; o++)
+    t[o - 1] = t[o];
+  t[i - 1] = 0;
+  let M = 0, b = 0, x = Number.EPSILON;
+  for (l = 0; l < i; l++) {
+    for (b = Math.max(b, Math.abs(n[l]) + Math.abs(t[l])), c = l; c < i && !(Math.abs(t[c]) <= x * b); )
+      c++;
+    if (c > l)
+      do {
+        for (s = n[l], a = (n[l + 1] - s) / (2 * t[l]), g = te(a, 1), a < 0 && (g = -g), n[l] = t[l] / (a + g), n[l + 1] = t[l] * (a + g), w = n[l + 1], r = s - n[l], o = l + 2; o < i; o++)
+          n[o] -= r;
+        for (M = M + r, a = n[c], m = 1, d = m, P = m, f = t[l + 1], p = 0, y = 0, o = c - 1; o >= l; o--)
+          for (P = d, d = m, y = p, s = m * t[o], r = m * a, g = te(a, t[o]), t[o + 1] = p * g, p = t[o] / g, m = a / g, a = m * n[o] - p * s, n[o + 1] = r + p * (m * s + p * n[o]), u = 0; u < i; u++)
+            r = e.get(u, o + 1), e.set(u, o + 1, p * e.get(u, o) + m * r), e.set(u, o, m * e.get(u, o) - p * r);
+        a = -p * y * P * f * t[l] / w, t[l] = p * a, n[l] = m * a;
+      } while (Math.abs(t[l]) > x * b);
+    n[l] = n[l] + M, t[l] = 0;
+  }
+  for (o = 0; o < i - 1; o++) {
+    for (u = o, a = n[o], h = o + 1; h < i; h++)
+      n[h] < a && (u = h, a = n[h]);
+    if (u !== o)
+      for (n[u] = n[o], n[o] = a, h = 0; h < i; h++)
+        a = e.get(h, o), e.set(h, o, e.get(h, u)), e.set(h, u, a);
+  }
+}
+function qn(i, t, n, e) {
+  let s = 0, r = i - 1, o, h, u, l, c, a, g;
+  for (a = s + 1; a <= r - 1; a++) {
+    for (g = 0, l = a; l <= r; l++)
+      g = g + Math.abs(t.get(l, a - 1));
+    if (g !== 0) {
+      for (u = 0, l = r; l >= a; l--)
+        n[l] = t.get(l, a - 1) / g, u += n[l] * n[l];
+      for (h = Math.sqrt(u), n[a] > 0 && (h = -h), u = u - n[a] * h, n[a] = n[a] - h, c = a; c < i; c++) {
+        for (o = 0, l = r; l >= a; l--)
+          o += n[l] * t.get(l, c);
+        for (o = o / u, l = a; l <= r; l++)
+          t.set(l, c, t.get(l, c) - o * n[l]);
+      }
+      for (l = 0; l <= r; l++) {
+        for (o = 0, c = r; c >= a; c--)
+          o += n[c] * t.get(l, c);
+        for (o = o / u, c = a; c <= r; c++)
+          t.set(l, c, t.get(l, c) - o * n[c]);
+      }
+      n[a] = g * n[a], t.set(a, a - 1, g * h);
+    }
+  }
+  for (l = 0; l < i; l++)
+    for (c = 0; c < i; c++)
+      e.set(l, c, l === c ? 1 : 0);
+  for (a = r - 1; a >= s + 1; a--)
+    if (t.get(a, a - 1) !== 0) {
+      for (l = a + 1; l <= r; l++)
+        n[l] = t.get(l, a - 1);
+      for (c = a; c <= r; c++) {
+        for (h = 0, l = a; l <= r; l++)
+          h += n[l] * e.get(l, c);
+        for (h = h / n[a] / t.get(a, a - 1), l = a; l <= r; l++)
+          e.set(l, c, e.get(l, c) + h * n[l]);
+      }
+    }
+}
+function Bn(i, t, n, e, s) {
+  let r = i - 1, o = 0, h = i - 1, u = Number.EPSILON, l = 0, c = 0, a = 0, g = 0, w = 0, m = 0, d = 0, P = 0, f, p, y, M, b, x, k, S, _, F, B, K, Q, Z, $;
+  for (f = 0; f < i; f++)
+    for ((f < o || f > h) && (n[f] = s.get(f, f), t[f] = 0), p = Math.max(f - 1, 0); p < i; p++)
+      c = c + Math.abs(s.get(f, p));
+  for (; r >= o; ) {
+    for (M = r; M > o && (m = Math.abs(s.get(M - 1, M - 1)) + Math.abs(s.get(M, M)), m === 0 && (m = c), !(Math.abs(s.get(M, M - 1)) < u * m)); )
+      M--;
+    if (M === r)
+      s.set(r, r, s.get(r, r) + l), n[r] = s.get(r, r), t[r] = 0, r--, P = 0;
+    else if (M === r - 1) {
+      if (k = s.get(r, r - 1) * s.get(r - 1, r), a = (s.get(r - 1, r - 1) - s.get(r, r)) / 2, g = a * a + k, d = Math.sqrt(Math.abs(g)), s.set(r, r, s.get(r, r) + l), s.set(r - 1, r - 1, s.get(r - 1, r - 1) + l), S = s.get(r, r), g >= 0) {
+        for (d = a >= 0 ? a + d : a - d, n[r - 1] = S + d, n[r] = n[r - 1], d !== 0 && (n[r] = S - k / d), t[r - 1] = 0, t[r] = 0, S = s.get(r, r - 1), m = Math.abs(S) + Math.abs(d), a = S / m, g = d / m, w = Math.sqrt(a * a + g * g), a = a / w, g = g / w, p = r - 1; p < i; p++)
+          d = s.get(r - 1, p), s.set(r - 1, p, g * d + a * s.get(r, p)), s.set(r, p, g * s.get(r, p) - a * d);
+        for (f = 0; f <= r; f++)
+          d = s.get(f, r - 1), s.set(f, r - 1, g * d + a * s.get(f, r)), s.set(f, r, g * s.get(f, r) - a * d);
+        for (f = o; f <= h; f++)
+          d = e.get(f, r - 1), e.set(f, r - 1, g * d + a * e.get(f, r)), e.set(f, r, g * e.get(f, r) - a * d);
+      } else
+        n[r - 1] = S + a, n[r] = S + a, t[r - 1] = d, t[r] = -d;
+      r = r - 2, P = 0;
+    } else {
+      if (S = s.get(r, r), _ = 0, k = 0, M < r && (_ = s.get(r - 1, r - 1), k = s.get(r, r - 1) * s.get(r - 1, r)), P === 10) {
+        for (l += S, f = o; f <= r; f++)
+          s.set(f, f, s.get(f, f) - S);
+        m = Math.abs(s.get(r, r - 1)) + Math.abs(s.get(r - 1, r - 2)), S = _ = 0.75 * m, k = -0.4375 * m * m;
+      }
+      if (P === 30 && (m = (_ - S) / 2, m = m * m + k, m > 0)) {
+        for (m = Math.sqrt(m), _ < S && (m = -m), m = S - k / ((_ - S) / 2 + m), f = o; f <= r; f++)
+          s.set(f, f, s.get(f, f) - m);
+        l += m, S = _ = k = 0.964;
+      }
+      for (P = P + 1, b = r - 2; b >= M && (d = s.get(b, b), w = S - d, m = _ - d, a = (w * m - k) / s.get(b + 1, b) + s.get(b, b + 1), g = s.get(b + 1, b + 1) - d - w - m, w = s.get(b + 2, b + 1), m = Math.abs(a) + Math.abs(g) + Math.abs(w), a = a / m, g = g / m, w = w / m, !(b === M || Math.abs(s.get(b, b - 1)) * (Math.abs(g) + Math.abs(w)) < u * (Math.abs(a) * (Math.abs(s.get(b - 1, b - 1)) + Math.abs(d) + Math.abs(s.get(b + 1, b + 1)))))); )
+        b--;
+      for (f = b + 2; f <= r; f++)
+        s.set(f, f - 2, 0), f > b + 2 && s.set(f, f - 3, 0);
+      for (y = b; y <= r - 1 && (Z = y !== r - 1, y !== b && (a = s.get(y, y - 1), g = s.get(y + 1, y - 1), w = Z ? s.get(y + 2, y - 1) : 0, S = Math.abs(a) + Math.abs(g) + Math.abs(w), S !== 0 && (a = a / S, g = g / S, w = w / S)), S !== 0); y++)
+        if (m = Math.sqrt(a * a + g * g + w * w), a < 0 && (m = -m), m !== 0) {
+          for (y !== b ? s.set(y, y - 1, -m * S) : M !== b && s.set(y, y - 1, -s.get(y, y - 1)), a = a + m, S = a / m, _ = g / m, d = w / m, g = g / a, w = w / a, p = y; p < i; p++)
+            a = s.get(y, p) + g * s.get(y + 1, p), Z && (a = a + w * s.get(y + 2, p), s.set(y + 2, p, s.get(y + 2, p) - a * d)), s.set(y, p, s.get(y, p) - a * S), s.set(y + 1, p, s.get(y + 1, p) - a * _);
+          for (f = 0; f <= Math.min(r, y + 3); f++)
+            a = S * s.get(f, y) + _ * s.get(f, y + 1), Z && (a = a + d * s.get(f, y + 2), s.set(f, y + 2, s.get(f, y + 2) - a * w)), s.set(f, y, s.get(f, y) - a), s.set(f, y + 1, s.get(f, y + 1) - a * g);
+          for (f = o; f <= h; f++)
+            a = S * e.get(f, y) + _ * e.get(f, y + 1), Z && (a = a + d * e.get(f, y + 2), e.set(f, y + 2, e.get(f, y + 2) - a * w)), e.set(f, y, e.get(f, y) - a), e.set(f, y + 1, e.get(f, y + 1) - a * g);
+        }
+    }
+  }
+  if (c !== 0) {
+    for (r = i - 1; r >= 0; r--)
+      if (a = n[r], g = t[r], g === 0)
+        for (M = r, s.set(r, r, 1), f = r - 1; f >= 0; f--) {
+          for (k = s.get(f, f) - a, w = 0, p = M; p <= r; p++)
+            w = w + s.get(f, p) * s.get(p, r);
+          if (t[f] < 0)
+            d = k, m = w;
+          else if (M = f, t[f] === 0 ? s.set(f, r, k !== 0 ? -w / k : -w / (u * c)) : (S = s.get(f, f + 1), _ = s.get(f + 1, f), g = (n[f] - a) * (n[f] - a) + t[f] * t[f], x = (S * m - d * w) / g, s.set(f, r, x), s.set(
+            f + 1,
+            r,
+            Math.abs(S) > Math.abs(d) ? (-w - k * x) / S : (-m - _ * x) / d
+          )), x = Math.abs(s.get(f, r)), u * x * x > 1)
+            for (p = f; p <= r; p++)
+              s.set(p, r, s.get(p, r) / x);
+        }
+      else if (g < 0)
+        for (M = r - 1, Math.abs(s.get(r, r - 1)) > Math.abs(s.get(r - 1, r)) ? (s.set(r - 1, r - 1, g / s.get(r, r - 1)), s.set(r - 1, r, -(s.get(r, r) - a) / s.get(r, r - 1))) : ($ = St(0, -s.get(r - 1, r), s.get(r - 1, r - 1) - a, g), s.set(r - 1, r - 1, $[0]), s.set(r - 1, r, $[1])), s.set(r, r - 1, 0), s.set(r, r, 1), f = r - 2; f >= 0; f--) {
+          for (F = 0, B = 0, p = M; p <= r; p++)
+            F = F + s.get(f, p) * s.get(p, r - 1), B = B + s.get(f, p) * s.get(p, r);
+          if (k = s.get(f, f) - a, t[f] < 0)
+            d = k, w = F, m = B;
+          else if (M = f, t[f] === 0 ? ($ = St(-F, -B, k, g), s.set(f, r - 1, $[0]), s.set(f, r, $[1])) : (S = s.get(f, f + 1), _ = s.get(f + 1, f), K = (n[f] - a) * (n[f] - a) + t[f] * t[f] - g * g, Q = (n[f] - a) * 2 * g, K === 0 && Q === 0 && (K = u * c * (Math.abs(k) + Math.abs(g) + Math.abs(S) + Math.abs(_) + Math.abs(d))), $ = St(
+            S * w - d * F + g * B,
+            S * m - d * B - g * F,
+            K,
+            Q
+          ), s.set(f, r - 1, $[0]), s.set(f, r, $[1]), Math.abs(S) > Math.abs(d) + Math.abs(g) ? (s.set(
+            f + 1,
+            r - 1,
+            (-F - k * s.get(f, r - 1) + g * s.get(f, r)) / S
+          ), s.set(
+            f + 1,
+            r,
+            (-B - k * s.get(f, r) - g * s.get(f, r - 1)) / S
+          )) : ($ = St(
+            -w - _ * s.get(f, r - 1),
+            -m - _ * s.get(f, r),
+            d,
+            g
+          ), s.set(f + 1, r - 1, $[0]), s.set(f + 1, r, $[1]))), x = Math.max(Math.abs(s.get(f, r - 1)), Math.abs(s.get(f, r))), u * x * x > 1)
+            for (p = f; p <= r; p++)
+              s.set(p, r - 1, s.get(p, r - 1) / x), s.set(p, r, s.get(p, r) / x);
+        }
+    for (f = 0; f < i; f++)
+      if (f < o || f > h)
+        for (p = f; p < i; p++)
+          e.set(f, p, s.get(f, p));
+    for (p = i - 1; p >= o; p--)
+      for (f = o; f <= h; f++) {
+        for (d = 0, y = o; y <= Math.min(p, h); y++)
+          d = d + e.get(f, y) * s.get(y, p);
+        e.set(f, p, d);
+      }
+  }
+}
+function St(i, t, n, e) {
+  let s, r;
+  return Math.abs(n) > Math.abs(e) ? (s = e / n, r = n + s * e, [(i + s * t) / r, (t - s * i) / r]) : (s = n / e, r = e + s * n, [(s * i + t) / r, (s * t - i) / r]);
+}
+function de(i, t = 1e-9) {
+  if (i[i.length - 1] == 0)
+    return de(i.slice(0, -1), t);
+  const n = new j([
+    i.slice(0, -1).reverse().map((s) => -s / i[i.length - 1]),
+    ...Array.from({ length: i.length - 2 }, (s, r) => {
+      const o = new Array(i.length - 1).fill(0);
+      return o[r] = 1, o;
+    })
+  ]), e = new Tn(n);
+  return e.realEigenvalues.filter(
+    (s, r) => Math.abs(e.imaginaryEigenvalues[r]) < t
+  );
+}
+function ye(i, t) {
+  const n = Math.max(i.precision, t.precision), e = i.coefficients, s = e.x2, r = e.xy, o = e.y2, h = e.x, u = e.y, l = e.c, c = t.coefficients, a = c.x2, g = c.xy, w = c.y2, m = c.x, d = c.y, P = c.c, f = {
+    z0: l * s * m * m + s * s * P * P - h * s * m * P + a * a * l * l - 2 * s * P * a * l - h * m * a * l + a * h * h * P,
+    z1: d * h * h * a - P * m * s * r - 2 * s * P * a * u - l * a * g * h + 2 * m * g * s * l + 2 * d * P * s * s + m * m * s * u - d * m * s * h - 2 * s * d * a * l - l * a * m * r + 2 * l * u * a * a - P * g * s * h - u * a * m * h + 2 * P * r * a * h,
+    z2: d * d * s * s + 2 * w * P * s * s - u * a * m * r + P * a * r * r - u * a * g * h - P * g * s * r - 2 * s * d * a * u + 2 * m * g * s * u - w * m * s * h - 2 * s * w * a * l + g * g * s * l + 2 * d * r * a * h + u * u * a * a - o * a * m * h - d * g * s * h + 2 * l * o * a * a - l * a * g * r + w * h * h * a + m * m * s * o - d * m * s * r - 2 * s * P * a * o,
+    z3: -2 * s * a * o * d + d * a * r * r + 2 * w * r * a * h - o * a * g * h + g * g * s * u - d * g * s * r - 2 * s * w * a * u - u * a * g * r - w * g * s * h + 2 * d * w * s * s + 2 * u * o * a * a - o * a * m * r + 2 * m * g * s * o - w * m * s * r,
+    z4: s * s * w * w - 2 * s * w * a * o + a * a * o * o - r * s * g * w - r * g * a * o + r * r * a * w + o * s * g * g
+  }, y = de(
+    [f.z0, f.z1, f.z2, f.z3, f.z4],
+    n
+  ).flatMap((M) => {
+    const b = s * g * M + s * m - a * r * M - a * h;
+    if (b)
+      return [[-(s * P + s * w * M * M - a * o * M * M + s * d * M - a * u * M - a * l) / b, M]];
+    const x = r * M + h, k = -x / (2 * s), S = o * M * M + u * M + l, _ = x * x / (4 * s * s) - S / s;
+    if (Math.abs(_) < n)
+      return [[k, M]];
+    if (_ > 0) {
+      const F = Math.sqrt(_);
+      return [[k + F, M], [k - F, M]];
+    }
+    return [];
+  });
+  return ct(y, n);
+}
+function ee(i, t) {
+  return ye(i, t).filter((e) => i.isOnSegment(e) && t.isOnSegment(e));
+}
+const On = (i) => {
+  const {
+    firstPoint: t,
+    lastPoint: n,
+    center: e,
+    majorRadius: s,
+    minorRadius: r,
+    tiltAngle: o,
+    clockwise: h
+  } = i;
+  return new C(
+    n,
+    t,
+    e,
+    s,
+    r,
+    o,
+    h,
+    {
+      ignoreChecks: !0,
+      angleUnits: "rad"
+    }
+  );
+}, An = (i, t) => {
+  if (i.isSame(t))
+    return [i];
+  const n = (s, r) => new C(
+    s,
+    r,
+    i.center,
+    i.majorRadius,
+    i.minorRadius,
+    i.tiltAngle,
+    i.clockwise,
+    { ignoreChecks: !0, angleUnits: "rad" }
+  ), e = ct(
+    [
+      t.isOnSegment(i.firstPoint) ? i.firstPoint : null,
+      t.isOnSegment(i.lastPoint) ? i.lastPoint : null,
+      i.isOnSegment(t.firstPoint) ? t.firstPoint : null,
+      i.isOnSegment(t.lastPoint) ? t.lastPoint : null
+    ].filter((s) => s !== null)
+    // We sort by the param value of the first arc. This means that the points
+    // will be sorted with the same orientation than arc1
+  ).sort((s, r) => i.pointToParam(s) - i.pointToParam(r));
+  if (e.length === 0)
+    return [];
+  if (e.length === 1)
+    return [];
+  if (e.length === 2)
+    return i.isSame(On(t)) ? [] : [n(e[0], e[1])];
+  if (e.length === 3) {
+    const s = E(e[0], t.lastPoint) || E(e[0], t.firstPoint) ? 1 : 0;
+    return [n(e[0 + s], e[1 + s])];
+  } else if (e.length === 4)
+    return [n(e[0], e[1]), n(e[2], e[3])];
+  throw new Error("Bug in the ellipse arc ellipse arc overlap algorithm");
+};
+function Nn(i, t, n = !1) {
+  const e = Math.max(i.precision, t.precision);
+  return E(i.center, t.center) && Math.abs(i.majorRadius - t.majorRadius) < e && Math.abs(i.minorRadius - t.minorRadius) < e && (Math.abs(i.tiltAngle - t.tiltAngle) < e || Math.abs(Math.abs(i.tiltAngle - t.tiltAngle) - Math.PI) < e) ? n ? An(i, t) : [] : ye(i, t).filter((o) => i.isOnSegment(o) && t.isOnSegment(o));
+}
+function gs(i, t, n) {
+  if (i instanceof R && t instanceof R) {
+    const e = ge(
+      i,
+      t,
+      !1,
+      n
+    );
+    return e === null ? [] : [e];
+  }
+  if (i instanceof R && t instanceof I)
+    return kt(i, t, n);
+  if (i instanceof I && t instanceof R)
+    return kt(t, i, n);
+  if (i instanceof I && t instanceof I)
+    return me(i, t, !1, n);
+  throw new Error("Not implemented");
+}
+function It(i, t, n) {
+  if (i instanceof R && t instanceof R) {
+    const e = ge(
+      i,
+      t,
+      !0,
+      n
+    );
+    return e === null ? { intersections: [], overlaps: [], count: 0 } : e instanceof R ? { intersections: [], overlaps: [e], count: 1 } : { intersections: [e], overlaps: [], count: 1 };
+  }
+  if (!i.boundingBox.overlaps(t.boundingBox))
+    return { intersections: [], overlaps: [], count: 0 };
+  if (i instanceof R && t instanceof I) {
+    const e = kt(i, t, n);
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (i instanceof I && t instanceof R) {
+    const e = kt(t, i, n);
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (i instanceof I && t instanceof I) {
+    const e = me(
+      i,
+      t,
+      !0,
+      n
+    );
+    return e.length ? e[0] instanceof I ? {
+      intersections: [],
+      overlaps: e,
+      count: e.length
+    } : {
+      intersections: e,
+      overlaps: [],
+      count: e.length
+    } : { intersections: [], overlaps: [], count: 0 };
+  }
+  if (i instanceof R && t instanceof C) {
+    const e = qt(
+      i,
+      t,
+      n
+    );
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (t instanceof R && i instanceof C) {
+    const e = qt(
+      t,
+      i,
+      n
+    );
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (i instanceof I && t instanceof C) {
+    const e = ee(i, t);
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (t instanceof I && i instanceof C) {
+    const e = ee(t, i);
+    return { intersections: e, overlaps: [], count: e.length };
+  }
+  if (i instanceof C && t instanceof C) {
+    const e = Nn(
+      i,
+      t,
+      !0
+    );
+    return e.length ? e[0] instanceof C ? {
+      intersections: [],
+      overlaps: e,
+      count: e.length
+    } : {
+      intersections: e,
+      overlaps: [],
+      count: e.length
+    } : { intersections: [], overlaps: [], count: 0 };
+  }
+  throw new Error("Not implemented");
+}
+function Pe(i) {
+  const t = [];
+  for (let n = 0; n < i; n++)
+    for (let e = 0; e <= n; e++)
+      t.push([n, e]);
+  return t;
+}
+function* Bt(i) {
+  for (const [t, n] of Pe(i.length))
+    t !== n && (yield [i[t], i[n]]);
+}
+class Me extends _t {
+  constructor(t, { ignoreChecks: n = !1 } = {}) {
+    super(), this._boundingBox = null, n || Ee(t), this.segments = t;
+  }
+  get repr() {
+    return this.segments.map((t) => t.repr).join(`
+`) + `
+`;
+  }
+  get info() {
+    return this.repr;
+  }
+  get firstPoint() {
+    return this.segments[0].firstPoint;
+  }
+  get lastPoint() {
+    return this.segments[this.segments.length - 1].lastPoint;
+  }
+  get segmentsCount() {
+    return this.segments.length;
+  }
+  onStroke(t) {
+    return this.segments.some((n) => n.isOnSegment(t));
+  }
+  intersects(t) {
+    return this.boundingBox.overlaps(t.boundingBox) ? this.segments.some(
+      (n) => t.segments.some(
+        (e) => It(n, e).count > 0
+      )
+    ) : !1;
+  }
+  overlappingSegments(t) {
+    return this.segments.flatMap((n) => t.segments.flatMap((e) => n.boundingBox.overlaps(e.boundingBox) ? It(n, e).overlaps : []));
+  }
+  get boundingBox() {
+    if (this._boundingBox === null) {
+      let t = this.segments[0].boundingBox;
+      this.segments.slice(1).forEach((n) => {
+        t = t.merge(n.boundingBox);
+      }), this._boundingBox = t;
+    }
+    return this._boundingBox;
+  }
+  [Symbol.for("nodejs.util.inspect.custom")]() {
+    return this.repr;
+  }
+}
+function Ln(i, t = "Stroke") {
+  Pe(i.length).forEach(
+    ([n, e]) => {
+      if (n === e)
+        return;
+      const s = i[n], r = i[e], o = It(s, r), h = Math.max(s.precision, r.precision);
+      if (o.count !== 0) {
+        if (o.count === 1 && !o.overlaps.length) {
+          const u = n - e, l = o.intersections[0];
+          if (u === 1 && E(s.firstPoint, l, h) || u === -1 && E(s.lastPoint, l, h) || u === i.length - 1 && E(s.lastPoint, l, h) && E(r.firstPoint, l, h) || -u === i.length - 1 && E(s.firstPoint, l, h) && E(r.lastPoint, l, h))
+            return;
+        }
+        if (!(o.count === 2 && i.length === 2 && (E(
+          s.firstPoint,
+          o.intersections[0],
+          h
+        ) && E(
+          s.lastPoint,
+          o.intersections[1],
+          h
+        ) || E(
+          s.firstPoint,
+          o.intersections[1],
+          h
+        ) && E(
+          s.lastPoint,
+          o.intersections[0],
+          h
+        ))))
+          throw new Error(
+            `${t} segments must not intersect, but segments ${s.info} and ${r.info} do at ${JSON.stringify(
+              o.intersections
+            )}`
+          );
+      }
+    }
+  );
+}
+function Ee(i, t = "Stroke") {
+  if (i.length === 0)
+    throw new Error(`${t} must have at least one segment`);
+  nt([i.slice(0, -1), i.slice(1)]).forEach(
+    ([n, e]) => {
+      if (!E(n.lastPoint, e.firstPoint))
+        throw new Error(
+          `${t} segments must be connected, but ${n.info} and ${e.info} are not`
+        );
+    }
+  ), Ln(i, t);
+}
+function ne(i, t) {
+  return !!(i instanceof R && t instanceof R && Vt(i.V, t.V) || i instanceof I && t instanceof I && E(i.center, t.center) && i.radius - t.radius < i.precision);
+}
+function se(i, t) {
+  if (i instanceof R && t instanceof R)
+    return new R(i.firstPoint, t.lastPoint);
+  if (i instanceof I && t instanceof I)
+    return new I(
+      i.firstPoint,
+      t.lastPoint,
+      i.center,
+      i.clockwise
+    );
+  throw new Error("Not implemented");
+}
+function Se(i) {
+  let t = !1;
+  const n = [];
+  for (const e of i.segments) {
+    if (n.length === 0) {
+      n.push(e);
+      continue;
+    }
+    const s = n[n.length - 1];
+    ne(s, e) ? (t = !0, n.pop(), n.push(se(s, e))) : n.push(e);
+  }
+  if (E(i.firstPoint, i.lastPoint) && ne(
+    n[0],
+    n[n.length - 1]
+  )) {
+    t = !0;
+    const e = n.pop();
+    n[0] = se(e, n[0]);
+  }
+  return t ? n : null;
+}
+class G extends Me {
+  constructor() {
+    super(...arguments), this.strokeType = "STRAND";
+  }
+  reverse() {
+    const t = this.segments.map((n) => n.reverse());
+    return t.reverse(), new G(t, { ignoreChecks: !0 });
+  }
+  clone() {
+    return new G(
+      this.segments.map((t) => t.clone()),
+      { ignoreChecks: !0 }
+    );
+  }
+  extend(t) {
+    if (!E(this.lastPoint, t.firstPoint))
+      throw console.error(this.repr, t.repr), new Error("Cannot extend strand: connection point is not the same");
+    return new G([...this.segments, ...t.segments]);
+  }
+  simplify() {
+    const t = Se(this);
+    return t ? new G(t, { ignoreChecks: !0 }) : this;
+  }
+  transform(t) {
+    return new G(
+      this.segments.map((n) => n.transform(t)),
+      { ignoreChecks: !0 }
+    );
+  }
+}
+const Vn = (i, t) => {
+  const n = Ft(t, {
+    V: [1, 0],
+    firstPoint: i,
+    precision: t.precision
+  });
+  if (n === "parallel")
+    return 0;
+  const { intersectionParam1: e, intersectionParam2: s } = n;
+  if (!t.isValidParameter(e) || s <= -t.precision)
+    return 0;
+  if (Math.abs(e) < t.precision || Math.abs(e - 1) < t.precision) {
+    const [, r] = t.midPoint;
+    return i[1] - r < 0 ? 1 : 0;
+  }
+  return 1;
+};
+class be {
+  constructor(t) {
+    this._count = 0, this.segment = t;
+  }
+  update(t, n = !1) {
+    !n && !this.segment.isOnSegment(t) || (E(t, this.segment.firstPoint) ? this._count += this.segment.tangentAtFirstPoint[1] > 0 ? 1 : 0 : E(t, this.segment.lastPoint) ? this._count += this.segment.tangentAtLastPoint[1] > 0 ? 0 : 1 : this._count += 1);
+  }
+  get count() {
+    return this._count;
+  }
+}
+const zn = (i, t) => {
+  const n = t.precision, e = Math.abs(i[1] - t.center[1]);
+  if (e > t.radius + n)
+    return 0;
+  const s = mt(i, t.center), r = t.radius * t.radius, o = n * n;
+  if (Math.abs(s - r) < o && t.isOnSegment(i))
+    return 0;
+  const h = s - r > o;
+  if (h && t.center[0] < i[0])
+    return 0;
+  const u = Math.sqrt(
+    t.radius * t.radius - e * e
+  ), l = new be(t);
+  return l.update([t.center[0] + u, i[1]]), h && l.update([t.center[0] - u, i[1]]), l.count;
+}, Dn = (i, t) => {
+  const n = t.boundingBox.xMax + t.boundingBox.width / 2, e = new R(i, [n, i[1]]), s = new be(t);
+  return qt(e, t).forEach((r) => {
+    s.update(r, !0);
+  }), s.count;
+};
+function Un(i, t) {
+  if (t instanceof R)
+    return Vn(i, t);
+  if (t instanceof I)
+    return zn(i, t);
+  if (t instanceof C)
+    return Dn(i, t);
+  throw new Error("Not implemented");
+}
+class et extends Me {
+  constructor(t, { ignoreChecks: n = !1 } = {}) {
+    super(t, { ignoreChecks: !0 }), this.strokeType = "LOOP", this._clockwise = null, n || Yn(t);
+  }
+  get clockwise() {
+    if (this._clockwise === null) {
+      const t = this.segments.flatMap((e) => e instanceof R ? [e.firstPoint] : [e.firstPoint, e.paramPoint(0.5)]), n = t.map((e, s) => {
+        const r = t[(s + 1) % t.length];
+        return (r[0] - e[0]) * (r[1] + e[1]);
+      }).reduce((e, s) => e + s, 0);
+      this._clockwise = n > 0;
+    }
+    return this._clockwise;
+  }
+  clone() {
+    return new et(
+      this.segments.map((t) => t.clone()),
+      { ignoreChecks: !0 }
+    );
+  }
+  reverse() {
+    const t = this.segments.map((n) => n.reverse());
+    return t.reverse(), new et(t, { ignoreChecks: !0 });
+  }
+  transform(t) {
+    return new et(
+      this.segments.map((n) => n.transform(t)),
+      { ignoreChecks: !0 }
+    );
+  }
+  contains(t) {
+    return this.onStroke(t) || !this.boundingBox.contains(t) ? !1 : this.segments.reduce((e, s) => e + Un(t, s), 0) % 2 === 1;
+  }
+  simplify() {
+    const t = Se(this);
+    return t ? new et(t, { ignoreChecks: !0 }) : this;
+  }
+}
+function Yn(i) {
+  if (Ee(i, "Loop"), !E(i[0].firstPoint, i[i.length - 1].lastPoint))
+    throw new Error("Loop segment must be closed");
+}
+function je(i) {
+  if (i instanceof R)
+    return {
+      type: i.segmentType,
+      firstPoint: i.firstPoint,
+      lastPoint: i.lastPoint
+    };
+  if (i instanceof I)
+    return {
+      type: i.segmentType,
+      firstPoint: i.firstPoint,
+      lastPoint: i.lastPoint,
+      center: i.center,
+      clockwise: i.clockwise
+    };
+  if (i instanceof C)
+    return {
+      type: i.segmentType,
+      firstPoint: i.firstPoint,
+      lastPoint: i.lastPoint,
+      center: i.center,
+      clockwise: i.clockwise,
+      majorRadius: i.majorRadius,
+      minorRadius: i.minorRadius,
+      tiltAngle: i.tiltAngle
+    };
+  throw new Error("Unknown segment type");
+}
+function Ot(i) {
+  return {
+    type: "LOOP",
+    segments: i.segments.map(je)
+  };
+}
+function ve(i) {
+  return {
+    type: "FIGURE",
+    contour: Ot(i.contour),
+    holes: i.holes.map(Ot)
+  };
+}
+function Xn(i) {
+  return {
+    type: "DIAGRAM",
+    figures: i.figures.map(ve)
+  };
+}
+function ie(i) {
+  if (i instanceof tt)
+    return Xn(i);
+  if (i instanceof U)
+    return ve(i);
+  if (i instanceof et)
+    return Ot(i);
+  if (i instanceof I || i instanceof R)
+    return je(i);
+  throw new Error("Unknown shape type");
+}
+class Gn {
+  constructor() {
+    this.ids = [], this.values = [], this.length = 0;
+  }
+  clear() {
+    this.length = 0;
+  }
+  push(t, n) {
+    let e = this.length++;
+    for (; e > 0; ) {
+      const s = e - 1 >> 1, r = this.values[s];
+      if (n >= r)
+        break;
+      this.ids[e] = this.ids[s], this.values[e] = r, e = s;
+    }
+    this.ids[e] = t, this.values[e] = n;
+  }
+  pop() {
+    if (this.length === 0)
+      return;
+    const t = this.ids[0];
+    if (this.length--, this.length > 0) {
+      const n = this.ids[0] = this.ids[this.length], e = this.values[0] = this.values[this.length], s = this.length >> 1;
+      let r = 0;
+      for (; r < s; ) {
+        let o = (r << 1) + 1;
+        const h = o + 1;
+        let u = this.ids[o], l = this.values[o];
+        const c = this.values[h];
+        if (h < this.length && c < l && (o = h, u = this.ids[h], l = c), l >= e)
+          break;
+        this.ids[r] = u, this.values[r] = l, r = o;
+      }
+      this.ids[r] = n, this.values[r] = e;
+    }
+    return t;
+  }
+  peek() {
+    if (this.length !== 0)
+      return this.ids[0];
+  }
+  peekValue() {
+    if (this.length !== 0)
+      return this.values[0];
+  }
+  shrink() {
+    this.ids.length = this.values.length = this.length;
+  }
+}
+const re = [
+  Int8Array,
+  Uint8Array,
+  Uint8ClampedArray,
+  Int16Array,
+  Uint16Array,
+  Int32Array,
+  Uint32Array,
+  Float32Array,
+  Float64Array
+], Ct = 3;
+class Dt {
+  static from(t) {
+    if (!t || t.byteLength === void 0 || t.buffer)
+      throw new Error("Data must be an instance of ArrayBuffer or SharedArrayBuffer.");
+    const [n, e] = new Uint8Array(t, 0, 2);
+    if (n !== 251)
+      throw new Error("Data does not appear to be in a Flatbush format.");
+    if (e >> 4 !== Ct)
+      throw new Error(`Got v${e >> 4} data when expected v${Ct}.`);
+    const [s] = new Uint16Array(t, 2, 1), [r] = new Uint32Array(t, 4, 1);
+    return new Dt(r, s, re[e & 15], void 0, t);
+  }
+  constructor(t, n = 16, e = Float64Array, s = ArrayBuffer, r) {
+    if (t === void 0)
+      throw new Error("Missing required argument: numItems.");
+    if (isNaN(t) || t <= 0)
+      throw new Error(`Unexpected numItems value: ${t}.`);
+    this.numItems = +t, this.nodeSize = Math.min(Math.max(+n, 2), 65535);
+    let o = t, h = o;
+    this._levelBounds = [o * 4];
+    do
+      o = Math.ceil(o / this.nodeSize), h += o, this._levelBounds.push(h * 4);
+    while (o !== 1);
+    this.ArrayType = e || Float64Array, this.IndexArrayType = h < 16384 ? Uint16Array : Uint32Array;
+    const u = re.indexOf(this.ArrayType), l = h * 4 * this.ArrayType.BYTES_PER_ELEMENT;
+    if (u < 0)
+      throw new Error(`Unexpected typed array class: ${e}.`);
+    r && r.byteLength !== void 0 && !r.buffer ? (this.data = r, this._boxes = new this.ArrayType(this.data, 8, h * 4), this._indices = new this.IndexArrayType(this.data, 8 + l, h), this._pos = h * 4, this.minX = this._boxes[this._pos - 4], this.minY = this._boxes[this._pos - 3], this.maxX = this._boxes[this._pos - 2], this.maxY = this._boxes[this._pos - 1]) : (this.data = new s(8 + l + h * this.IndexArrayType.BYTES_PER_ELEMENT), this._boxes = new this.ArrayType(this.data, 8, h * 4), this._indices = new this.IndexArrayType(this.data, 8 + l, h), this._pos = 0, this.minX = 1 / 0, this.minY = 1 / 0, this.maxX = -1 / 0, this.maxY = -1 / 0, new Uint8Array(this.data, 0, 2).set([251, (Ct << 4) + u]), new Uint16Array(this.data, 2, 1)[0] = n, new Uint32Array(this.data, 4, 1)[0] = t), this._queue = new Gn();
+  }
+  add(t, n, e, s) {
+    const r = this._pos >> 2, o = this._boxes;
+    return this._indices[r] = r, o[this._pos++] = t, o[this._pos++] = n, o[this._pos++] = e, o[this._pos++] = s, t < this.minX && (this.minX = t), n < this.minY && (this.minY = n), e > this.maxX && (this.maxX = e), s > this.maxY && (this.maxY = s), r;
+  }
+  finish() {
+    if (this._pos >> 2 !== this.numItems)
+      throw new Error(`Added ${this._pos >> 2} items when expected ${this.numItems}.`);
+    const t = this._boxes;
+    if (this.numItems <= this.nodeSize) {
+      t[this._pos++] = this.minX, t[this._pos++] = this.minY, t[this._pos++] = this.maxX, t[this._pos++] = this.maxY;
+      return;
+    }
+    const n = this.maxX - this.minX || 1, e = this.maxY - this.minY || 1, s = new Uint32Array(this.numItems), r = (1 << 16) - 1;
+    for (let o = 0, h = 0; o < this.numItems; o++) {
+      const u = t[h++], l = t[h++], c = t[h++], a = t[h++], g = Math.floor(r * ((u + c) / 2 - this.minX) / n), w = Math.floor(r * ((l + a) / 2 - this.minY) / e);
+      s[o] = Jn(g, w);
+    }
+    At(s, t, this._indices, 0, this.numItems - 1, this.nodeSize);
+    for (let o = 0, h = 0; o < this._levelBounds.length - 1; o++) {
+      const u = this._levelBounds[o];
+      for (; h < u; ) {
+        const l = h;
+        let c = t[h++], a = t[h++], g = t[h++], w = t[h++];
+        for (let m = 1; m < this.nodeSize && h < u; m++)
+          c = Math.min(c, t[h++]), a = Math.min(a, t[h++]), g = Math.max(g, t[h++]), w = Math.max(w, t[h++]);
+        this._indices[this._pos >> 2] = l, t[this._pos++] = c, t[this._pos++] = a, t[this._pos++] = g, t[this._pos++] = w;
+      }
+    }
+  }
+  search(t, n, e, s, r) {
+    if (this._pos !== this._boxes.length)
+      throw new Error("Data not yet indexed - call index.finish().");
+    let o = this._boxes.length - 4;
+    const h = [], u = [];
+    for (; o !== void 0; ) {
+      const l = Math.min(o + this.nodeSize * 4, he(o, this._levelBounds));
+      for (let c = o; c < l; c += 4) {
+        if (e < this._boxes[c] || s < this._boxes[c + 1] || t > this._boxes[c + 2] || n > this._boxes[c + 3])
+          continue;
+        const a = this._indices[c >> 2] | 0;
+        o >= this.numItems * 4 ? h.push(a) : (r === void 0 || r(a)) && u.push(a);
+      }
+      o = h.pop();
+    }
+    return u;
+  }
+  neighbors(t, n, e = 1 / 0, s = 1 / 0, r) {
+    if (this._pos !== this._boxes.length)
+      throw new Error("Data not yet indexed - call index.finish().");
+    let o = this._boxes.length - 4;
+    const h = this._queue, u = [], l = s * s;
+    for (; o !== void 0; ) {
+      const c = Math.min(o + this.nodeSize * 4, he(o, this._levelBounds));
+      for (let a = o; a < c; a += 4) {
+        const g = this._indices[a >> 2] | 0, w = oe(t, this._boxes[a], this._boxes[a + 2]), m = oe(n, this._boxes[a + 1], this._boxes[a + 3]), d = w * w + m * m;
+        o >= this.numItems * 4 ? h.push(g << 1, d) : (r === void 0 || r(g)) && h.push((g << 1) + 1, d);
+      }
+      for (; h.length && h.peek() & 1; )
+        if (h.peekValue() > l || (u.push(h.pop() >> 1), u.length === e))
+          return h.clear(), u;
+      o = h.pop() >> 1;
+    }
+    return h.clear(), u;
+  }
+}
+function oe(i, t, n) {
+  return i < t ? t - i : i <= n ? 0 : i - n;
+}
+function he(i, t) {
+  let n = 0, e = t.length - 1;
+  for (; n < e; ) {
+    const s = n + e >> 1;
+    t[s] > i ? e = s : n = s + 1;
+  }
+  return t[n];
+}
+function At(i, t, n, e, s, r) {
+  if (Math.floor(e / r) >= Math.floor(s / r))
+    return;
+  const o = i[e + s >> 1];
+  let h = e - 1, u = s + 1;
+  for (; ; ) {
+    do
+      h++;
+    while (i[h] < o);
+    do
+      u--;
+    while (i[u] > o);
+    if (h >= u)
+      break;
+    Wn(i, t, n, h, u);
+  }
+  At(i, t, n, e, u, r), At(i, t, n, u + 1, s, r);
+}
+function Wn(i, t, n, e, s) {
+  const r = i[e];
+  i[e] = i[s], i[s] = r;
+  const o = 4 * e, h = 4 * s, u = t[o], l = t[o + 1], c = t[o + 2], a = t[o + 3];
+  t[o] = t[h], t[o + 1] = t[h + 1], t[o + 2] = t[h + 2], t[o + 3] = t[h + 3], t[h] = u, t[h + 1] = l, t[h + 2] = c, t[h + 3] = a;
+  const g = n[e];
+  n[e] = n[s], n[s] = g;
+}
+function Jn(i, t) {
+  let n = i ^ t, e = 65535 ^ n, s = 65535 ^ (i | t), r = i & (t ^ 65535), o = n | e >> 1, h = n >> 1 ^ n, u = s >> 1 ^ e & r >> 1 ^ s, l = n & s >> 1 ^ r >> 1 ^ r;
+  n = o, e = h, s = u, r = l, o = n & n >> 2 ^ e & e >> 2, h = n & e >> 2 ^ e & (n ^ e) >> 2, u ^= n & s >> 2 ^ e & r >> 2, l ^= e & s >> 2 ^ (n ^ e) & r >> 2, n = o, e = h, s = u, r = l, o = n & n >> 4 ^ e & e >> 4, h = n & e >> 4 ^ e & (n ^ e) >> 4, u ^= n & s >> 4 ^ e & r >> 4, l ^= e & s >> 4 ^ (n ^ e) & r >> 4, n = o, e = h, s = u, r = l, u ^= n & s >> 8 ^ e & r >> 8, l ^= e & s >> 8 ^ (n ^ e) & r >> 8, n = u ^ u >> 1, e = l ^ l >> 1;
+  let c = i ^ t, a = e | 65535 ^ (c | n);
+  return c = (c | c << 8) & 16711935, c = (c | c << 4) & 252645135, c = (c | c << 2) & 858993459, c = (c | c << 1) & 1431655765, a = (a | a << 8) & 16711935, a = (a | a << 4) & 252645135, a = (a | a << 2) & 858993459, a = (a | a << 1) & 1431655765, (a << 1 | c) >>> 0;
+}
+function ke(i, t = 1e-7) {
+  if (i.length === 0)
+    return [];
+  if (i.length === 1)
+    return [i];
+  const n = new Dt(i.length);
+  i.forEach((r) => {
+    const [o, h] = r.firstPoint;
+    n.add(o - t, h - t, o + t, h + t);
+  }), n.finish();
+  const e = [], s = /* @__PURE__ */ new Set();
+  return i.forEach((r, o) => {
+    if (s.has(o))
+      return;
+    const h = [r];
+    let u = o;
+    s.add(o);
+    let l = i.length;
+    for (; ; ) {
+      if (l-- < 0)
+        throw new Error("Infinite loop detected");
+      const c = h[h.length - 1].lastPoint, [a, g] = c, w = n.search(
+        a - t,
+        g - t,
+        a + t,
+        g + t
+      ), m = (p) => Math.abs((u - p) % i.length), d = w.filter((p) => !s.has(p)).map((p) => [
+        i[p],
+        p,
+        m(p)
+      ]).sort(([, , p], [, , y]) => m(p) - m(y));
+      if (d.length === 0) {
+        e.push(h);
+        break;
+      }
+      const [P, f] = d[0];
+      h.push(P), s.add(f), u = f;
+    }
+  }), e;
+}
+class U extends _t {
+  constructor(t, n = [], { ignoreChecks: e = !1 } = {}) {
+    super(), e || Kn(t, n), this.contour = t, this.holes = n;
+  }
+  get boundingBox() {
+    return this.contour.boundingBox;
+  }
+  get isFull() {
+    return this.holes.length === 0;
+  }
+  get allLoops() {
+    return [this.contour, ...this.holes];
+  }
+  clone() {
+    return new U(
+      this.contour.clone(),
+      this.holes.map((t) => t.clone())
+    );
+  }
+  transform(t) {
+    return new U(
+      this.contour.transform(t),
+      this.holes.map((n) => n.transform(t))
+    );
+  }
+  contains(t) {
+    return this.contour.contains(t) && !this.holes.some((n) => n.contains(t));
+  }
+  intersects(t) {
+    return this.allLoops.some(
+      (n) => t.allLoops.some((e) => n.intersects(e))
+    );
+  }
+  overlappingStrands(t) {
+    const n = t instanceof U ? t.allLoops : [t], e = this.allLoops.flatMap((s) => n.flatMap((r) => s.overlappingSegments(r)));
+    return ke(e).map((s) => new G(s));
+  }
+}
+function Kn(i, t = []) {
+  if (!i)
+    throw new Error("Figure must have a contour");
+  for (const [n, e] of Bt([i, ...t]))
+    if (n.intersects(e))
+      throw new Error("Loops in a figure must not intersect");
+  if (t.some(
+    (n) => !i.contains(n.firstPoint) && !i.onStroke(n.firstPoint)
+  ))
+    throw new Error("Holes must be inside the contour");
+  for (const [n, e] of Bt(t))
+    if (n.contains(e.firstPoint))
+      throw console.error(ie(n), ie(e)), new Error("Holes must not be inside other holes");
+}
+const Qn = (i) => {
+  const t = i.map((s, r) => i.slice(r + 1).map((o, h) => [h + r + 1, o]).filter(([, o]) => s.boundingBox.overlaps(o.boundingBox)).map(([o]) => o)), n = [], e = Array(t.length);
+  return t.forEach((s, r) => {
+    let o = e[r];
+    o || (o = [], n.push(o)), o.push(i[r]), s.length && s.forEach((h) => {
+      e[h] = o;
+    });
+  }), n;
+}, Ie = (i) => i.map((t, n) => {
+  const s = t.segments[0].midPoint, r = i.filter((o, h) => n === h ? !1 : o.contains(s));
+  return {
+    loop: t,
+    isIn: r
+  };
+}), Zn = (i, t) => i.flatMap(({ loop: n }) => Ut(
+  t.filter(
+    ({ loop: e, isIn: s }) => e === n || s.indexOf(n) !== -1
+  )
+)), Hn = (i, t) => {
+  const n = t.filter(({ isIn: s }) => s.length <= 1), e = Ut(
+    Ie(i.map(({ loop: s }) => s))
+  );
+  return [n, ...e];
+}, Ut = (i) => {
+  if (!i.length)
+    return [];
+  const t = i.filter(({ isIn: e }) => !e.length), n = i.filter(({ isIn: e }) => e.length > 1);
+  return t.length === 1 && n.length === 0 ? [i] : t.length > 1 ? Zn(t, i) : Hn(n, i);
+};
+function gt(i) {
+  return Qn(i).map(Ie).flatMap(Ut).map((n) => {
+    if (n.length === 1)
+      return new U(n[0].loop);
+    n.sort((r, o) => r.isIn.length - o.isIn.length);
+    const [e, ...s] = n.map(({ loop: r }) => r);
+    return new U(e, s);
+  });
+}
+function ts(i, t) {
+  const n = [];
+  for (const e of i)
+    for (const s of t)
+      n.push([e, s]);
+  return n;
+}
+function* le(i, t, n) {
+  const e = (o) => t.some((h) => E(h, o.lastPoint)), s = (o) => n.some((h) => o.isSame(h));
+  let r = [];
+  for (const o of i)
+    e(o) ? (r.push(o), yield new G(r, { ignoreChecks: !0 }), r = []) : s(o) ? (r.length && (yield new G(r, { ignoreChecks: !0 }), r = []), yield new G([o], { ignoreChecks: !0 })) : r.push(o);
+  r.length && (yield new G(r, { ignoreChecks: !0 }));
+}
+const ue = (i, t) => {
+  const n = i.findIndex((r) => E(t, r.firstPoint)), e = i.slice(0, n);
+  return i.slice(n).concat(e);
+}, ae = (i, t) => {
+  let n = i;
+  const e = (h) => E(h.firstPoint, t.firstPoint) && E(h.lastPoint, t.lastPoint);
+  let s = i.findIndex(e);
+  if (s === -1) {
+    const h = i.map((u) => u.reverse());
+    if (h.reverse(), s = h.findIndex(e), s === -1)
+      throw console.error(
+        h.map((u) => u.repr),
+        t.repr
+      ), new Error("Failed to rotate to segment start");
+    n = h;
+  }
+  const r = n.slice(0, s);
+  return n.slice(s).concat(r);
+};
+function es(i, t, n) {
+  return i.filter((e) => {
+    const s = t.filter((h) => E(h.firstPoint, e) || E(h.lastPoint, e));
+    if (s.length % 2)
+      throw new Error("Bug in the intersection algo on non crossing point");
+    const r = s.map((h) => n.contains(h.midPoint));
+    return !(r.every((h) => h) || !r.some((h) => h));
+  });
+}
+function ns(i, t, n) {
+  let e = [];
+  const s = [], r = new Array(i.segments.length).fill(0).map(() => []), o = new Array(t.segments.length).fill(0).map(() => []);
+  if (i.segments.forEach((g, w) => {
+    t.segments.forEach((m, d) => {
+      const { intersections: P, overlaps: f } = It(
+        g,
+        m,
+        n
+      );
+      e.push(...P), r[w].push(...P), o[d].push(...P), s.push(...f);
+      const p = f.flatMap((y) => [
+        y.firstPoint,
+        y.lastPoint
+      ]);
+      e.push(...p), r[w].push(...p), o[d].push(...p);
+    });
+  }), e = ct(e, n), !e.length || e.length === 1)
+    return null;
+  const h = ([g, w]) => w.length ? g.splitAt(w) : [g];
+  let u = nt([i.segments, r]).flatMap(h), l = nt([t.segments, o]).flatMap(h);
+  if (e = es(
+    e,
+    u,
+    t
+  ), !e.length && !s.length)
+    return null;
+  if (s.length) {
+    const g = s[0];
+    u = ae(
+      u,
+      g
+    ), l = ae(
+      l,
+      g
+    );
+  } else {
+    const g = e[0];
+    u = ue(u, g), l = ue(l, g);
+  }
+  let c = Array.from(
+    le(
+      u,
+      e,
+      s
+    )
+  ), a = Array.from(
+    le(
+      l,
+      e,
+      s
+    )
+  );
+  return (!E(
+    a[0].lastPoint,
+    c[0].lastPoint
+  ) || s.length > 0 && a[0].segmentsCount !== 1) && (a = a.map((g) => g.reverse()).reverse(), E(a[0].lastPoint, c[0].lastPoint) || (c = c.map((g) => g.reverse()).reverse())), nt([c, a]).map(([g, w]) => g.segmentsCount === 1 && s.some((m) => g.segments[0].isSame(m)) ? [g, "same"] : [g, w]);
+}
+function ce(i) {
+  let t = i[0];
+  for (const n of i.slice(1))
+    t = t.extend(n);
+  if (!E(t.firstPoint, t.lastPoint))
+    throw console.error(
+      T(t.firstPoint),
+      T(t.lastPoint)
+    ), new Error("Bug in the intersection algo on non closing strand");
+  return new et(t.segments);
+}
+function ss(i, t) {
+  const n = nt([
+    t.slice(0, -1),
+    t.slice(1)
+  ]).map(([s, r]) => ce(i.slice(s, r)));
+  let e = i.slice(
+    t[t.length - 1]
+  );
+  return t[0] !== 0 && (e = e.concat(i.slice(0, t[0]))), n.push(ce(e)), n;
+}
+function is(i) {
+  if (!i.length)
+    return [];
+  const t = i.map((s) => s.firstPoint);
+  let n = i.map((s) => s.lastPoint);
+  n = n.slice(-1).concat(n.slice(0, -1));
+  const e = nt([t, n]).flatMap(
+    ([s, r], o) => E(s, r) ? [] : o
+  );
+  try {
+    return ss(i, e);
+  } catch {
+    return ke(i.flatMap((r) => r.segments)).filter((r) => r.length > 1).filter((r) => E(r[0].firstPoint, r.at(-1).lastPoint)).map((r) => new et(r));
+  }
+}
+const fe = (i, t) => {
+  if (i.length === 0)
+    return [t];
+  const n = i.at(-1);
+  return E(n.lastPoint, t.firstPoint) ? i.slice(0, -1).concat([n.extend(t)]) : E(n.lastPoint, t.lastPoint) ? i.slice(0, -1).concat([n.extend(t.reverse())]) : i.concat([t]);
+}, rs = (i, t) => i.length === 0 ? [t] : E(i[0].firstPoint, t.lastPoint) ? [t.extend(i[0])].concat(i.slice(1)) : [t].concat(i);
+function Yt(i, t, {
+  firstInside: n,
+  secondInside: e
+}) {
+  const s = ns(i, t);
+  if (!s) {
+    const u = i.segments[0].midPoint, l = t.contains(u), c = t.segments[0].midPoint, a = i.contains(c);
+    return {
+      identical: !1,
+      firstCurveInSecond: l,
+      secondCurveInFirst: a
+    };
+  }
+  if (s.every(([, u]) => u === "same"))
+    return { identical: !0 };
+  let r = null, o = null;
+  const h = s.flatMap(([u, l]) => {
+    let c = [], a = 0;
+    if (l === "same")
+      return o === 1 ? (o = 1, u) : o === 2 || o === 0 ? (o = null, []) : o === null ? (r ? r = r.extend(u) : r = u, []) : (console.error("weird situation"), []);
+    const g = u.segments[0].midPoint, w = t.contains(g);
+    (n === "keep" && w || n === "remove" && !w) && (a += 1, c = fe(c, u));
+    const m = l.segments[0].midPoint, d = i.contains(m);
+    if (e === "keep" && d || e === "remove" && !d) {
+      const P = l;
+      a += 1, a === 2 && c.length ? (c = fe(c, P), r = null) : c = [P];
+    }
+    return o === null && a === 1 && r && (c = rs(c, r)), a === 1 && (o = a, r = null), c.length ? c : (r = null, []);
+  });
+  return is(h);
+}
+const os = (i, t) => {
+  const n = Yt(i, t, {
+    firstInside: "remove",
+    secondInside: "remove"
+  });
+  return Array.isArray(n) ? n : n.identical ? [i] : n.firstCurveInSecond ? [t] : n.secondCurveInFirst ? [i] : [i, t];
+}, xt = (i, t) => {
+  const n = Yt(i, t, {
+    firstInside: "remove",
+    secondInside: "keep"
+  });
+  return Array.isArray(n) ? n : n.identical ? [] : n.firstCurveInSecond ? [] : n.secondCurveInFirst ? [i, t] : [i];
+}, Xt = (i, t) => {
+  const n = Yt(i, t, {
+    firstInside: "keep",
+    secondInside: "keep"
+  });
+  return Array.isArray(n) ? n : n.identical ? [i] : n.firstCurveInSecond ? [i] : n.secondCurveInFirst ? [t] : [];
+};
+function hs(i) {
+  const t = /* @__PURE__ */ new Map(), n = [];
+  return i.forEach((e, s) => {
+    let r;
+    t.has(s) ? r = t.get(s) : (r = { current: [e], fusedWith: /* @__PURE__ */ new Set([s]) }, n.push(r)), i.slice(s + 1).forEach((o, h) => {
+      const u = r.current, l = s + h + 1;
+      if (r.fusedWith.has(l))
+        return;
+      let c = [o], a = !1;
+      if (t.has(l) && (c = t.get(l).current, a = !0), !u.some(
+        (m) => c.some((d) => m.intersects(d))
+      ))
+        return;
+      let w;
+      u.length > 1 || c.length > 1 ? w = pt(u, c) : w = xe(u[0], c[0]), r.fusedWith.add(l), r.current = w, a || t.set(l, r);
+    });
+  }), n.flatMap(({ current: e }) => e);
+}
+function xe(i, t) {
+  const n = os(i.contour, t.contour), e = t.holes.flatMap((o) => xt(o, i.contour)), s = i.holes.flatMap((o) => xt(o, t.contour)), r = ts(i.holes, t.holes).flatMap(
+    ([o, h]) => Xt(o, h)
+  );
+  return gt([
+    ...n,
+    ...e,
+    ...s,
+    ...r
+  ]);
+}
+function Rt(i, t) {
+  if (i.isFull && t.isFull)
+    return gt(xt(i.contour, t.contour));
+  if (i.isFull) {
+    const e = xt(i.contour, t.contour), s = t.holes.flatMap(
+      (r) => Xt(r, i.contour)
+    );
+    return gt([...e, ...s]);
+  } else if (t.isFull && !i.contour.intersects(t.contour))
+    if (i.contour.contains(t.contour.firstPoint)) {
+      const e = pt(
+        i.holes.map((s) => new U(s)),
+        [t]
+      );
+      return gt([
+        i.contour,
+        ...e.flatMap((s) => s.allLoops)
+      ]);
+    } else
+      return [i];
+  let n = Rt(new U(i.contour), t);
+  return i.holes.forEach((e) => {
+    n = n.flatMap((s) => Rt(s, new U(e)));
+  }), n;
+}
+function ls(i, t) {
+  const n = Xt(i.contour, t.contour);
+  if (!n.length)
+    return [];
+  let e = gt(n);
+  return e = wt(
+    e,
+    i.holes.map((s) => new U(s))
+  ), wt(
+    e,
+    t.holes.map((s) => new U(s))
+  );
+}
+function pt(i, t) {
+  if (!i.length)
+    return t;
+  if (!t.length)
+    return i;
+  if (i.length === 1 && t.length > 1 || t.length === 1 && i.length > 1)
+    return hs([...i, ...t]);
+  if (i.length > 1 && t.length > 1) {
+    let n = pt([i[0]], t);
+    return i.slice(1).forEach((e) => {
+      n = pt([e], n);
+    }), n;
+  }
+  return i.length === 1 && t.length === 1 ? xe(i[0], t[0]) : [];
+}
+function wt(i, t) {
+  if (!i.length)
+    return [];
+  if (!t.length)
+    return i;
+  if (i.length === 1 && t.length === 1)
+    return Rt(i[0], t[0]);
+  if (i.length > 1)
+    return i.flatMap((e) => wt([e], t));
+  let n = Rt(i[0], t[0]);
+  return t.slice(1).forEach((e) => {
+    n = wt(n, [e]);
+  }), n;
+}
+function Nt(i, t) {
+  return !i.length || !t.length ? [] : i.length === 1 && t.length === 1 ? ls(i[0], t[0]) : i.length > 1 ? i.flatMap((n) => Nt([n], t)) : t.flatMap((n) => Nt(i, [n]));
+}
+class tt extends _t {
+  constructor(t = [], { ignoreChecks: n = !1 } = {}) {
+    super(), this._boundingBox = null, n || us(t), this.figures = t;
+  }
+  get isEmpty() {
+    return this.figures.length === 0;
+  }
+  get boundingBox() {
+    if (this.isEmpty)
+      return new it();
+    if (this._boundingBox === null) {
+      let t = this.figures[0].boundingBox;
+      for (const n of this.figures.slice(1))
+        t = t.merge(n.boundingBox);
+      this._boundingBox = t;
+    }
+    return this._boundingBox;
+  }
+  clone() {
+    return new tt(this.figures.map((t) => t.clone()));
+  }
+  transform(t) {
+    return new tt(this.figures.map((n) => n.transform(t)));
+  }
+  contains(t) {
+    return this.figures.some((n) => n.contains(t));
+  }
+  intersects(t) {
+    return this.figures.some(
+      (n) => t.figures.some((e) => n.intersects(e))
+    );
+  }
+  overlappingStrands(t) {
+    return this.figures.flatMap((n) => t instanceof tt ? t.figures.flatMap(
+      (e) => n.overlappingStrands(e)
+    ) : n.overlappingStrands(t));
+  }
+  fuse(t) {
+    return new tt(pt(this.figures, t.figures));
+  }
+  cut(t) {
+    return new tt(wt(this.figures, t.figures));
+  }
+  intersect(t) {
+    return new tt(Nt(this.figures, t.figures));
+  }
+}
+function us(i) {
+  for (const [t, n] of Bt(i))
+    if (t.intersects(n))
+      throw new Error("Diagram figures must not intersect");
+}
+export {
+  I as A,
+  it as B,
+  gs as C,
+  tt as D,
+  C as E,
+  U as F,
+  bt as G,
+  as as H,
+  fs as I,
+  Xe as J,
+  ft as K,
+  et as L,
+  Lt as M,
+  Ue as R,
+  G as S,
+  V as T,
+  pt as a,
+  kt as b,
+  wt as c,
+  q as d,
+  A as e,
+  It as f,
+  O as g,
+  D as h,
+  Nt as i,
+  me as j,
+  Ye as k,
+  Ft as l,
+  mt as m,
+  W as n,
+  R as o,
+  Qe as p,
+  ke as q,
+  ct as r,
+  le as s,
+  E as t,
+  cs as u,
+  at as v,
+  ie as w,
+  _t as x,
+  J as y,
+  nt as z
+};
+//# sourceMappingURL=Diagram-d848c815.js.map