mathjs 10.6.2 → 11.0.0

package/lib/cjs/function/matrix/matrixFromRows.js

@@ -7,7 +7,7 @@ exports.createMatrixFromRows = void 0;
 
 var _factory = require("../../utils/factory.js");
 
-function _createForOfIteratorHelper(o, allowArrayLike) { var it = typeof Symbol !== "undefined" && o[Symbol.iterator] || o["@@iterator"]; if (!it) { if (Array.isArray(o) || (it = _unsupportedIterableToArray(o)) || allowArrayLike && o && typeof o.length === "number") { if (it) o = it; var i = 0; var F = function F() {}; return { s: F, n: function n() { if (i >= o.length) return { done: true }; return { done: false, value: o[i++] }; }, e: function e(_e) { throw _e; }, f: F }; } throw new TypeError("Invalid attempt to iterate non-iterable instance.\nIn order to be iterable, non-array objects must have a [Symbol.iterator]() method."); } var normalCompletion = true, didErr = false, err; return { s: function s() { it = it.call(o); }, n: function n() { var step = it.next(); normalCompletion = step.done; return step; }, e: function e(_e2) { didErr = true; err = _e2; }, f: function f() { try { if (!normalCompletion && it.return != null) it.return(); } finally { if (didErr) throw err; } } }; }
+function _createForOfIteratorHelper(o, allowArrayLike) { var it = typeof Symbol !== "undefined" && o[Symbol.iterator] || o["@@iterator"]; if (!it) { if (Array.isArray(o) || (it = _unsupportedIterableToArray(o)) || allowArrayLike && o && typeof o.length === "number") { if (it) o = it; var i = 0; var F = function F() {}; return { s: F, n: function n() { if (i >= o.length) return { done: true }; return { done: false, value: o[i++] }; }, e: function e(_e) { throw _e; }, f: F }; } throw new TypeError("Invalid attempt to iterate non-iterable instance.\nIn order to be iterable, non-array objects must have a [Symbol.iterator]() method."); } var normalCompletion = true, didErr = false, err; return { s: function s() { it = it.call(o); }, n: function n() { var step = it.next(); normalCompletion = step.done; return step; }, e: function e(_e2) { didErr = true; err = _e2; }, f: function f() { try { if (!normalCompletion && it["return"] != null) it["return"](); } finally { if (didErr) throw err; } } }; }
 
 function _unsupportedIterableToArray(o, minLen) { if (!o) return; if (typeof o === "string") return _arrayLikeToArray(o, minLen); var n = Object.prototype.toString.call(o).slice(8, -1); if (n === "Object" && o.constructor) n = o.constructor.name; if (n === "Map" || n === "Set") return Array.from(o); if (n === "Arguments" || /^(?:Ui|I)nt(?:8|16|32)(?:Clamped)?Array$/.test(n)) return _arrayLikeToArray(o, minLen); }
 

package/lib/cjs/function/matrix/ones.js

@@ -38,6 +38,7 @@ var createOnes = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functi
    *
    * Examples:
    *
+   *    math.ones()                    // returns []
    *    math.ones(3)                   // returns [1, 1, 1]
    *    math.ones(3, 2)                // returns [[1, 1], [1, 1], [1, 1]]
    *    math.ones(3, 2, 'dense')       // returns Dense Matrix [[1, 1], [1, 1], [1, 1]]
@@ -49,7 +50,7 @@ var createOnes = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functi
    *
    *    zeros, identity, size, range
    *
-   * @param {...number | Array} size    The size of each dimension of the matrix
+   * @param {...(number|BigNumber) | Array} size    The size of each dimension of the matrix
    * @param {string} [format]           The Matrix storage format
    *
    * @return {Array | Matrix | number}  A matrix filled with ones

package/lib/cjs/function/matrix/sqrtm.js

@@ -14,12 +14,13 @@ var _array = require("../../utils/array.js");
 var _factory = require("../../utils/factory.js");
 
 var name = 'sqrtm';
-var dependencies = ['typed', 'abs', 'add', 'multiply', 'sqrt', 'subtract', 'inv', 'size', 'max', 'identity'];
+var dependencies = ['typed', 'abs', 'add', 'multiply', 'map', 'sqrt', 'subtract', 'inv', 'size', 'max', 'identity'];
 var createSqrtm = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
       abs = _ref.abs,
       add = _ref.add,
       multiply = _ref.multiply,
+      map = _ref.map,
       sqrt = _ref.sqrt,
       subtract = _ref.subtract,
       inv = _ref.inv,
