mathjs 8.1.0 → 12.2.1

package/lib/cjs/plain/number/logical.js

@@ -3,33 +3,25 @@
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
+exports.andNumber = andNumber;
 exports.notNumber = notNumber;
 exports.orNumber = orNumber;
 exports.xorNumber = xorNumber;
-exports.andNumber = andNumber;
 var n1 = 'number';
 var n2 = 'number, number';
-
 function notNumber(x) {
   return !x;
 }
-
 notNumber.signature = n1;
-
 function orNumber(x, y) {
   return !!(x || y);
 }
-
 orNumber.signature = n2;
-
 function xorNumber(x, y) {
   return !!x !== !!y;
 }
-
 xorNumber.signature = n2;
-
 function andNumber(x, y) {
   return !!(x && y);
 }
-
 andNumber.signature = n2;

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

@@ -3,37 +3,32 @@
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
+exports.gammaG = void 0;
 exports.gammaNumber = gammaNumber;
-exports.gammaP = exports.gammaG = void 0;
-
+exports.lgammaN = exports.lgammaG = exports.gammaP = void 0;
+exports.lgammaNumber = lgammaNumber;
+exports.lnSqrt2PI = exports.lgammaSeries = void 0;
 var _number = require("../../utils/number.js");
-
 var _product = require("../../utils/product.js");
-
 /* eslint-disable no-loss-of-precision */
+
 function gammaNumber(n) {
   var x;
-
   if ((0, _number.isInteger)(n)) {
     if (n <= 0) {
       return isFinite(n) ? Infinity : NaN;
     }
-
     if (n > 171) {
       return Infinity; // Will overflow
     }
-
     return (0, _product.product)(1, n - 1);
   }
-
   if (n < 0.5) {
     return Math.PI / (Math.sin(Math.PI * n) * gammaNumber(1 - n));
   }
-
   if (n >= 171.35) {
     return Infinity; // will overflow
   }
-
   if (n > 85.0) {
     // Extended Stirling Approx
     var twoN = n * n;
@@ -42,21 +37,49 @@ function gammaNumber(n) {
     var fiveN = fourN * n;
     return Math.sqrt(2 * Math.PI / n) * Math.pow(n / Math.E, n) * (1 + 1 / (12 * n) + 1 / (288 * twoN) - 139 / (51840 * threeN) - 571 / (2488320 * fourN) + 163879 / (209018880 * fiveN) + 5246819 / (75246796800 * fiveN * n));
   }
-
   --n;
   x = gammaP[0];
-
   for (var i = 1; i < gammaP.length; ++i) {
     x += gammaP[i] / (n + i);
   }
-
   var t = n + gammaG + 0.5;
   return Math.sqrt(2 * Math.PI) * Math.pow(t, n + 0.5) * Math.exp(-t) * x;
 }
+gammaNumber.signature = 'number';
+
+// TODO: comment on the variables g and p
+
+var gammaG = exports.gammaG = 4.7421875;
+var gammaP = exports.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];
 
-gammaNumber.signature = 'number'; // TODO: comment on the variables g and p
+// lgamma implementation ref: https://mrob.com/pub/ries/lanczos-gamma.html#code
 
-var gammaG = 4.7421875;
-exports.gammaG = gammaG;
-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];
-exports.gammaP = gammaP;
+// log(2 * pi) / 2
+var lnSqrt2PI = exports.lnSqrt2PI = 0.91893853320467274178;
+var lgammaG = exports.lgammaG = 5; // Lanczos parameter "g"
+var lgammaN = exports.lgammaN = 7; // Range of coefficients "n"
+
+var lgammaSeries = exports.lgammaSeries = [1.000000000190015, 76.18009172947146, -86.50532032941677, 24.01409824083091, -1.231739572450155, 0.1208650973866179e-2, -0.5395239384953e-5];
+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/cjs/plain/number/trigonometry.js

@@ -13,8 +13,8 @@ exports.asecNumber = asecNumber;
 exports.asechNumber = asechNumber;
 exports.asinNumber = asinNumber;
 exports.asinhNumber = asinhNumber;
-exports.atanNumber = atanNumber;
 exports.atan2Number = atan2Number;
+exports.atanNumber = atanNumber;
 exports.atanhNumber = atanhNumber;
 exports.cosNumber = cosNumber;
 exports.coshNumber = coshNumber;
@@ -28,124 +28,85 @@ exports.sinNumber = sinNumber;
 exports.sinhNumber = sinhNumber;
 exports.tanNumber = tanNumber;
 exports.tanhNumber = tanhNumber;
-
 var _number = require("../../utils/number.js");
-
 var n1 = 'number';
 var n2 = 'number, number';
-
 function acosNumber(x) {
   return Math.acos(x);
 }
-
 acosNumber.signature = n1;
-
 function acoshNumber(x) {
   return (0, _number.acosh)(x);
 }
-
 acoshNumber.signature = n1;
-
 function acotNumber(x) {
   return Math.atan(1 / x);
 }
-
 acotNumber.signature = n1;
-
 function acothNumber(x) {
   return isFinite(x) ? (Math.log((x + 1) / x) + Math.log(x / (x - 1))) / 2 : 0;
 }
