mathjs 10.0.0 → 10.6.4

package/lib/esm/plain/number/arithmetic.js

@@ -1,4 +1,4 @@
-import { isInteger, log2, log10, cbrt, expm1, sign, toFixed, log1p } from '../../utils/number.js';
+import { cbrt, expm1, isInteger, log10, log1p, log2, sign, toFixed } from '../../utils/number.js';
 var n1 = 'number';
 var n2 = 'number, number';
 export function absNumber(a) {
@@ -33,10 +33,6 @@ export function cbrtNumber(x) {
   return cbrt(x);
 }
 cbrtNumber.signature = n1;
-export function ceilNumber(x) {
-  return Math.ceil(x);
-}
-ceilNumber.signature = n1;
 export function cubeNumber(x) {
   return x * x * x;
 }
@@ -49,14 +45,6 @@ export function expm1Number(x) {
   return expm1(x);
 }
 expm1Number.signature = n1;
-export function fixNumber(x) {
-  return x > 0 ? Math.floor(x) : Math.ceil(x);
-}
-fixNumber.signature = n1;
-export function floorNumber(x) {
-  return Math.floor(x);
-}
-floorNumber.signature = n1;
 /**
  * Calculate gcd for numbers
  * @param {number} a
@@ -112,15 +100,19 @@ export function lcmNumber(a, b) {
 }
 lcmNumber.signature = n2;
 /**
- * Calculate the logarithm of a value.
+ * Calculate the logarithm of a value, optionally to a given base.
  * @param {number} x
+ * @param {number | null | undefined} base
  * @return {number}
  */
 
-export function logNumber(x) {
+export function logNumber(x, y) {
+  if (y) {
+    return Math.log(x) / Math.log(y);
+  }
+
   return Math.log(x);
 }
-logNumber.signature = n1;
 /**
  * Calculate the 10-base logarithm of a number
  * @param {number} x
@@ -178,11 +170,12 @@ modNumber.signature = n2;
  * Calculate the nth root of a, solve x^root == a
  * http://rosettacode.org/wiki/Nth_root#JavaScript
  * @param {number} a
- * @param {number} root
+ * @param {number} [2] root
  * @private
  */
 
-export function nthRootNumber(a, root) {
+export function nthRootNumber(a) {
+  var root = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 2;
   var inv = root < 0;
 
   if (inv) {
@@ -232,7 +225,6 @@ export function nthRootNumber(a, root) {
    return inv ? 1 / x : x
   */
 }
-nthRootNumber.signature = n2;
 export function signNumber(x) {
   return sign(x);
 }
@@ -321,9 +313,13 @@ powNumber.signature = n2;
 
 export function roundNumber(value) {
   var decimals = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+
+  if (!isInteger(decimals) || decimals < 0 || decimals > 15) {
+    throw new Error('Number of decimals in function round must be an integer from 0 to 15 inclusive');
+  }
+
   return parseFloat(toFixed(value, decimals));
 }
-roundNumber.signature = n2;
 /**
  * Calculate the norm of a number, the absolute value.
  * @param {number} x

package/lib/esm/plain/number/combinations.js

@@ -14,14 +14,26 @@ export function combinationsNumber(n, k) {
   }
 
   var nMinusk = n - k;
-  var prodrange;
+  var answer = 1;
+  var firstnumerator = k < nMinusk ? nMinusk + 1 : k + 1;
+  var nextdivisor = 2;
+  var lastdivisor = k < nMinusk ? k : nMinusk; // balance multiplications and divisions to try to keep intermediate values
+  // in exact-integer range as long as possible
 
-  if (k < nMinusk) {
-    prodrange = product(nMinusk + 1, n);
-    return prodrange / product(1, k);
+  for (var nextnumerator = firstnumerator; nextnumerator <= n; ++nextnumerator) {
+    answer *= nextnumerator;
+
+    while (nextdivisor <= lastdivisor && answer % nextdivisor === 0) {
+      answer /= nextdivisor;
+      ++nextdivisor;
+    }
+  } // for big n, k, floating point may have caused weirdness in remainder
+
+
+  if (nextdivisor <= lastdivisor) {
+    answer /= product(nextdivisor, lastdivisor);
   }
 
-  prodrange = product(k + 1, n);
-  return prodrange / product(1, nMinusk);
+  return answer;
 }
 combinationsNumber.signature = 'number, number';

package/lib/esm/plain/number/probability.js

@@ -46,4 +46,36 @@ export function gammaNumber(n) {
 gammaNumber.signature = 'number'; // TODO: comment on the variables g and p
 
 export var gammaG = 4.7421875;
-export var gammaP = [0.99999999999999709182, 57.156235665862923517, -59.597960355475491248, 14.136097974741747174, -0.49191381609762019978, 0.33994649984811888699e-4, 0.46523628927048575665e-4, -0.98374475304879564677e-4, 0.15808870322491248884e-3, -0.21026444172410488319e-3, 0.21743961811521264320e-3, -0.16431810653676389022e-3, 0.84418223983852743293e-4, -0.26190838401581408670e-4, 0.36899182659531622704e-5];
+export var gammaP = [0.99999999999999709182, 57.156235665862923517, -59.597960355475491248, 14.136097974741747174, -0.49191381609762019978, 0.33994649984811888699e-4, 0.46523628927048575665e-4, -0.98374475304879564677e-4, 0.15808870322491248884e-3, -0.21026444172410488319e-3, 0.21743961811521264320e-3, -0.16431810653676389022e-3, 0.84418223983852743293e-4, -0.26190838401581408670e-4, 0.36899182659531622704e-5]; // lgamma implementation ref: https://mrob.com/pub/ries/lanczos-gamma.html#code
+// log(2 * pi) / 2
+
+export var lnSqrt2PI = 0.91893853320467274178;
+export var lgammaG = 5; // Lanczos parameter "g"
+
+export var lgammaN = 7; // Range of coefficients "n"
+
+export var lgammaSeries = [1.000000000190015, 76.18009172947146, -86.50532032941677, 24.01409824083091, -1.231739572450155, 0.1208650973866179e-2, -0.5395239384953e-5];
+export function lgammaNumber(n) {
+  if (n < 0) return NaN;
+  if (n === 0) return Infinity;
+  if (!isFinite(n)) return n;
+
+  if (n < 0.5) {
+    // Use Euler's reflection formula:
+    // gamma(z) = PI / (sin(PI * z) * gamma(1 - z))
+    return Math.log(Math.PI / Math.sin(Math.PI * n)) - lgammaNumber(1 - n);
+  } // Compute the logarithm of the Gamma function using the Lanczos method
+
+
+  n = n - 1;
+  var base = n + lgammaG + 0.5; // Base of the Lanczos exponential
+
+  var sum = lgammaSeries[0]; // We start with the terms that have the smallest coefficients and largest denominator
+
+  for (var i = lgammaN - 1; i >= 1; i--) {
+    sum += lgammaSeries[i] / (n + i);
+  }
+
+  return lnSqrt2PI + (n + 0.5) * Math.log(base) - base + Math.log(sum);
+}
+lgammaNumber.signature = 'number';

package/lib/esm/type/complex/Complex.js

@@ -138,8 +138,8 @@ export var createComplexClass = /* #__PURE__ */factory(name, dependencies, () =>
 
             if (isNumber(phi)) {
               return new Complex({
-                r: r,
-                phi: phi
+                r,
+                phi
               });
             }
 

package/lib/esm/type/fraction/function/fraction.js

@@ -9,27 +9,39 @@ export var createFraction = /* #__PURE__ */factory(name, dependencies, _ref => {
   } = _ref;
 
   /**
-   * Create a fraction convert a value to a fraction.
+   * Create a fraction or convert a value to a fraction.
+   *
+   * With one numeric argument, produces the closest rational approximation to the
+   * input.
+   * With two arguments, the first is the numerator and the second is the denominator,
+   * and creates the corresponding fraction. Both numerator and denominator must be
+   * integers.
+   * With one object argument, looks for the integer numerator as the value of property
+   * 'n' and the integer denominator as the value of property 'd'.
+   * With a matrix argument, creates a matrix of the same shape with entries
+   * converted into fractions.
    *
    * Syntax:
+   *     math.fraction(value)
    *     math.fraction(numerator, denominator)
    *     math.fraction({n: numerator, d: denominator})
-   *     math.fraction(matrix: Array | Matrix)         Turn all matrix entries
-   *                                                   into fractions
+   *     math.fraction(matrix: Array | Matrix)
    *
    * Examples:
    *
-   *     math.fraction(1, 3)
-   *     math.fraction('2/3')
-   *     math.fraction({n: 2, d: 3})
-   *     math.fraction([0.2, 0.25, 1.25])
+   *     math.fraction(6.283)             // returns Fraction 6283/1000
+   *     math.fraction(1, 3)              // returns Fraction 1/3
+   *     math.fraction('2/3')             // returns Fraction 2/3
+   *     math.fraction({n: 2, d: 3})      // returns Fraction 2/3
+   *     math.fraction([0.2, 0.25, 1.25]) // returns Array [1/5, 1/4, 5/4]
+   *     math.fraction(4, 5.1)            // throws Error: Parameters must be integer
    *
    * See also:
    *
    *    bignumber, number, string, unit
    *
    * @param {number | string | Fraction | BigNumber | Array | Matrix} [args]
-   *            Arguments specifying the numerator and denominator of
+   *            Arguments specifying the value, or numerator and denominator of
    *            the fraction
    * @return {Fraction | Array | Matrix} Returns a fraction
    */

package/lib/esm/type/matrix/DenseMatrix.js

@@ -793,7 +793,7 @@ export var createDenseMatrixClass = /* #__PURE__ */factory(name, dependencies, _
 
 
     return new DenseMatrix({
-      data: data,
+      data,
       size: [n],
       datatype: this._datatype
     });
@@ -918,7 +918,7 @@ export var createDenseMatrixClass = /* #__PURE__ */factory(name, dependencies, _
 
 
     return new DenseMatrix({
-      data: data,
+      data,
       size: [rows, columns]
     });
   };

package/lib/esm/type/matrix/FibonacciHeap.js

@@ -39,8 +39,8 @@ export var createFibonacciHeapClass = /* #__PURE__ */factory(name, dependencies,
   FibonacciHeap.prototype.insert = function (key, value) {
     // create node
     var node = {
-      key: key,
-      value: value,
+      key,
+      value,
       degree: 0
     }; // check we have a node in the minimum
 

package/lib/esm/type/matrix/SparseMatrix.js

@@ -350,10 +350,10 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     ptr.push(index.length); // return matrix
 
     return new SparseMatrix({
-      values: values,
-      index: index,
-      ptr: ptr,
-      size: size,
+      values,
+      index,
+      ptr,
+      size,
       datatype: matrix._datatype
     });
   }
@@ -417,23 +417,27 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
 
       if (!deepStrictEqual(iSize, sSize)) {
         throw new DimensionError(iSize, sSize, '>');
-      } // offsets
+      } // insert the sub matrix
 
 
-      var x0 = index.min()[0];
-      var y0 = index.min()[1]; // submatrix rows and columns
-
-      var m = sSize[0];
-      var n = sSize[1]; // loop submatrix
-
-      for (var x = 0; x < m; x++) {
-        // loop columns
-        for (var y = 0; y < n; y++) {
-          // value at i, j
-          var v = submatrix[x][y]; // invoke set (zero value will remove entry from matrix)
-
-          matrix.set([x + x0, y + y0], v, defaultValue);
-        }
+      if (iSize.length === 1) {
+        // if the replacement index only has 1 dimension, go trough each one and set its value
+        var range = index.dimension(0);
+        range.forEach(function (dataIndex, subIndex) {
+          validateIndex(dataIndex);
+          matrix.set([dataIndex, 0], submatrix[subIndex[0]], defaultValue);
+        });
+      } else {
+        // if the replacement index has 2 dimensions, go through each one and set the value in the correct index
+        var firstDimensionRange = index.dimension(0);
+        var secondDimensionRange = index.dimension(1);
+        firstDimensionRange.forEach(function (firstDataIndex, firstSubIndex) {
+          validateIndex(firstDataIndex);
+          secondDimensionRange.forEach(function (secondDataIndex, secondSubIndex) {
+            validateIndex(secondDataIndex);
+            matrix.set([firstDataIndex, secondDataIndex], submatrix[firstSubIndex[0]][secondSubIndex[0]], defaultValue);
+          });
+        });
       }
     }
 
@@ -1000,9 +1004,9 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     ptr.push(values.length); // return sparse matrix
 
     return new SparseMatrix({
-      values: values,
-      index: index,
-      ptr: ptr,
+      values,
+      index,
+      ptr,
       size: [maxRow - minRow + 1, maxColumn - minColumn + 1]
     });
   }
@@ -1269,9 +1273,9 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     ptr.push(values.length); // return matrix
 
     return new SparseMatrix({
-      values: values,
-      index: index,
-      ptr: ptr,
+      values,
+      index,
+      ptr,
       size: [n, 1]
     });
   };
@@ -1428,9 +1432,9 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     ptr.push(values.length); // create SparseMatrix
 
     return new SparseMatrix({
-      values: values,
-      index: index,
-      ptr: ptr,
+      values,
+      index,
+      ptr,
       size: [rows, columns]
     });
   };

package/lib/esm/type/matrix/function/sparse.js

@@ -11,6 +11,9 @@ export var createSparse = /* #__PURE__ */factory(name, dependencies, _ref => {
    * Create a Sparse Matrix. The function creates a new `math.Matrix` object from
    * an `Array`. A Matrix has utility functions to manipulate the data in the
    * matrix, like getting the size and getting or setting values in the matrix.
+   * Note that a Sparse Matrix is always 2-dimensional, so for example if
+   * you create one from a plain array of _n_ numbers, you get an _n_ by 1
+   * Sparse "column vector".
    *
    * Syntax:
    *
@@ -25,6 +28,9 @@ export var createSparse = /* #__PURE__ */factory(name, dependencies, _ref => {
    *    m.resize([3, 2], 5)
    *    m.valueOf()                     // Array [[1, 2], [3, 4], [5, 5]]
    *    m.get([1, 0])                    // number 3
+   *    let v = math.sparse([0, 0, 1])
+   *    v.size()                        // Array [3, 1]
+   *    v.get([2, 0])                   // number 1
    *
    * See also:
    *

package/lib/esm/type/unit/Unit.js

@@ -1,9 +1,9 @@
 import _extends from "@babel/runtime/helpers/extends";
 import _defineProperty from "@babel/runtime/helpers/defineProperty";
 
-function ownKeys(object, enumerableOnly) { var keys = Object.keys(object); if (Object.getOwnPropertySymbols) { var symbols = Object.getOwnPropertySymbols(object); if (enumerableOnly) { symbols = symbols.filter(function (sym) { return Object.getOwnPropertyDescriptor(object, sym).enumerable; }); } keys.push.apply(keys, symbols); } return keys; }
+function ownKeys(object, enumerableOnly) { var keys = Object.keys(object); if (Object.getOwnPropertySymbols) { var symbols = Object.getOwnPropertySymbols(object); enumerableOnly && (symbols = symbols.filter(function (sym) { return Object.getOwnPropertyDescriptor(object, sym).enumerable; })), keys.push.apply(keys, symbols); } return keys; }
 
-function _objectSpread(target) { for (var i = 1; i < arguments.length; i++) { var source = arguments[i] != null ? arguments[i] : {}; if (i % 2) { ownKeys(Object(source), true).forEach(function (key) { _defineProperty(target, key, source[key]); }); } else if (Object.getOwnPropertyDescriptors) { Object.defineProperties(target, Object.getOwnPropertyDescriptors(source)); } else { ownKeys(Object(source)).forEach(function (key) { Object.defineProperty(target, key, Object.getOwnPropertyDescriptor(source, key)); }); } } return target; }
+function _objectSpread(target) { for (var i = 1; i < arguments.length; i++) { var source = null != arguments[i] ? arguments[i] : {}; i % 2 ? ownKeys(Object(source), !0).forEach(function (key) { _defineProperty(target, key, source[key]); }) : Object.getOwnPropertyDescriptors ? Object.defineProperties(target, Object.getOwnPropertyDescriptors(source)) : ownKeys(Object(source)).forEach(function (key) { Object.defineProperty(target, key, Object.getOwnPropertyDescriptor(source, key)); }); } return target; }
 
 import { isComplex, isUnit, typeOf } from '../../utils/is.js';
 import { factory } from '../../utils/factory.js';
@@ -361,7 +361,7 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
       unit.units.push({
         unit: res.unit,
         prefix: res.prefix,
-        power: power
+        power
       });
 
       for (var i = 0; i < BASE_DIMENSIONS.length; i++) {
@@ -485,34 +485,23 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
 
 
   Unit.prototype._normalize = function (value) {
-    var unitValue, unitOffset, unitPower, unitPrefixValue;
-    var convert;
-
     if (value === null || value === undefined || this.units.length === 0) {
       return value;
-    } else if (this._isDerived()) {
-      // This is a derived unit, so do not apply offsets.
-      // For example, with J kg^-1 degC^-1 you would NOT want to apply the offset.
-      var res = value;
-      convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
-
-      for (var i = 0; i < this.units.length; i++) {
-        unitValue = convert(this.units[i].unit.value);
-        unitPrefixValue = convert(this.units[i].prefix.value);
-        unitPower = convert(this.units[i].power);
-        res = multiplyScalar(res, pow(multiplyScalar(unitValue, unitPrefixValue), unitPower));
-      }
+    }
+
+    var res = value;
+
+    var convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
 
-      return res;
-    } else {
-      // This is a single unit of power 1, like kg or degC
-      convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
 
-      unitValue = convert(this.units[0].unit.value);
-      unitOffset = convert(this.units[0].unit.offset);
-      unitPrefixValue = convert(this.units[0].prefix.value);
-      return multiplyScalar(addScalar(value, unitOffset), multiplyScalar(unitValue, unitPrefixValue));
+    for (var i = 0; i < this.units.length; i++) {
+      var unitValue = convert(this.units[i].unit.value);
+      var unitPrefixValue = convert(this.units[i].prefix.value);
+      var unitPower = convert(this.units[i].power);
+      res = multiplyScalar(res, pow(multiplyScalar(unitValue, unitPrefixValue), unitPower));
     }
+
+    return res;
   };
   /**
    * Denormalize a value, based on its currently set unit(s)
@@ -525,40 +514,23 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
 
 
   Unit.prototype._denormalize = function (value, prefixValue) {
-    var unitValue, unitOffset, unitPower, unitPrefixValue;
-    var convert;
-
     if (value === null || value === undefined || this.units.length === 0) {
       return value;
-    } else if (this._isDerived()) {
-      // This is a derived unit, so do not apply offsets.
-      // For example, with J kg^-1 degC^-1 you would NOT want to apply the offset.
-      // Also, prefixValue is ignored--but we will still use the prefix value stored in each unit, since kg is usually preferable to g unless the user decides otherwise.
-      var res = value;
-      convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
-
-      for (var i = 0; i < this.units.length; i++) {
-        unitValue = convert(this.units[i].unit.value);
-        unitPrefixValue = convert(this.units[i].prefix.value);
-        unitPower = convert(this.units[i].power);
-        res = divideScalar(res, pow(multiplyScalar(unitValue, unitPrefixValue), unitPower));
-      }
+    }
 
-      return res;
-    } else {
-      // This is a single unit of power 1, like kg or degC
-      convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
+    var res = value;
 
-      unitValue = convert(this.units[0].unit.value);
-      unitPrefixValue = convert(this.units[0].prefix.value);
-      unitOffset = convert(this.units[0].unit.offset);
+    var convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
 
-      if (prefixValue === undefined || prefixValue === null) {
-        return subtract(divideScalar(divideScalar(value, unitValue), unitPrefixValue), unitOffset);
-      } else {
-        return subtract(divideScalar(divideScalar(value, unitValue), prefixValue), unitOffset);
-      }
+
+    for (var i = 0; i < this.units.length; i++) {
+      var unitValue = convert(this.units[i].unit.value);
+      var unitPrefixValue = convert(this.units[i].prefix.value);
+      var unitPower = convert(this.units[i].power);
+      res = divideScalar(res, pow(multiplyScalar(unitValue, unitPrefixValue), unitPower));
     }
+
+    return res;
   };
   /**
    * Find a unit from a string
@@ -806,10 +778,24 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
 
 
   Unit.prototype.abs = function () {
-    // This gives correct, but unexpected, results for units with an offset.
-    // For example, abs(-283.15 degC) = -263.15 degC !!!
     var ret = this.clone();
-    ret.value = ret.value !== null ? abs(ret.value) : null;
+
+    if (ret.value !== null) {
+      if (ret._isDerived() || ret.units[0].unit.offset === 0) {
+        ret.value = abs(ret.value);
+      } else {
+        // To give the correct, but unexpected, results for units with an offset.
+        // For example, abs(-283.15 degC) = -263.15 degC !!!
+        // We must take the offset into consideration here
+        var convert = Unit._getNumberConverter(typeOf(ret.value)); // convert to Fraction or BigNumber if needed
+
+
+        var unitValue = convert(ret.units[0].unit.value);
+        var nominalOffset = convert(ret.units[0].unit.offset);
+        var unitOffset = multiplyScalar(unitValue, nominalOffset);
+        ret.value = subtract(abs(addScalar(ret.value, unitOffset)), unitOffset);
+      }
+    }
 
     for (var i in ret.units) {
       if (ret.units[i].unit.name === 'VA' || ret.units[i].unit.name === 'VAR') {
@@ -828,42 +814,44 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
 
 
   Unit.prototype.to = function (valuelessUnit) {
-    var other;
     var value = this.value === null ? this._normalize(1) : this.value;
+    var other;
 
     if (typeof valuelessUnit === 'string') {
-      // other = new Unit(null, valuelessUnit)
       other = Unit.parse(valuelessUnit);
+    } else if (isUnit(valuelessUnit)) {
+      other = valuelessUnit.clone();
+    } else {
+      throw new Error('String or Unit expected as parameter');
+    }
 
-      if (!this.equalBase(other)) {
-        throw new Error("Units do not match ('".concat(other.toString(), "' != '").concat(this.toString(), "')"));
-      }
+    if (!this.equalBase(other)) {
+      throw new Error("Units do not match ('".concat(other.toString(), "' != '").concat(this.toString(), "')"));
+    }
 
-      if (other.value !== null) {
-        throw new Error('Cannot convert to a unit with a value');
-      }
+    if (other.value !== null) {
+      throw new Error('Cannot convert to a unit with a value');
+    }
 
+    if (this.value === null || this._isDerived() || this.units[0].unit.offset === other.units[0].unit.offset) {
       other.value = clone(value);
-      other.fixPrefix = true;
-      other.skipAutomaticSimplification = true;
-      return other;
-    } else if (isUnit(valuelessUnit)) {
-      if (!this.equalBase(valuelessUnit)) {
-        throw new Error("Units do not match ('".concat(valuelessUnit.toString(), "' != '").concat(this.toString(), "')"));
-      }
+    } else {
+      /* Need to adjust value by difference in offset to convert */
+      var convert = Unit._getNumberConverter(typeOf(value)); // convert to Fraction or BigNumber if needed
 
-      if (valuelessUnit.value !== null) {
-        throw new Error('Cannot convert to a unit with a value');
-      }
 
-      other = valuelessUnit.clone();
-      other.value = clone(value);
-      other.fixPrefix = true;
-      other.skipAutomaticSimplification = true;
-      return other;
-    } else {
-      throw new Error('String or Unit expected as parameter');
+      var thisUnitValue = convert(this.units[0].unit.value);
+      var thisNominalOffset = convert(this.units[0].unit.offset);
+      var thisUnitOffset = multiplyScalar(thisUnitValue, thisNominalOffset);
+      var otherUnitValue = convert(other.units[0].unit.value);
+      var otherNominalOffset = convert(other.units[0].unit.offset);
+      var otherUnitOffset = multiplyScalar(otherUnitValue, otherNominalOffset);
+      other.value = subtract(addScalar(value, thisUnitOffset), otherUnitOffset);
     }
+
+    other.fixPrefix = true;
+    other.skipAutomaticSimplification = true;
+    return other;
   };
   /**
    * Return the value of the unit when represented with given valueless unit
@@ -3891,11 +3879,11 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
       newBaseUnit.key = baseName;
       BASE_UNITS[baseName] = newBaseUnit;
       newUnit = {
-        name: name,
+        name,
         value: 1,
         dimensions: BASE_UNITS[baseName].dimensions.slice(0),
-        prefixes: prefixes,
-        offset: offset,
+        prefixes,
+        offset,
         base: BASE_UNITS[baseName]
       };
       currentUnitSystem[baseName] = {
@@ -3904,11 +3892,11 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
       };
     } else {
       newUnit = {
-        name: name,
+        name,
         value: defUnit.value,
         dimensions: defUnit.dimensions.slice(0),
-        prefixes: prefixes,
-        offset: offset
+        prefixes,
+        offset
       }; // Create a new base if no matching base exists
 
       var anyMatch = false;

package/lib/esm/utils/collection.js

@@ -1,6 +1,7 @@
 import { isCollection, isMatrix } from './is.js';
 import { IndexError } from '../error/IndexError.js';
 import { arraySize } from './array.js';
+import { _switch } from './switch.js';
 /**
  * Test whether an array contains collections
  * @param {Array} array
@@ -127,32 +128,6 @@ function _reduce(mat, dim, callback) {
 
     return ret;
   }
-}
-/**
- * Transpose a matrix
- * @param {Array} mat
- * @returns {Array} ret
- * @private
- */
-
-
-function _switch(mat) {
-  var I = mat.length;
-  var J = mat[0].length;
-  var i, j;
-  var ret = [];
-
-  for (j = 0; j < J; j++) {
-    var tmp = [];
-
-    for (i = 0; i < I; i++) {
-      tmp.push(mat[i][j]);
-    }
-
-    ret.push(tmp);
-  }
-
-  return ret;
 } // TODO: document function scatter