@@ -88,7 +89,7 @@ var createSqrtm = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
         case 1:
           // Single element Array | Matrix
           if (size[0] === 1) {
-            return sqrt(A);
+            return map(A, sqrt);
           } else {
             throw new RangeError('Matrix must be square ' + '(size: ' + (0, _string.format)(size) + ')');
           }

package/lib/cjs/function/matrix/transpose.js

@@ -39,59 +39,58 @@ var createTranspose = /* #__PURE__ */(0, _factory.factory)(name, dependencies, f
    * @param {Array | Matrix} x  Matrix to be transposed
    * @return {Array | Matrix}   The transposed matrix
    */
-  return typed('transpose', {
+  return typed(name, {
     Array: function Array(x) {
-      // use dense matrix implementation
-      return this(matrix(x)).valueOf();
+      return transposeMatrix(matrix(x)).valueOf();
     },
-    Matrix: function Matrix(x) {
-      // matrix size
-      var size = x.size(); // result
-
-      var c; // process dimensions
-
-      switch (size.length) {
-        case 1:
-          // vector
-          c = x.clone();
-          break;
-
-        case 2:
-          {
-            // rows and columns
-            var rows = size[0];
-            var columns = size[1]; // check columns
-
-            if (columns === 0) {
-              // throw exception
-              throw new RangeError('Cannot transpose a 2D matrix with no columns (size: ' + (0, _string.format)(size) + ')');
-            } // process storage format
-
-
-            switch (x.storage()) {
-              case 'dense':
-                c = _denseTranspose(x, rows, columns);
-                break;
-
-              case 'sparse':
-                c = _sparseTranspose(x, rows, columns);
-                break;
-            }
-          }
-          break;
+    Matrix: transposeMatrix,
+    any: _object.clone // scalars
 
-        default:
-          // multi dimensional
-          throw new RangeError('Matrix must be a vector or two dimensional (size: ' + (0, _string.format)(this._size) + ')');
-      }
+  });
 
-      return c;
-    },
-    // scalars
-    any: function any(x) {
-      return (0, _object.clone)(x);
+  function transposeMatrix(x) {
+    // matrix size
+    var size = x.size(); // result
+
+    var c; // process dimensions
+
+    switch (size.length) {
+      case 1:
+        // vector
+        c = x.clone();
+        break;
+
+      case 2:
+        {
+          // rows and columns
+          var rows = size[0];
+          var columns = size[1]; // check columns
+
+          if (columns === 0) {
+            // throw exception
+            throw new RangeError('Cannot transpose a 2D matrix with no columns (size: ' + (0, _string.format)(size) + ')');
+          } // process storage format
+
+
+          switch (x.storage()) {
+            case 'dense':
+              c = _denseTranspose(x, rows, columns);
+              break;
+
+            case 'sparse':
+              c = _sparseTranspose(x, rows, columns);
+              break;
+          }
+        }
+        break;
+
+      default:
+        // multi dimensional
+        throw new RangeError('Matrix must be a vector or two dimensional (size: ' + (0, _string.format)(size) + ')');
     }
-  });
+
+    return c;
+  }
 
   function _denseTranspose(m, rows, columns) {
     // matrix array

package/lib/cjs/function/matrix/zeros.js

@@ -36,6 +36,7 @@ var createZeros = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
    *
    * Examples:
    *
+   *    math.zeros()                   // returns []
    *    math.zeros(3)                  // returns [0, 0, 0]
    *    math.zeros(3, 2)               // returns [[0, 0], [0, 0], [0, 0]]
    *    math.zeros(3, 'dense')         // returns [0, 0, 0]
@@ -47,7 +48,7 @@ var createZeros = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
    *
    *    ones, identity, size, range
    *
-   * @param {...number | Array} size    The size of each dimension of the matrix
+   * @param {...(number|BigNumber) | Array} size    The size of each dimension of the matrix
    * @param {string} [format]           The Matrix storage format
    *
    * @return {Array | Matrix}           A matrix filled with zeros

package/lib/cjs/function/probability/factorial.js

@@ -52,9 +52,11 @@ var createFactorial = /* #__PURE__ */(0, _factory.factory)(name, dependencies, f
 
       return gamma(n.plus(1));
     },
-    'Array | Matrix': function ArrayMatrix(n) {
-      return (0, _collection.deepMap)(n, this);
-    }
+    'Array | Matrix': typed.referToSelf(function (self) {
+      return function (n) {
+        return (0, _collection.deepMap)(n, self);
+      };
+    })
   });
 });
 exports.createFactorial = createFactorial;

