mathjs 11.11.2 → 12.0.0

package/lib/esm/expression/parse.js

@@ -155,6 +155,19 @@ export var createParse = /* #__PURE__ */factory(name, dependencies, _ref => {
     undefined
   };
   var NUMERIC_CONSTANTS = ['NaN', 'Infinity'];
+  var ESCAPE_CHARACTERS = {
+    '"': '"',
+    "'": "'",
+    '\\': '\\',
+    '/': '/',
+    b: '\b',
+    f: '\f',
+    n: '\n',
+    r: '\r',
+    t: '\t'
+    // note that \u is handled separately in parseStringToken()
+  };
+
   function initialState() {
     return {
       extraNodes: {},
@@ -1227,7 +1240,7 @@ export var createParse = /* #__PURE__ */factory(name, dependencies, _ref => {
       node = parseAccessors(state, node);
       return node;
     }
-    return parseDoubleQuotesString(state);
+    return parseString(state);
   }
 
   /**
@@ -1313,61 +1326,14 @@ export var createParse = /* #__PURE__ */factory(name, dependencies, _ref => {
   }
 
   /**
-   * Parse a double quotes string.
+   * Parse a single or double quoted string.
    * @return {Node} node
    * @private
    */
-  function parseDoubleQuotesString(state) {
+  function parseString(state) {
     var node, str;
-    if (state.token === '"') {
-      str = parseDoubleQuotesStringToken(state);
-
-      // create constant
-      node = new ConstantNode(str);
-
-      // parse index parameters
-      node = parseAccessors(state, node);
-      return node;
-    }
-    return parseSingleQuotesString(state);
-  }
-
-  /**
-   * Parse a string surrounded by double quotes "..."
-   * @return {string}
-   */
-  function parseDoubleQuotesStringToken(state) {
-    var str = '';
-    while (currentCharacter(state) !== '' && currentCharacter(state) !== '"') {
-      if (currentCharacter(state) === '\\') {
-        // escape character, immediately process the next
-        // character to prevent stopping at a next '\"'
-        var cNext = nextCharacter(state);
-        if (cNext !== "'") {
-          str += currentCharacter(state);
-        }
-        next(state);
-      }
-      str += currentCharacter(state);
-      next(state);
-    }
-    getToken(state);
-    if (state.token !== '"') {
-      throw createSyntaxError(state, 'End of string " expected');
-    }
-    getToken(state);
-    return JSON.parse('"' + str + '"'); // unescape escaped characters
-  }
-
-  /**
-   * Parse a single quotes string.
-   * @return {Node} node
-   * @private
-   */
-  function parseSingleQuotesString(state) {
-    var node, str;
-    if (state.token === '\'') {
-      str = parseSingleQuotesStringToken(state);
+    if (state.token === '"' || state.token === "'") {
+      str = parseStringToken(state, state.token);
 
       // create constant
       node = new ConstantNode(str);
@@ -1380,30 +1346,47 @@ export var createParse = /* #__PURE__ */factory(name, dependencies, _ref => {
   }
 
   /**
-   * Parse a string surrounded by single quotes '...'
+   * Parse a string surrounded by single or double quotes
+   * @param {Object} state
+   * @param {"'" | "\""} quote
    * @return {string}
    */
-  function parseSingleQuotesStringToken(state) {
+  function parseStringToken(state, quote) {
     var str = '';
-    while (currentCharacter(state) !== '' && currentCharacter(state) !== '\'') {
+    while (currentCharacter(state) !== '' && currentCharacter(state) !== quote) {
       if (currentCharacter(state) === '\\') {
-        // escape character, immediately process the next
-        // character to prevent stopping at a next '\''
-        var cNext = nextCharacter(state);
-        if (cNext !== "'" && cNext !== '"') {
-          str += currentCharacter(state);
+        next(state);
+        var char = currentCharacter(state);
+        var escapeChar = ESCAPE_CHARACTERS[char];
+        if (escapeChar !== undefined) {
+          // an escaped control character like \" or \n
+          str += escapeChar;
+          state.index += 1;
+        } else if (char === 'u') {
+          // escaped unicode character
+          var unicode = state.expression.slice(state.index + 1, state.index + 5);
+          if (/^[0-9A-Fa-f]{4}$/.test(unicode)) {
+            // test whether the string holds four hexadecimal values
+            str += String.fromCharCode(parseInt(unicode, 16));
+            state.index += 5;
+          } else {
+            throw createSyntaxError(state, "Invalid unicode character \\u".concat(unicode));
+          }
+        } else {
+          throw createSyntaxError(state, "Bad escape character \\".concat(char));
         }
+      } else {
+        // any regular character
+        str += currentCharacter(state);
         next(state);
       }
-      str += currentCharacter(state);
-      next(state);
     }
     getToken(state);
-    if (state.token !== '\'') {
-      throw createSyntaxError(state, 'End of string \' expected');
+    if (state.token !== quote) {
+      throw createSyntaxError(state, "End of string ".concat(quote, " expected"));
     }
     getToken(state);
-    return JSON.parse('"' + str.replace(/"/g, '\\"') + '"'); // unescape escaped characters
+    return str;
   }
 
   /**
@@ -1498,10 +1481,8 @@ export var createParse = /* #__PURE__ */factory(name, dependencies, _ref => {
         getToken(state);
         if (state.token !== '}') {
           // parse key
-          if (state.token === '"') {
-            key = parseDoubleQuotesStringToken(state);
-          } else if (state.token === '\'') {
-            key = parseSingleQuotesStringToken(state);
+          if (state.token === '"' || state.token === "'") {
+            key = parseStringToken(state, state.token);
           } else if (state.tokenType === TOKENTYPE.SYMBOL || state.tokenType === TOKENTYPE.DELIMITER && state.token in NAMED_DELIMITERS) {
             key = state.token;
             getToken(state);

package/lib/esm/factoriesAny.js

@@ -33,6 +33,7 @@ export { createUnaryPlus } from './function/arithmetic/unaryPlus.js';
 export { createAbs } from './function/arithmetic/abs.js';
 export { createApply } from './function/matrix/apply.js';
 export { createAddScalar } from './function/arithmetic/addScalar.js';
+export { createSubtractScalar } from './function/arithmetic/subtractScalar.js';
 export { createCbrt } from './function/arithmetic/cbrt.js';
 export { createCeil } from './function/arithmetic/ceil.js';
 export { createCube } from './function/arithmetic/cube.js';

package/lib/esm/factoriesNumber.js

@@ -29,6 +29,7 @@ export var createUnaryMinus = /* #__PURE__ */createNumberFactory('unaryMinus', u
 export var createUnaryPlus = /* #__PURE__ */createNumberFactory('unaryPlus', unaryPlusNumber);
 export var createAbs = /* #__PURE__ */createNumberFactory('abs', absNumber);
 export var createAddScalar = /* #__PURE__ */createNumberFactory('addScalar', addNumber);
+export var createSubtractScalar = /* #__PURE__ */createNumberFactory('subtractScalar', subtractNumber);
 export var createCbrt = /* #__PURE__ */createNumberFactory('cbrt', cbrtNumber);
 export { createCeilNumber as createCeil } from './function/arithmetic/ceil.js';
 export var createCube = /* #__PURE__ */createNumberFactory('cube', cubeNumber);

package/lib/esm/function/algebra/decomposition/lup.js

@@ -1,7 +1,7 @@
 import { clone } from '../../../utils/object.js';
 import { factory } from '../../../utils/factory.js';
 var name = 'lup';
-var dependencies = ['typed', 'matrix', 'abs', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtract', 'larger', 'equalScalar', 'unaryMinus', 'DenseMatrix', 'SparseMatrix', 'Spa'];
+var dependencies = ['typed', 'matrix', 'abs', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'larger', 'equalScalar', 'unaryMinus', 'DenseMatrix', 'SparseMatrix', 'Spa'];
 export var createLup = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
@@ -10,7 +10,7 @@ export var createLup = /* #__PURE__ */factory(name, dependencies, _ref => {
     addScalar,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     larger,
     equalScalar,
     unaryMinus,
@@ -100,7 +100,7 @@ export var createLup = /* #__PURE__ */factory(name, dependencies, _ref => {
             // s = l[i, k] - data[k, j]
             s = addScalar(s, multiplyScalar(data[i][k], data[k][j]));
           }
-          data[i][j] = subtract(data[i][j], s);
+          data[i][j] = subtractScalar(data[i][j], s);
         }
       }
       // row with larger value in cvector, row >= j

package/lib/esm/function/algebra/decomposition/qr.js

@@ -1,7 +1,7 @@
 import _extends from "@babel/runtime/helpers/extends";
 import { factory } from '../../../utils/factory.js';
 var name = 'qr';
-var dependencies = ['typed', 'matrix', 'zeros', 'identity', 'isZero', 'equal', 'sign', 'sqrt', 'conj', 'unaryMinus', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtract', 'complex'];
+var dependencies = ['typed', 'matrix', 'zeros', 'identity', 'isZero', 'equal', 'sign', 'sqrt', 'conj', 'unaryMinus', 'addScalar', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'complex'];
 export var createQr = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
@@ -17,7 +17,7 @@ export var createQr = /* #__PURE__ */factory(name, dependencies, _ref => {
     addScalar,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     complex
   } = _ref;
   /**
@@ -132,7 +132,7 @@ export var createQr = /* #__PURE__ */factory(name, dependencies, _ref => {
       var alpha = multiplyScalar(sgn, sqrt(alphaSquared));
       if (!isZero(alpha)) {
         // first element in vector u
-        var u1 = subtract(pivot, alpha);
+        var u1 = subtractScalar(pivot, alpha);
 
         // w = v * u1 / |u|    (only elements k to (rows-1) are used)
         w[k] = 1;
@@ -167,7 +167,7 @@ export var createQr = /* #__PURE__ */factory(name, dependencies, _ref => {
           // calculate the jth element of [tau * transpose(w) * R]
           s = multiplyScalar(s, tau);
           for (i = k; i < rows; i++) {
-            Rdata[i][j] = multiplyScalar(subtract(Rdata[i][j], multiplyScalar(w[i], s)), conjSgn);
+            Rdata[i][j] = multiplyScalar(subtractScalar(Rdata[i][j], multiplyScalar(w[i], s)), conjSgn);
           }
         }
         /*
@@ -188,7 +188,7 @@ export var createQr = /* #__PURE__ */factory(name, dependencies, _ref => {
           // calculate the ith element of [tau * Q * w]
           s = multiplyScalar(s, tau);
           for (j = k; j < rows; ++j) {
-            Qdata[i][j] = divideScalar(subtract(Qdata[i][j], multiplyScalar(s, conj(w[j]))), conjSgn);
+            Qdata[i][j] = divideScalar(subtractScalar(Qdata[i][j], multiplyScalar(s, conj(w[j]))), conjSgn);
           }
         }
       }

package/lib/esm/function/algebra/solver/lsolve.js

@@ -1,14 +1,14 @@
 import { factory } from '../../../utils/factory.js';
 import { createSolveValidation } from './utils/solveValidation.js';
 var name = 'lsolve';
-var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtract', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
 export var createLsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     equalScalar,
     DenseMatrix
   } = _ref;
@@ -79,7 +79,7 @@ export var createLsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
 
         // loop rows
         for (var i = j + 1; i < rows; i++) {
-          bdata[i] = [subtract(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
+          bdata[i] = [subtractScalar(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
         }
       } else {
         // degenerate row, we can choose any value
@@ -141,7 +141,7 @@ export var createLsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
         var xj = divideScalar(bj, vjj);
         for (var _k = 0, l = jIndices.length; _k < l; _k++) {
           var _i = jIndices[_k];
-          bdata[_i] = [subtract(bdata[_i][0] || 0, multiplyScalar(xj, jValues[_k]))];
+          bdata[_i] = [subtractScalar(bdata[_i][0] || 0, multiplyScalar(xj, jValues[_k]))];
         }
         x[j] = [xj];
       } else {

package/lib/esm/function/algebra/solver/lsolveAll.js

@@ -1,14 +1,14 @@
 import { factory } from '../../../utils/factory.js';
 import { createSolveValidation } from './utils/solveValidation.js';
 var name = 'lsolveAll';
-var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtract', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
 export var createLsolveAll = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     equalScalar,
     DenseMatrix
   } = _ref;
@@ -76,7 +76,7 @@ export var createLsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           b[i] = divideScalar(b[i], M[i][i]);
           for (var j = i + 1; j < columns; j++) {
             // b[j] -= b[i] * M[j,i]
-            b[j] = subtract(b[j], multiplyScalar(b[i], M[j][i]));
+            b[j] = subtractScalar(b[j], multiplyScalar(b[i], M[j][i]));
           }
         } else if (!equalScalar(b[i], 0)) {
           // singular row, nonzero RHS
@@ -96,7 +96,7 @@ export var createLsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           var bNew = [...b];
           bNew[i] = 1;
           for (var _j = i + 1; _j < columns; _j++) {
-            bNew[_j] = subtract(bNew[_j], M[_j][i]);
+            bNew[_j] = subtractScalar(bNew[_j], M[_j][i]);
           }
           B.push(bNew);
         }
@@ -151,7 +151,7 @@ export var createLsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           b[i] = divideScalar(b[i], Mii);
           for (var _j2 = 0, _lastIndex = iIndices.length; _j2 < _lastIndex; _j2++) {
             var _J = iIndices[_j2];
-            b[_J] = subtract(b[_J], multiplyScalar(b[i], iValues[_j2]));
+            b[_J] = subtractScalar(b[_J], multiplyScalar(b[i], iValues[_j2]));
           }
         } else if (!equalScalar(b[i], 0)) {
           // singular row, nonzero RHS
@@ -172,7 +172,7 @@ export var createLsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           bNew[i] = 1;
           for (var _j3 = 0, _lastIndex2 = iIndices.length; _j3 < _lastIndex2; _j3++) {
             var _J2 = iIndices[_j3];
-            bNew[_J2] = subtract(bNew[_J2], iValues[_j3]);
+            bNew[_J2] = subtractScalar(bNew[_J2], iValues[_j3]);
           }
           B.push(bNew);
         }

package/lib/esm/function/algebra/solver/usolve.js

@@ -1,14 +1,14 @@
 import { factory } from '../../../utils/factory.js';
 import { createSolveValidation } from './utils/solveValidation.js';
 var name = 'usolve';
-var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtract', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
 export var createUsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     equalScalar,
     DenseMatrix
   } = _ref;
@@ -81,7 +81,7 @@ export var createUsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
         // loop rows
         for (var i = j - 1; i >= 0; i--) {
           // update copy of b
-          bdata[i] = [subtract(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
+          bdata[i] = [subtractScalar(bdata[i][0] || 0, multiplyScalar(xj, mdata[i][j]))];
         }
       } else {
         // zero value at j
@@ -145,7 +145,7 @@ export var createUsolve = /* #__PURE__ */factory(name, dependencies, _ref => {
         var xj = divideScalar(bj, vjj);
         for (var _k = 0, _lastIndex = jIndices.length; _k < _lastIndex; _k++) {
           var _i = jIndices[_k];
-          bdata[_i] = [subtract(bdata[_i][0], multiplyScalar(xj, jValues[_k]))];
+          bdata[_i] = [subtractScalar(bdata[_i][0], multiplyScalar(xj, jValues[_k]))];
         }
         x[j] = [xj];
       } else {

package/lib/esm/function/algebra/solver/usolveAll.js

@@ -1,14 +1,14 @@
 import { factory } from '../../../utils/factory.js';
 import { createSolveValidation } from './utils/solveValidation.js';
 var name = 'usolveAll';
-var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtract', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'divideScalar', 'multiplyScalar', 'subtractScalar', 'equalScalar', 'DenseMatrix'];
 export var createUsolveAll = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     divideScalar,
     multiplyScalar,
-    subtract,
+    subtractScalar,
     equalScalar,
     DenseMatrix
   } = _ref;
@@ -76,7 +76,7 @@ export var createUsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           b[i] = divideScalar(b[i], M[i][i]);
           for (var j = i - 1; j >= 0; j--) {
             // b[j] -= b[i] * M[j,i]
-            b[j] = subtract(b[j], multiplyScalar(b[i], M[j][i]));
+            b[j] = subtractScalar(b[j], multiplyScalar(b[i], M[j][i]));
           }
         } else if (!equalScalar(b[i], 0)) {
           // singular row, nonzero RHS
@@ -96,7 +96,7 @@ export var createUsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           var bNew = [...b];
           bNew[i] = 1;
           for (var _j = i - 1; _j >= 0; _j--) {
-            bNew[_j] = subtract(bNew[_j], M[_j][i]);
+            bNew[_j] = subtractScalar(bNew[_j], M[_j][i]);
           }
           B.push(bNew);
         }
@@ -153,7 +153,7 @@ export var createUsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           // loop upper triangular
           for (var _j2 = 0, _lastIndex = iIndices.length; _j2 < _lastIndex; _j2++) {
             var _J = iIndices[_j2];
-            b[_J] = subtract(b[_J], multiplyScalar(b[i], iValues[_j2]));
+            b[_J] = subtractScalar(b[_J], multiplyScalar(b[i], iValues[_j2]));
           }
         } else if (!equalScalar(b[i], 0)) {
           // singular row, nonzero RHS
@@ -176,7 +176,7 @@ export var createUsolveAll = /* #__PURE__ */factory(name, dependencies, _ref =>
           // loop upper triangular
           for (var _j3 = 0, _lastIndex2 = iIndices.length; _j3 < _lastIndex2; _j3++) {
             var _J2 = iIndices[_j3];
-            bNew[_J2] = subtract(bNew[_J2], iValues[_j3]);
+            bNew[_J2] = subtractScalar(bNew[_J2], iValues[_j3]);
           }
           B.push(bNew);
         }

package/lib/esm/function/arithmetic/subtract.js

@@ -6,13 +6,13 @@ import { createMatAlgo10xSids } from '../../type/matrix/utils/matAlgo10xSids.js'
 import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'subtract';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'addScalar', 'unaryMinus', 'DenseMatrix', 'concat'];
+var dependencies = ['typed', 'matrix', 'equalScalar', 'subtractScalar', 'unaryMinus', 'DenseMatrix', 'concat'];
 export var createSubtract = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     equalScalar,
-    addScalar,
+    subtractScalar,
     unaryMinus,
     DenseMatrix,
     concat
@@ -69,34 +69,14 @@ export var createSubtract = /* #__PURE__ */factory(name, dependencies, _ref => {
    *
    *    add
    *
-   * @param  {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} x
-   *            Initial value
-   * @param  {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} y
-   *            Value to subtract from `x`
-   * @return {number | BigNumber | Fraction | Complex | Unit | Array | Matrix}
-   *            Subtraction of `x` and `y`
+   * @param  {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} x Initial value
+   * @param  {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} y Value to subtract from `x`
+   * @return {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} Subtraction of `x` and `y`
    */
   return typed(name, {
-    'number, number': (x, y) => x - y,
-    'Complex, Complex': (x, y) => x.sub(y),
-    'BigNumber, BigNumber': (x, y) => x.minus(y),
-    'Fraction, Fraction': (x, y) => x.sub(y),
-    'Unit, Unit': typed.referToSelf(self => (x, y) => {
-      if (x.value === null) {
-        throw new Error('Parameter x contains a unit with undefined value');
-      }
-      if (y.value === null) {
-        throw new Error('Parameter y contains a unit with undefined value');
-      }
-      if (!x.equalBase(y)) {
-        throw new Error('Units do not match');
-      }
-      var res = x.clone();
-      res.value = typed.find(self, [res.valueType(), y.valueType()])(res.value, y.value);
-      res.fixPrefix = false;
-      return res;
-    })
+    'any, any': subtractScalar
   }, matrixAlgorithmSuite({
+    elop: subtractScalar,
     SS: matAlgo05xSfSf,
     DS: matAlgo01xDSid,
     SD: matAlgo03xDSf,

package/lib/esm/function/arithmetic/subtractScalar.js

@@ -0,0 +1,46 @@
+import { factory } from '../../utils/factory.js';
+import { subtractNumber } from '../../plain/number/index.js';
+var name = 'subtractScalar';
+var dependencies = ['typed'];
+export var createSubtractScalar = /* #__PURE__ */factory(name, dependencies, _ref => {
+  var {
+    typed
+  } = _ref;
+  /**
+   * Subtract two scalar values, `x - y`.
+   * This function is meant for internal use: it is used by the public function
+   * `subtract`
+   *
+   * This function does not support collections (Array or Matrix).
+   *
+   * @param  {number | BigNumber | Fraction | Complex | Unit} x   First value
+   * @param  {number | BigNumber | Fraction | Complex} y          Second value to be subtracted from `x`
+   * @return {number | BigNumber | Fraction | Complex | Unit}     Difference of `x` and `y`
+   * @private
+   */
+  return typed(name, {
+    'number, number': subtractNumber,
+    'Complex, Complex': function ComplexComplex(x, y) {
+      return x.sub(y);
+    },
+    'BigNumber, BigNumber': function BigNumberBigNumber(x, y) {
+      return x.minus(y);
+    },
+    'Fraction, Fraction': function FractionFraction(x, y) {
+      return x.sub(y);
+    },
+    'Unit, Unit': typed.referToSelf(self => (x, y) => {
+      if (x.value === null || x.value === undefined) {
+        throw new Error('Parameter x contains a unit with undefined value');
+      }
+      if (y.value === null || y.value === undefined) {
+        throw new Error('Parameter y contains a unit with undefined value');
+      }
+      if (!x.equalBase(y)) throw new Error('Units do not match');
+      var res = x.clone();
+      res.value = typed.find(self, [res.valueType(), y.valueType()])(res.value, y.value);
+      res.fixPrefix = false;
+      return res;
+    })
+  });
+});

package/lib/esm/function/combinatorics/stirlingS2.js

@@ -1,12 +1,12 @@
 import { factory } from '../../utils/factory.js';
 import { isNumber } from '../../utils/is.js';
 var name = 'stirlingS2';
-var dependencies = ['typed', 'addScalar', 'subtract', 'multiplyScalar', 'divideScalar', 'pow', 'factorial', 'combinations', 'isNegative', 'isInteger', 'number', '?bignumber', 'larger'];
+var dependencies = ['typed', 'addScalar', 'subtractScalar', 'multiplyScalar', 'divideScalar', 'pow', 'factorial', 'combinations', 'isNegative', 'isInteger', 'number', '?bignumber', 'larger'];
 export var createStirlingS2 = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     addScalar,
-    subtract,
+    subtractScalar,
     multiplyScalar,
     divideScalar,
     pow,

package/lib/esm/function/geometry/distance.js

@@ -1,12 +1,12 @@
 import { isBigNumber } from '../../utils/is.js';
 import { factory } from '../../utils/factory.js';
 var name = 'distance';
-var dependencies = ['typed', 'addScalar', 'subtract', 'divideScalar', 'multiplyScalar', 'deepEqual', 'sqrt', 'abs'];
+var dependencies = ['typed', 'addScalar', 'subtractScalar', 'divideScalar', 'multiplyScalar', 'deepEqual', 'sqrt', 'abs'];
 export var createDistance = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     addScalar,
-    subtract,
+    subtractScalar,
     multiplyScalar,
     divideScalar,
     deepEqual,
@@ -80,9 +80,9 @@ export var createDistance = /* #__PURE__ */factory(name, dependencies, _ref => {
         if (deepEqual(y, z)) {
           throw new TypeError('LinePoint1 should not be same with LinePoint2');
         }
-        var xCoeff = subtract(z[1], y[1]);
-        var yCoeff = subtract(y[0], z[0]);
-        var constant = subtract(multiplyScalar(z[0], y[1]), multiplyScalar(y[0], z[1]));
+        var xCoeff = subtractScalar(z[1], y[1]);
+        var yCoeff = subtractScalar(y[0], z[0]);
+        var constant = subtractScalar(multiplyScalar(z[0], y[1]), multiplyScalar(y[0], z[1]));
         return _distancePointLine2D(x[0], x[1], xCoeff, yCoeff, constant);
       } else {
         throw new TypeError('Invalid Arguments: Try again');
@@ -103,9 +103,9 @@ export var createDistance = /* #__PURE__ */factory(name, dependencies, _ref => {
           throw new TypeError('LinePoint1 should not be same with LinePoint2');
         }
         if ('pointX' in x && 'pointY' in x && 'lineOnePtX' in y && 'lineOnePtY' in y && 'lineTwoPtX' in z && 'lineTwoPtY' in z) {
-          var xCoeff = subtract(z.lineTwoPtY, y.lineOnePtY);
-          var yCoeff = subtract(y.lineOnePtX, z.lineTwoPtX);
-          var constant = subtract(multiplyScalar(z.lineTwoPtX, y.lineOnePtY), multiplyScalar(y.lineOnePtX, z.lineTwoPtY));
+          var xCoeff = subtractScalar(z.lineTwoPtY, y.lineOnePtY);
+          var yCoeff = subtractScalar(y.lineOnePtX, z.lineTwoPtX);
+          var constant = subtractScalar(multiplyScalar(z.lineTwoPtX, y.lineOnePtY), multiplyScalar(y.lineOnePtX, z.lineTwoPtY));
           return _distancePointLine2D(x.pointX, x.pointY, xCoeff, yCoeff, constant);
         } else {
           throw new TypeError('Key names do not match');
@@ -269,7 +269,7 @@ export var createDistance = /* #__PURE__ */factory(name, dependencies, _ref => {
     return divideScalar(num, den);
   }
   function _distancePointLine3D(x, y, z, x0, y0, z0, a, b, c) {
-    var num = [subtract(multiplyScalar(subtract(y0, y), c), multiplyScalar(subtract(z0, z), b)), subtract(multiplyScalar(subtract(z0, z), a), multiplyScalar(subtract(x0, x), c)), subtract(multiplyScalar(subtract(x0, x), b), multiplyScalar(subtract(y0, y), a))];
+    var num = [subtractScalar(multiplyScalar(subtractScalar(y0, y), c), multiplyScalar(subtractScalar(z0, z), b)), subtractScalar(multiplyScalar(subtractScalar(z0, z), a), multiplyScalar(subtractScalar(x0, x), c)), subtractScalar(multiplyScalar(subtractScalar(x0, x), b), multiplyScalar(subtractScalar(y0, y), a))];
     num = sqrt(addScalar(addScalar(multiplyScalar(num[0], num[0]), multiplyScalar(num[1], num[1])), multiplyScalar(num[2], num[2])));
     var den = sqrt(addScalar(addScalar(multiplyScalar(a, a), multiplyScalar(b, b)), multiplyScalar(c, c)));
     return divideScalar(num, den);
@@ -279,7 +279,7 @@ export var createDistance = /* #__PURE__ */factory(name, dependencies, _ref => {
     var result = 0;
     var diff = 0;
     for (var i = 0; i < vectorSize; i++) {
-      diff = subtract(x[i], y[i]);
+      diff = subtractScalar(x[i], y[i]);
       result = addScalar(multiplyScalar(diff, diff), result);
     }
     return sqrt(result);

package/lib/esm/function/matrix/det.js

@@ -3,12 +3,12 @@ import { clone } from '../../utils/object.js';
 import { format } from '../../utils/string.js';
 import { factory } from '../../utils/factory.js';
 var name = 'det';
-var dependencies = ['typed', 'matrix', 'subtract', 'multiply', 'divideScalar', 'isZero', 'unaryMinus'];
+var dependencies = ['typed', 'matrix', 'subtractScalar', 'multiply', 'divideScalar', 'isZero', 'unaryMinus'];
 export var createDet = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
-    subtract,
+    subtractScalar,
     multiply,
     divideScalar,
     isZero,
@@ -104,7 +104,7 @@ export var createDet = /* #__PURE__ */factory(name, dependencies, _ref => {
     } else if (rows === 2) {
       // this is a 2 x 2 matrix
       // the determinant of [a11,a12;a21,a22] is det = a11*a22-a21*a12
-      return subtract(multiply(matrix[0][0], matrix[1][1]), multiply(matrix[1][0], matrix[0][1]));
+      return subtractScalar(multiply(matrix[0][0], matrix[1][1]), multiply(matrix[1][0], matrix[0][1]));
     } else {
       // Bareiss algorithm
       // this algorithm have same complexity as LUP decomposition (O(n^3))
@@ -132,7 +132,7 @@ export var createDet = /* #__PURE__ */factory(name, dependencies, _ref => {
         for (var i = k + 1; i < rows; i++) {
           var i_ = rowIndices[i];
           for (var j = k + 1; j < rows; j++) {
-            matrix[i_][j] = divideScalar(subtract(multiply(matrix[i_][j], piv), multiply(matrix[i_][k], matrix[k_][j])), piv_);
+            matrix[i_][j] = divideScalar(subtractScalar(multiply(matrix[i_][j], piv), multiply(matrix[i_][k], matrix[k_][j])), piv_);
           }
         }
       }