-
 acothNumber.signature = n1;
-
 function acscNumber(x) {
   return Math.asin(1 / x);
 }
-
 acscNumber.signature = n1;
-
 function acschNumber(x) {
   var xInv = 1 / x;
   return Math.log(xInv + Math.sqrt(xInv * xInv + 1));
 }
-
 acschNumber.signature = n1;
-
 function asecNumber(x) {
   return Math.acos(1 / x);
 }
-
 asecNumber.signature = n1;
-
 function asechNumber(x) {
   var xInv = 1 / x;
   var ret = Math.sqrt(xInv * xInv - 1);
   return Math.log(ret + xInv);
 }
-
 asechNumber.signature = n1;
-
 function asinNumber(x) {
   return Math.asin(x);
 }
-
 asinNumber.signature = n1;
-
 function asinhNumber(x) {
   return (0, _number.asinh)(x);
 }
-
 asinhNumber.signature = n1;
-
 function atanNumber(x) {
   return Math.atan(x);
 }
-
 atanNumber.signature = n1;
-
 function atan2Number(y, x) {
   return Math.atan2(y, x);
 }
-
 atan2Number.signature = n2;
-
 function atanhNumber(x) {
   return (0, _number.atanh)(x);
 }
-
 atanhNumber.signature = n1;
-
 function cosNumber(x) {
   return Math.cos(x);
 }
-
 cosNumber.signature = n1;
-
 function coshNumber(x) {
   return (0, _number.cosh)(x);
 }
-
 coshNumber.signature = n1;
-
 function cotNumber(x) {
   return 1 / Math.tan(x);
 }
-
 cotNumber.signature = n1;
-
 function cothNumber(x) {
   var e = Math.exp(2 * x);
   return (e + 1) / (e - 1);
 }
-
 cothNumber.signature = n1;
-
 function cscNumber(x) {
   return 1 / Math.sin(x);
 }
-
 cscNumber.signature = n1;
-
 function cschNumber(x) {
   // consider values close to zero (+/-)
   if (x === 0) {
@@ -154,41 +115,28 @@ function cschNumber(x) {
     return Math.abs(2 / (Math.exp(x) - Math.exp(-x))) * (0, _number.sign)(x);
   }
 }
-
 cschNumber.signature = n1;
-
 function secNumber(x) {
   return 1 / Math.cos(x);
 }
-
 secNumber.signature = n1;
-
 function sechNumber(x) {
   return 2 / (Math.exp(x) + Math.exp(-x));
 }
-
 sechNumber.signature = n1;
-
 function sinNumber(x) {
   return Math.sin(x);
 }
-
 sinNumber.signature = n1;
-
 function sinhNumber(x) {
   return (0, _number.sinh)(x);
 }
-
 sinhNumber.signature = n1;
-
 function tanNumber(x) {
   return Math.tan(x);
 }
-
 tanNumber.signature = n1;
-
 function tanhNumber(x) {
   return (0, _number.tanh)(x);
 }
-
 tanhNumber.signature = n1;

package/lib/cjs/plain/number/utils.js

@@ -4,41 +4,29 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.isIntegerNumber = isIntegerNumber;
+exports.isNaNNumber = isNaNNumber;
 exports.isNegativeNumber = isNegativeNumber;
 exports.isPositiveNumber = isPositiveNumber;
 exports.isZeroNumber = isZeroNumber;
-exports.isNaNNumber = isNaNNumber;
-
 var _number = require("../../utils/number.js");
-
 var n1 = 'number';
-
 function isIntegerNumber(x) {
   return (0, _number.isInteger)(x);
 }
-
 isIntegerNumber.signature = n1;
-
 function isNegativeNumber(x) {
   return x < 0;
 }
-
 isNegativeNumber.signature = n1;
-
 function isPositiveNumber(x) {
   return x > 0;
 }
-
 isPositiveNumber.signature = n1;
-
 function isZeroNumber(x) {
   return x === 0;
 }
-
 isZeroNumber.signature = n1;
-
 function isNaNNumber(x) {
   return Number.isNaN(x);
 }
-
 isNaNNumber.signature = n1;

package/lib/cjs/type/bignumber/BigNumber.js

@@ -1,59 +1,51 @@
 "use strict";
 
+var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.createBigNumberClass = void 0;
-
 var _decimal = _interopRequireDefault(require("decimal.js"));
-
 var _factory = require("../../utils/factory.js");
-
-function _interopRequireDefault(obj) { return obj && obj.__esModule ? obj : { "default": obj }; }
-
 var name = 'BigNumber';
 var dependencies = ['?on', 'config'];
-var createBigNumberClass = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var createBigNumberClass = exports.createBigNumberClass = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var on = _ref.on,
-      config = _ref.config;
-  var EUCLID = 9; // Use euclidian division for mod calculation
-
+    config = _ref.config;
   var BigNumber = _decimal["default"].clone({
     precision: config.precision,
-    modulo: EUCLID
+    modulo: _decimal["default"].EUCLID
   });
+  BigNumber.prototype = Object.create(BigNumber.prototype);
+
   /**
    * Attach type information
    */
-
-
   BigNumber.prototype.type = 'BigNumber';
   BigNumber.prototype.isBigNumber = true;
+
   /**
    * Get a JSON representation of a BigNumber containing
    * type information
    * @returns {Object} Returns a JSON object structured as:
    *                   `{"mathjs": "BigNumber", "value": "0.2"}`
    */
-
   BigNumber.prototype.toJSON = function () {
     return {
       mathjs: 'BigNumber',
       value: this.toString()
     };
   };
+
   /**
    * Instantiate a BigNumber from a JSON object
    * @param {Object} json  a JSON object structured as:
    *                       `{"mathjs": "BigNumber", "value": "0.2"}`
    * @return {BigNumber}
    */
-
-
   BigNumber.fromJSON = function (json) {
     return new BigNumber(json.value);
   };
-
   if (on) {
     // listen for changed in the configuration, automatically apply changed precision
     on('config', function (curr, prev) {
@@ -64,9 +56,7 @@ var createBigNumberClass = /* #__PURE__ */(0, _factory.factory)(name, dependenci
       }
     });
   }
-
   return BigNumber;
 }, {
   isClass: true
-});
-exports.createBigNumberClass = createBigNumberClass;
+});