package/lib/cjs/function/probability/gamma.js

@@ -5,8 +5,6 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.createGamma = void 0;
 
-var _collection = require("../../utils/collection.js");
-
 var _factory = require("../../utils/factory.js");
 
 var _index = require("../../plain/number/index.js");
@@ -19,13 +17,14 @@ var createGamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
       multiplyScalar = _ref.multiplyScalar,
       pow = _ref.pow,
       _BigNumber = _ref.BigNumber,
-      _Complex = _ref.Complex;
+      Complex = _ref.Complex;
 
   /**
    * Compute the gamma function of a value using Lanczos approximation for
    * small values, and an extended Stirling approximation for large values.
    *
-   * For matrices, the function is evaluated element wise.
+   * To avoid confusion with the matrix Gamma function, this function does
+   * not apply to matrices.
    *
    * Syntax:
    *
@@ -41,51 +40,53 @@ var createGamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
    *
    *    combinations, factorial, permutations
    *
-   * @param {number | Array | Matrix} n   A real or complex number
-   * @return {number | Array | Matrix}    The gamma of `n`
+   * @param {number | BigNumber | Complex} n   A real or complex number
+   * @return {number | BigNumber | Complex}    The gamma of `n`
    */
-  return typed(name, {
-    number: _index.gammaNumber,
-    Complex: function Complex(n) {
-      if (n.im === 0) {
-        return this(n.re);
-      } // Lanczos approximation doesn't work well with real part lower than 0.5
-      // So reflection formula is required
+  function gammaComplex(n) {
+    if (n.im === 0) {
+      return (0, _index.gammaNumber)(n.re);
+    } // Lanczos approximation doesn't work well with real part lower than 0.5
+    // So reflection formula is required
 
 
-      if (n.re < 0.5) {
-        // Euler's reflection formula
-        // gamma(1-z) * gamma(z) = PI / sin(PI * z)
-        // real part of Z should not be integer [sin(PI) == 0 -> 1/0 - undefined]
-        // thanks to imperfect sin implementation sin(PI * n) != 0
-        // we can safely use it anyway
-        var _t = new _Complex(1 - n.re, -n.im);
+    if (n.re < 0.5) {
+      // Euler's reflection formula
+      // gamma(1-z) * gamma(z) = PI / sin(PI * z)
+      // real part of Z should not be integer [sin(PI) == 0 -> 1/0 - undefined]
+      // thanks to imperfect sin implementation sin(PI * n) != 0
+      // we can safely use it anyway
+      var _t = new Complex(1 - n.re, -n.im);
 
-        var r = new _Complex(Math.PI * n.re, Math.PI * n.im);
-        return new _Complex(Math.PI).div(r.sin()).div(this(_t));
-      } // Lanczos approximation
-      // z -= 1
+      var r = new Complex(Math.PI * n.re, Math.PI * n.im);
+      return new Complex(Math.PI).div(r.sin()).div(gammaComplex(_t));
+    } // Lanczos approximation
+    // z -= 1
 
 
-      n = new _Complex(n.re - 1, n.im); // x = gammaPval[0]
+    n = new Complex(n.re - 1, n.im); // x = gammaPval[0]
 
-      var x = new _Complex(_index.gammaP[0], 0); // for (i, gammaPval) in enumerate(gammaP):
+    var x = new Complex(_index.gammaP[0], 0); // for (i, gammaPval) in enumerate(gammaP):
 
-      for (var i = 1; i < _index.gammaP.length; ++i) {
-        // x += gammaPval / (z + i)
-        var gammaPval = new _Complex(_index.gammaP[i], 0);
-        x = x.add(gammaPval.div(n.add(i)));
-      } // t = z + gammaG + 0.5
+    for (var i = 1; i < _index.gammaP.length; ++i) {
+      // x += gammaPval / (z + i)
+      var gammaPval = new Complex(_index.gammaP[i], 0);
+      x = x.add(gammaPval.div(n.add(i)));
+    } // t = z + gammaG + 0.5
 
 
-      var t = new _Complex(n.re + _index.gammaG + 0.5, n.im); // y = sqrt(2 * pi) * t ** (z + 0.5) * exp(-t) * x
+    var t = new Complex(n.re + _index.gammaG + 0.5, n.im); // y = sqrt(2 * pi) * t ** (z + 0.5) * exp(-t) * x
 
-      var twoPiSqrt = Math.sqrt(2 * Math.PI);
-      var tpow = t.pow(n.add(0.5));
-      var expt = t.neg().exp(); // y = [x] * [sqrt(2 * pi)] * [t ** (z + 0.5)] * [exp(-t)]
+    var twoPiSqrt = Math.sqrt(2 * Math.PI);
+    var tpow = t.pow(n.add(0.5));
+    var expt = t.neg().exp(); // y = [x] * [sqrt(2 * pi)] * [t ** (z + 0.5)] * [exp(-t)]
 
-      return x.mul(twoPiSqrt).mul(tpow).mul(expt);
-    },
+    return x.mul(twoPiSqrt).mul(tpow).mul(expt);
+  }
+
+  return typed(name, {
+    number: _index.gammaNumber,
+    Complex: gammaComplex,
     BigNumber: function BigNumber(n) {
       if (n.isInteger()) {
         return n.isNegative() || n.isZero() ? new _BigNumber(Infinity) : bigFactorial(n.minus(1));
@@ -96,9 +97,6 @@ var createGamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
       }
 
       throw new Error('Integer BigNumber expected');
-    },
-    'Array | Matrix': function ArrayMatrix(n) {
-      return (0, _collection.deepMap)(n, this);
     }
   });
   /**

package/lib/cjs/function/probability/kldivergence.js

@@ -8,13 +8,14 @@ exports.createKldivergence = void 0;
 var _factory = require("../../utils/factory.js");
 
 var name = 'kldivergence';
-var dependencies = ['typed', 'matrix', 'divide', 'sum', 'multiply', 'dotDivide', 'log', 'isNumeric'];
+var dependencies = ['typed', 'matrix', 'divide', 'sum', 'multiply', 'map', 'dotDivide', 'log', 'isNumeric'];
 var createKldivergence = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
       matrix = _ref.matrix,
       divide = _ref.divide,
       sum = _ref.sum,
       multiply = _ref.multiply,
+      map = _ref.map,
       dotDivide = _ref.dotDivide,
       log = _ref.log,
       isNumeric = _ref.isNumeric;
@@ -81,7 +82,9 @@ var createKldivergence = /* #__PURE__ */(0, _factory.factory)(name, dependencies
 
     var qnorm = divide(q, sum(q));
     var pnorm = divide(p, sum(p));
-    var result = sum(multiply(qnorm, log(dotDivide(qnorm, pnorm))));
+    var result = sum(multiply(qnorm, map(dotDivide(qnorm, pnorm), function (x) {
+      return log(x);
+    })));
 
     if (isNumeric(result)) {
       return result;

package/lib/cjs/function/probability/lgamma.js

@@ -19,7 +19,7 @@ var _number = require("../../utils/number.js");
 var name = 'lgamma';
 var dependencies = ['Complex', 'typed'];
 var createLgamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
-  var _Complex = _ref.Complex,
+  var Complex = _ref.Complex,
       typed = _ref.typed;
   // Stirling series is non-convergent, we need to use the recurrence `lgamma(z) = lgamma(z+1) - log z` to get
   // sufficient accuracy.
@@ -65,38 +65,38 @@ var createLgamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, func
 
   return typed(name, {
     number: _index.lgammaNumber,
-    Complex: function Complex(n) {
-      var TWOPI = 6.2831853071795864769252842; // 2*pi
-
-      var LOGPI = 1.1447298858494001741434262; // log(pi)
-
-      var REFLECTION = 0.1;
-
-      if (n.isNaN()) {
-        return new _Complex(NaN, NaN);
-      } else if (n.im === 0) {
-        return new _Complex((0, _index.lgammaNumber)(n.re), 0);
-      } else if (n.re >= SMALL_RE || Math.abs(n.im) >= SMALL_IM) {
-        return lgammaStirling(n);
-      } else if (n.re <= REFLECTION) {
-        // Reflection formula. see Proposition 3.1 in [1]
-        var tmp = (0, _number.copysign)(TWOPI, n.im) * Math.floor(0.5 * n.re + 0.25); // TODO: `complex.js sin` doesn't have extremely high precision, so this value `a` may lose a little precision,
-        // causing the computation results to be less accurate than the lgamma of real numbers
-
-        var a = n.mul(Math.PI).sin().log();
-        var b = this(new _Complex(1 - n.re, -n.im));
-        return new _Complex(LOGPI, tmp).sub(a).sub(b);
-      } else if (n.im >= 0) {
-        return lgammaRecurrence(n);
-      } else {
-        return lgammaRecurrence(n.conjugate()).conjugate();
-      }
-    },
+    Complex: lgammaComplex,
     BigNumber: function BigNumber() {
       throw new Error("mathjs doesn't yet provide an implementation of the algorithm lgamma for BigNumber");
     }
   });
 
+  function lgammaComplex(n) {
+    var TWOPI = 6.2831853071795864769252842; // 2*pi
+
+    var LOGPI = 1.1447298858494001741434262; // log(pi)
+
+    var REFLECTION = 0.1;
+
+    if (n.isNaN()) {
+      return new Complex(NaN, NaN);
+    } else if (n.im === 0) {
+      return new Complex((0, _index.lgammaNumber)(n.re), 0);
+    } else if (n.re >= SMALL_RE || Math.abs(n.im) >= SMALL_IM) {
+      return lgammaStirling(n);
+    } else if (n.re <= REFLECTION) {
+      // Reflection formula. see Proposition 3.1 in [1]
+      var tmp = (0, _number.copysign)(TWOPI, n.im) * Math.floor(0.5 * n.re + 0.25);
+      var a = n.mul(Math.PI).sin().log();
+      var b = lgammaComplex(new Complex(1 - n.re, -n.im));
+      return new Complex(LOGPI, tmp).sub(a).sub(b);
+    } else if (n.im >= 0) {
+      return lgammaRecurrence(n);
+    } else {
+      return lgammaRecurrence(n.conjugate()).conjugate();
+    }
+  }
+
   function lgammaStirling(z) {
     // formula ref in [2]
     // computation ref:
@@ -106,7 +106,7 @@ var createLgamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, func
     // => (x - 0.5) * log(x) - x + log(2PI) / 2
     var leftPart = z.sub(0.5).mul(z.log()).sub(z).add(_index.lnSqrt2PI); // right part
 
-    var rz = new _Complex(1, 0).div(z);
+    var rz = new Complex(1, 0).div(z);
     var rzz = rz.div(z);
     var a = coeffs[0];
     var b = coeffs[1];
@@ -140,7 +140,7 @@ var createLgamma = /* #__PURE__ */(0, _factory.factory)(name, dependencies, func
       z = z.add(1);
     }
 
-    return lgammaStirling(z).sub(shiftprod.log()).sub(new _Complex(0, signflips * 2 * Math.PI * 1));
+    return lgammaStirling(z).sub(shiftprod.log()).sub(new Complex(0, signflips * 2 * Math.PI * 1));
   }
 });
 exports.createLgamma = createLgamma;

package/lib/cjs/function/probability/util/seededRNG.js

@@ -9,13 +9,13 @@ exports.createRng = createRng;
 
 var _seedrandom = _interopRequireDefault(require("seedrandom"));
 
-var singletonRandom = /* #__PURE__ */(0, _seedrandom.default)(Date.now());
+var singletonRandom = /* #__PURE__ */(0, _seedrandom["default"])(Date.now());
 
 function createRng(randomSeed) {
   var random; // create a new random generator with given seed
 
   function setSeed(seed) {
-    random = seed === null ? singletonRandom : (0, _seedrandom.default)(String(seed));
+    random = seed === null ? singletonRandom : (0, _seedrandom["default"])(String(seed));
   } // initialize a seeded pseudo random number generator with config's random seed
 
 

package/lib/cjs/function/relational/compare.js

@@ -11,15 +11,15 @@ var _number = require("../../utils/number.js");
 
 var _factory = require("../../utils/factory.js");
 
-var _algorithm = require("../../type/matrix/utils/algorithm03.js");
+var _matAlgo03xDSf = require("../../type/matrix/utils/matAlgo03xDSf.js");
 
-var _algorithm2 = require("../../type/matrix/utils/algorithm12.js");
+var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 
-var _algorithm3 = require("../../type/matrix/utils/algorithm14.js");
+var _matAlgo05xSfSf = require("../../type/matrix/utils/matAlgo05xSfSf.js");
 
-var _algorithm4 = require("../../type/matrix/utils/algorithm13.js");
+var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 
-var _algorithm5 = require("../../type/matrix/utils/algorithm05.js");
+var _compareUnits = require("./compareUnits.js");
 
 var name = 'compare';
 var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'BigNumber', 'Fraction', 'DenseMatrix'];
@@ -31,21 +31,22 @@ var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
       BigNumber = _ref.BigNumber,
       Fraction = _ref.Fraction,
       DenseMatrix = _ref.DenseMatrix;
-  var algorithm03 = (0, _algorithm.createAlgorithm03)({
+  var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
-  var algorithm05 = (0, _algorithm5.createAlgorithm05)({
+  var matAlgo05xSfSf = (0, _matAlgo05xSfSf.createMatAlgo05xSfSf)({
     typed: typed,
     equalScalar: equalScalar
   });
-  var algorithm12 = (0, _algorithm2.createAlgorithm12)({
+  var matAlgo12xSfs = (0, _matAlgo12xSfs.createMatAlgo12xSfs)({
     typed: typed,
     DenseMatrix: DenseMatrix
   });
-  var algorithm13 = (0, _algorithm4.createAlgorithm13)({
-    typed: typed
+  var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
+    typed: typed,
+    matrix: matrix
   });
-  var algorithm14 = (0, _algorithm3.createAlgorithm14)({
+  var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed
   });
   /**
@@ -86,13 +87,13 @@ var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
    *                                                          1 when x > y, -1 when x < y, and 0 when x == y.
    */
 
-  return typed(name, {
+  return typed(name, createCompareNumber({
+    typed: typed,
+    config: config
+  }), {
     'boolean, boolean': function booleanBoolean(x, y) {
       return x === y ? 0 : x > y ? 1 : -1;
     },
-    'number, number': function numberNumber(x, y) {
-      return (0, _number.nearlyEqual)(x, y, config.epsilon) ? 0 : x > y ? 1 : -1;
-    },
     'BigNumber, BigNumber': function BigNumberBigNumber(x, y) {
       return (0, _nearlyEqual.nearlyEqual)(x, y, config.epsilon) ? new BigNumber(0) : new BigNumber(x.cmp(y));
     },
@@ -101,59 +102,12 @@ var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
     },
     'Complex, Complex': function ComplexComplex() {
       throw new TypeError('No ordering relation is defined for complex numbers');
-    },
-    'Unit, Unit': function UnitUnit(x, y) {
-      if (!x.equalBase(y)) {
-        throw new Error('Cannot compare units with different base');
-      }
-
-      return this(x.value, y.value);
-    },
-    'SparseMatrix, SparseMatrix': function SparseMatrixSparseMatrix(x, y) {
-      return algorithm05(x, y, this);
-    },
-    'SparseMatrix, DenseMatrix': function SparseMatrixDenseMatrix(x, y) {
-      return algorithm03(y, x, this, true);
-    },
-    'DenseMatrix, SparseMatrix': function DenseMatrixSparseMatrix(x, y) {
-      return algorithm03(x, y, this, false);
-    },
-    'DenseMatrix, DenseMatrix': function DenseMatrixDenseMatrix(x, y) {
-      return algorithm13(x, y, this);
-    },
-    'Array, Array': function ArrayArray(x, y) {
-      // use matrix implementation
-      return this(matrix(x), matrix(y)).valueOf();
-    },
-    'Array, Matrix': function ArrayMatrix(x, y) {
-      // use matrix implementation
-      return this(matrix(x), y);
-    },
-    'Matrix, Array': function MatrixArray(x, y) {
-      // use matrix implementation
-      return this(x, matrix(y));
-    },
-    'SparseMatrix, any': function SparseMatrixAny(x, y) {
-      return algorithm12(x, y, this, false);
-    },
-    'DenseMatrix, any': function DenseMatrixAny(x, y) {
-      return algorithm14(x, y, this, false);
-    },
-    'any, SparseMatrix': function anySparseMatrix(x, y) {
-      return algorithm12(y, x, this, true);
-    },
-    'any, DenseMatrix': function anyDenseMatrix(x, y) {
-      return algorithm14(y, x, this, true);
-    },
-    'Array, any': function ArrayAny(x, y) {
-      // use matrix implementation
-      return algorithm14(matrix(x), y, this, false).valueOf();
-    },
-    'any, Array': function anyArray(x, y) {
-      // use matrix implementation
-      return algorithm14(matrix(y), x, this, true).valueOf();
     }
-  });
+  }, compareUnits, matrixAlgorithmSuite({
+    SS: matAlgo05xSfSf,
+    DS: matAlgo03xDSf,
+    Ss: matAlgo12xSfs
+  }));
 });
 exports.createCompare = createCompare;
 var createCompareNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'config'], function (_ref2) {