universal-physics-tensor 10.0.0 → 12.0.0

package/dist/dimensional/errors.js

@@ -11,7 +11,6 @@ export class UPTError extends Error {
         super(message, options);
         this.name = 'UPTError';
         // Restore prototype chain (needed for instanceof after ES5 transpilation).
-        Object.setPrototypeOf(this, UPTError.prototype);
     }
 }
 /**
@@ -27,7 +26,6 @@ export class DimensionMismatchError extends UPTError {
         this.name = 'DimensionMismatchError';
         this.expected = expected;
         this.actual = actual;
-        Object.setPrototypeOf(this, DimensionMismatchError.prototype);
     }
 }
 /**
@@ -50,7 +48,6 @@ export class DuplicateIndexLabelError extends UPTError {
         this.name = 'DuplicateIndexLabelError';
         this.tensorName = tensorName;
         this.label = label;
-        Object.setPrototypeOf(this, DuplicateIndexLabelError.prototype);
     }
 }
 /**
@@ -78,7 +75,6 @@ export class IndexLabelCollisionError extends UPTError {
         this.label = label;
         this.totalCount = totalCount;
         this.sources = sources;
-        Object.setPrototypeOf(this, IndexLabelCollisionError.prototype);
     }
 }
 /**
@@ -99,7 +95,6 @@ export class VarianceMismatchError extends UPTError {
         this.name = 'VarianceMismatchError';
         this.label = label;
         this.variance = variance;
-        Object.setPrototypeOf(this, VarianceMismatchError.prototype);
     }
 }
 /**
@@ -115,7 +110,6 @@ export class TensorInScalarOpError extends UPTError {
             `non-empty freeIndices. Use 'tensor-product' for tensor multiplication.`);
         this.name = 'TensorInScalarOpError';
         this.op = op;
-        Object.setPrototypeOf(this, TensorInScalarOpError.prototype);
     }
 }
 /**
@@ -128,7 +122,6 @@ export class FreeIndexMismatchError extends UPTError {
     constructor(message) {
         super(message);
         this.name = 'FreeIndexMismatchError';
-        Object.setPrototypeOf(this, FreeIndexMismatchError.prototype);
     }
 }
 /**
@@ -142,7 +135,6 @@ export class TensorProductChildInferenceError extends UPTError {
     constructor(message) {
         super(message);
         this.name = 'TensorProductChildInferenceError';
-        Object.setPrototypeOf(this, TensorProductChildInferenceError.prototype);
     }
 }
 /**
@@ -160,7 +152,6 @@ export class InvalidMetricRankError extends UPTError {
         this.name = 'InvalidMetricRankError';
         this.tensorName = tensorName;
         this.actualRank = actualRank;
-        Object.setPrototypeOf(this, InvalidMetricRankError.prototype);
     }
 }
 /**
@@ -177,7 +168,6 @@ export class MetricSignatureError extends UPTError {
         this.name = 'MetricSignatureError';
         this.tensorName = tensorName;
         this.reason = reason;
-        Object.setPrototypeOf(this, MetricSignatureError.prototype);
     }
 }
 /**
@@ -191,7 +181,6 @@ export class InvalidKroneckerRankError extends UPTError {
             `(exactly one upper + one lower index). Per Part-VIII §VIII.3.`);
         this.name = 'InvalidKroneckerRankError';
         this.actualRank = actualRank;
-        Object.setPrototypeOf(this, InvalidKroneckerRankError.prototype);
     }
 }
 /**
@@ -207,7 +196,6 @@ export class KroneckerVarianceError extends UPTError {
             `Per Part-VIII §VIII.3.`);
         this.name = 'KroneckerVarianceError';
         this.bothVariance = bothVariance;
-        Object.setPrototypeOf(this, KroneckerVarianceError.prototype);
     }
 }
 /**
@@ -224,7 +212,6 @@ export class PartialDerivativeIndexVarianceError extends UPTError {
             `be 'lower'. Per Part-VIII §VIII.4.`);
         this.name = 'PartialDerivativeIndexVarianceError';
         this.label = label;
-        Object.setPrototypeOf(this, PartialDerivativeIndexVarianceError.prototype);
     }
 }
 /**
@@ -238,7 +225,7 @@ export class PartialDerivativeIndexVarianceError extends UPTError {
  * a dimensional → numerical import cycle; re-exported from numerical/index.ts
  * for the public API surface per v0.4.0-Implementation-Plan Task 13.
  *
- * Uses `Object.setPrototypeOf` for correct `instanceof` after ES5 transpilation.
+ * A plain `class extends Error` under the ES2022 target: `instanceof` holds without a prototype reset.
  *
  * @public
  */
@@ -248,7 +235,6 @@ export class DuplicateCoordinateWarning extends Error {
             `'${conflictingIndex}' of the operand; result will silently misbehave. ` +
             `Rename the operand's free index or pick a different wrt label.`);
         this.name = 'DuplicateCoordinateWarning';
-        Object.setPrototypeOf(this, DuplicateCoordinateWarning.prototype);
     }
 }
 //# sourceMappingURL=errors.js.map

package/dist/dimensional/errors.js.map

