mathjs 10.4.2 → 10.5.1

package/lib/cjs/function/arithmetic/ceil.js

@@ -5,7 +5,7 @@ var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefau
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.createCeil = void 0;
+exports.createCeilNumber = exports.createCeil = void 0;
 
 var _slicedToArray2 = _interopRequireDefault(require("@babel/runtime/helpers/slicedToArray"));
 
@@ -19,27 +19,73 @@ var _number = require("../../utils/number.js");
 
 var _nearlyEqual = require("../../utils/bignumber/nearlyEqual.js");
 
-var _index = require("../../plain/number/index.js");
-
 var _algorithm = require("../../type/matrix/utils/algorithm11.js");
 
-var _algorithm2 = require("../../type/matrix/utils/algorithm14.js");
+var _algorithm2 = require("../../type/matrix/utils/algorithm12.js");
+
+var _algorithm3 = require("../../type/matrix/utils/algorithm14.js");
 
 var name = 'ceil';
-var dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar'];
-var createCeil = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar', 'zeros', 'DenseMatrix'];
+var createCeilNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'config', 'round'], function (_ref) {
   var typed = _ref.typed,
       config = _ref.config,
-      round = _ref.round,
-      matrix = _ref.matrix,
-      equalScalar = _ref.equalScalar;
+      round = _ref.round;
+  return typed(name, {
+    number: function number(x) {
+      if ((0, _number.nearlyEqual)(x, round(x), config.epsilon)) {
+        return round(x);
+      } else {
+        return Math.ceil(x);
+      }
+    },
+    'number, number': function numberNumber(x, n) {
+      if ((0, _number.nearlyEqual)(x, round(x, n), config.epsilon)) {
+        return round(x, n);
+      } else {
+        var _$split = "".concat(x, "e").split('e'),
+            _$split2 = (0, _slicedToArray2.default)(_$split, 2),
+            number = _$split2[0],
+            exponent = _$split2[1];
+
+        var result = Math.ceil(Number("".concat(number, "e").concat(Number(exponent) + n)));
+
+        var _$split3 = "".concat(result, "e").split('e');
+
+        var _$split4 = (0, _slicedToArray2.default)(_$split3, 2);
+
+        number = _$split4[0];
+        exponent = _$split4[1];
+        return Number("".concat(number, "e").concat(Number(exponent) - n));
+      }
+    }
+  });
+});
+exports.createCeilNumber = createCeilNumber;
+var createCeil = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref2) {
+  var typed = _ref2.typed,
+      config = _ref2.config,
+      round = _ref2.round,
+      matrix = _ref2.matrix,
+      equalScalar = _ref2.equalScalar,
+      zeros = _ref2.zeros,
+      DenseMatrix = _ref2.DenseMatrix;
   var algorithm11 = (0, _algorithm.createAlgorithm11)({
     typed: typed,
     equalScalar: equalScalar
   });
-  var algorithm14 = (0, _algorithm2.createAlgorithm14)({
+  var algorithm12 = (0, _algorithm2.createAlgorithm12)({
+    typed: typed,
+    DenseMatrix: DenseMatrix
+  });
+  var algorithm14 = (0, _algorithm3.createAlgorithm14)({
     typed: typed
   });
+  var ceilNumber = createCeilNumber({
+    typed: typed,
+    config: config,
+    round: round
+  });
   /**
    * Round a value towards plus infinity
    * If `x` is complex, both real and imaginary part are rounded towards plus infinity.
@@ -79,39 +125,17 @@ var createCeil = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functi
    */
 
   return typed('ceil', {
-    number: function number(x) {
-      if ((0, _number.nearlyEqual)(x, round(x), config.epsilon)) {
-        return round(x);
-      } else {
-        return (0, _index.ceilNumber)(x);
-      }
-    },
-    'number, number': function numberNumber(x, n) {
-      if ((0, _number.nearlyEqual)(x, round(x, n), config.epsilon)) {
-        return round(x, n);
-      } else {
-        var _$split = "".concat(x, "e").split('e'),
-            _$split2 = (0, _slicedToArray2.default)(_$split, 2),
-            number = _$split2[0],
-            exponent = _$split2[1];
-
-        var result = Math.ceil(Number("".concat(number, "e").concat(Number(exponent) + n)));
-
-        var _$split3 = "".concat(result, "e").split('e');
-
-        var _$split4 = (0, _slicedToArray2.default)(_$split3, 2);
-
-        number = _$split4[0];
-        exponent = _$split4[1];
-        return Number("".concat(number, "e").concat(Number(exponent) - n));
-      }
-    },
+    number: ceilNumber.signatures.number,
+    'number,number': ceilNumber.signatures['number,number'],
     Complex: function Complex(x) {
       return x.ceil();
     },
     'Complex, number': function ComplexNumber(x, n) {
       return x.ceil(n);
     },
+    'Complex, BigNumber': function ComplexBigNumber(x, n) {
+      return x.ceil(n.toNumber());
+    },
     BigNumber: function BigNumber(x) {
       if ((0, _nearlyEqual.nearlyEqual)(x, round(x), config.epsilon)) {
         return round(x);
@@ -132,11 +156,14 @@ var createCeil = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functi
     'Fraction, number': function FractionNumber(x, n) {
       return x.ceil(n);
     },
+    'Fraction, BigNumber': function FractionBigNumber(x, n) {
+      return x.ceil(n.toNumber());
+    },
     'Array | Matrix': function ArrayMatrix(x) {
       // deep map collection, skip zeros since ceil(0) = 0
       return (0, _collection.deepMap)(x, this, true);
     },
-    'Array | Matrix, number': function ArrayMatrixNumber(x, n) {
+    'Array, number | BigNumber': function ArrayNumberBigNumber(x, n) {
       var _this = this;
 
       // deep map collection, skip zeros since ceil(0) = 0
@@ -150,9 +177,18 @@ var createCeil = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functi
     'DenseMatrix, number | BigNumber': function DenseMatrixNumberBigNumber(x, y) {
       return algorithm14(x, y, this, false);
     },
-    'number | Complex | BigNumber, Array': function numberComplexBigNumberArray(x, y) {
+    'number | Complex | Fraction | BigNumber, Array': function numberComplexFractionBigNumberArray(x, y) {
       // use matrix implementation
       return algorithm14(matrix(y), x, this, true).valueOf();
+    },
+    'number | Complex | Fraction | BigNumber, Matrix': function numberComplexFractionBigNumberMatrix(x, y) {
+      if (equalScalar(x, 0)) return zeros(y.size(), y.storage());
+
+      if (y.storage() === 'dense') {
+        return algorithm14(y, x, this, true);
+      }
+
+      return algorithm12(y, x, this, true);
     }
   });
 });

package/lib/cjs/function/arithmetic/fix.js

@@ -3,25 +3,53 @@
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.createFix = void 0;
+exports.createFixNumber = exports.createFix = void 0;
 
 var _factory = require("../../utils/factory.js");
 
 var _collection = require("../../utils/collection.js");
 
-var _algorithm = require("../../type/matrix/utils/algorithm14.js");
+var _algorithm = require("../../type/matrix/utils/algorithm12.js");
+
+var _algorithm2 = require("../../type/matrix/utils/algorithm14.js");
 
 var name = 'fix';
-var dependencies = ['typed', 'Complex', 'matrix', 'ceil', 'floor'];
-var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var dependencies = ['typed', 'Complex', 'matrix', 'ceil', 'floor', 'equalScalar', 'zeros', 'DenseMatrix'];
+var createFixNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'ceil', 'floor'], function (_ref) {
   var typed = _ref.typed,
-      _Complex = _ref.Complex,
-      matrix = _ref.matrix,
       ceil = _ref.ceil,
       floor = _ref.floor;
-  var algorithm14 = (0, _algorithm.createAlgorithm14)({
+  return typed(name, {
+    number: function number(x) {
+      return x > 0 ? floor(x) : ceil(x);
+    },
+    'number, number': function numberNumber(x, n) {
+      return x > 0 ? floor(x, n) : ceil(x, n);
+    }
+  });
+});
+exports.createFixNumber = createFixNumber;
+var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref2) {
+  var typed = _ref2.typed,
+      _Complex = _ref2.Complex,
+      matrix = _ref2.matrix,
+      ceil = _ref2.ceil,
+      floor = _ref2.floor,
+      equalScalar = _ref2.equalScalar,
+      zeros = _ref2.zeros,
+      DenseMatrix = _ref2.DenseMatrix;
+  var algorithm12 = (0, _algorithm.createAlgorithm12)({
+    typed: typed,
+    DenseMatrix: DenseMatrix
+  });
+  var algorithm14 = (0, _algorithm2.createAlgorithm14)({
     typed: typed
   });
+  var fixNumber = createFixNumber({
+    typed: typed,
+    ceil: ceil,
+    floor: floor
+  });
   /**
    * Round a value towards zero.
    * For matrices, the function is evaluated element wise.
@@ -29,6 +57,7 @@ var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
    * Syntax:
    *
    *    math.fix(x)
+   *    math.fix(x,n)
    *
    * Examples:
    *
@@ -59,16 +88,16 @@ var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
    */
 
   return typed('fix', {
-    number: function number(x) {
-      return x > 0 ? floor(x) : ceil(x);
-    },
-    'number, number | BigNumber': function numberNumberBigNumber(x, n) {
-      return x > 0 ? floor(x, n) : ceil(x, n);
-    },
+    number: fixNumber.signatures.number,
+    'number, number | BigNumber': fixNumber.signatures['number,number'],
     Complex: function Complex(x) {
       return new _Complex(x.re > 0 ? Math.floor(x.re) : Math.ceil(x.re), x.im > 0 ? Math.floor(x.im) : Math.ceil(x.im));
     },
-    'Complex, number | BigNumber': function ComplexNumberBigNumber(x, n) {
+    'Complex, number': function ComplexNumber(x, n) {
+      return new _Complex(x.re > 0 ? floor(x.re, n) : ceil(x.re, n), x.im > 0 ? floor(x.im, n) : ceil(x.im, n));
+    },
+    'Complex, BigNumber': function ComplexBigNumber(x, bn) {
+      var n = bn.toNumber();
       return new _Complex(x.re > 0 ? floor(x.re, n) : ceil(x.re, n), x.im > 0 ? floor(x.im, n) : ceil(x.im, n));
     },
     BigNumber: function BigNumber(x) {
@@ -81,7 +110,7 @@ var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
       return x.s < 0 ? x.ceil() : x.floor();
     },
     'Fraction, number | BigNumber': function FractionNumberBigNumber(x, n) {
-      return x.s < 0 ? x.ceil(n) : x.floor(n);
+      return x.s < 0 ? ceil(x, n) : floor(x, n);
     },
     'Array | Matrix': function ArrayMatrix(x) {
       // deep map collection, skip zeros since fix(0) = 0
@@ -95,9 +124,18 @@ var createFix = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
         return _this(i, n);
       }, true);
     },
-    'number | Complex | BigNumber, Array': function numberComplexBigNumberArray(x, y) {
+    'number | Complex | Fraction | BigNumber, Array': function numberComplexFractionBigNumberArray(x, y) {
       // use matrix implementation
       return algorithm14(matrix(y), x, this, true).valueOf();
+    },
+    'number | Complex | Fraction | BigNumber, Matrix': function numberComplexFractionBigNumberMatrix(x, y) {
+      if (equalScalar(x, 0)) return zeros(y.size(), y.storage());
+
+      if (y.storage() === 'dense') {
+        return algorithm14(y, x, this, true);
+      }
+
+      return algorithm12(y, x, this, true);
     }
   });
 });

package/lib/cjs/function/arithmetic/floor.js

@@ -5,7 +5,7 @@ var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefau
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.createFloor = void 0;
+exports.createFloorNumber = exports.createFloor = void 0;
 
 var _slicedToArray2 = _interopRequireDefault(require("@babel/runtime/helpers/slicedToArray"));
 
@@ -21,23 +21,71 @@ var _nearlyEqual = require("../../utils/bignumber/nearlyEqual.js");
 
 var _algorithm = require("../../type/matrix/utils/algorithm11.js");
 
-var _algorithm2 = require("../../type/matrix/utils/algorithm14.js");
+var _algorithm2 = require("../../type/matrix/utils/algorithm12.js");
+
+var _algorithm3 = require("../../type/matrix/utils/algorithm14.js");
 
 var name = 'floor';
-var dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar'];
-var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar', 'zeros', 'DenseMatrix'];
+var createFloorNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'config', 'round'], function (_ref) {
   var typed = _ref.typed,
       config = _ref.config,
-      round = _ref.round,
-      matrix = _ref.matrix,
-      equalScalar = _ref.equalScalar;
+      round = _ref.round;
+  return typed(name, {
+    number: function number(x) {
+      if ((0, _number.nearlyEqual)(x, round(x), config.epsilon)) {
+        return round(x);
+      } else {
+        return Math.floor(x);
+      }
+    },
+    'number, number': function numberNumber(x, n) {
+      if ((0, _number.nearlyEqual)(x, round(x, n), config.epsilon)) {
+        return round(x, n);
+      } else {
+        var _$split = "".concat(x, "e").split('e'),
+            _$split2 = (0, _slicedToArray2.default)(_$split, 2),
+            number = _$split2[0],
+            exponent = _$split2[1];
+
+        var result = Math.floor(Number("".concat(number, "e").concat(Number(exponent) + n)));
+
+        var _$split3 = "".concat(result, "e").split('e');
+
+        var _$split4 = (0, _slicedToArray2.default)(_$split3, 2);
+
+        number = _$split4[0];
+        exponent = _$split4[1];
+        return Number("".concat(number, "e").concat(Number(exponent) - n));
+      }
+    }
+  });
+});
+exports.createFloorNumber = createFloorNumber;
+var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref2) {
+  var typed = _ref2.typed,
+      config = _ref2.config,
+      round = _ref2.round,
+      matrix = _ref2.matrix,
+      equalScalar = _ref2.equalScalar,
+      zeros = _ref2.zeros,
+      DenseMatrix = _ref2.DenseMatrix;
   var algorithm11 = (0, _algorithm.createAlgorithm11)({
     typed: typed,
     equalScalar: equalScalar
   });
-  var algorithm14 = (0, _algorithm2.createAlgorithm14)({
+  var algorithm12 = (0, _algorithm2.createAlgorithm12)({
+    typed: typed,
+    DenseMatrix: DenseMatrix
+  });
+  var algorithm14 = (0, _algorithm3.createAlgorithm14)({
     typed: typed
   });
+  var floorNumber = createFloorNumber({
+    typed: typed,
+    config: config,
+    round: round
+  });
   /**
    * Round a value towards minus infinity.
    * For matrices, the function is evaluated element wise.
@@ -66,6 +114,10 @@ var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
    *    math.floor([3.2, 3.8, -4.7])       // returns Array [3, 3, -5]
    *    math.floor([3.21, 3.82, -4.71], 1)  // returns Array [3.2, 3.8, -4.8]
    *
+   *    math.floor(math.tau, [2, 3])  // returns Array [6.28, 6.283]
+   *
+   *    // Note that floor(array, array) currently not implemented.
+   *
    * See also:
    *
    *    ceil, fix, round
@@ -76,39 +128,17 @@ var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
    */
 
   return typed('floor', {
-    number: function number(x) {
-      if ((0, _number.nearlyEqual)(x, round(x), config.epsilon)) {
-        return round(x);
-      } else {
-        return Math.floor(x);
-      }
-    },
-    'number, number': function numberNumber(x, n) {
-      if ((0, _number.nearlyEqual)(x, round(x, n), config.epsilon)) {
-        return round(x, n);
-      } else {
-        var _$split = "".concat(x, "e").split('e'),
-            _$split2 = (0, _slicedToArray2.default)(_$split, 2),
-            number = _$split2[0],
-            exponent = _$split2[1];
-
-        var result = Math.floor(Number("".concat(number, "e").concat(Number(exponent) + n)));
-
-        var _$split3 = "".concat(result, "e").split('e');
-
-        var _$split4 = (0, _slicedToArray2.default)(_$split3, 2);
-
-        number = _$split4[0];
-        exponent = _$split4[1];
-        return Number("".concat(number, "e").concat(Number(exponent) - n));
-      }
-    },
+    number: floorNumber.signatures.number,
+    'number,number': floorNumber.signatures['number,number'],
     Complex: function Complex(x) {
       return x.floor();
     },
     'Complex, number': function ComplexNumber(x, n) {
       return x.floor(n);
     },
+    'Complex, BigNumber': function ComplexBigNumber(x, n) {
+      return x.floor(n.toNumber());
+    },
     BigNumber: function BigNumber(x) {
       if ((0, _nearlyEqual.nearlyEqual)(x, round(x), config.epsilon)) {
         return round(x);
@@ -129,11 +159,14 @@ var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
     'Fraction, number': function FractionNumber(x, n) {
       return x.floor(n);
     },
+    'Fraction, BigNumber': function FractionBigNumber(x, n) {
+      return x.floor(n.toNumber());
+    },
     'Array | Matrix': function ArrayMatrix(x) {
       // deep map collection, skip zeros since floor(0) = 0
       return (0, _collection.deepMap)(x, this, true);
     },
-    'Array | Matrix, number': function ArrayMatrixNumber(x, n) {
+    'Array, number | BigNumber': function ArrayNumberBigNumber(x, n) {
       var _this = this;
 
       // deep map collection, skip zeros since ceil(0) = 0
@@ -147,9 +180,18 @@ var createFloor = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
     'DenseMatrix, number | BigNumber': function DenseMatrixNumberBigNumber(x, y) {
       return algorithm14(x, y, this, false);
     },
-    'number | Complex | BigNumber, Array': function numberComplexBigNumberArray(x, y) {
+    'number | Complex | Fraction | BigNumber, Array': function numberComplexFractionBigNumberArray(x, y) {
       // use matrix implementation
       return algorithm14(matrix(y), x, this, true).valueOf();
+    },
+    'number | Complex | Fraction | BigNumber, Matrix': function numberComplexFractionBigNumberMatrix(x, y) {
+      if (equalScalar(x, 0)) return zeros(y.size(), y.storage());
+
+      if (y.storage() === 'dense') {
+        return algorithm14(y, x, this, true);
+      }
+
+      return algorithm12(y, x, this, true);
     }
   });
 });

package/lib/cjs/function/arithmetic/nthRoot.js

@@ -81,9 +81,7 @@ var createNthRoot = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
 
   var complexErr = '' + 'Complex number not supported in function nthRoot. ' + 'Use nthRoots instead.';
   return typed(name, {
-    number: function number(x) {
-      return (0, _index.nthRootNumber)(x, 2);
-    },
+    number: _index.nthRootNumber,
     'number, number': _index.nthRootNumber,
     BigNumber: function BigNumber(x) {
       return _bigNthRoot(x, new _BigNumber(2));

package/lib/cjs/function/arithmetic/pow.js

@@ -14,21 +14,25 @@ var _array = require("../../utils/array.js");
 var _index = require("../../plain/number/index.js");
 
 var name = 'pow';
-var dependencies = ['typed', 'config', 'identity', 'multiply', 'matrix', 'fraction', 'number', 'Complex'];
+var dependencies = ['typed', 'config', 'identity', 'multiply', 'matrix', 'inv', 'fraction', 'number', 'Complex'];
 var createPow = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
       config = _ref.config,
       identity = _ref.identity,
       multiply = _ref.multiply,
       matrix = _ref.matrix,
+      inv = _ref.inv,
       number = _ref.number,
       fraction = _ref.fraction,
       Complex = _ref.Complex;
 
   /**
    * Calculates the power of x to y, `x ^ y`.
-   * Matrix exponentiation is supported for square matrices `x`, and positive
-   * integer exponents `y`.
+   *
+   * Matrix exponentiation is supported for square matrices `x` and integers `y`:
+   * when `y` is nonnegative, `x` may be any square matrix; and when `y` is
+   * negative, `x` must be invertible, and then this function returns
+   * inv(x)^(-y).
    *
    * For cubic roots of negative numbers, the function returns the principal
    * root by default. In order to let the function return the real root,
@@ -49,6 +53,9 @@ var createPow = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
    *    const b = [[1, 2], [4, 3]]
    *    math.pow(b, 2)               // returns Array [[9, 8], [16, 17]]
    *
+   *    const c = [[1, 2], [4, 3]]
+   *    math.pow(c, -1)               // returns Array [[-0.6, 0.4], [0.8, -0.2]]
+   *
    * See also:
    *
    *    multiply, sqrt, cbrt, nthRoot
@@ -144,15 +151,15 @@ var createPow = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
   /**
    * Calculate the power of a 2d array
    * @param {Array} x     must be a 2 dimensional, square matrix
-   * @param {number} y    a positive, integer value
+   * @param {number} y    a integer value (positive if `x` is not invertible)
    * @returns {Array}
    * @private
    */
 
 
   function _powArray(x, y) {
-    if (!(0, _number.isInteger)(y) || y < 0) {
-      throw new TypeError('For A^b, b must be a positive integer (value is ' + y + ')');
+    if (!(0, _number.isInteger)(y)) {
+      throw new TypeError('For A^b, b must be an integer (value is ' + y + ')');
     } // verify that A is a 2 dimensional square matrix
 
 
@@ -166,6 +173,18 @@ var createPow = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
       throw new Error('For A^b, A must be square (size is ' + s[0] + 'x' + s[1] + ')');
     }
 
+    if (y < 0) {
+      try {
+        return _powArray(inv(x), -y);
+      } catch (error) {
+        if (error.message === 'Cannot calculate inverse, determinant is zero') {
+          throw new TypeError('For A^b, when A is not invertible, b must be a positive integer (value is ' + y + ')');
+        }
+
+        throw error;
+      }
+    }
+
     var res = identity(s[0]).valueOf();
     var px = x;