package/lib/cjs/type/bignumber/function/bignumber.js

@@ -4,17 +4,13 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.createBignumber = void 0;
-
 var _factory = require("../../../utils/factory.js");
-
 var _collection = require("../../../utils/collection.js");
-
 var name = 'bignumber';
 var dependencies = ['typed', 'BigNumber'];
-var createBignumber = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var createBignumber = exports.createBignumber = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
-      BigNumber = _ref.BigNumber;
-
+    BigNumber = _ref.BigNumber;
   /**
    * Create a BigNumber, which can store numbers with arbitrary precision.
    * When a matrix is provided, all elements will be converted to BigNumber.
@@ -49,21 +45,45 @@ var createBignumber = /* #__PURE__ */(0, _factory.factory)(name, dependencies, f
       return new BigNumber(x + '');
     },
     string: function string(x) {
+      var wordSizeSuffixMatch = x.match(/(0[box][0-9a-fA-F]*)i([0-9]*)/);
+      if (wordSizeSuffixMatch) {
+        // x has a word size suffix
+        var size = wordSizeSuffixMatch[2];
+        var n = BigNumber(wordSizeSuffixMatch[1]);
+        var twoPowSize = new BigNumber(2).pow(Number(size));
+        if (n.gt(twoPowSize.sub(1))) {
+          throw new SyntaxError("String \"".concat(x, "\" is out of range"));
+        }
+        var twoPowSizeSubOne = new BigNumber(2).pow(Number(size) - 1);
+        if (n.gte(twoPowSizeSubOne)) {
+          return n.sub(twoPowSize);
+        } else {
+          return n;
+        }
+      }
       return new BigNumber(x);
     },
     BigNumber: function BigNumber(x) {
       // we assume a BigNumber is immutable
       return x;
     },
+    Unit: typed.referToSelf(function (self) {
+      return function (x) {
+        var clone = x.clone();
+        clone.value = self(x.value);
+        return clone;
+      };
+    }),
     Fraction: function Fraction(x) {
       return new BigNumber(x.n).div(x.d).times(x.s);
     },
     "null": function _null(x) {
       return new BigNumber(0);
     },
-    'Array | Matrix': function ArrayMatrix(x) {
-      return (0, _collection.deepMap)(x, this);
-    }
+    'Array | Matrix': typed.referToSelf(function (self) {
+      return function (x) {
+        return (0, _collection.deepMap)(x, self);
+      };
+    })
   });
-});
-exports.createBignumber = createBignumber;
+});

package/lib/cjs/type/boolean.js

@@ -4,16 +4,12 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.createBoolean = void 0;
-
 var _factory = require("../utils/factory.js");
-
 var _collection = require("../utils/collection.js");
-
 var name = 'boolean';
 var dependencies = ['typed'];
-var createBoolean = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+var createBoolean = exports.createBoolean = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed;
-
   /**
    * Create a boolean or convert a string or number to a boolean.
    * In case of a number, `true` is returned for non-zero numbers, and `false` in
@@ -60,25 +56,23 @@ var createBoolean = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
     string: function string(x) {
       // try case insensitive
       var lcase = x.toLowerCase();
-
       if (lcase === 'true') {
         return true;
       } else if (lcase === 'false') {
         return false;
-      } // test whether value is a valid number
-
+      }
 
+      // test whether value is a valid number
       var num = Number(x);
-
       if (x !== '' && !isNaN(num)) {
         return !!num;
       }
-
       throw new Error('Cannot convert "' + x + '" to a boolean');
     },
-    'Array | Matrix': function ArrayMatrix(x) {
-      return (0, _collection.deepMap)(x, this);
-    }
+    'Array | Matrix': typed.referToSelf(function (self) {
+      return function (x) {
+        return (0, _collection.deepMap)(x, self);
+      };
+    })
   });
-});
-exports.createBoolean = createBoolean;
+});