@@ -1 +1 @@
-{"version":3,"file":"errors.js","sourceRoot":"","sources":["../../src/dimensional/errors.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAIH,4CAA4C;AAC5C,MAAM,OAAO,QAAS,SAAQ,KAAK;IACjC,YAAY,OAAe,EAAE,OAAsC;QACjE,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACxB,IAAI,CAAC,IAAI,GAAG,UAAU,CAAC;QACvB,2EAA2E;QAC3E,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,QAAQ,CAAC,SAAS,CAAC,CAAC;IAClD,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,QAAQ,CAAY;IACpB,MAAM,CAAY;IAClC,YAAY,OAAe,EAAE,QAAmB,EAAE,MAAiB;QACjE,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,sBAAsB,CAAC,SAAS,CAAC,CAAC;IAChE,CAAC;CACF;AAED;;;;;;;;;GASG;AACH,MAAM,OAAO,wBAAyB,SAAQ,QAAQ;IACpC,UAAU,CAAS;IACnB,KAAK,CAAS;IAC9B,YAAY,UAAkB,EAAE,KAAa;QAC3C,KAAK,CACH,WAAW,UAAU,+BAA+B,KAAK,KAAK;YAC9D,mEAAmE;YACnE,wCAAwC,CACzC,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;QACvC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,wBAAwB,CAAC,SAAS,CAAC,CAAC;IAClE,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,wBAAyB,SAAQ,QAAQ;IACpC,KAAK,CAAS;IACd,UAAU,CAAS;IACnC;;;;OAIG;IACa,OAAO,CAAyB;IAChD,YACE,KAAa,EACb,UAAkB,EAClB,OAA+B;QAE/B,MAAM,YAAY,GAChB,OAAO,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC;YAC3B,CAAC,CAAC,cAAc,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG;YACrC,CAAC,CAAC,EAAE,CAAC;QACT,KAAK,CACH,gBAAgB,KAAK,aAAa,UAAU,oBAAoB;YAChE,0BAA0B,YAAY,gCAAgC;YACtE,oDAAoD,CACrD,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,wBAAwB,CAAC,SAAS,CAAC,CAAC;IAClE,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjC,KAAK,CAAS;IACd,QAAQ,CAAoB;IAC5C,YAAY,KAAa,EAAE,QAA2B;QACpD,KAAK,CACH,gBAAgB,KAAK,2CAA2C,QAAQ,KAAK;YAC7E,yDAAyD;YACzD,yDAAyD,KAAK,QAAQ;YACtE,sBAAsB,KAAK,yCAAyC;YACpE,qBAAqB,CACtB,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjC,EAAE,CAAS;IAC3B,YAAY,EAAU;QACpB,KAAK,CACH,aAAa,EAAE,oDAAoD;YACnE,wEAAwE,CACzE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,qBAAqB,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClD,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,sBAAsB,CAAC,SAAS,CAAC,CAAC;IAChE,CAAC;CACF;AAED;;;;;;GAMG;AACH,MAAM,OAAO,gCAAiC,SAAQ,QAAQ;IAC5D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,kCAAkC,CAAC;QAC/C,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,gCAAgC,CAAC,SAAS,CAAC,CAAC;IAC1E,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,UAAU,CAAS;IACnB,UAAU,CAAS;IACnC,YAAY,UAAkB,EAAE,UAAkB;QAChD,KAAK,CACH,kBAAkB,UAAU,cAAc,UAAU,UAAU;YAC9D,qEAAqE,CACtE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,sBAAsB,CAAC,SAAS,CAAC,CAAC;IAChE,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,oBAAqB,SAAQ,QAAQ;IAChC,UAAU,CAAS;IACnB,MAAM,CAAS;IAC/B,YAAY,UAAkB,EAAE,MAAc;QAC5C,KAAK,CACH,kBAAkB,UAAU,MAAM,MAAM,0BAA0B,CACnE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,sBAAsB,CAAC;QACnC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,oBAAoB,CAAC,SAAS,CAAC,CAAC;IAC9D,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,OAAO,yBAA0B,SAAQ,QAAQ;IACrC,UAAU,CAAS;IACnC,YAAY,UAAkB;QAC5B,KAAK,CACH,4BAA4B,UAAU,uBAAuB;YAC7D,+DAA+D,CAChE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,2BAA2B,CAAC;QACxC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,yBAAyB,CAAC,SAAS,CAAC,CAAC;IACnE,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,YAAY,CAAoB;IAChD,YAAY,YAA+B;QACzC,KAAK,CACH,+CAA+C,YAAY,KAAK;YAChE,wDAAwD;YACxD,wBAAwB,CACzB,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,YAAY,GAAG,YAAY,CAAC;QACjC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,sBAAsB,CAAC,SAAS,CAAC,CAAC;IAChE,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,mCAAoC,SAAQ,QAAQ;IAC/C,KAAK,CAAS;IAC9B,YAAY,KAAa;QACvB,KAAK,CACH,6BAA6B,KAAK,8BAA8B;YAChE,iEAAiE;YACjE,oCAAoC,CACrC,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,qCAAqC,CAAC;QAClD,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,mCAAmC,CAAC,SAAS,CAAC,CAAC;IAC7E,CAAC;CACF;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,OAAO,0BAA2B,SAAQ,KAAK;IACnD,YAAY,KAAa,EAAE,gBAAwB;QACjD,KAAK,CACH,6BAA6B,KAAK,sCAAsC;YACxE,IAAI,gBAAgB,oDAAoD;YACxE,gEAAgE,CACjE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,4BAA4B,CAAC;QACzC,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,0BAA0B,CAAC,SAAS,CAAC,CAAC;IACpE,CAAC;CACF","sourcesContent":["/**\n * UPT error hierarchy. All UPT-source errors subclass UPTError so\n * downstream consumers (mathjs, threejs, TensorJS) can discriminate\n * with `err instanceof UPTError`.\n *\n * @module dimensional/errors\n */\n\nimport { Dimension } from './types.js';\n\n/** Base class for all UPT-source errors. */\nexport class UPTError extends Error {\n  constructor(message: string, options?: { readonly cause?: unknown }) {\n    super(message, options);\n    this.name = 'UPTError';\n    // Restore prototype chain (needed for instanceof after ES5 transpilation).\n    Object.setPrototypeOf(this, UPTError.prototype);\n  }\n}\n\n/**\n * Thrown by add/sub and equation-level checks when dimensions disagree.\n * (Moved from src/dimensional/algebra.ts as part of the v0.2.0 UPTError\n * refactor. Behavior unchanged.)\n */\nexport class DimensionMismatchError extends UPTError {\n  public readonly expected: Dimension;\n  public readonly actual: Dimension;\n  constructor(message: string, expected: Dimension, actual: Dimension) {\n    super(message);\n    this.name = 'DimensionMismatchError';\n    this.expected = expected;\n    this.actual = actual;\n    Object.setPrototypeOf(this, DimensionMismatchError.prototype);\n  }\n}\n\n/**\n * Thrown when a tensor-symbol's indices list contains the same label\n * more than once (e.g., T^μ_μ_μ). Per Part-VII §VII.4.\n *\n * NOTE: this error fires on DECLARATION-TIME duplicates within a single\n * tensor-symbol's indices list. Those indices are *free*, not \"dummy\"\n * (a dummy index in tensor calculus is one that is summed over). The\n * earlier name `RepeatedDummyLabelError` was a misnomer; that deprecated\n * alias was removed in v0.4.5.\n */\nexport class DuplicateIndexLabelError extends UPTError {\n  public readonly tensorName: string;\n  public readonly label: string;\n  constructor(tensorName: string, label: string) {\n    super(\n      `Tensor '${tensorName}' has repeated index label '${label}'. ` +\n      `A label may appear at most twice (once upper, once lower) within ` +\n      `a single tensor-symbol's indices list.`,\n    );\n    this.name = 'DuplicateIndexLabelError';\n    this.tensorName = tensorName;\n    this.label = label;\n    Object.setPrototypeOf(this, DuplicateIndexLabelError.prototype);\n  }\n}\n\n/**\n * Thrown when an index label appears more than twice across a\n * tensor-product's operands. Einstein contraction requires exactly two\n * occurrences (one upper + one lower) per dummy label. Per Part-VII §VII.4.\n */\nexport class IndexLabelCollisionError extends UPTError {\n  public readonly label: string;\n  public readonly totalCount: number;\n  /**\n   * Tensor names that contributed occurrences of the colliding label,\n   * if the caller has them available. May be undefined for callers that\n   * don't track per-operand provenance.\n   */\n  public readonly sources?: ReadonlyArray<string>;\n  constructor(\n    label: string,\n    totalCount: number,\n    sources?: ReadonlyArray<string>,\n  ) {\n    const sourceClause =\n      sources && sources.length > 0\n        ? ` (sources: ${sources.join(', ')})`\n        : '';\n    super(\n      `Index label '${label}' appears ${totalCount} times across the ` +\n      `tensor-product operands${sourceClause} (max allowed: 2 for Einstein ` +\n      `contraction). Rename one of the offending indices.`,\n    );\n    this.name = 'IndexLabelCollisionError';\n    this.label = label;\n    this.totalCount = totalCount;\n    this.sources = sources;\n    Object.setPrototypeOf(this, IndexLabelCollisionError.prototype);\n  }\n}\n\n/**\n * Thrown when an index label appears twice in a tensor-product but both\n * occurrences share the same variance (both upper or both lower). v0.2.0\n * has no metric tensor available to raise/lower indices, so such a\n * contraction is rejected at validation time. Per Part-VII §VII.4.\n */\nexport class VarianceMismatchError extends UPTError {\n  public readonly label: string;\n  public readonly variance: 'upper' | 'lower';\n  constructor(label: string, variance: 'upper' | 'lower') {\n    super(\n      `Index label '${label}' appears twice but both with variance '${variance}'. ` +\n      `Einstein contraction requires one upper and one lower. ` +\n      `If this is intentional, use raise(operand, gInverse, '${label}') or ` +\n      `lower(operand, g, '${label}') to traverse variance via the metric ` +\n      `before contracting.`,\n    );\n    this.name = 'VarianceMismatchError';\n    this.label = label;\n    this.variance = variance;\n    Object.setPrototypeOf(this, VarianceMismatchError.prototype);\n  }\n}\n\n/**\n * Thrown by `op '*'` / `'/'` / `'^'` when an operand carries non-empty\n * freeIndices. Scalar operators in v0.2.0 are strict: any tensor-valued\n * operand must instead pass through `tensor-product`. Per Part-VII §VII.5\n * (op-tensor boundary rules).\n */\nexport class TensorInScalarOpError extends UPTError {\n  public readonly op: string;\n  constructor(op: string) {\n    super(\n      `Operator '${op}' is scalar-only; received a tensor argument with ` +\n      `non-empty freeIndices. Use 'tensor-product' for tensor multiplication.`,\n    );\n    this.name = 'TensorInScalarOpError';\n    this.op = op;\n    Object.setPrototypeOf(this, TensorInScalarOpError.prototype);\n  }\n}\n\n/**\n * Thrown by `op '+'` / `'-'` when the operand freeIndices maps differ.\n * Addition / subtraction across tensors requires identical free-index\n * signatures (same labels, same per-label upper/lower counts). Per\n * Part-VII §VII.5.\n */\nexport class FreeIndexMismatchError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'FreeIndexMismatchError';\n    Object.setPrototypeOf(this, FreeIndexMismatchError.prototype);\n  }\n}\n\n/**\n * Thrown by the validator's tensor-product child resolver when a non-tensor\n * operand to a `tensor-product` fails dimension inference. Subclassing\n * UPTError (rather than throwing a plain `Error`) preserves the §14.7\n * forward-compat invariant that downstream consumers can discriminate\n * UPT-source errors with `err instanceof UPTError`.\n */\nexport class TensorProductChildInferenceError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'TensorProductChildInferenceError';\n    Object.setPrototypeOf(this, TensorProductChildInferenceError.prototype);\n  }\n}\n\n/**\n * Thrown when a `metric-tensor` node has rank other than 2. Per\n * Part-VIII §VIII.2 (metric-tensor-rank-2-only TENSOR-RULE). The metric\n * is fundamentally a rank-2 object in v0.3.0; higher-rank generalizations\n * are out of scope.\n */\nexport class InvalidMetricRankError extends UPTError {\n  public readonly tensorName: string;\n  public readonly actualRank: number;\n  constructor(tensorName: string, actualRank: number) {\n    super(\n      `Metric tensor '${tensorName}' has rank ${actualRank}, but a ` +\n      `metric must be rank-2 (exactly two indices). Per Part-VIII §VIII.2.`,\n    );\n    this.name = 'InvalidMetricRankError';\n    this.tensorName = tensorName;\n    this.actualRank = actualRank;\n    Object.setPrototypeOf(this, InvalidMetricRankError.prototype);\n  }\n}\n\n/**\n * Thrown when a `metric-tensor` has mixed-variance indices, an empty\n * signature, an invalid signature string, or other structural problem\n * orthogonal to rank. Per Part-VIII §VIII.2 and §VIII.5 (raise/lower\n * variance requirements).\n */\nexport class MetricSignatureError extends UPTError {\n  public readonly tensorName: string;\n  public readonly reason: string;\n  constructor(tensorName: string, reason: string) {\n    super(\n      `Metric tensor '${tensorName}': ${reason}. Per Part-VIII §VIII.2.`,\n    );\n    this.name = 'MetricSignatureError';\n    this.tensorName = tensorName;\n    this.reason = reason;\n    Object.setPrototypeOf(this, MetricSignatureError.prototype);\n  }\n}\n\n/**\n * Thrown when a `kronecker-delta` node has rank other than 2.\n * Per Part-VIII §VIII.3.\n */\nexport class InvalidKroneckerRankError extends UPTError {\n  public readonly actualRank: number;\n  constructor(actualRank: number) {\n    super(\n      `Kronecker delta has rank ${actualRank}, but must be rank-2 ` +\n      `(exactly one upper + one lower index). Per Part-VIII §VIII.3.`,\n    );\n    this.name = 'InvalidKroneckerRankError';\n    this.actualRank = actualRank;\n    Object.setPrototypeOf(this, InvalidKroneckerRankError.prototype);\n  }\n}\n\n/**\n * Thrown when a `kronecker-delta`'s two indices share the same variance\n * (both upper or both lower). Per Part-VIII §VIII.3, the canonical\n * δ^μ_ν requires mixed variance.\n */\nexport class KroneckerVarianceError extends UPTError {\n  public readonly bothVariance: 'upper' | 'lower';\n  constructor(bothVariance: 'upper' | 'lower') {\n    super(\n      `Kronecker delta indices both have variance '${bothVariance}'; ` +\n      `must be one upper + one lower (canonical δ^μ_ν form). ` +\n      `Per Part-VIII §VIII.3.`,\n    );\n    this.name = 'KroneckerVarianceError';\n    this.bothVariance = bothVariance;\n    Object.setPrototypeOf(this, KroneckerVarianceError.prototype);\n  }\n}\n\n/**\n * Thrown when a `tensor-partial-derivative` node's `wrtIndex.variance`\n * is `'upper'`. The differentiation operator ∂/∂x^μ is fundamentally\n * covariant (lower) regardless of the variance of the coordinate. Per\n * Part-VIII §VIII.4 (pderiv-wrtIndex-always-lower TENSOR-RULE).\n */\nexport class PartialDerivativeIndexVarianceError extends UPTError {\n  public readonly label: string;\n  constructor(label: string) {\n    super(\n      `Partial-derivative index '${label}' has variance 'upper', but ` +\n      `∂_μ is fundamentally covariant — wrtIndex.variance must always ` +\n      `be 'lower'. Per Part-VIII §VIII.4.`,\n    );\n    this.name = 'PartialDerivativeIndexVarianceError';\n    this.label = label;\n    Object.setPrototypeOf(this, PartialDerivativeIndexVarianceError.prototype);\n  }\n}\n\n/**\n * Emitted (not thrown) via `process.emitWarning` when a covariant-derivative's\n * `wrt` coordinate label collides with an existing free index of the operand\n * and the env var `UPT_ALLOW_COORD_SHADOW=1` has opted in to downgraded\n * handling. By default the collision throws `MetricSignatureError` (soundness\n * over friendliness — silent wrong results are the worst failure mode).\n *\n * Lives here (dimensional/errors.ts) rather than numerical/index.ts to avoid\n * a dimensional → numerical import cycle; re-exported from numerical/index.ts\n * for the public API surface per v0.4.0-Implementation-Plan Task 13.\n *\n * Uses `Object.setPrototypeOf` for correct `instanceof` after ES5 transpilation.\n *\n * @public\n */\nexport class DuplicateCoordinateWarning extends Error {\n  constructor(coord: string, conflictingIndex: string) {\n    super(\n      `Covariant-derivative wrt='${coord}' collides with existing free index ` +\n      `'${conflictingIndex}' of the operand; result will silently misbehave. ` +\n      `Rename the operand's free index or pick a different wrt label.`,\n    );\n    this.name = 'DuplicateCoordinateWarning';\n    Object.setPrototypeOf(this, DuplicateCoordinateWarning.prototype);\n  }\n}\n"]}
+{"version":3,"file":"errors.js","sourceRoot":"","sources":["../../src/dimensional/errors.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAIH,4CAA4C;AAC5C,MAAM,OAAO,QAAS,SAAQ,KAAK;IACjC,YAAY,OAAe,EAAE,OAAsC;QACjE,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACxB,IAAI,CAAC,IAAI,GAAG,UAAU,CAAC;QACvB,2EAA2E;IAC7E,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,QAAQ,CAAY;IACpB,MAAM,CAAY;IAClC,YAAY,OAAe,EAAE,QAAmB,EAAE,MAAiB;QACjE,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;CACF;AAED;;;;;;;;;GASG;AACH,MAAM,OAAO,wBAAyB,SAAQ,QAAQ;IACpC,UAAU,CAAS;IACnB,KAAK,CAAS;IAC9B,YAAY,UAAkB,EAAE,KAAa;QAC3C,KAAK,CACH,WAAW,UAAU,+BAA+B,KAAK,KAAK;YAC9D,mEAAmE;YACnE,wCAAwC,CACzC,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;QACvC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACrB,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,wBAAyB,SAAQ,QAAQ;IACpC,KAAK,CAAS;IACd,UAAU,CAAS;IACnC;;;;OAIG;IACa,OAAO,CAAyB;IAChD,YACE,KAAa,EACb,UAAkB,EAClB,OAA+B;QAE/B,MAAM,YAAY,GAChB,OAAO,IAAI,OAAO,CAAC,MAAM,GAAG,CAAC;YAC3B,CAAC,CAAC,cAAc,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG;YACrC,CAAC,CAAC,EAAE,CAAC;QACT,KAAK,CACH,gBAAgB,KAAK,aAAa,UAAU,oBAAoB;YAChE,0BAA0B,YAAY,gCAAgC;YACtE,oDAAoD,CACrD,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;IACzB,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjC,KAAK,CAAS;IACd,QAAQ,CAAoB;IAC5C,YAAY,KAAa,EAAE,QAA2B;QACpD,KAAK,CACH,gBAAgB,KAAK,2CAA2C,QAAQ,KAAK;YAC7E,yDAAyD;YACzD,yDAAyD,KAAK,QAAQ;YACtE,sBAAsB,KAAK,yCAAyC;YACpE,qBAAqB,CACtB,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;IAC3B,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,qBAAsB,SAAQ,QAAQ;IACjC,EAAE,CAAS;IAC3B,YAAY,EAAU;QACpB,KAAK,CACH,aAAa,EAAE,oDAAoD;YACnE,wEAAwE,CACzE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;QACpC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClD,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;IACvC,CAAC;CACF;AAED;;;;;;GAMG;AACH,MAAM,OAAO,gCAAiC,SAAQ,QAAQ;IAC5D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,kCAAkC,CAAC;IACjD,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,UAAU,CAAS;IACnB,UAAU,CAAS;IACnC,YAAY,UAAkB,EAAE,UAAkB;QAChD,KAAK,CACH,kBAAkB,UAAU,cAAc,UAAU,UAAU;YAC9D,qEAAqE,CACtE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;IAC/B,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,oBAAqB,SAAQ,QAAQ;IAChC,UAAU,CAAS;IACnB,MAAM,CAAS;IAC/B,YAAY,UAAkB,EAAE,MAAc;QAC5C,KAAK,CACH,kBAAkB,UAAU,MAAM,MAAM,0BAA0B,CACnE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,sBAAsB,CAAC;QACnC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;QAC7B,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,OAAO,yBAA0B,SAAQ,QAAQ;IACrC,UAAU,CAAS;IACnC,YAAY,UAAkB;QAC5B,KAAK,CACH,4BAA4B,UAAU,uBAAuB;YAC7D,+DAA+D,CAChE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,2BAA2B,CAAC;QACxC,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;IAC/B,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,sBAAuB,SAAQ,QAAQ;IAClC,YAAY,CAAoB;IAChD,YAAY,YAA+B;QACzC,KAAK,CACH,+CAA+C,YAAY,KAAK;YAChE,wDAAwD;YACxD,wBAAwB,CACzB,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;QACrC,IAAI,CAAC,YAAY,GAAG,YAAY,CAAC;IACnC,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAM,OAAO,mCAAoC,SAAQ,QAAQ;IAC/C,KAAK,CAAS;IAC9B,YAAY,KAAa;QACvB,KAAK,CACH,6BAA6B,KAAK,8BAA8B;YAChE,iEAAiE;YACjE,oCAAoC,CACrC,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,qCAAqC,CAAC;QAClD,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACrB,CAAC;CACF;AAED;;;;;;;;;;;;;;GAcG;AACH,MAAM,OAAO,0BAA2B,SAAQ,KAAK;IACnD,YAAY,KAAa,EAAE,gBAAwB;QACjD,KAAK,CACH,6BAA6B,KAAK,sCAAsC;YACxE,IAAI,gBAAgB,oDAAoD;YACxE,gEAAgE,CACjE,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,4BAA4B,CAAC;IAC3C,CAAC;CACF","sourcesContent":["/**\n * UPT error hierarchy. All UPT-source errors subclass UPTError so\n * downstream consumers (mathjs, threejs, TensorJS) can discriminate\n * with `err instanceof UPTError`.\n *\n * @module dimensional/errors\n */\n\nimport { Dimension } from './types.js';\n\n/** Base class for all UPT-source errors. */\nexport class UPTError extends Error {\n  constructor(message: string, options?: { readonly cause?: unknown }) {\n    super(message, options);\n    this.name = 'UPTError';\n    // Restore prototype chain (needed for instanceof after ES5 transpilation).\n  }\n}\n\n/**\n * Thrown by add/sub and equation-level checks when dimensions disagree.\n * (Moved from src/dimensional/algebra.ts as part of the v0.2.0 UPTError\n * refactor. Behavior unchanged.)\n */\nexport class DimensionMismatchError extends UPTError {\n  public readonly expected: Dimension;\n  public readonly actual: Dimension;\n  constructor(message: string, expected: Dimension, actual: Dimension) {\n    super(message);\n    this.name = 'DimensionMismatchError';\n    this.expected = expected;\n    this.actual = actual;\n  }\n}\n\n/**\n * Thrown when a tensor-symbol's indices list contains the same label\n * more than once (e.g., T^μ_μ_μ). Per Part-VII §VII.4.\n *\n * NOTE: this error fires on DECLARATION-TIME duplicates within a single\n * tensor-symbol's indices list. Those indices are *free*, not \"dummy\"\n * (a dummy index in tensor calculus is one that is summed over). The\n * earlier name `RepeatedDummyLabelError` was a misnomer; that deprecated\n * alias was removed in v0.4.5.\n */\nexport class DuplicateIndexLabelError extends UPTError {\n  public readonly tensorName: string;\n  public readonly label: string;\n  constructor(tensorName: string, label: string) {\n    super(\n      `Tensor '${tensorName}' has repeated index label '${label}'. ` +\n      `A label may appear at most twice (once upper, once lower) within ` +\n      `a single tensor-symbol's indices list.`,\n    );\n    this.name = 'DuplicateIndexLabelError';\n    this.tensorName = tensorName;\n    this.label = label;\n  }\n}\n\n/**\n * Thrown when an index label appears more than twice across a\n * tensor-product's operands. Einstein contraction requires exactly two\n * occurrences (one upper + one lower) per dummy label. Per Part-VII §VII.4.\n */\nexport class IndexLabelCollisionError extends UPTError {\n  public readonly label: string;\n  public readonly totalCount: number;\n  /**\n   * Tensor names that contributed occurrences of the colliding label,\n   * if the caller has them available. May be undefined for callers that\n   * don't track per-operand provenance.\n   */\n  public readonly sources?: ReadonlyArray<string>;\n  constructor(\n    label: string,\n    totalCount: number,\n    sources?: ReadonlyArray<string>,\n  ) {\n    const sourceClause =\n      sources && sources.length > 0\n        ? ` (sources: ${sources.join(', ')})`\n        : '';\n    super(\n      `Index label '${label}' appears ${totalCount} times across the ` +\n      `tensor-product operands${sourceClause} (max allowed: 2 for Einstein ` +\n      `contraction). Rename one of the offending indices.`,\n    );\n    this.name = 'IndexLabelCollisionError';\n    this.label = label;\n    this.totalCount = totalCount;\n    this.sources = sources;\n  }\n}\n\n/**\n * Thrown when an index label appears twice in a tensor-product but both\n * occurrences share the same variance (both upper or both lower). v0.2.0\n * has no metric tensor available to raise/lower indices, so such a\n * contraction is rejected at validation time. Per Part-VII §VII.4.\n */\nexport class VarianceMismatchError extends UPTError {\n  public readonly label: string;\n  public readonly variance: 'upper' | 'lower';\n  constructor(label: string, variance: 'upper' | 'lower') {\n    super(\n      `Index label '${label}' appears twice but both with variance '${variance}'. ` +\n      `Einstein contraction requires one upper and one lower. ` +\n      `If this is intentional, use raise(operand, gInverse, '${label}') or ` +\n      `lower(operand, g, '${label}') to traverse variance via the metric ` +\n      `before contracting.`,\n    );\n    this.name = 'VarianceMismatchError';\n    this.label = label;\n    this.variance = variance;\n  }\n}\n\n/**\n * Thrown by `op '*'` / `'/'` / `'^'` when an operand carries non-empty\n * freeIndices. Scalar operators in v0.2.0 are strict: any tensor-valued\n * operand must instead pass through `tensor-product`. Per Part-VII §VII.5\n * (op-tensor boundary rules).\n */\nexport class TensorInScalarOpError extends UPTError {\n  public readonly op: string;\n  constructor(op: string) {\n    super(\n      `Operator '${op}' is scalar-only; received a tensor argument with ` +\n      `non-empty freeIndices. Use 'tensor-product' for tensor multiplication.`,\n    );\n    this.name = 'TensorInScalarOpError';\n    this.op = op;\n  }\n}\n\n/**\n * Thrown by `op '+'` / `'-'` when the operand freeIndices maps differ.\n * Addition / subtraction across tensors requires identical free-index\n * signatures (same labels, same per-label upper/lower counts). Per\n * Part-VII §VII.5.\n */\nexport class FreeIndexMismatchError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'FreeIndexMismatchError';\n  }\n}\n\n/**\n * Thrown by the validator's tensor-product child resolver when a non-tensor\n * operand to a `tensor-product` fails dimension inference. Subclassing\n * UPTError (rather than throwing a plain `Error`) preserves the §14.7\n * forward-compat invariant that downstream consumers can discriminate\n * UPT-source errors with `err instanceof UPTError`.\n */\nexport class TensorProductChildInferenceError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'TensorProductChildInferenceError';\n  }\n}\n\n/**\n * Thrown when a `metric-tensor` node has rank other than 2. Per\n * Part-VIII §VIII.2 (metric-tensor-rank-2-only TENSOR-RULE). The metric\n * is fundamentally a rank-2 object in v0.3.0; higher-rank generalizations\n * are out of scope.\n */\nexport class InvalidMetricRankError extends UPTError {\n  public readonly tensorName: string;\n  public readonly actualRank: number;\n  constructor(tensorName: string, actualRank: number) {\n    super(\n      `Metric tensor '${tensorName}' has rank ${actualRank}, but a ` +\n      `metric must be rank-2 (exactly two indices). Per Part-VIII §VIII.2.`,\n    );\n    this.name = 'InvalidMetricRankError';\n    this.tensorName = tensorName;\n    this.actualRank = actualRank;\n  }\n}\n\n/**\n * Thrown when a `metric-tensor` has mixed-variance indices, an empty\n * signature, an invalid signature string, or other structural problem\n * orthogonal to rank. Per Part-VIII §VIII.2 and §VIII.5 (raise/lower\n * variance requirements).\n */\nexport class MetricSignatureError extends UPTError {\n  public readonly tensorName: string;\n  public readonly reason: string;\n  constructor(tensorName: string, reason: string) {\n    super(\n      `Metric tensor '${tensorName}': ${reason}. Per Part-VIII §VIII.2.`,\n    );\n    this.name = 'MetricSignatureError';\n    this.tensorName = tensorName;\n    this.reason = reason;\n  }\n}\n\n/**\n * Thrown when a `kronecker-delta` node has rank other than 2.\n * Per Part-VIII §VIII.3.\n */\nexport class InvalidKroneckerRankError extends UPTError {\n  public readonly actualRank: number;\n  constructor(actualRank: number) {\n    super(\n      `Kronecker delta has rank ${actualRank}, but must be rank-2 ` +\n      `(exactly one upper + one lower index). Per Part-VIII §VIII.3.`,\n    );\n    this.name = 'InvalidKroneckerRankError';\n    this.actualRank = actualRank;\n  }\n}\n\n/**\n * Thrown when a `kronecker-delta`'s two indices share the same variance\n * (both upper or both lower). Per Part-VIII §VIII.3, the canonical\n * δ^μ_ν requires mixed variance.\n */\nexport class KroneckerVarianceError extends UPTError {\n  public readonly bothVariance: 'upper' | 'lower';\n  constructor(bothVariance: 'upper' | 'lower') {\n    super(\n      `Kronecker delta indices both have variance '${bothVariance}'; ` +\n      `must be one upper + one lower (canonical δ^μ_ν form). ` +\n      `Per Part-VIII §VIII.3.`,\n    );\n    this.name = 'KroneckerVarianceError';\n    this.bothVariance = bothVariance;\n  }\n}\n\n/**\n * Thrown when a `tensor-partial-derivative` node's `wrtIndex.variance`\n * is `'upper'`. The differentiation operator ∂/∂x^μ is fundamentally\n * covariant (lower) regardless of the variance of the coordinate. Per\n * Part-VIII §VIII.4 (pderiv-wrtIndex-always-lower TENSOR-RULE).\n */\nexport class PartialDerivativeIndexVarianceError extends UPTError {\n  public readonly label: string;\n  constructor(label: string) {\n    super(\n      `Partial-derivative index '${label}' has variance 'upper', but ` +\n      `∂_μ is fundamentally covariant — wrtIndex.variance must always ` +\n      `be 'lower'. Per Part-VIII §VIII.4.`,\n    );\n    this.name = 'PartialDerivativeIndexVarianceError';\n    this.label = label;\n  }\n}\n\n/**\n * Emitted (not thrown) via `process.emitWarning` when a covariant-derivative's\n * `wrt` coordinate label collides with an existing free index of the operand\n * and the env var `UPT_ALLOW_COORD_SHADOW=1` has opted in to downgraded\n * handling. By default the collision throws `MetricSignatureError` (soundness\n * over friendliness — silent wrong results are the worst failure mode).\n *\n * Lives here (dimensional/errors.ts) rather than numerical/index.ts to avoid\n * a dimensional → numerical import cycle; re-exported from numerical/index.ts\n * for the public API surface per v0.4.0-Implementation-Plan Task 13.\n *\n * A plain `class extends Error` under the ES2022 target: `instanceof` holds without a prototype reset.\n *\n * @public\n */\nexport class DuplicateCoordinateWarning extends Error {\n  constructor(coord: string, conflictingIndex: string) {\n    super(\n      `Covariant-derivative wrt='${coord}' collides with existing free index ` +\n      `'${conflictingIndex}' of the operand; result will silently misbehave. ` +\n      `Rename the operand's free index or pick a different wrt label.`,\n    );\n    this.name = 'DuplicateCoordinateWarning';\n  }\n}\n"]}

package/dist/dimensional/exact-scale.d.ts

@@ -12,9 +12,15 @@
  * through its shortest round-trip decimal (`String(x)`), which is the
  * literal a constant was written with.
  *
+ * The arithmetic is MathTS's (`@danielsimonjr/mathts-core`); the UPT
+ * {@link ExactScale} shape and function names are the thin layer over it, and
+ * a scale expression (`data/units.json`) reads through MathTS
+ * `readScaleExpression` under this module's UPT-facing signature.
+ *
  * @module dimensional/exact-scale
  * @internal
  */
+import { type ExactScale as CoreExactScale } from '@danielsimonjr/mathts-core';
 /** `num / den × irrational`, with `den > 0` and the fraction in lowest terms. @internal */
 export interface ExactScale {
     readonly num: bigint;
@@ -22,6 +28,10 @@ export interface ExactScale {
     /** The factor no rational states (π, ln 2). 1 for a rational scale. */
     readonly irrational: number;
 }
+/** UPT {@link ExactScale} to the MathTS core shape. @internal */
+export declare function coreExactScaleFromUpt(s: ExactScale): CoreExactScale;
+/** MathTS core exact scale to UPT {@link ExactScale}. @internal */
+export declare function uptExactScaleFromCore(s: CoreExactScale): ExactScale;
 /** The scale 1. @internal */
 export declare const UNIT_SCALE: ExactScale;
 /**
@@ -42,7 +52,7 @@ export declare function multiplyScales(...factors: readonly ExactScale[]): Exact
 export declare function divideScales(a: ExactScale, b: ExactScale): ExactScale;
 /**
  * `a ^ n`. An integer power stays exact. A fractional power has no exact
- * rational in general, so it becomes one floating factor.
+ * rational in general and throws.
  * @internal
  */
 export declare function powerScale(a: ExactScale, n: number): ExactScale;
@@ -51,13 +61,13 @@ export declare function addScales(a: ExactScale, b: ExactScale): ExactScale;
 /** The scale as a double: the rational rounded once, times the irrational factor. @internal */
 export declare function scaleToNumber(s: ExactScale): number;
 /**
- * The two sides of a unit text or a scale expression under the single-solidus
- * convention of ISO 80000-1: everything after the first `/` is the
- * denominator, so a chain is one denominator (`km/s/Mpc` is km/(s·Mpc), and
- * `W/m*K` is W/(m·K), not (W/m)·K). The denominator comes back as one
- * `*`-joined product, each `/`-part's outer parentheses dropped. `units.ts`
- * and {@link readScaleExpression} both split here, so a unit text and a scale
- * expression in `data/units.json` read a `/` the same way.
+ * The two sides of a unit text under the single-solidus convention of
+ * ISO 80000-1: everything after the first `/` is the denominator, so a chain
+ * is one denominator (`km/s/Mpc` is km/(s·Mpc), and `W/m*K` is W/(m·K), not
+ * (W/m)·K). The denominator comes back as one `*`-joined product, each
+ * `/`-part's outer parentheses dropped. `units.ts` splits here, so a unit text
+ * reads a `/` the same way a scale expression does. MathTS states the same
+ * rule inside its `readScaleExpression` and does not export the split.
  * @internal
  */
 export declare function solidusSides(text: string): {
@@ -70,9 +80,12 @@ export declare function solidusSides(text: string): {
  * each a decimal literal or a name with an optional integer power
  * (`4.1868*453.59237/1.8`, `pi/180`, `lbf/inch^2`). A `/` follows the unit
  * grammar ({@link solidusSides}): `a/b*c` is a/(b·c). `resolve` gives a name
- * its scale and returns undefined for a name it does not know, which throws
- * here. A decimal literal is read digit for digit, so the result is the
- * rational the text states.
+ * its scale and returns undefined for a name it does not know, which throws.
+ * A decimal literal is read digit for digit, so the result is the rational the
+ * text states.
+ *
+ * MathTS's `readScaleExpression` is the reader; this is the UPT-facing
+ * signature over it, the one adapter between its exact scale and UPT's.
  * @internal
  */
 export declare function readScaleExpression(text: string, resolve: (name: string) => ExactScale | undefined): ExactScale;

package/dist/dimensional/exact-scale.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"exact-scale.d.ts","sourceRoot":"","sources":["../../src/dimensional/exact-scale.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,2FAA2F;AAC3F,MAAM,WAAW,UAAU;IACzB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,uEAAuE;IACvE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAmBD,6BAA6B;AAC7B,eAAO,MAAM,UAAU,EAAE,UAAgD,CAAC;AAI1E;;;;GAIG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,UAAU,GAAG,IAAI,CAW5D;AAED,2EAA2E;AAC3E,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU,CAG7C;AAED,sCAAsC;AACtC,wBAAgB,UAAU,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,GAAG,UAAU,CAEjF;AAED,gEAAgE;AAChE,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU,CAGrD;AAED,6BAA6B;AAC7B,wBAAgB,cAAc,CAAC,GAAG,OAAO,EAAE,SAAS,UAAU,EAAE,GAAG,UAAU,CAU5E;AAED,yBAAyB;AACzB,wBAAgB,YAAY,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU,GAAG,UAAU,CAGrE;AAED;;;;GAIG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,GAAG,UAAU,CAO/D;AAED,+FAA+F;AAC/F,wBAAgB,SAAS,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU,GAAG,UAAU,CAGlE;AAwDD,+FAA+F;AAC/F,wBAAgB,aAAa,CAAC,CAAC,EAAE,UAAU,GAAG,MAAM,CAGnD;AAID;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG;IAAE,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAA;CAAE,CASvG;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,UAAU,GAAG,SAAS,GAAG,UAAU,CAiB/G"}
+{"version":3,"file":"exact-scale.d.ts","sourceRoot":"","sources":["../../src/dimensional/exact-scale.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,OAAO,EAYL,KAAK,UAAU,IAAI,cAAc,EAClC,MAAM,4BAA4B,CAAC;AAKpC,2FAA2F;AAC3F,MAAM,WAAW,UAAU;IACzB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,uEAAuE;IACvE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AA+BD,iEAAiE;AACjE,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,UAAU,GAAG,cAAc,CAEnE;AAED,mEAAmE;AACnE,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,cAAc,GAAG,UAAU,CAEnE;AAgDD,6BAA6B;AAC7B,eAAO,MAAM,UAAU,EAAE,UAAgD,CAAC;AAE1E;;;;GAIG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,UAAU,GAAG,IAAI,CAG5D;AAED,2EAA2E;AAC3E,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU,CAG7C;AAED,sCAAsC;AACtC,wBAAgB,UAAU,CAAC,GAAG,EAAE,MAAM,GAAG,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,GAAG,UAAU,CAEjF;AAED,gEAAgE;AAChE,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU,CAGrD;AAED,6BAA6B;AAC7B,wBAAgB,cAAc,CAAC,GAAG,OAAO,EAAE,SAAS,UAAU,EAAE,GAAG,UAAU,CAG5E;AAED,yBAAyB;AACzB,wBAAgB,YAAY,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU,GAAG,UAAU,CAErE;AAED;;;;GAIG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,GAAG,UAAU,CAG/D;AAED,+FAA+F;AAC/F,wBAAgB,SAAS,CAAC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU,GAAG,UAAU,CAIlE;AAED,+FAA+F;AAC/F,wBAAgB,aAAa,CAAC,CAAC,EAAE,UAAU,GAAG,MAAM,CAEnD;AAED;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG;IAAE,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAA;CAAE,CASvG;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,UAAU,GAAG,SAAS,GAAG,UAAU,CAO/G"}

package/dist/dimensional/exact-scale.js

@@ -12,15 +12,21 @@
  * through its shortest round-trip decimal (`String(x)`), which is the
  * literal a constant was written with.
  *
+ * The arithmetic is MathTS's (`@danielsimonjr/mathts-core`); the UPT
+ * {@link ExactScale} shape and function names are the thin layer over it, and
+ * a scale expression (`data/units.json`) reads through MathTS
+ * `readScaleExpression` under this module's UPT-facing signature.
+ *
  * @module dimensional/exact-scale
  * @internal
  */
-function gcd(a, b) {
-    let x = a < 0n ? -a : a;
-    let y = b < 0n ? -b : b;
-    while (y !== 0n)
-        [x, y] = [y, x % y];
-    return x;
+import { addExactScales, decimalExactScale, divideExactScales, exactScaleOf, exactScaleToNumber, LN2_EXACT_SCALE, multiplyExactScales, PI_EXACT_SCALE, powerExactScale, ratioExactScale, readScaleExpression as mathtsReadScaleExpression, } from '@danielsimonjr/mathts-core';
+import { gcd } from '@danielsimonjr/mathts-functions';
+import { bigintFrom } from '../core/mathts-read.js';
+/** The bigint gcd, MathTS `gcd` on its `bigint, bigint` signature, read back as the bigint it is. */
+function bigGcd(a, b) {
+    const g = bigintFrom('gcd', gcd(a, b));
+    return g < 0n ? -g : g;
 }
 function make(num, den, irrational) {
     if (den === 0n)
@@ -29,42 +35,104 @@ function make(num, den, irrational) {
         num = -num;
         den = -den;
     }
-    const g = gcd(num, den);
+    const g = bigGcd(num, den);
     return g > 1n ? { num: num / g, den: den / g, irrational } : { num, den, irrational };
 }
+function irrationalFactor(irrationals) {
+    let factor = 1;
+    const pi = irrationals.pi ?? 0;
+    const ln2 = irrationals.ln2 ?? 0;
+    if (pi !== 0)
+        factor *= Math.PI ** pi;
+    if (ln2 !== 0)
+        factor *= Math.LN2 ** ln2;
+    return factor;
+}
+function fromCore(s) {
+    return make(s.num, s.den, irrationalFactor(s.irrationals));
+}
+/** UPT {@link ExactScale} to the MathTS core shape. @internal */
+export function coreExactScaleFromUpt(s) {
+    return toCore(s);
+}
+/** MathTS core exact scale to UPT {@link ExactScale}. @internal */
+export function uptExactScaleFromCore(s) {
+    return fromCore(s);
+}
+/** When `x / base` is a nonzero integer, that integer; otherwise null. */
+function integerMultiple(x, base) {
+    if (!Number.isFinite(x) || !Number.isFinite(base) || base === 0 || x === 0)
+        return null;
+    const q = x / base;
+    if (!Number.isFinite(q))
+        return null;
+    const n = Math.trunc(q);
+    if (n === 0 || q !== n)
+        return null;
+    if (Math.abs(n * base - x) > Math.max(1e-12, Math.abs(x) * 1e-14))
+        return null;
+    return BigInt(n);
+}
+/**
+ * When `x` is `base` raised to a small integer power (π⁻¹ for a parsec-scale,
+ * ln 2 for a bit), that integer; otherwise null. The factor is compared
+ * exactly: `irrationalFactor` built it as `base ** n`, so the same expression
+ * reads it back.
+ */
+function irrationalExponent(x, base) {
+    if (!Number.isFinite(x) || x <= 0 || !Number.isFinite(base) || base <= 0)
+        return null;
+    for (let n = -8; n <= 8; n++) {
+        if (x === base ** n)
+            return n;
+    }
+    return null;
+}
+function toCore(s) {
+    let core = ratioExactScale(s.num, s.den);
+    const ir = s.irrational;
+    if (ir === 1)
+        return core;
+    const piPower = irrationalExponent(ir, Math.PI);
+    if (piPower !== null)
+        return multiplyExactScales(core, powerExactScale(PI_EXACT_SCALE, piPower));
+    const ln2Power = irrationalExponent(ir, Math.LN2);
+    if (ln2Power !== null)
+        return multiplyExactScales(core, powerExactScale(LN2_EXACT_SCALE, ln2Power));
+    const piMultiple = integerMultiple(ir, Math.PI);
+    if (piMultiple !== null) {
+        if (piMultiple !== 1n)
+            core = multiplyExactScales(core, ratioExactScale(piMultiple, 1n));
+        return multiplyExactScales(core, PI_EXACT_SCALE);
+    }
+    const ln2Multiple = integerMultiple(ir, Math.LN2);
+    if (ln2Multiple !== null) {
+        if (ln2Multiple !== 1n)
+            core = multiplyExactScales(core, ratioExactScale(ln2Multiple, 1n));
+        return multiplyExactScales(core, LN2_EXACT_SCALE);
+    }
+    return multiplyExactScales(core, exactScaleOf(ir));
+}
 /** The scale 1. @internal */
 export const UNIT_SCALE = { num: 1n, den: 1n, irrational: 1 };
-const DECIMAL = /^([+-]?)(\d*)(?:\.(\d*))?(?:[eE]([+-]?\d+))?$/;
 /**
  * The exact value of a decimal literal (`72`, `0.45359237`, `1.602176634e-19`),
  * or null when `text` is not one.
  * @internal
  */
 export function decimalScale(text) {
-    const m = DECIMAL.exec(text.trim());
-    if (m === null)
-        return null;
-    const whole = m[2] ?? '';
-    const fraction = m[3] ?? '';
-    if (whole === '' && fraction === '')
-        return null;
-    const exponent = Number(m[4] ?? '0') - fraction.length;
-    if (!Number.isSafeInteger(exponent) || Math.abs(exponent) > 10_000)
-        return null;
-    let num = BigInt(whole + fraction);
-    if (m[1] === '-')
-        num = -num;
-    return exponent >= 0 ? make(num * 10n ** BigInt(exponent), 1n, 1) : make(num, 10n ** BigInt(-exponent), 1);
+    const scaled = decimalExactScale(text.trim());
+    return scaled === null ? null : fromCore(scaled);
 }
 /** The exact scale a finite JavaScript number was written as. @internal */
 export function scaleOf(x) {
     if (!Number.isFinite(x))
         throw new RangeError(`exact scale: ${x} is not finite`);
-    return decimalScale(String(x));
+    return fromCore(exactScaleOf(x));
 }
 /** `num / den`, exactly. @internal */
 export function ratioScale(num, den) {
-    return make(BigInt(num), BigInt(den), 1);
+    return fromCore(ratioExactScale(num, den));
 }
 /** A factor no rational states, such as π or ln 2. @internal */
 export function irrationalScale(x) {
@@ -74,108 +142,43 @@ export function irrationalScale(x) {
 }
 /** `a × b × …`. @internal */
 export function multiplyScales(...factors) {
-    let num = 1n;
-    let den = 1n;
-    let irrational = 1;
-    for (const f of factors) {
-        num *= f.num;
-        den *= f.den;
-        irrational *= f.irrational;
-    }
-    return make(num, den, irrational);
+    if (factors.length === 0)
+        return UNIT_SCALE;
+    return fromCore(multiplyExactScales(...factors.map(toCore)));
 }
 /** `a / b`. @internal */
 export function divideScales(a, b) {
-    if (b.num === 0n)
-        throw new RangeError('exact scale: division by zero');
-    return make(a.num * b.den, a.den * b.num, a.irrational / b.irrational);
+    return fromCore(divideExactScales(toCore(a), toCore(b)));
 }
 /**
  * `a ^ n`. An integer power stays exact. A fractional power has no exact
- * rational in general, so it becomes one floating factor.
+ * rational in general and throws.
  * @internal
  */
 export function powerScale(a, n) {
-    if (Number.isInteger(n)) {
-        const k = BigInt(Math.abs(n));
-        const up = make(a.num ** k, a.den ** k, a.irrational ** Math.abs(n));
-        return n >= 0 ? up : divideScales(UNIT_SCALE, up);
-    }
-    return irrationalScale(scaleToNumber(a) ** n);
+    if (!Number.isInteger(n))
+        throw new RangeError(`exact scale: power ${n} is not an integer`);
+    return fromCore(powerExactScale(toCore(a), n));
 }
 /** `a + b`: exact when both are rational multiples of the same irrational factor. @internal */
 export function addScales(a, b) {
-    if (a.irrational === b.irrational)
-        return make(a.num * b.den + b.num * a.den, a.den * b.den, a.irrational);
-    return irrationalScale(scaleToNumber(a) + scaleToNumber(b));
-}
-/** The fewest significant digits kept before the sticky digit; far beyond a double's 17. */
-const MIN_DIGITS = 40;
-/** Decimal places of `1 / den` when it terminates (`den` = 2^a·5^b): max(a, b). Null otherwise. */
-function terminatingPlaces(den) {
-    let d = den;
-    let twos = 0;
-    let fives = 0;
-    while (d % 2n === 0n) {
-        d /= 2n;
-        twos++;
-    }
-    while (d % 5n === 0n) {
-        d /= 5n;
-        fives++;
-    }
-    return d === 1n ? Math.max(twos, fives) : null;
-}
-/**
- * The rational `num / den` (`den > 0`), reduced, rounded once to the nearest double.
- *
- * A terminating quotient (`den` = 2^a·5^b) is written out exactly, every digit,
- * and parsed: a tie between two doubles is a dyadic rational, so it always
- * lands here and rounds half to even. Any other quotient is not a tie, and it
- * is at least `1 / (den·2^54)` of its magnitude away from every tie; truncating
- * it to `max(40, digits(den) + 20)` significant digits stays closer than that,
- * and a sticky `1` puts the string strictly inside the truncation interval, so
- * the string rounds the way the quotient does for any denominator.
- */
-function rationalToNumber(num, den) {
-    if (num === 0n)
-        return 0;
-    const negative = num < 0n;
-    const g = gcd(num, den);
-    const a = (negative ? -num : num) / g;
-    den /= g;
-    const places = terminatingPlaces(den);
-    if (places !== null) {
-        const value = Number(`${(a * 10n ** BigInt(places)) / den}e-${places}`);
-        return negative ? -value : value;
-    }
-    const digits = Math.max(MIN_DIGITS, den.toString().length + 20);
-    const shift = digits - (a.toString().length - den.toString().length);
-    let q;
-    if (shift >= 0) {
-        q = (a * 10n ** BigInt(shift)) / den;
-    }
-    else {
-        q = a / (den * 10n ** BigInt(-shift));
-    }
-    // Not terminating, so the remainder is never zero: the sticky digit is always due.
-    const value = Number(`${q.toString()}1e${-(shift + 1)}`);
-    return negative ? -value : value;
+    const sum = addExactScales(toCore(a), toCore(b));
+    if (sum === null)
+        throw new RangeError('exact scale: unlike irrational factors');
+    return fromCore(sum);
 }
 /** The scale as a double: the rational rounded once, times the irrational factor. @internal */
 export function scaleToNumber(s) {
-    const rational = rationalToNumber(s.num, s.den);
-    return s.irrational === 1 ? rational : rational * s.irrational;
+    return exactScaleToNumber(toCore(s));
 }
-const FACTOR = /^([^*/^\s]+)(?:\^([+-]?\d+))?$/;
 /**
- * The two sides of a unit text or a scale expression under the single-solidus
- * convention of ISO 80000-1: everything after the first `/` is the
- * denominator, so a chain is one denominator (`km/s/Mpc` is km/(s·Mpc), and
- * `W/m*K` is W/(m·K), not (W/m)·K). The denominator comes back as one
- * `*`-joined product, each `/`-part's outer parentheses dropped. `units.ts`
- * and {@link readScaleExpression} both split here, so a unit text and a scale
- * expression in `data/units.json` read a `/` the same way.
+ * The two sides of a unit text under the single-solidus convention of
+ * ISO 80000-1: everything after the first `/` is the denominator, so a chain
+ * is one denominator (`km/s/Mpc` is km/(s·Mpc), and `W/m*K` is W/(m·K), not
+ * (W/m)·K). The denominator comes back as one `*`-joined product, each
+ * `/`-part's outer parentheses dropped. `units.ts` splits here, so a unit text
+ * reads a `/` the same way a scale expression does. MathTS states the same
+ * rule inside its `readScaleExpression` and does not export the split.
  * @internal
  */
 export function solidusSides(text) {
@@ -194,29 +197,18 @@ export function solidusSides(text) {
  * each a decimal literal or a name with an optional integer power
  * (`4.1868*453.59237/1.8`, `pi/180`, `lbf/inch^2`). A `/` follows the unit
  * grammar ({@link solidusSides}): `a/b*c` is a/(b·c). `resolve` gives a name
- * its scale and returns undefined for a name it does not know, which throws
- * here. A decimal literal is read digit for digit, so the result is the
- * rational the text states.
+ * its scale and returns undefined for a name it does not know, which throws.
+ * A decimal literal is read digit for digit, so the result is the rational the
+ * text states.
+ *
+ * MathTS's `readScaleExpression` is the reader; this is the UPT-facing
+ * signature over it, the one adapter between its exact scale and UPT's.
  * @internal
  */
 export function readScaleExpression(text, resolve) {
-    const product = (side) => {
-        let scale = UNIT_SCALE;
-        for (const part of side.split('*')) {
-            const m = FACTOR.exec(part);
-            if (m === null)
-                throw new RangeError(`exact scale: '${text}' is not a product of factors`);
-            const atom = m[1];
-            const base = decimalScale(atom) ?? resolve(atom);
-            if (base === undefined)
-                throw new RangeError(`exact scale: '${text}' names '${atom}', which is not a known scale`);
-            scale = multiplyScales(scale, m[2] === undefined ? base : powerScale(base, Number(m[2])));
-        }
-        return scale;
-    };
-    const trimmed = text.trim();
-    const { numerator, denominator } = solidusSides(trimmed);
-    const top = product(numerator);
-    return trimmed.includes('/') ? divideScales(top, product(denominator)) : top;
+    return uptExactScaleFromCore(mathtsReadScaleExpression(text.trim(), (name) => {
+        const scale = resolve(name);
+        return scale === undefined ? undefined : coreExactScaleFromUpt(scale);
+    }));
 }
 //# sourceMappingURL=exact-scale.js.map