@danielsimonjr/mathts-core 0.14.1 → 0.15.0

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAOH,YAAY,EACV,WAAW,EACX,MAAM,EACN,WAAW,EACX,WAAW,EACX,aAAa,EACb,OAAO,EACP,QAAQ,EACR,SAAS,EACT,UAAU,EACV,gBAAgB,GACjB,MAAM,uBAAuB,CAAC;AAG/B,OAAO,EACL,OAAO,EACP,SAAS,EACT,CAAC,EACD,YAAY,EACZ,WAAW,EACX,eAAe,GAChB,MAAM,oBAAoB,CAAC;AAG5B,OAAO,EAAE,gBAAgB,EAAE,KAAK,EAAE,MAAM,yBAAyB,CAAC;AAClE,YAAY,EACV,oBAAoB,EACpB,gBAAgB,EAChB,kBAAkB,EAClB,SAAS,GACV,MAAM,yBAAyB,CAAC;AAGjC,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAE/C,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,YAAY,EAAE,aAAa,EAAE,iBAAiB,EAAE,MAAM,uBAAuB,CAAC;AAG9E,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,gBAAgB,CAAC;AAG3F,OAAO,EACL,QAAQ,EACR,UAAU,EACV,aAAa,EACb,YAAY,EACZ,gBAAgB,EAChB,aAAa,EACb,cAAc,EACd,gBAAgB,GACjB,MAAM,qBAAqB,CAAC;AAG7B,OAAO,EACL,SAAS,EACT,WAAW,EACX,cAAc,EACd,aAAa,EACb,iBAAiB,EACjB,aAAa,EACb,YAAY,EACZ,WAAW,EACX,aAAa,EACb,cAAc,GACf,MAAM,sBAAsB,CAAC;AAC9B,YAAY,EAAE,eAAe,EAAE,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAG1E,OAAO,EACL,SAAS,EACT,cAAc,EACd,cAAc,EACd,YAAY,EACZ,GAAG,EACH,GAAG,EACH,GAAG,EACH,KAAK,EACL,KAAK,EACL,SAAS,EACT,MAAM,GACP,MAAM,wBAAwB,CAAC;AAChC,YAAY,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAG5D,OAAO,EACL,IAAI,EACJ,MAAM,IAAI,WAAW,EACrB,sBAAsB,EACtB,cAAc,EACd,aAAa,EACb,GAAG,GACJ,MAAM,iBAAiB,CAAC;AACzB,YAAY,EAAE,UAAU,EAAE,OAAO,EAAE,YAAY,EAAE,MAAM,iBAAiB,CAAC;AACzE,OAAO,EACL,UAAU,EACV,aAAa,EACb,SAAS,EACT,YAAY,EACZ,UAAU,GACX,MAAM,6BAA6B,CAAC;AACrC,OAAO,EAAE,WAAW,EAAE,aAAa,EAAE,SAAS,EAAE,MAAM,0BAA0B,CAAC;AAMjF,OAAO,EAEL,SAAS,EACT,iBAAiB,EACjB,KAAK,EACL,MAAM,EACN,mBAAmB,EACnB,YAAY,EAGZ,YAAY,EACZ,kBAAkB,EAGlB,QAAQ,EACR,SAAS,EACT,QAAQ,EACR,QAAQ,EACR,OAAO,EACP,UAAU,EACV,QAAQ,EACR,MAAM,EACN,WAAW,EAGX,QAAQ,EACR,aAAa,EACb,cAAc,EAGd,MAAM,EAGN,aAAa,EACb,oBAAoB,EAGpB,mBAAmB,GACpB,MAAM,kBAAkB,CAAC;AAE1B,YAAY,EACV,aAAa,EACb,aAAa,EACb,OAAO,EACP,aAAa,EACb,iBAAiB,EACjB,OAAO,EACP,WAAW,GACZ,MAAM,kBAAkB,CAAC;AAM1B,OAAO,EACL,gBAAgB,EAChB,aAAa,EACb,QAAQ,EACR,IAAI,EACJ,cAAc,GACf,MAAM,oBAAoB,CAAC;AAE5B,YAAY,EACV,YAAY,EACZ,eAAe,EACf,mBAAmB,EACnB,aAAa,GACd,MAAM,oBAAoB,CAAC;AAW5B,eAAO,MAAM,OAAO,EAAE,MAAwB,CAAC;AAI/C,OAAO,EACL,WAAW,EACX,WAAW,EACX,KAAK,EACL,WAAW,EACX,cAAc,EACd,cAAc,EACd,oBAAoB,GACrB,MAAM,qBAAqB,CAAC"}

package/dist/index.js

@@ -44,7 +44,7 @@ import {
   round,
   sign,
   subtractScalar
-} from "./chunk-VPK5YTKG.js";
+} from "./chunk-L5MWGOLB.js";
 
 // src/types/matrix/Range.ts
 var name = "Range";
