@danielsimonjr/mathts-core 0.8.0 → 0.10.0

package/dist/{chunk-AMOHH5MP.js → chunk-2WW76EEF.js}

@@ -2410,23 +2410,28 @@ var BigNumber = class _BigNumber {
   // ============================================================
   // Comparison Methods
   // ============================================================
+  // All comparisons coerce their argument (number/string → BigNumber) exactly like the arithmetic
+  // methods, so `bignumber(8).lessThanOrEqual(3)` behaves like `.add(3)` does. Without the
+  // coercion the raw number reached the static comparator, which read `undefined._sign` and
+  // returned garbage — silently breaking every factory-layer consumer that compares against a
+  // number literal (isPrime, quantiles, …).
   equals(other) {
-    return _BigNumber.compare(this, other) === 0;
+    return _BigNumber.compare(this, this.ensureBigNumber(other)) === 0;
   }
   lessThan(other) {
-    return _BigNumber.compare(this, other) < 0;
+    return _BigNumber.compare(this, this.ensureBigNumber(other)) < 0;
   }
   lessThanOrEqual(other) {
-    return _BigNumber.compare(this, other) <= 0;
+    return _BigNumber.compare(this, this.ensureBigNumber(other)) <= 0;
   }
   greaterThan(other) {
-    return _BigNumber.compare(this, other) > 0;
+    return _BigNumber.compare(this, this.ensureBigNumber(other)) > 0;
   }
   greaterThanOrEqual(other) {
-    return _BigNumber.compare(this, other) >= 0;
+    return _BigNumber.compare(this, this.ensureBigNumber(other)) >= 0;
   }
   compareTo(other) {
-    return _BigNumber.compare(this, other);
+    return _BigNumber.compare(this, this.ensureBigNumber(other));
   }
   compare(other) {
     return this.compareTo(other);
@@ -2437,7 +2442,7 @@ var BigNumber = class _BigNumber {
    * Named `gt` for Decimal.js / mathjs formatter compatibility.
    */
   gt(other) {
-    return this.greaterThan(this.ensureBigNumber(other));
+    return this.greaterThan(other);
   }
   // ============================================================
   // Rounding Methods

package/dist/index.d.ts

@@ -1720,13 +1720,13 @@ declare class BigNumber implements MathTSValue {
      * Square root
      */
     sqrt(): BigNumber;
-    equals(other: BigNumber): boolean;
-    lessThan(other: BigNumber): boolean;
-    lessThanOrEqual(other: BigNumber): boolean;
-    greaterThan(other: BigNumber): boolean;
-    greaterThanOrEqual(other: BigNumber): boolean;
-    compareTo(other: BigNumber): number;
-    compare(other: BigNumber): number;
+    equals(other: BigNumber | number | string): boolean;
+    lessThan(other: BigNumber | number | string): boolean;
+    lessThanOrEqual(other: BigNumber | number | string): boolean;
+    greaterThan(other: BigNumber | number | string): boolean;
+    greaterThanOrEqual(other: BigNumber | number | string): boolean;
+    compareTo(other: BigNumber | number | string): number;
+    compare(other: BigNumber | number | string): number;
     /**
      * Returns true iff this > other.
      * Accepts the same argument types as `add`/`multiply` (BigNumber, number, or string).
@@ -2538,6 +2538,23 @@ declare function pairwiseDot(a: ArrayLike<number>, b: ArrayLike<number>, start?:
  * The two ranges are assumed the same length; callers guard that.
  */
 declare function scaledDistance(a: ArrayLike<number>, b: ArrayLike<number>): number;
+/**
+ * Sum of squared deviations from the mean, `Σ(xᵢ − x̄)²` — the **corrected two-pass** form, the
+ * numerator of a numerically stable variance.
+ *
+ * Variance is where large means bite: the deviations `xᵢ − x̄` can be O(1) while the values sit on
+ * a huge pedestal, so any error in `x̄` rides straight into every deviation. Two things fix it:
+ *   1. compute the mean with {@link pairwiseSum}, not a naive running total; and
+ *   2. subtract the residual mean-bias term `(Σd)²/n` — `Σd` is zero in exact arithmetic but not in
+ *      floating point, and that leftover is exactly the systematic error a plain `Σd²` carries.
+ *
+ * Measured on 1e9-pedestal data: the plain naive two-pass (what shipped) lands ~1e-7 relative error;
+ * this lands ~1e-16 — better than `np.var`, which uses the uncorrected two-pass (~1e-13).
+ *
+ * Divide the result by `n` (uncorrected), `n − 1` (unbiased/sample), or `n + 1` (biased) for the
+ * corresponding variance. Returns 0 for fewer than two elements.
+ */
+declare function sumSquaredDeviations(xs: ArrayLike<number>): number;
 /**
  * Cumulative sum with **Neumaier compensation** — each prefix total is written to `out[i]` carrying
  * the low-order bits a naive running sum throws away.
@@ -2562,4 +2579,4 @@ declare function neumaierCumsum(xs: ArrayLike<number>, out: {
 
 declare const VERSION = "0.1.0";
 
-export { ALL_UNITS, BASE_UNITS, BEST_PREFIXES, BIGNUMBER_E, BIGNUMBER_LN10, BIGNUMBER_LN2, BIGNUMBER_NEG_ONE, BIGNUMBER_ONE, BIGNUMBER_PI, BIGNUMBER_TEN, BIGNUMBER_ZERO, type BackendType, BigNumber, type BigNumberConfig, COMPLEX_NEG_ONE, COMPLEX_ONE, COMPLEX_ZERO, Complex, type ConversionDef, DEFAULT_CONFIG, DERIVED_UNITS, DIMENSIONLESS, DUAL_UNARY_RULES, DimensionMismatchError, type Dimensions, Dual, type DualUnaryRule, type DualUnaryRuleName, E, FRACTION_HALF, FRACTION_NEG_ONE, FRACTION_ONE, FRACTION_QUARTER, FRACTION_THIRD, FRACTION_ZERO, type FactoryDependencies, type FactoryFunction, type FactoryImport, Fraction, FunctionRegistry, I, type IBigNumber, type IComplex, type IFraction, type IMatrix, LN10, LN2, LOG10E, LOG2E, MATHTS_CONVERSIONS, MATHTS_TYPES, type MathTSConfig, type MathTSValue, type MatrixBackend, type MatrixDimensions, type NumericScalar, type NumericType, PHI, PI, Range, type RangeForEachCallback, type RangeFormatOptions, type RangeJSON, type RangeMapCallback, type RoundingMode, SI_PREFIXES, SQRT1_2, SQRT2, type Scalar, TAU, type TypeDef, TypeRegistry, UNIT_ALIASES, type UnitDef, UnitParseError, VERSION, abs, addScalar, createFactory, createMathTSTyped, createRangeClass, createTypedFunction, dim, divideScalar, equal, fix, getPrefix, getUnitDef, initTypedWasm, isArray, isBigInt, isBigNumber, isBoolean, isComplex, isDenseMatrix, isDual, isFraction, isFunction, isMatrix, isNull, isNumber, isNumeric, isObject, isSparseMatrix, isString, isTypedWasmAvailable, isUndefined, isUnit, math, mathTyped, multiplyScalar, neumaierCumsum, neumaierSum, norm2, number, pairwiseDot, pairwiseSum, pow, registerNativeTypes, registry, round, scaledDistance, subtractScalar };
+export { ALL_UNITS, BASE_UNITS, BEST_PREFIXES, BIGNUMBER_E, BIGNUMBER_LN10, BIGNUMBER_LN2, BIGNUMBER_NEG_ONE, BIGNUMBER_ONE, BIGNUMBER_PI, BIGNUMBER_TEN, BIGNUMBER_ZERO, type BackendType, BigNumber, type BigNumberConfig, COMPLEX_NEG_ONE, COMPLEX_ONE, COMPLEX_ZERO, Complex, type ConversionDef, DEFAULT_CONFIG, DERIVED_UNITS, DIMENSIONLESS, DUAL_UNARY_RULES, DimensionMismatchError, type Dimensions, Dual, type DualUnaryRule, type DualUnaryRuleName, E, FRACTION_HALF, FRACTION_NEG_ONE, FRACTION_ONE, FRACTION_QUARTER, FRACTION_THIRD, FRACTION_ZERO, type FactoryDependencies, type FactoryFunction, type FactoryImport, Fraction, FunctionRegistry, I, type IBigNumber, type IComplex, type IFraction, type IMatrix, LN10, LN2, LOG10E, LOG2E, MATHTS_CONVERSIONS, MATHTS_TYPES, type MathTSConfig, type MathTSValue, type MatrixBackend, type MatrixDimensions, type NumericScalar, type NumericType, PHI, PI, Range, type RangeForEachCallback, type RangeFormatOptions, type RangeJSON, type RangeMapCallback, type RoundingMode, SI_PREFIXES, SQRT1_2, SQRT2, type Scalar, TAU, type TypeDef, TypeRegistry, UNIT_ALIASES, type UnitDef, UnitParseError, VERSION, abs, addScalar, createFactory, createMathTSTyped, createRangeClass, createTypedFunction, dim, divideScalar, equal, fix, getPrefix, getUnitDef, initTypedWasm, isArray, isBigInt, isBigNumber, isBoolean, isComplex, isDenseMatrix, isDual, isFraction, isFunction, isMatrix, isNull, isNumber, isNumeric, isObject, isSparseMatrix, isString, isTypedWasmAvailable, isUndefined, isUnit, math, mathTyped, multiplyScalar, neumaierCumsum, neumaierSum, norm2, number, pairwiseDot, pairwiseSum, pow, registerNativeTypes, registry, round, scaledDistance, subtractScalar, sumSquaredDeviations };

package/dist/index.js

@@ -44,7 +44,7 @@ import {
   round,
   sign,
   subtractScalar
-} from "./chunk-AMOHH5MP.js";
+} from "./chunk-2WW76EEF.js";
 
 // src/types/matrix/Range.ts
 var name = "Range";
@@ -1297,6 +1297,17 @@ function scaledDistance(a, b) {
   }
   return scale === 0 ? 0 : scale * Math.sqrt(ssq);
 }
+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 corrected = sumDD - sumD * sumD / n;
+  return corrected < 0 ? 0 : corrected;
+}
 function neumaierCumsum(xs, out) {
   let sum = 0;
   let c = 0;
@@ -1411,5 +1422,6 @@ export {
   round,
   scaledDistance,
   subtractScalar,
+  sumSquaredDeviations,
   typed
 };

package/dist/internal.js

@@ -69,7 +69,7 @@ import {
   typeOf,
   unitDependencies,
   warnOnce
-} from "./chunk-AMOHH5MP.js";
+} from "./chunk-2WW76EEF.js";
 export {
   Unit as CoreUnit,
   acosh,

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-core",
-  "version": "0.8.0",
+  "version": "0.10.0",
   "description": "Core types, utilities, and configuration for MathTS",
   "author": "Daniel Simon Jr.",
   "license": "MIT",