rf-touchstone 0.0.0 → 0.0.1

package/dist/Touchstone.es.js

@@ -159,13 +159,13 @@ function Wu(e) {
 function Vu(e) {
   return e instanceof Date;
 }
-function Ju(e) {
+function Gu(e) {
   return e instanceof RegExp;
 }
 function Kt(e) {
   return !!(e && typeof e == "object" && e.constructor === Object && !ii(e) && !ui(e));
 }
-function Gu(e) {
+function Ju(e) {
   return e ? e instanceof Map || e instanceof $u || typeof e.set == "function" && typeof e.get == "function" && typeof e.keys == "function" && typeof e.has == "function" : !1;
 }
 function Yu(e) {
@@ -781,7 +781,7 @@ function si() {
   function Y(F, _) {
     let N = _;
     if (F.some(($) => $.hasConversion)) {
-      const $ = y(F), I = F.map(G);
+      const $ = y(F), I = F.map(J);
       N = function() {
         const U = [], Z = $ ? arguments.length - 1 : arguments.length;
         for (let P = 0; P < Z; P++)
@@ -798,7 +798,7 @@ function si() {
     }
     return T;
   }
-  function G(F) {
+  function J(F) {
     let _, N, T, $;
     const I = [], H = [];
     switch (F.types.forEach(function(U) {
@@ -874,7 +874,7 @@ function si() {
     return I ? H ? T === $ : $ >= T : H ? T >= $ : T === $;
   }
   function ie(F) {
-    return F.map((_) => sr(_) ? Ge(_.referToSelf.callback) : ar(_) ? Me(_.referTo.references, _.referTo.callback) : _);
+    return F.map((_) => sr(_) ? Je(_.referToSelf.callback) : ar(_) ? Me(_.referTo.references, _.referTo.callback) : _);
   }
   function fe(F, _, N) {
     const T = [];
@@ -992,7 +992,7 @@ function si() {
   function ae(F, _, N) {
     return Array.prototype.slice.call(F, _, N);
   }
-  function Je(F, _) {
+  function Ge(F, _) {
     for (let N = 0; N < F.length; N++)
       if (_(F[N]))
         return F[N];
@@ -1014,7 +1014,7 @@ function si() {
       }
     };
   }
-  function Ge(F) {
+  function Je(F) {
     if (typeof F != "function")
       throw new TypeError("Callback function expected as first argument");
     return {
@@ -1080,7 +1080,7 @@ function si() {
       _ || (T = Ar(T, Z)), zt($, U);
     }
     return Ae(T || "", $);
-  }, o.create = si, o.createCount = ru.createCount, o.onMismatch = we, o.throwMismatchError = we, o.createError = R, o.clear = l, o.clearConversions = f, o.addTypes = c, o._findType = a, o.referTo = le, o.referToSelf = Ge, o.convert = s, o.findSignature = p, o.find = m, o.isTypedFunction = D, o.warnAgainstDeprecatedThis = !0, o.addType = function(F, _) {
+  }, o.create = si, o.createCount = ru.createCount, o.onMismatch = we, o.throwMismatchError = we, o.createError = R, o.clear = l, o.clearConversions = f, o.addTypes = c, o._findType = a, o.referTo = le, o.referToSelf = Je, o.convert = s, o.findSignature = p, o.find = m, o.isTypedFunction = D, o.warnAgainstDeprecatedThis = !0, o.addType = function(F, _) {
     let N = "any";
     _ !== !1 && n.has("Object") && (N = "Object"), o.addTypes([F], N);
   };
@@ -1114,7 +1114,7 @@ function si() {
     F.forEach((N) => o.addConversion(N, _));
   }, o.removeConversion = function(F) {
     on(F);
-    const _ = a(F.to), N = Je(_.conversionsTo, ($) => $.from === F.from);
+    const _ = a(F.to), N = Ge(_.conversionsTo, ($) => $.from === F.from);
     if (!N)
       throw new Error("Attempt to remove nonexistent conversion from " + F.from + " to " + F.to);
     if (N.convert !== F.convert)
@@ -1132,7 +1132,7 @@ function si() {
   }, o;
 }
 const ot = si();
-function J(e, r, t, n) {
+function G(e, r, t, n) {
   function i(u) {
     var o = mo(u, r.map(Fo));
     return wo(e, r, u), t(o);
@@ -1337,7 +1337,7 @@ function cn(e, r) {
 }
 var hi = function() {
   return hi = ot.create, ot;
-}, _o = ["?BigNumber", "?Complex", "?DenseMatrix", "?Fraction"], Bo = /* @__PURE__ */ J("typed", _o, function(r) {
+}, _o = ["?BigNumber", "?Complex", "?DenseMatrix", "?Fraction"], Bo = /* @__PURE__ */ G("typed", _o, function(r) {
   var {
     BigNumber: t,
     Complex: n,
@@ -1432,7 +1432,7 @@ var hi = function() {
     },
     {
       name: "RegExp",
-      test: Ju
+      test: Gu
     },
     {
       name: "null",
@@ -1508,7 +1508,7 @@ var hi = function() {
     },
     {
       name: "Map",
-      test: Gu
+      test: Ju
     },
     {
       name: "Object",
@@ -1773,7 +1773,7 @@ var Sr = 9e15, vr = 1e9, Wt = "0123456789abcdef", at = "2.3025850929940456840179
   // Whether to use cryptographically-secure random number generation, if available.
   crypto: !1
   // true/false
-}, pi, cr, j = !0, Et = "[DecimalError] ", dr = Et + "Invalid argument: ", Di = Et + "Precision limit exceeded", mi = Et + "crypto unavailable", di = "[object Decimal]", $e = Math.floor, Ne = Math.pow, Mo = /^0b([01]+(\.[01]*)?|\.[01]+)(p[+-]?\d+)?$/i, So = /^0x([0-9a-f]+(\.[0-9a-f]*)?|\.[0-9a-f]+)(p[+-]?\d+)?$/i, xo = /^0o([0-7]+(\.[0-7]*)?|\.[0-7]+)(p[+-]?\d+)?$/i, gi = /^(\d+(\.\d*)?|\.\d+)(e[+-]?\d+)?$/i, rr = 1e7, K = 7, To = 9007199254740991, Io = at.length - 1, Jt = st.length - 1, z = { toStringTag: di };
+}, pi, cr, j = !0, Et = "[DecimalError] ", dr = Et + "Invalid argument: ", Di = Et + "Precision limit exceeded", mi = Et + "crypto unavailable", di = "[object Decimal]", $e = Math.floor, Ne = Math.pow, Mo = /^0b([01]+(\.[01]*)?|\.[01]+)(p[+-]?\d+)?$/i, So = /^0x([0-9a-f]+(\.[0-9a-f]*)?|\.[0-9a-f]+)(p[+-]?\d+)?$/i, xo = /^0o([0-7]+(\.[0-7]*)?|\.[0-7]+)(p[+-]?\d+)?$/i, gi = /^(\d+(\.\d*)?|\.\d+)(e[+-]?\d+)?$/i, rr = 1e7, K = 7, To = 9007199254740991, Io = at.length - 1, Gt = st.length - 1, z = { toStringTag: di };
 z.absoluteValue = z.abs = function() {
   var e = new this.constructor(this);
   return e.s < 0 && (e.s = 1), X(e);
@@ -1897,11 +1897,11 @@ z.inverseTangent = z.atan = function() {
   if (l.isFinite()) {
     if (l.isZero())
       return new f(l);
-    if (l.abs().eq(1) && h + 4 <= Jt)
+    if (l.abs().eq(1) && h + 4 <= Gt)
       return o = nr(f, h + 4, D).times(0.25), o.s = l.s, o;
   } else {
     if (!l.s) return new f(NaN);
-    if (h + 4 <= Jt)
+    if (h + 4 <= Gt)
       return o = nr(f, h + 4, D).times(0.5), o.s = l.s, o;
   }
   for (f.precision = a = h + 10, f.rounding = 1, t = Math.min(28, a / K + 2 | 0), e = t; e; --e) l = l.div(l.times(l).plus(1).sqrt().plus(1));
@@ -1999,7 +1999,7 @@ z.modulo = z.mod = function(e) {
   return e = new n(e), !t.d || !e.s || e.d && !e.d[0] ? new n(NaN) : !e.d || t.d && !t.d[0] ? X(new n(t), n.precision, n.rounding) : (j = !1, n.modulo == 9 ? (r = he(t, e.abs(), 0, 3, 1), r.s *= e.s) : r = he(t, e, 0, n.modulo, 1), r = r.times(e), j = !0, t.minus(r));
 };
 z.naturalExponential = z.exp = function() {
-  return Gt(this);
+  return Jt(this);
 };
 z.naturalLogarithm = z.ln = function() {
   return Dr(this);
@@ -2134,7 +2134,7 @@ z.toPower = z.pow = function(e) {
     if ((e.d[r] & 1) == 0 && (o = 1), a.e == 0 && a.d[0] == 1 && a.d.length == 1)
       return a.s = o, a;
   }
-  return t = Ne(+a, l), r = t == 0 || !isFinite(t) ? $e(l * (Math.log("0." + Ie(a.d)) / Math.LN10 + a.e + 1)) : new c(t + "").e, r > c.maxE + 1 || r < c.minE - 1 ? new c(r > 0 ? o / 0 : 0) : (j = !1, c.rounding = a.s = 1, t = Math.min(12, (r + "").length), i = Gt(e.times(Dr(a, n + t)), n), i.d && (i = X(i, n + 5, 1), Lr(i.d, n, u) && (r = n + 10, i = X(Gt(e.times(Dr(a, r + t)), r), r + 5, 1), +Ie(i.d).slice(n + 1, n + 15) + 1 == 1e14 && (i = X(i, n + 1, 0)))), i.s = o, j = !0, c.rounding = u, X(i, n, u));
+  return t = Ne(+a, l), r = t == 0 || !isFinite(t) ? $e(l * (Math.log("0." + Ie(a.d)) / Math.LN10 + a.e + 1)) : new c(t + "").e, r > c.maxE + 1 || r < c.minE - 1 ? new c(r > 0 ? o / 0 : 0) : (j = !1, c.rounding = a.s = 1, t = Math.min(12, (r + "").length), i = Jt(e.times(Dr(a, n + t)), n), i.d && (i = X(i, n + 5, 1), Lr(i.d, n, u) && (r = n + 10, i = X(Jt(e.times(Dr(a, r + t)), r), r + 5, 1), +Ie(i.d).slice(n + 1, n + 15) + 1 == 1e14 && (i = X(i, n + 1, 0)))), i.s = o, j = !0, c.rounding = u, X(i, n, u));
 };
 z.toPrecision = function(e, r) {
   var t, n = this, i = n.constructor;
@@ -2305,7 +2305,7 @@ function ft(e, r, t) {
   return X(new e(at), r, 1, !0);
 }
 function nr(e, r, t) {
-  if (r > Jt) throw Error(Di);
+  if (r > Gt) throw Error(Di);
   return X(new e(st), r, t, !0);
 }
 function vi(e) {
@@ -2344,7 +2344,7 @@ function yi(e, r, t) {
   }
   return u;
 }
-function Gt(e, r) {
+function Jt(e, r) {
   var t, n, i, u, o, a, c, l = 0, f = 0, h = 0, D = e.constructor, p = D.rounding, m = D.precision;
   if (!e.d || !e.d[0] || e.e > 17)
     return new D(e.d ? e.d[0] ? e.s < 0 ? 0 : 1 / 0 : 1 : e.s ? e.s < 0 ? 0 : e : NaN);
@@ -2525,10 +2525,10 @@ function Vo(e, r) {
   var t, n = this.precision, i = this.rounding, u = n + 4;
   return !e.s || !r.s ? t = new this(NaN) : !e.d && !r.d ? (t = nr(this, u, 1).times(r.s > 0 ? 0.25 : 0.75), t.s = e.s) : !r.d || e.isZero() ? (t = r.s < 0 ? nr(this, n, i) : new this(0), t.s = e.s) : !e.d || r.isZero() ? (t = nr(this, u, 1).times(0.5), t.s = e.s) : r.s < 0 ? (this.precision = u, this.rounding = 1, t = this.atan(he(e, r, u, 1)), r = nr(this, u, 1), this.precision = n, this.rounding = i, t = e.s < 0 ? t.minus(r) : t.plus(r)) : t = this.atan(he(e, r, u, 1)), t;
 }
-function Jo(e) {
+function Go(e) {
   return new this(e).cbrt();
 }
-function Go(e) {
+function Jo(e) {
   return X(e = new this(e), e.e + 1, 2);
 }
 function Yo(e, r, t) {
@@ -2613,7 +2613,7 @@ function Ai(e) {
       return u < 0 ? (u = -u, l.s = -1) : l.s = 1, tt(l, u.toString());
     throw Error(dr + u);
   }
-  if (i.prototype = z, i.ROUND_UP = 0, i.ROUND_DOWN = 1, i.ROUND_CEIL = 2, i.ROUND_FLOOR = 3, i.ROUND_HALF_UP = 4, i.ROUND_HALF_DOWN = 5, i.ROUND_HALF_EVEN = 6, i.ROUND_HALF_CEIL = 7, i.ROUND_HALF_FLOOR = 8, i.EUCLID = 9, i.config = i.set = ko, i.clone = Ai, i.isDecimal = pn, i.abs = Po, i.acos = Uo, i.acosh = qo, i.add = Ro, i.asin = Lo, i.asinh = Ho, i.atan = Zo, i.atanh = Wo, i.atan2 = Vo, i.cbrt = Jo, i.ceil = Go, i.clamp = Yo, i.cos = Xo, i.cosh = Qo, i.div = Ko, i.exp = jo, i.floor = ea, i.hypot = ra, i.ln = ta, i.log = na, i.log10 = ua, i.log2 = ia, i.max = oa, i.min = aa, i.mod = sa, i.mul = fa, i.pow = ca, i.random = la, i.round = ha, i.sign = pa, i.sin = Da, i.sinh = ma, i.sqrt = da, i.sub = ga, i.sum = va, i.tan = wa, i.tanh = ya, i.trunc = Fa, e === void 0 && (e = {}), e && e.defaults !== !0)
+  if (i.prototype = z, i.ROUND_UP = 0, i.ROUND_DOWN = 1, i.ROUND_CEIL = 2, i.ROUND_FLOOR = 3, i.ROUND_HALF_UP = 4, i.ROUND_HALF_DOWN = 5, i.ROUND_HALF_EVEN = 6, i.ROUND_HALF_CEIL = 7, i.ROUND_HALF_FLOOR = 8, i.EUCLID = 9, i.config = i.set = ko, i.clone = Ai, i.isDecimal = pn, i.abs = Po, i.acos = Uo, i.acosh = qo, i.add = Ro, i.asin = Lo, i.asinh = Ho, i.atan = Zo, i.atanh = Wo, i.atan2 = Vo, i.cbrt = Go, i.ceil = Jo, i.clamp = Yo, i.cos = Xo, i.cosh = Qo, i.div = Ko, i.exp = jo, i.floor = ea, i.hypot = ra, i.ln = ta, i.log = na, i.log10 = ua, i.log2 = ia, i.max = oa, i.min = aa, i.mod = sa, i.mul = fa, i.pow = ca, i.random = la, i.round = ha, i.sign = pa, i.sin = Da, i.sinh = ma, i.sqrt = da, i.sub = ga, i.sum = va, i.tan = wa, i.tanh = ya, i.trunc = Fa, e === void 0 && (e = {}), e && e.defaults !== !0)
     for (n = ["precision", "rounding", "toExpNeg", "toExpPos", "maxE", "minE", "modulo", "crypto"], r = 0; r < n.length; ) e.hasOwnProperty(t = n[r++]) || (e[t] = this[t]);
   return i.config(e), i;
 }
@@ -2728,7 +2728,7 @@ z[Symbol.toStringTag] = "Decimal";
 var zr = z.constructor = Ai(Vt);
 at = new zr(at);
 st = new zr(st);
-var Ea = "BigNumber", Aa = ["?on", "config"], ba = /* @__PURE__ */ J(Ea, Aa, (e) => {
+var Ea = "BigNumber", Aa = ["?on", "config"], ba = /* @__PURE__ */ G(Ea, Aa, (e) => {
   var {
     on: r,
     config: t
@@ -3495,7 +3495,7 @@ x.E = new x(Math.E, 0);
 x.INFINITY = new x(1 / 0, 1 / 0);
 x.NAN = new x(NaN, NaN);
 x.EPSILON = 1e-15;
-var Ba = "Complex", Na = [], Ma = /* @__PURE__ */ J(Ba, Na, () => (Object.defineProperty(x, "name", {
+var Ba = "Complex", Na = [], Ma = /* @__PURE__ */ G(Ba, Na, () => (Object.defineProperty(x, "name", {
   value: "Complex"
 }), x.prototype.constructor = x, x.prototype.type = "Complex", x.prototype.isComplex = !0, x.prototype.toJSON = function() {
   return {
@@ -4068,7 +4068,7 @@ Xe.prototype = {
     return this;
   }
 };
-var za = "Fraction", Oa = [], $a = /* @__PURE__ */ J(za, Oa, () => (Object.defineProperty(Xe, "name", {
+var za = "Fraction", Oa = [], $a = /* @__PURE__ */ G(za, Oa, () => (Object.defineProperty(Xe, "name", {
   value: "Fraction"
 }), Xe.prototype.constructor = Xe, Xe.prototype.type = "Fraction", Xe.prototype.isFraction = !0, Xe.prototype.toJSON = function() {
   return {
@@ -4080,7 +4080,7 @@ var za = "Fraction", Oa = [], $a = /* @__PURE__ */ J(za, Oa, () => (Object.defin
   return new Xe(e);
 }, Xe), {
   isClass: !0
-}), Pa = "Matrix", Ua = [], qa = /* @__PURE__ */ J(Pa, Ua, () => {
+}), Pa = "Matrix", Ua = [], qa = /* @__PURE__ */ G(Pa, Ua, () => {
   function e() {
     if (!(this instanceof e))
       throw new SyntaxError("Constructor must be called with the new operator");
@@ -4360,7 +4360,7 @@ function kt(e, r, t, n) {
   }
 }
 function _i(e, r) {
-  var t = Ja(e, !0), n = t.length;
+  var t = Ga(e, !0), n = t.length;
   if (!Array.isArray(e) || !Array.isArray(r))
     throw new TypeError("Array expected");
   if (r.length === 0)
@@ -4420,7 +4420,7 @@ function Mi(e, r, t) {
       e = [e];
   return e;
 }
-function Ja(e) {
+function Ga(e) {
   var r = arguments.length > 1 && arguments[1] !== void 0 ? arguments[1] : !1;
   if (!Array.isArray(e))
     return e;
@@ -4503,10 +4503,10 @@ function Xt(e, r) {
   var n = Ti(t, r), i = n.length, u = [...Array(i - t.length).fill(1), ...t], o = Ya(e);
   t.length < i && (o = _i(o, u), t = ze(o));
   for (var a = 0; a < i; a++)
-    t[a] < n[a] && (o = Ga(o, n[a], a), t = ze(o));
+    t[a] < n[a] && (o = Ja(o, n[a], a), t = ze(o));
   return o;
 }
-function Ga(e, r, t) {
+function Ja(e, r, t) {
   return xi(...Array(r).fill(e), t);
 }
 function zi(e, r) {
@@ -4624,7 +4624,7 @@ function Ka(e, r, t, n) {
   } else
     throw new TypeError("Function ".concat(t, " cannot apply callback arguments ") + "to function ".concat(n, ": ").concat(e.message));
 }
-var ja = "DenseMatrix", es = ["Matrix"], rs = /* @__PURE__ */ J(ja, es, (e) => {
+var ja = "DenseMatrix", es = ["Matrix"], rs = /* @__PURE__ */ G(ja, es, (e) => {
   var {
     Matrix: r
   } = e;
@@ -4993,7 +4993,7 @@ function Ve(e, r, t) {
   var n = (i) => i === 0 ? i : r(i);
   return ve(e) ? e.map((i) => n(i), !1, !0) : yn(e, n, !0);
 }
-var En = "isInteger", ts = ["typed"], ns = /* @__PURE__ */ J(En, ts, (e) => {
+var En = "isInteger", ts = ["typed"], ns = /* @__PURE__ */ G(En, ts, (e) => {
   var {
     typed: r
   } = e;
@@ -5059,7 +5059,7 @@ function gr(e, r) {
     throw new Error("Absolute tolerance must be at least 0");
   return e.isNaN() || r.isNaN() ? !1 : !e.isFinite() || !r.isFinite() ? e.eq(r) : e.eq(r) ? !0 : e.minus(r).abs().lte(e.constructor.max(e.constructor.max(e.abs(), r.abs()).mul(t), n));
 }
-var An = "isPositive", os = ["typed", "config"], as = /* @__PURE__ */ J(An, os, (e) => {
+var An = "isPositive", os = ["typed", "config"], as = /* @__PURE__ */ G(An, os, (e) => {
   var {
     typed: r,
     config: t
@@ -5072,7 +5072,7 @@ var An = "isPositive", os = ["typed", "config"], as = /* @__PURE__ */ J(An, os,
     Unit: r.referToSelf((n) => (i) => r.find(n, i.valueType())(i.value)),
     "Array | Matrix": r.referToSelf((n) => (i) => Ve(i, n))
   });
-}), bn = "isZero", ss = ["typed", "equalScalar"], fs = /* @__PURE__ */ J(bn, ss, (e) => {
+}), bn = "isZero", ss = ["typed", "equalScalar"], fs = /* @__PURE__ */ G(bn, ss, (e) => {
   var {
     typed: r,
     equalScalar: t
@@ -5087,7 +5087,7 @@ var An = "isPositive", os = ["typed", "config"], as = /* @__PURE__ */ J(An, os,
 function cs(e, r, t, n) {
   return tr(e.re, r.re, t, n) && tr(e.im, r.im, t, n);
 }
-var Wr = /* @__PURE__ */ J("compareUnits", ["typed"], (e) => {
+var Wr = /* @__PURE__ */ G("compareUnits", ["typed"], (e) => {
   var {
     typed: r
   } = e;
@@ -5098,7 +5098,7 @@ var Wr = /* @__PURE__ */ J("compareUnits", ["typed"], (e) => {
       return r.find(t, [n.valueType(), i.valueType()])(n.value, i.value);
     })
   };
-}), Dt = "equalScalar", ls = ["typed", "config"], hs = /* @__PURE__ */ J(Dt, ls, (e) => {
+}), Dt = "equalScalar", ls = ["typed", "config"], hs = /* @__PURE__ */ G(Dt, ls, (e) => {
   var {
     typed: r,
     config: t
@@ -5126,7 +5126,7 @@ var Wr = /* @__PURE__ */ J("compareUnits", ["typed"], (e) => {
     }
   }, n);
 });
-J(Dt, ["typed", "config"], (e) => {
+G(Dt, ["typed", "config"], (e) => {
   var {
     typed: r,
     config: t
@@ -5137,7 +5137,7 @@ J(Dt, ["typed", "config"], (e) => {
     }
   });
 });
-var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PURE__ */ J(ps, Ds, (e) => {
+var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PURE__ */ G(ps, Ds, (e) => {
   var {
     typed: r,
     equalScalar: t,
@@ -5227,12 +5227,12 @@ var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PUR
     for (v = 0, A = s._size.length; v < A; v++)
       De(C[v], s._size[v]), De(b[v], s._size[v]);
     var B = s._values, M = s._index, S = s._ptr, R = d.dimension(0), V = d.dimension(1), q = [], W = [];
-    R.forEach(function(G, ee) {
-      W[G] = ee[0], q[G] = !0;
+    R.forEach(function(J, ee) {
+      W[J] = ee[0], q[J] = !0;
     });
     var O = B ? [] : void 0, Q = [], Y = [];
-    return V.forEach(function(G) {
-      for (Y.push(Q.length), y = S[G], w = S[G + 1]; y < w; y++)
+    return V.forEach(function(J) {
+      for (Y.push(Q.length), y = S[J], w = S[J + 1]; y < w; y++)
         v = M[y], q[v] === !0 && (Q.push(W[v]), O && O.push(B[y]));
     }), Y.push(Q.length), new i({
       values: O,
@@ -5371,8 +5371,8 @@ var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PUR
       throw new TypeError("Array expected");
     if (s.length !== 2)
       throw new Error("Sparse matrices can only be reshaped in two dimensions");
-    s.forEach(function(G) {
-      if (!de(G) || !ge(G) || G <= -2 || G === 0)
+    s.forEach(function(J) {
+      if (!de(J) || !ge(J) || J <= -2 || J === 0)
         throw new TypeError("Invalid size, must contain positive integers or -1 (size: " + Be(s) + ")");
     });
     var g = this._size[0] * this._size[1];
@@ -5433,8 +5433,8 @@ var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PUR
         }
       else {
         for (var Q = {}, Y = V; Y < q; Y++) {
-          var G = s._index[Y];
-          Q[G] = s._values[Y];
+          var J = s._index[Y];
+          Q[J] = s._values[Y];
         }
         for (var ee = d; ee <= g; ee++) {
           var oe = ee in Q ? Q[ee] : 0;
@@ -5545,10 +5545,10 @@ var ps = "SparseMatrix", Ds = ["typed", "equalScalar", "Matrix"], ms = /* @__PUR
       throw new TypeError("Array expected, size parameter");
     if (s.length !== 2)
       throw new Error("Only two dimensions matrix are supported");
-    if (s = s.map(function(G) {
-      if (Ee(G) && (G = G.toNumber()), !de(G) || !ge(G) || G < 1)
+    if (s = s.map(function(J) {
+      if (Ee(J) && (J = J.toNumber()), !de(J) || !ge(J) || J < 1)
         throw new Error("Size values must be positive integers");
-      return G;
+      return J;
     }), g) {
       if (Ee(g) && (g = g.toNumber()), !de(g) || !ge(g))
         throw new TypeError("The parameter k must be an integer number");
@@ -5648,7 +5648,7 @@ function ws(e) {
     throw new SyntaxError('String "' + e.input + '" is not a valid number');
   return u;
 }
-var ys = /* @__PURE__ */ J(ds, gs, (e) => {
+var ys = /* @__PURE__ */ G(ds, gs, (e) => {
   var {
     typed: r
   } = e, t = r("number", {
@@ -5699,7 +5699,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
   return t.fromJSON = function(n) {
     return parseFloat(n.value);
   }, t;
-}), Fs = "bignumber", Es = ["typed", "BigNumber"], As = /* @__PURE__ */ J(Fs, Es, (e) => {
+}), Fs = "bignumber", Es = ["typed", "BigNumber"], As = /* @__PURE__ */ G(Fs, Es, (e) => {
   var {
     typed: r,
     BigNumber: t
@@ -5740,7 +5740,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
     },
     "Array | Matrix": r.referToSelf((n) => (i) => Ve(i, n))
   });
-}), bs = "complex", Cs = ["typed", "Complex"], _s = /* @__PURE__ */ J(bs, Cs, (e) => {
+}), bs = "complex", Cs = ["typed", "Complex"], _s = /* @__PURE__ */ G(bs, Cs, (e) => {
   var {
     typed: r,
     Complex: t
@@ -5780,7 +5780,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
     },
     "Array | Matrix": r.referToSelf((n) => (i) => Ve(i, n))
   });
-}), Bs = "fraction", Ns = ["typed", "Fraction"], Ms = /* @__PURE__ */ J(Bs, Ns, (e) => {
+}), Bs = "fraction", Ns = ["typed", "Fraction"], Ms = /* @__PURE__ */ G(Bs, Ns, (e) => {
   var {
     typed: r,
     Fraction: t
@@ -5821,7 +5821,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
     },
     "Array | Matrix": r.referToSelf((n) => (i) => Ve(i, n))
   });
-}), Cn = "matrix", Ss = ["typed", "Matrix", "DenseMatrix", "SparseMatrix"], xs = /* @__PURE__ */ J(Cn, Ss, (e) => {
+}), Cn = "matrix", Ss = ["typed", "Matrix", "DenseMatrix", "SparseMatrix"], xs = /* @__PURE__ */ G(Cn, Ss, (e) => {
   var {
     typed: r,
     Matrix: t,
@@ -5854,7 +5854,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       return new i(o, c);
     throw new TypeError("Unknown matrix type " + JSON.stringify(a) + ".");
   }
-}), _n = "unaryMinus", Ts = ["typed"], Is = /* @__PURE__ */ J(_n, Ts, (e) => {
+}), _n = "unaryMinus", Ts = ["typed"], Is = /* @__PURE__ */ G(_n, Ts, (e) => {
   var {
     typed: r
   } = e;
@@ -5870,7 +5870,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
     "Array | Matrix": r.referToSelf((t) => (n) => Ve(n, t, !0))
     // TODO: add support for string
   });
-}), Bn = "abs", zs = ["typed"], Os = /* @__PURE__ */ J(Bn, zs, (e) => {
+}), Bn = "abs", zs = ["typed"], Os = /* @__PURE__ */ G(Bn, zs, (e) => {
   var {
     typed: r
   } = e;
@@ -5881,7 +5881,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
     // deep map collection, skip zeros since abs(0) = 0
     "Array | Matrix": r.referToSelf((t) => (n) => Ve(n, t, !0))
   });
-}), Nn = "addScalar", $s = ["typed"], Ps = /* @__PURE__ */ J(Nn, $s, (e) => {
+}), Nn = "addScalar", $s = ["typed"], Ps = /* @__PURE__ */ G(Nn, $s, (e) => {
   var {
     typed: r
   } = e;
@@ -5909,7 +5909,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       return u.value = r.find(t, [u.valueType(), i.valueType()])(u.value, i.value), u.fixPrefix = !1, u;
     })
   });
-}), Mn = "subtractScalar", Us = ["typed"], qs = /* @__PURE__ */ J(Mn, Us, (e) => {
+}), Mn = "subtractScalar", Us = ["typed"], qs = /* @__PURE__ */ G(Mn, Us, (e) => {
   var {
     typed: r
   } = e;
@@ -5937,7 +5937,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       return u.value = r.find(t, [u.valueType(), i.valueType()])(u.value, i.value), u.fixPrefix = !1, u;
     })
   });
-}), Rs = "matAlgo11xS0s", Ls = ["typed", "equalScalar"], Hi = /* @__PURE__ */ J(Rs, Ls, (e) => {
+}), Rs = "matAlgo11xS0s", Ls = ["typed", "equalScalar"], Hi = /* @__PURE__ */ G(Rs, Ls, (e) => {
   var {
     typed: r,
     equalScalar: t
@@ -5963,7 +5963,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       datatype: s
     });
   };
-}), Hs = "matAlgo12xSfs", Zs = ["typed", "DenseMatrix"], Vr = /* @__PURE__ */ J(Hs, Zs, (e) => {
+}), Hs = "matAlgo12xSfs", Zs = ["typed", "DenseMatrix"], Vr = /* @__PURE__ */ G(Hs, Zs, (e) => {
   var {
     typed: r,
     DenseMatrix: t
@@ -5988,7 +5988,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       datatype: s
     });
   };
-}), Ws = "matAlgo14xDs", Vs = ["typed"], nn = /* @__PURE__ */ J(Ws, Vs, (e) => {
+}), Ws = "matAlgo14xDs", Vs = ["typed"], nn = /* @__PURE__ */ G(Ws, Vs, (e) => {
   var {
     typed: r
   } = e;
@@ -6012,7 +6012,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
         f[D] = t(n, i + 1, u, u[i + 1], a[D], c, l);
     return f;
   }
-}), Js = "matAlgo03xDSf", Gs = ["typed"], Bt = /* @__PURE__ */ J(Js, Gs, (e) => {
+}), Gs = "matAlgo03xDSf", Js = ["typed"], Bt = /* @__PURE__ */ G(Gs, Js, (e) => {
   var {
     typed: r
   } = e;
@@ -6042,7 +6042,7 @@ var ys = /* @__PURE__ */ J(ds, gs, (e) => {
       datatype: l === n._datatype && m === i._datatype ? g : void 0
     });
   };
-}), Ys = "matAlgo13xDD", ks = ["typed"], Xs = /* @__PURE__ */ J(Ys, ks, (e) => {
+}), Ys = "matAlgo13xDD", ks = ["typed"], Xs = /* @__PURE__ */ G(Ys, ks, (e) => {
   var {
     typed: r
   } = e;
@@ -6084,7 +6084,7 @@ function Se(e, r) {
 function Qs(e, r) {
   return mr(e.size(), r) ? e : e.create(Xt(e.valueOf(), r), e.datatype());
 }
-var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Jr = /* @__PURE__ */ J(Ks, js, (e) => {
+var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Gr = /* @__PURE__ */ G(Ks, js, (e) => {
   var {
     typed: r,
     matrix: t
@@ -6111,7 +6111,7 @@ var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Jr = /* @__PURE__ */
     var D = o.sS !== void 0 ? o.sS : o.Ss;
     return a ? (o.Ss && (l["SparseMatrix," + f] = (p, m) => o.Ss(p, m, a, !1)), D && (l[f + ", SparseMatrix"] = (p, m) => D(m, p, a, !0))) : (o.Ss && (l["SparseMatrix," + f] = r.referToSelf((p) => (m, s) => o.Ss(m, s, p, !1))), D && (l[f + ", SparseMatrix"] = r.referToSelf((p) => (m, s) => D(s, m, p, !0)))), a && a.signatures && Do(l, a.signatures), l;
   };
-}), ef = "matAlgo01xDSid", rf = ["typed"], tf = /* @__PURE__ */ J(ef, rf, (e) => {
+}), ef = "matAlgo01xDSid", rf = ["typed"], tf = /* @__PURE__ */ G(ef, rf, (e) => {
   var {
     typed: r
   } = e;
@@ -6139,7 +6139,7 @@ var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Jr = /* @__PURE__ */
       datatype: l === n._datatype && m === i._datatype ? g : void 0
     });
   };
-}), nf = "matAlgo04xSidSid", uf = ["typed", "equalScalar"], of = /* @__PURE__ */ J(nf, uf, (e) => {
+}), nf = "matAlgo04xSidSid", uf = ["typed", "equalScalar"], of = /* @__PURE__ */ G(nf, uf, (e) => {
   var {
     typed: r,
     equalScalar: t
@@ -6155,20 +6155,20 @@ var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Jr = /* @__PURE__ */
     var C = a && D ? [] : void 0, b = [], B = [], M = a && D ? [] : void 0, S = a && D ? [] : void 0, R = [], V = [], q, W, O, Q, Y;
     for (W = 0; W < E; W++) {
       B[W] = b.length;
-      var G = W + 1;
+      var J = W + 1;
       for (Q = l[W], Y = l[W + 1], O = Q; O < Y; O++)
-        q = c[O], b.push(q), R[q] = G, M && (M[q] = a[O]);
+        q = c[O], b.push(q), R[q] = J, M && (M[q] = a[O]);
       for (Q = m[W], Y = m[W + 1], O = Q; O < Y; O++)
-        if (q = p[O], R[q] === G) {
+        if (q = p[O], R[q] === J) {
           if (M) {
             var ee = w(M[q], D[O]);
             A(ee, y) ? R[q] = null : M[q] = ee;
           }
         } else
-          b.push(q), V[q] = G, S && (S[q] = D[O]);
+          b.push(q), V[q] = J, S && (S[q] = D[O]);
       if (M && S)
         for (O = B[W]; O < b.length; )
-          q = b[O], R[q] === G ? (C[O] = M[q], O++) : V[q] === G ? (C[O] = S[q], O++) : b.splice(O, 1);
+          q = b[O], R[q] === J ? (C[O] = M[q], O++) : V[q] === J ? (C[O] = S[q], O++) : b.splice(O, 1);
     }
     return B[E] = b.length, i.createSparseMatrix({
       values: C,
@@ -6178,7 +6178,7 @@ var Ks = "matrixAlgorithmSuite", js = ["typed", "matrix"], Jr = /* @__PURE__ */
       datatype: h === i._datatype && d === u._datatype ? v : void 0
     });
   };
-}), af = "matAlgo10xSids", sf = ["typed", "DenseMatrix"], ff = /* @__PURE__ */ J(af, sf, (e) => {
+}), af = "matAlgo10xSids", sf = ["typed", "DenseMatrix"], ff = /* @__PURE__ */ G(af, sf, (e) => {
   var {
     typed: r,
     DenseMatrix: t
@@ -6214,7 +6214,7 @@ function cf(e, r, t, n) {
     return n(i.toNumber());
   };
 }
-var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__PURE__ */ J(Sn, lf, (e) => {
+var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__PURE__ */ G(Sn, lf, (e) => {
   var {
     typed: r,
     config: t,
@@ -6237,7 +6237,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     },
     "Array | Matrix": r.referToSelf((o) => (a) => Ve(a, o))
   });
-}), Df = "multiplyScalar", mf = ["typed"], df = /* @__PURE__ */ J(Df, mf, (e) => {
+}), Df = "multiplyScalar", mf = ["typed"], df = /* @__PURE__ */ G(Df, mf, (e) => {
   var {
     typed: r
   } = e;
@@ -6258,7 +6258,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     "number | Fraction | BigNumber | Complex, Unit": (t, n) => n.multiply(t),
     "Unit, number | Fraction | BigNumber | Complex | Unit": (t, n) => t.multiply(n)
   });
-}), xn = "multiply", gf = ["typed", "matrix", "addScalar", "multiplyScalar", "equalScalar", "dot"], vf = /* @__PURE__ */ J(xn, gf, (e) => {
+}), xn = "multiply", gf = ["typed", "matrix", "addScalar", "multiplyScalar", "equalScalar", "dot"], vf = /* @__PURE__ */ G(xn, gf, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -6319,10 +6319,10 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
   function D(y, w) {
     var C = y._data, b = y._size, B = y._datatype || y.getDataType(), M = w._data, S = w._size, R = w._datatype || w.getDataType(), V = b[0], q = S[1], W, O = n, Q = i;
     B && R && B === R && typeof B == "string" && B !== "mixed" && (W = B, O = r.find(n, [W, W]), Q = r.find(i, [W, W]));
-    for (var Y = [], G = 0; G < q; G++) {
-      for (var ee = Q(C[0], M[0][G]), oe = 1; oe < V; oe++)
-        ee = O(ee, Q(C[oe], M[oe][G]));
-      Y[G] = ee;
+    for (var Y = [], J = 0; J < q; J++) {
+      for (var ee = Q(C[0], M[0][J]), oe = 1; oe < V; oe++)
+        ee = O(ee, Q(C[oe], M[oe][J]));
+      Y[J] = ee;
     }
     return y.createDenseMatrix({
       data: Y,
@@ -6343,8 +6343,8 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     var C = y._data, b = y._size, B = y._datatype || y.getDataType(), M = w._data, S = w._datatype || w.getDataType(), R = b[0], V = b[1], q, W = n, O = i;
     B && S && B === S && typeof B == "string" && B !== "mixed" && (q = B, W = r.find(n, [q, q]), O = r.find(i, [q, q]));
     for (var Q = [], Y = 0; Y < R; Y++) {
-      for (var G = C[Y], ee = O(G[0], M[0]), oe = 1; oe < V; oe++)
-        ee = W(ee, O(G[oe], M[oe]));
+      for (var J = C[Y], ee = O(J[0], M[0]), oe = 1; oe < V; oe++)
+        ee = W(ee, O(J[oe], M[oe]));
       Q[Y] = ee;
     }
     return y.createDenseMatrix({
@@ -6356,17 +6356,17 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
   function d(y, w) {
     var C = y._data, b = y._size, B = y._datatype || y.getDataType(), M = w._data, S = w._size, R = w._datatype || w.getDataType(), V = b[0], q = b[1], W = S[1], O, Q = n, Y = i;
     B && R && B === R && typeof B == "string" && B !== "mixed" && B !== "mixed" && (O = B, Q = r.find(n, [O, O]), Y = r.find(i, [O, O]));
-    for (var G = [], ee = 0; ee < V; ee++) {
+    for (var J = [], ee = 0; ee < V; ee++) {
       var oe = C[ee];
-      G[ee] = [];
+      J[ee] = [];
       for (var ie = 0; ie < W; ie++) {
         for (var fe = Y(oe[0], M[0][ie]), Ce = 1; Ce < q; Ce++)
           fe = Q(fe, Y(oe[Ce], M[Ce][ie]));
-        G[ee][ie] = fe;
+        J[ee][ie] = fe;
       }
     }
     return y.createDenseMatrix({
-      data: G,
+      data: J,
       size: [V, W],
       datatype: B === y._datatype && R === w._datatype ? O : void 0
     });
@@ -6375,8 +6375,8 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     var C = y._data, b = y._size, B = y._datatype || y.getDataType(), M = w._values, S = w._index, R = w._ptr, V = w._size, q = w._datatype || w._data === void 0 ? w._datatype : w.getDataType();
     if (!M)
       throw new Error("Cannot multiply Dense Matrix times Pattern only Matrix");
-    var W = b[0], O = V[1], Q, Y = n, G = i, ee = u, oe = 0;
-    B && q && B === q && typeof B == "string" && B !== "mixed" && (Q = B, Y = r.find(n, [Q, Q]), G = r.find(i, [Q, Q]), ee = r.find(u, [Q, Q]), oe = r.convert(0, Q));
+    var W = b[0], O = V[1], Q, Y = n, J = i, ee = u, oe = 0;
+    B && q && B === q && typeof B == "string" && B !== "mixed" && (Q = B, Y = r.find(n, [Q, Q]), J = r.find(i, [Q, Q]), ee = r.find(u, [Q, Q]), oe = r.convert(0, Q));
     for (var ie = [], fe = [], Ce = [], be = w.createSparseMatrix({
       values: ie,
       index: fe,
@@ -6388,11 +6388,11 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       var we = R[Ae], qe = R[Ae + 1];
       if (qe > we)
         for (var ce = 0, ae = 0; ae < W; ae++) {
-          for (var Je = ae + 1, Te = void 0, le = we; le < qe; le++) {
+          for (var Ge = ae + 1, Te = void 0, le = we; le < qe; le++) {
             var Me = S[le];
-            ce !== Je ? (Te = G(C[ae][Me], M[le]), ce = Je) : Te = Y(Te, G(C[ae][Me], M[le]));
+            ce !== Ge ? (Te = J(C[ae][Me], M[le]), ce = Ge) : Te = Y(Te, J(C[ae][Me], M[le]));
           }
-          ce === Je && !ee(Te, oe) && (fe.push(ae), ie.push(Te));
+          ce === Ge && !ee(Te, oe) && (fe.push(ae), ie.push(Te));
         }
     }
     return Ce[O] = fe.length, be;
@@ -6401,8 +6401,8 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     var C = y._values, b = y._index, B = y._ptr, M = y._datatype || y._data === void 0 ? y._datatype : y.getDataType();
     if (!C)
       throw new Error("Cannot multiply Pattern only Matrix times Dense Matrix");
-    var S = w._data, R = w._datatype || w.getDataType(), V = y._size[0], q = w._size[0], W = [], O = [], Q = [], Y, G = n, ee = i, oe = u, ie = 0;
-    M && R && M === R && typeof M == "string" && M !== "mixed" && (Y = M, G = r.find(n, [Y, Y]), ee = r.find(i, [Y, Y]), oe = r.find(u, [Y, Y]), ie = r.convert(0, Y));
+    var S = w._data, R = w._datatype || w.getDataType(), V = y._size[0], q = w._size[0], W = [], O = [], Q = [], Y, J = n, ee = i, oe = u, ie = 0;
+    M && R && M === R && typeof M == "string" && M !== "mixed" && (Y = M, J = r.find(n, [Y, Y]), ee = r.find(i, [Y, Y]), oe = r.find(u, [Y, Y]), ie = r.convert(0, Y));
     var fe = [], Ce = [];
     Q[0] = 0;
     for (var be = 0; be < q; be++) {
@@ -6410,10 +6410,10 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       if (!oe(Ae, ie))
         for (var we = B[be], qe = B[be + 1], ce = we; ce < qe; ce++) {
           var ae = b[ce];
-          Ce[ae] ? fe[ae] = G(fe[ae], ee(Ae, C[ce])) : (Ce[ae] = !0, O.push(ae), fe[ae] = ee(Ae, C[ce]));
+          Ce[ae] ? fe[ae] = J(fe[ae], ee(Ae, C[ce])) : (Ce[ae] = !0, O.push(ae), fe[ae] = ee(Ae, C[ce]));
         }
     }
-    for (var Je = O.length, Te = 0; Te < Je; Te++) {
+    for (var Ge = O.length, Te = 0; Te < Ge; Te++) {
       var le = O[Te];
       W[Te] = fe[le];
     }
@@ -6429,8 +6429,8 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     var C = y._values, b = y._index, B = y._ptr, M = y._datatype || y._data === void 0 ? y._datatype : y.getDataType();
     if (!C)
       throw new Error("Cannot multiply Pattern only Matrix times Dense Matrix");
-    var S = w._data, R = w._datatype || w.getDataType(), V = y._size[0], q = w._size[0], W = w._size[1], O, Q = n, Y = i, G = u, ee = 0;
-    M && R && M === R && typeof M == "string" && M !== "mixed" && (O = M, Q = r.find(n, [O, O]), Y = r.find(i, [O, O]), G = r.find(u, [O, O]), ee = r.convert(0, O));
+    var S = w._data, R = w._datatype || w.getDataType(), V = y._size[0], q = w._size[0], W = w._size[1], O, Q = n, Y = i, J = u, ee = 0;
+    M && R && M === R && typeof M == "string" && M !== "mixed" && (O = M, Q = r.find(n, [O, O]), Y = r.find(i, [O, O]), J = r.find(u, [O, O]), ee = r.convert(0, O));
     for (var oe = [], ie = [], fe = [], Ce = y.createSparseMatrix({
       values: oe,
       index: ie,
@@ -6441,13 +6441,13 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       fe[we] = ie.length;
       for (var qe = we + 1, ce = 0; ce < q; ce++) {
         var ae = S[ce][we];
-        if (!G(ae, ee))
-          for (var Je = B[ce], Te = B[ce + 1], le = Je; le < Te; le++) {
+        if (!J(ae, ee))
+          for (var Ge = B[ce], Te = B[ce + 1], le = Ge; le < Te; le++) {
             var Me = b[le];
             Ae[Me] !== qe ? (Ae[Me] = qe, ie.push(Me), be[Me] = Y(ae, C[le])) : be[Me] = Q(be[Me], Y(ae, C[le]));
           }
       }
-      for (var Ge = fe[we], ar = ie.length, sr = Ge; sr < ar; sr++) {
+      for (var Je = fe[we], ar = ie.length, sr = Je; sr < ar; sr++) {
         var Ar = ie[sr];
         oe[sr] = be[Ar];
       }
@@ -6455,26 +6455,26 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     return fe[W] = ie.length, Ce;
   }
   function A(y, w) {
-    var C = y._values, b = y._index, B = y._ptr, M = y._datatype || y._data === void 0 ? y._datatype : y.getDataType(), S = w._values, R = w._index, V = w._ptr, q = w._datatype || w._data === void 0 ? w._datatype : w.getDataType(), W = y._size[0], O = w._size[1], Q = C && S, Y, G = n, ee = i;
-    M && q && M === q && typeof M == "string" && M !== "mixed" && (Y = M, G = r.find(n, [Y, Y]), ee = r.find(i, [Y, Y]));
+    var C = y._values, b = y._index, B = y._ptr, M = y._datatype || y._data === void 0 ? y._datatype : y.getDataType(), S = w._values, R = w._index, V = w._ptr, q = w._datatype || w._data === void 0 ? w._datatype : w.getDataType(), W = y._size[0], O = w._size[1], Q = C && S, Y, J = n, ee = i;
+    M && q && M === q && typeof M == "string" && M !== "mixed" && (Y = M, J = r.find(n, [Y, Y]), ee = r.find(i, [Y, Y]));
     for (var oe = Q ? [] : void 0, ie = [], fe = [], Ce = y.createSparseMatrix({
       values: oe,
       index: ie,
       ptr: fe,
       size: [W, O],
       datatype: M === y._datatype && q === w._datatype ? Y : void 0
-    }), be = Q ? [] : void 0, Ae = [], we, qe, ce, ae, Je, Te, le, Me, Ge = 0; Ge < O; Ge++) {
-      fe[Ge] = ie.length;
-      var ar = Ge + 1;
-      for (Je = V[Ge], Te = V[Ge + 1], ae = Je; ae < Te; ae++)
+    }), be = Q ? [] : void 0, Ae = [], we, qe, ce, ae, Ge, Te, le, Me, Je = 0; Je < O; Je++) {
+      fe[Je] = ie.length;
+      var ar = Je + 1;
+      for (Ge = V[Je], Te = V[Je + 1], ae = Ge; ae < Te; ae++)
         if (Me = R[ae], Q)
           for (qe = B[Me], ce = B[Me + 1], we = qe; we < ce; we++)
-            le = b[we], Ae[le] !== ar ? (Ae[le] = ar, ie.push(le), be[le] = ee(S[ae], C[we])) : be[le] = G(be[le], ee(S[ae], C[we]));
+            le = b[we], Ae[le] !== ar ? (Ae[le] = ar, ie.push(le), be[le] = ee(S[ae], C[we])) : be[le] = J(be[le], ee(S[ae], C[we]));
         else
           for (qe = B[Me], ce = B[Me + 1], we = qe; we < ce; we++)
             le = b[we], Ae[le] !== ar && (Ae[le] = ar, ie.push(le));
       if (Q)
-        for (var sr = fe[Ge], Ar = ie.length, br = sr; br < Ar; br++) {
+        for (var sr = fe[Je], Ar = ie.length, br = sr; br < Ar; br++) {
           var zt = ie[br];
           oe[br] = be[zt];
         }
@@ -6519,7 +6519,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       return B;
     })
   });
-}), wf = "matAlgo07xSSf", yf = ["typed", "SparseMatrix"], Nt = /* @__PURE__ */ J(wf, yf, (e) => {
+}), wf = "matAlgo07xSSf", yf = ["typed", "SparseMatrix"], Nt = /* @__PURE__ */ G(wf, yf, (e) => {
   var {
     typed: r,
     SparseMatrix: t
@@ -6555,7 +6555,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       o[m] = c, a[m] = l[D];
     }
   }
-}), Tn = "arg", Ff = ["typed"], Ef = /* @__PURE__ */ J(Tn, Ff, (e) => {
+}), Tn = "arg", Ff = ["typed"], Ef = /* @__PURE__ */ G(Tn, Ff, (e) => {
   var {
     typed: r
   } = e;
@@ -6572,7 +6572,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     // TODO: implement BigNumber support for function arg
     "Array | Matrix": r.referToSelf((t) => (n) => Ve(n, t))
   });
-}), In = "conj", Af = ["typed"], bf = /* @__PURE__ */ J(In, Af, (e) => {
+}), In = "conj", Af = ["typed"], bf = /* @__PURE__ */ G(In, Af, (e) => {
   var {
     typed: r
   } = e;
@@ -6582,7 +6582,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
     Unit: r.referToSelf((t) => (n) => new n.constructor(t(n.toNumeric()), n.formatUnits())),
     "Array | Matrix": r.referToSelf((t) => (n) => Ve(n, t))
   });
-}), zn = "concat", Cf = ["typed", "matrix", "isInteger"], _f = /* @__PURE__ */ J(zn, Cf, (e) => {
+}), zn = "concat", Cf = ["typed", "matrix", "isInteger"], _f = /* @__PURE__ */ G(zn, Cf, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -6617,7 +6617,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       return u.join("");
     }
   });
-}), On = "getMatrixDataType", Bf = ["typed"], Nf = /* @__PURE__ */ J(On, Bf, (e) => {
+}), On = "getMatrixDataType", Bf = ["typed"], Nf = /* @__PURE__ */ G(On, Bf, (e) => {
   var {
     typed: r
   } = e;
@@ -6629,7 +6629,7 @@ var Sn = "log10", lf = ["typed", "config", "Complex"], hf = pt(16), pf = /* @__P
       return n.getDataType();
     }
   });
-}), $n = "identity", Mf = ["typed", "config", "matrix", "BigNumber", "DenseMatrix", "SparseMatrix"], Sf = /* @__PURE__ */ J($n, Mf, (e) => {
+}), $n = "identity", Mf = ["typed", "config", "matrix", "BigNumber", "DenseMatrix", "SparseMatrix"], Sf = /* @__PURE__ */ G($n, Mf, (e) => {
   var {
     typed: r,
     config: t,
@@ -6710,7 +6710,7 @@ function xf() {
 function Wi() {
   throw new Error('No "matrix" implementation available');
 }
-var Pn = "range", Tf = ["typed", "config", "?matrix", "?bignumber", "smaller", "smallerEq", "larger", "largerEq", "add", "isPositive"], If = /* @__PURE__ */ J(Pn, Tf, (e) => {
+var Pn = "range", Tf = ["typed", "config", "?matrix", "?bignumber", "smaller", "smallerEq", "larger", "largerEq", "add", "isPositive"], If = /* @__PURE__ */ G(Pn, Tf, (e) => {
   var {
     typed: r,
     config: t,
@@ -6842,7 +6842,7 @@ var Pn = "range", Tf = ["typed", "config", "?matrix", "?bignumber", "smaller", "
         return null;
     }
   }
-}), Un = "size", zf = ["typed", "config", "?matrix"], Of = /* @__PURE__ */ J(Un, zf, (e) => {
+}), Un = "size", zf = ["typed", "config", "?matrix"], Of = /* @__PURE__ */ G(Un, zf, (e) => {
   var {
     typed: r,
     config: t,
@@ -6860,7 +6860,7 @@ var Pn = "range", Tf = ["typed", "config", "?matrix", "?bignumber", "smaller", "
       return t.matrix === "Array" ? [] : n ? n([], "dense", "number") : Wi();
     }
   });
-}), qn = "subset", $f = ["typed", "matrix", "zeros", "add"], Pf = /* @__PURE__ */ J(qn, $f, (e) => {
+}), qn = "subset", $f = ["typed", "matrix", "zeros", "add"], Pf = /* @__PURE__ */ G(qn, $f, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -6981,7 +6981,7 @@ function Rf(e, r, t) {
   var i = me(e);
   return ni(i, n, t), i;
 }
-var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PURE__ */ J(Ln, Lf, (e) => {
+var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PURE__ */ G(Ln, Lf, (e) => {
   var {
     typed: r,
     config: t,
@@ -7032,7 +7032,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
         throw new Error("Parameters in function zeros must be positive integers");
     });
   }
-}), Zf = "numeric", Wf = ["number", "?bignumber", "?fraction"], Vf = /* @__PURE__ */ J(Zf, Wf, (e) => {
+}), Zf = "numeric", Wf = ["number", "?bignumber", "?fraction"], Vf = /* @__PURE__ */ G(Zf, Wf, (e) => {
   var {
     number: r,
     bignumber: t,
@@ -7059,7 +7059,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
       throw new TypeError("Cannot convert " + a + ' to type "' + c + '"; valid output types are ' + Object.keys(u).join(", "));
     return c === f ? a : u[c](a);
   };
-}), Hn = "divideScalar", Jf = ["typed", "numeric"], Gf = /* @__PURE__ */ J(Hn, Jf, (e) => {
+}), Hn = "divideScalar", Gf = ["typed", "numeric"], Jf = /* @__PURE__ */ G(Hn, Gf, (e) => {
   var {
     typed: r,
     numeric: t
@@ -7083,7 +7083,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     "Unit, number | Complex | Fraction | BigNumber | Unit": (n, i) => n.divide(i),
     "number | Fraction | Complex | BigNumber, Unit": (n, i) => i.divideInto(n)
   });
-}), Zn = "pow", Yf = ["typed", "config", "identity", "multiply", "matrix", "inv", "fraction", "number", "Complex"], kf = /* @__PURE__ */ J(Zn, Yf, (e) => {
+}), Zn = "pow", Yf = ["typed", "config", "identity", "multiply", "matrix", "inv", "fraction", "number", "Complex"], kf = /* @__PURE__ */ G(Zn, Yf, (e) => {
   var {
     typed: r,
     config: t,
@@ -7155,7 +7155,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   function D(p, m) {
     return u(h(p.valueOf(), m));
   }
-}), Mr = "Number of decimals in function round must be an integer", Wn = "round", Xf = ["typed", "config", "matrix", "equalScalar", "zeros", "BigNumber", "DenseMatrix"], Qf = /* @__PURE__ */ J(Wn, Xf, (e) => {
+}), Mr = "Number of decimals in function round must be an integer", Wn = "round", Xf = ["typed", "config", "matrix", "equalScalar", "zeros", "BigNumber", "DenseMatrix"], Qf = /* @__PURE__ */ G(Wn, Xf, (e) => {
   var {
     typed: r,
     config: t,
@@ -7252,7 +7252,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     "number | Complex | BigNumber | Fraction, DenseMatrix": r.referToSelf((D) => (p, m) => i(p, 0) ? u(m.size(), m.storage()) : f(m, p, D, !0)),
     "number | Complex | BigNumber | Fraction, Array": r.referToSelf((D) => (p, m) => f(n(m), p, D, !0).valueOf())
   });
-}), mt = "smaller", Kf = ["typed", "config", "bignumber", "matrix", "DenseMatrix", "concat", "SparseMatrix"], jf = /* @__PURE__ */ J(mt, Kf, (e) => {
+}), mt = "smaller", Kf = ["typed", "config", "bignumber", "matrix", "DenseMatrix", "concat", "SparseMatrix"], jf = /* @__PURE__ */ G(mt, Kf, (e) => {
   var {
     typed: r,
     config: t,
@@ -7269,7 +7269,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   }), f = Vr({
     typed: r,
     DenseMatrix: u
-  }), h = Jr({
+  }), h = Gr({
     typed: r,
     matrix: i,
     concat: o
@@ -7301,7 +7301,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     DS: c,
     Ss: f
   }));
-}), ec = /* @__PURE__ */ J(mt, ["typed", "config"], (e) => {
+}), ec = /* @__PURE__ */ G(mt, ["typed", "config"], (e) => {
   var {
     typed: r,
     config: t
@@ -7311,7 +7311,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
       return i < u && !tr(i, u, t.relTol, t.absTol);
     }
   });
-}), dt = "smallerEq", rc = ["typed", "config", "matrix", "DenseMatrix", "concat", "SparseMatrix"], tc = /* @__PURE__ */ J(dt, rc, (e) => {
+}), dt = "smallerEq", rc = ["typed", "config", "matrix", "DenseMatrix", "concat", "SparseMatrix"], tc = /* @__PURE__ */ G(dt, rc, (e) => {
   var {
     typed: r,
     config: t,
@@ -7327,7 +7327,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   }), l = Vr({
     typed: r,
     DenseMatrix: i
-  }), f = Jr({
+  }), f = Gr({
     typed: r,
     matrix: n,
     concat: u
@@ -7352,7 +7352,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     DS: a,
     Ss: l
   }));
-}), nc = /* @__PURE__ */ J(dt, ["typed", "config"], (e) => {
+}), nc = /* @__PURE__ */ G(dt, ["typed", "config"], (e) => {
   var {
     typed: r,
     config: t
@@ -7362,7 +7362,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
       return i <= u || tr(i, u, t.relTol, t.absTol);
     }
   });
-}), gt = "larger", ic = ["typed", "config", "bignumber", "matrix", "DenseMatrix", "concat", "SparseMatrix"], uc = /* @__PURE__ */ J(gt, ic, (e) => {
+}), gt = "larger", ic = ["typed", "config", "bignumber", "matrix", "DenseMatrix", "concat", "SparseMatrix"], uc = /* @__PURE__ */ G(gt, ic, (e) => {
   var {
     typed: r,
     config: t,
@@ -7379,7 +7379,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   }), f = Vr({
     typed: r,
     DenseMatrix: u
-  }), h = Jr({
+  }), h = Gr({
     typed: r,
     matrix: i,
     concat: o
@@ -7411,7 +7411,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     DS: c,
     Ss: f
   }));
-}), oc = /* @__PURE__ */ J(gt, ["typed", "config"], (e) => {
+}), oc = /* @__PURE__ */ G(gt, ["typed", "config"], (e) => {
   var {
     typed: r,
     config: t
@@ -7421,7 +7421,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
       return i > u && !tr(i, u, t.relTol, t.absTol);
     }
   });
-}), vt = "largerEq", ac = ["typed", "config", "matrix", "DenseMatrix", "concat", "SparseMatrix"], sc = /* @__PURE__ */ J(vt, ac, (e) => {
+}), vt = "largerEq", ac = ["typed", "config", "matrix", "DenseMatrix", "concat", "SparseMatrix"], sc = /* @__PURE__ */ G(vt, ac, (e) => {
   var {
     typed: r,
     config: t,
@@ -7437,7 +7437,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   }), l = Vr({
     typed: r,
     DenseMatrix: i
-  }), f = Jr({
+  }), f = Gr({
     typed: r,
     matrix: n,
     concat: u
@@ -7464,7 +7464,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
     DS: a,
     Ss: l
   }));
-}), fc = /* @__PURE__ */ J(vt, ["typed", "config"], (e) => {
+}), fc = /* @__PURE__ */ G(vt, ["typed", "config"], (e) => {
   var {
     typed: r,
     config: t
@@ -7474,7 +7474,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
       return i >= u || tr(i, u, t.relTol, t.absTol);
     }
   });
-}), cc = "ImmutableDenseMatrix", lc = ["smaller", "DenseMatrix"], hc = /* @__PURE__ */ J(cc, lc, (e) => {
+}), cc = "ImmutableDenseMatrix", lc = ["smaller", "DenseMatrix"], hc = /* @__PURE__ */ G(cc, lc, (e) => {
   var {
     smaller: r,
     DenseMatrix: t
@@ -7554,7 +7554,7 @@ var Ln = "zeros", Lf = ["typed", "config", "matrix", "BigNumber"], Hf = /* @__PU
   }, n;
 }, {
   isClass: !0
-}), pc = "Index", Dc = ["ImmutableDenseMatrix", "getMatrixDataType"], mc = /* @__PURE__ */ J(pc, Dc, (e) => {
+}), pc = "Index", Dc = ["ImmutableDenseMatrix", "getMatrixDataType"], mc = /* @__PURE__ */ G(pc, Dc, (e) => {
   var {
     ImmutableDenseMatrix: r,
     getMatrixDataType: t
@@ -7723,7 +7723,7 @@ Mt(function(e) {
 function St(e) {
   return e[0].precision;
 }
-var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatrix", "SparseMatrix", "concat"], vc = /* @__PURE__ */ J(Jn, gc, (e) => {
+var Gn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatrix", "SparseMatrix", "concat"], vc = /* @__PURE__ */ G(Gn, gc, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -7740,12 +7740,12 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
   }), f = ff({
     typed: r,
     DenseMatrix: u
-  }), h = Jr({
+  }), h = Gr({
     typed: r,
     matrix: t,
     concat: a
   });
-  return r(Jn, {
+  return r(Gn, {
     "any, any": n,
     "any, any, ...any": r.referToSelf((D) => (p, m, s) => {
       for (var d = D(p, m), g = 0; g < s.length; g++)
@@ -7758,7 +7758,7 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
     SS: l,
     Ss: f
   }));
-}), Gn = "dot", wc = ["typed", "addScalar", "multiplyScalar", "conj", "size"], yc = /* @__PURE__ */ J(Gn, wc, (e) => {
+}), Jn = "dot", wc = ["typed", "addScalar", "multiplyScalar", "conj", "size"], yc = /* @__PURE__ */ G(Jn, wc, (e) => {
   var {
     typed: r,
     addScalar: t,
@@ -7766,7 +7766,7 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
     conj: i,
     size: u
   } = e;
-  return r(Gn, {
+  return r(Jn, {
     "Array | DenseMatrix, Array | DenseMatrix": a,
     "SparseMatrix, SparseMatrix": c
   });
@@ -7834,7 +7834,7 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
   function l(f) {
     return ve(f) ? f.size() : u(f);
   }
-}), Yn = "index", Fc = ["typed", "Index"], Ec = /* @__PURE__ */ J(Yn, Fc, (e) => {
+}), Yn = "index", Fc = ["typed", "Index"], Ec = /* @__PURE__ */ G(Yn, Fc, (e) => {
   var {
     typed: r,
     Index: t
@@ -7849,7 +7849,7 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
       return t.apply(o, u), o;
     }
   });
-}), kn = "det", Ac = ["typed", "matrix", "subtractScalar", "multiply", "divideScalar", "isZero", "unaryMinus"], bc = /* @__PURE__ */ J(kn, Ac, (e) => {
+}), kn = "det", Ac = ["typed", "matrix", "subtractScalar", "multiply", "divideScalar", "isZero", "unaryMinus"], bc = /* @__PURE__ */ G(kn, Ac, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -7910,7 +7910,7 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
     var w = l[p[f - 1]][f - 1];
     return D ? a(w) : w;
   }
-}), Xn = "inv", Cc = ["typed", "matrix", "divideScalar", "addScalar", "multiply", "unaryMinus", "det", "identity", "abs"], _c = /* @__PURE__ */ J(Xn, Cc, (e) => {
+}), Xn = "inv", Cc = ["typed", "matrix", "divideScalar", "addScalar", "multiply", "unaryMinus", "det", "identity", "abs"], _c = /* @__PURE__ */ G(Xn, Cc, (e) => {
   var {
     typed: r,
     matrix: t,
@@ -7995,13 +7995,13 @@ var Jn = "add", gc = ["typed", "matrix", "addScalar", "equalScalar", "DenseMatri
   return r.number === "BigNumber" ? Vi(t) : is;
 });
 function Nc(e, r, t) {
-  return J(e, r, t, {
+  return G(e, r, t, {
     recreateOnConfigChange: !0
   });
 }
 var $r = /* @__PURE__ */ ba({
   config: Pe
-}), xt = /* @__PURE__ */ Ma({}), Ji = /* @__PURE__ */ $a({}), un = /* @__PURE__ */ qa({}), or = /* @__PURE__ */ rs({
+}), xt = /* @__PURE__ */ Ma({}), Gi = /* @__PURE__ */ $a({}), un = /* @__PURE__ */ qa({}), or = /* @__PURE__ */ rs({
   Matrix: un
 }), et = /* @__PURE__ */ Bc({
   BigNumber: $r,
@@ -8010,7 +8010,7 @@ var $r = /* @__PURE__ */ ba({
   BigNumber: $r,
   Complex: xt,
   DenseMatrix: or,
-  Fraction: Ji
+  Fraction: Gi
 }), Qt = /* @__PURE__ */ Os({
   typed: te
 }), Tt = /* @__PURE__ */ Ps({
@@ -8022,7 +8022,7 @@ var $r = /* @__PURE__ */ ba({
   typed: te
 }), Mc = /* @__PURE__ */ bf({
   typed: te
-}), Gr = /* @__PURE__ */ hs({
+}), Jr = /* @__PURE__ */ hs({
   config: Pe,
   typed: te
 }), Sc = /* @__PURE__ */ Nf({
@@ -8033,19 +8033,19 @@ var $r = /* @__PURE__ */ ba({
   config: Pe,
   typed: te
 }), Ic = /* @__PURE__ */ fs({
-  equalScalar: Gr,
+  equalScalar: Jr,
   typed: te
 }), zc = /* @__PURE__ */ pf({
   Complex: xt,
   config: Pe,
   typed: te
-}), Gi = /* @__PURE__ */ df({
+}), Ji = /* @__PURE__ */ df({
   typed: te
 }), Yi = /* @__PURE__ */ ys({
   typed: te
 }), Er = /* @__PURE__ */ ms({
   Matrix: un,
-  equalScalar: Gr,
+  equalScalar: Jr,
   typed: te
 }), Oc = /* @__PURE__ */ qs({
   typed: te
@@ -8067,7 +8067,7 @@ var $r = /* @__PURE__ */ ba({
   matrix: Ue,
   typed: te
 }), Xi = /* @__PURE__ */ Ms({
-  Fraction: Ji,
+  Fraction: Gi,
   typed: te
 }), Qi = /* @__PURE__ */ Sf({
   BigNumber: $r,
@@ -8091,7 +8091,7 @@ var $r = /* @__PURE__ */ ba({
   BigNumber: $r,
   DenseMatrix: or,
   config: Pe,
-  equalScalar: Gr,
+  equalScalar: Jr,
   matrix: Ue,
   typed: te,
   zeros: ki
@@ -8114,10 +8114,10 @@ var $r = /* @__PURE__ */ ba({
   SparseMatrix: Er,
   addScalar: Tt,
   concat: Yr,
-  equalScalar: Gr,
+  equalScalar: Jr,
   matrix: Ue,
   typed: te
-}), eu = /* @__PURE__ */ Gf({
+}), eu = /* @__PURE__ */ Jf({
   numeric: Pc,
   typed: te
 }), qc = /* @__PURE__ */ hc({
@@ -8149,7 +8149,7 @@ var $r = /* @__PURE__ */ ba({
 }), Zc = /* @__PURE__ */ yc({
   addScalar: Tt,
   conj: Mc,
-  multiplyScalar: Gi,
+  multiplyScalar: Ji,
   size: Uc,
   typed: te
 }), Lt = /* @__PURE__ */ Ec({
@@ -8158,9 +8158,9 @@ var $r = /* @__PURE__ */ ba({
 }), Tr = /* @__PURE__ */ vf({
   addScalar: Tt,
   dot: Zc,
-  equalScalar: Gr,
+  equalScalar: Jr,
   matrix: Ue,
-  multiplyScalar: Gi,
+  multiplyScalar: Ji,
   typed: te
 }), Ht = /* @__PURE__ */ If({
   bignumber: It,
@@ -8191,7 +8191,7 @@ var $r = /* @__PURE__ */ ba({
   multiply: Tr,
   typed: te,
   unaryMinus: ji
-}), Jc = /* @__PURE__ */ kf({
+}), Gc = /* @__PURE__ */ kf({
   Complex: xt,
   config: Pe,
   fraction: Xi,
@@ -8202,7 +8202,7 @@ var $r = /* @__PURE__ */ ba({
   number: Yi,
   typed: te
 });
-const jn = 299792458, lr = {
+const Yc = ["Hz", "kHz", "MHz", "GHz"], jn = 299792458, lr = {
   Hz: 1,
   kHz: 1e3,
   MHz: 1e6,
@@ -8215,7 +8215,7 @@ const jn = 299792458, lr = {
   um: 1e-6,
   nm: 1e-9
 };
-class Gc {
+class Jc {
   /**
    * Internal storage for frequency unit
    * Defaults to 'Hz' as the base SI unit for frequency
@@ -8734,8 +8734,8 @@ class Gc {
     this._setWavelengthFromTargetUnit(r, "nm");
   }
 }
-const Yc = ["RI", "MA", "DB"], kc = ["S", "Y", "Z", "G", "H"];
-class Xc {
+const kc = ["RI", "MA", "DB"], Xc = ["S", "Y", "Z", "G", "H"];
+class Qc {
   /**
    * Comments in the file header with "!" symbol at the beginning of each row
    */
@@ -8957,12 +8957,36 @@ class Xc {
    * 2. Extracts frequency points
    * 3. Converts raw data into complex numbers based on format
    * 4. Stores the results in the matrix property
+   *
+   * @example
+   * ```typescript
+   * import { Touchstone, Frequency } from 'rf-touchstone';
+   *
+   * const s1pString = `
+   * ! This is a 1-port S-parameter file
+   * # MHz S MA R 50
+   * 100 0.99 -4
+   * 200 0.80 -22
+   * 300 0.707 -45
+   * `;
+   *
+   * const touchstone = new Touchstone();
+   * touchstone.readContent(s1pString, 1);
+   *
+   * console.log(touchstone.comments); // Outputs: [ 'This is a 1-port S-parameter file' ]
+   * console.log(touchstone.format); // Outputs: 'MA'
+   * console.log(touchstone.parameter); // Outputs: 'S'
+   * console.log(touchstone.impedance); // Outputs: 50
+   * console.log(touchstone.nports); // Outputs: 1
+   * console.log(touchstone.frequency?.f_scaled); // Outputs: [ 100, 200, 300 ]
+   * console.log(touchstone.matrix); // Outputs: the parsed matrix data
+   * ```
    */
   readContent(r, t) {
     this.nports = t;
     const n = r.split(`
 `).map((f) => f.trim()).filter((f) => f !== "");
-    this.comments = n.filter((f) => f.startsWith("!")).map((f) => f.slice(1).trim()), this.frequency = new Gc();
+    this.comments = n.filter((f) => f.startsWith("!")).map((f) => f.slice(1).trim()), this.frequency = new Jc();
     const i = n.filter((f) => f.startsWith("#"));
     if (i.length === 0)
       throw new Error('Unable to find the option line starting with "#"');
@@ -9041,7 +9065,7 @@ class Xc {
               break;
             case "DB":
               s = qt({
-                r: Jc(10, D[m] / 20),
+                r: Gc(10, D[m] / 20),
                 phi: p[m] / 180 * et
               });
               break;
@@ -9120,6 +9144,37 @@ class Xc {
    * 1. Comments (if any)
    * 2. Option line with format, parameter type, and impedance
    * 3. Network parameter data in the specified format
+   *
+   * @example
+   * ```typescript
+   * import { Touchstone, Frequency } from 'rf-touchstone';
+   * import { complex } from 'mathjs';
+   *
+   * const touchstone = new Touchstone();
+   *
+   * // Set properties and matrix data (using simplified complex number representation for example)
+   * touchstone.comments = ['Generated by rf-touchstone'];
+   * touchstone.nports = 1;
+   * touchstone.frequency = new Frequency();
+   * touchstone.frequency.unit = 'GHz';
+   * touchstone.frequency.f_scaled = [1.0, 2.0];
+   * touchstone.parameter = 'S';
+   * touchstone.format = 'MA';
+   * touchstone.impedance = 50;
+   * touchstone.matrix = [
+   *   [ // Output port 1
+   *     [complex(0.5, 0.1), complex(0.4, 0.2)] // S11 at each frequency
+   *   ]
+   * ];
+   *
+   * const s1pString = touchstone.writeContent();
+   * console.log(s1pString);
+   * // Expected output (approximately, due to floating point precision):
+   * // ! Generated by rf-touchstone
+   * // # GHz S MA R 50
+   * // 1 0.5099 11.3099
+   * // 2 0.4472 26.5651
+   * ```
    */
   writeContent() {
     this.validate();
@@ -9155,7 +9210,22 @@ class Xc {
   }
 }
 export {
-  Xc as Touchstone,
-  Yc as TouchstoneFormats,
-  kc as TouchstoneParameters
+  Jc as Frequency,
+  Yc as FrequencyUnits,
+  jn as SPEED_OF_LIGHT,
+  Qc as Touchstone,
+  kc as TouchstoneFormats,
+  Xc as TouchstoneParameters,
+  Qt as abs,
+  wt as add,
+  Qn as arg,
+  qt as complex,
+  Lt as index,
+  zc as log10,
+  Tr as multiply,
+  et as pi,
+  Gc as pow,
+  Ht as range,
+  Kn as round,
+  Rt as subset
 };