@@ -75,12 +75,21 @@ var createRangeClass = /* @__PURE__ */ factory(
        */
       step;
       /**
-       * Cached primitive array representation
+       * Cached primitive array representation.
+       *
+       * @internal Not part of the public API.
+       *
+       * These four were ES-private (`#cache`) until TypeScript 7. This class is
+       * returned from a `factory()` call, so declaration emit describes it as an
+       * ANONYMOUS class type -- and TS4094 forbids private or protected members
+       * there, because an anonymous type has no name to refer them to. The
+       * underscore prefix keeps the "do not touch" signal without the emit
+       * failure. Restore `#` if the class is ever hoisted out of the factory.
        */
-      #cache = null;
-      #cacheStart;
-      #cacheEnd;
-      #cacheStep;
+      _cache = null;
+      _cacheStart;
+      _cacheEnd;
+      _cacheStep;
       constructor(start, end, step) {
         if (!(this instanceof Range2)) {
           throw new SyntaxError("Constructor must be called with the new operator");
@@ -280,14 +289,14 @@ var createRangeClass = /* @__PURE__ */ factory(
        * @returns {Array} array
        */
       valueOf() {
-        if (this.#cache !== null && this.start === this.#cacheStart && this.end === this.#cacheEnd && this.step === this.#cacheStep) {
-          return this.#cache;
+        if (this._cache !== null && this.start === this._cacheStart && this.end === this._cacheEnd && this.step === this._cacheStep) {
+          return this._cache;
         }
-        this.#cache = this.toArray();
-        this.#cacheStart = this.start;
-        this.#cacheEnd = this.end;
-        this.#cacheStep = this.step;
-        return this.#cache;
+        this._cache = this.toArray();
+        this._cacheStart = this.start;
+        this._cacheEnd = this.end;
+        this._cacheStep = this.step;
+        return this._cache;
       }
       /**
        * Get a string representation of the range, with optional formatting options.
@@ -1301,13 +1310,60 @@ function sumSquaredDeviations(xs) {
   const n = xs.length;
   if (n < 2) return 0;
   const mean = pairwiseSum(xs) / n;
-  const d = new Float64Array(n);
-  for (let i = 0; i < n; i++) d[i] = xs[i] - mean;
-  const sumD = pairwiseSum(d);
-  const sumDD = pairwiseDot(d, d);
+  const { sumD, sumDD } = pairwiseDeviationSums(xs, mean, 0, n);
   const corrected = sumDD - sumD * sumD / n;
   return corrected < 0 ? 0 : corrected;
 }
+function pairwiseDeviationSums(xs, mean, start, end) {
+  const n = end - start;
+  if (n <= PAIRWISE_BLOCK) {
+    let d0 = 0, d1 = 0, d2 = 0, d3 = 0, d4 = 0, d5 = 0, d6 = 0, d7 = 0;
+    let q0 = 0, q1 = 0, q2 = 0, q3 = 0, q4 = 0, q5 = 0, q6 = 0, q7 = 0;
+    let i = start;
+    const limit = start + (n - n % 8);
+    for (; i < limit; i += 8) {
+      const a0 = xs[i] - mean;
+      const a1 = xs[i + 1] - mean;
+      const a2 = xs[i + 2] - mean;
+      const a3 = xs[i + 3] - mean;
+      const a4 = xs[i + 4] - mean;
+      const a5 = xs[i + 5] - mean;
+      const a6 = xs[i + 6] - mean;
+      const a7 = xs[i + 7] - mean;
+      d0 += a0;
+      q0 += a0 * a0;
+      d1 += a1;
+      q1 += a1 * a1;
+      d2 += a2;
+      q2 += a2 * a2;
+      d3 += a3;
+      q3 += a3 * a3;
+      d4 += a4;
+      q4 += a4 * a4;
+      d5 += a5;
+      q5 += a5 * a5;
+      d6 += a6;
+      q6 += a6 * a6;
+      d7 += a7;
+      q7 += a7 * a7;
+    }
+    let tailD = 0;
+    let tailDD = 0;
+    for (; i < end; i++) {
+      const a = xs[i] - mean;
+      tailD += a;
+      tailDD += a * a;
+    }
+    return {
+      sumD: d0 + d1 + (d2 + d3) + (d4 + d5 + (d6 + d7)) + tailD,
+      sumDD: q0 + q1 + (q2 + q3) + (q4 + q5 + (q6 + q7)) + tailDD
+    };
+  }
+  const half = (n >> 1) - (n >> 1) % 8;
+  const left = pairwiseDeviationSums(xs, mean, start, start + half);
+  const right = pairwiseDeviationSums(xs, mean, start + half, end);
+  return { sumD: left.sumD + right.sumD, sumDD: left.sumDD + right.sumDD };
+}
 function neumaierCumsum(xs, out) {
   let sum = 0;
   let c = 0;
@@ -1321,7 +1377,7 @@ function neumaierCumsum(xs, out) {
 }
 
 // src/index.ts
-var VERSION = "0.14.1";
+var VERSION = "0.15.0";
 export {
   ALL_UNITS,
   BASE_UNITS,