mathjs 14.1.0 → 14.2.0

package/lib/cjs/entry/pureFunctionsNumber.generated.js

@@ -3,8 +3,8 @@
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.median = exports.mean = exports.max = exports.matrix = exports.map = exports.mad = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.index = exports.hypot = exports.hasNumericValue = exports.gcd = exports.gamma = exports.format = exports.forEach = exports.floor = exports.fix = exports.filter = exports.factorial = exports.expm1 = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.e = exports.divideScalar = exports.divide = exports.deepEqual = exports.cumsum = exports.cube = exports.csch = exports.csc = exports.coth = exports.cot = exports.cosh = exports.cos = exports.corr = exports.composition = exports.compareText = exports.compareNatural = exports.compare = exports.combinationsWithRep = exports.combinations = exports.clone = exports.ceil = exports.cbrt = exports.catalan = exports.boolean = exports.bitXor = exports.bitOr = exports.bitNot = exports.bitAnd = exports.bigint = exports.bellNumbers = exports.atanh = exports.atan2 = exports.atan = exports.asinh = exports.asin = exports.asech = exports.asec = exports.apply = exports.and = exports.addScalar = exports.add = exports.acsch = exports.acsc = exports.acoth = exports.acot = exports.acosh = exports.acos = exports.abs = exports._true = exports._null = exports._false = exports._NaN = exports._Infinity = exports.SQRT2 = exports.SQRT1_2 = exports.ResultSet = exports.Range = exports.LOG2E = exports.LOG10E = exports.LN2 = exports.LN10 = void 0;
-exports.zeta = exports.xor = exports.xgcd = exports.version = exports.variance = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.tau = exports.tanh = exports.tan = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.std = exports.square = exports.sqrt = exports.smallerEq = exports.smaller = exports.size = exports.sinh = exports.sin = exports.sign = exports.sech = exports.sec = exports.round = exports.rightLogShift = exports.rightArithShift = exports.replacer = exports.range = exports.randomInt = exports.random = exports.quantileSeq = exports.prod = exports.print = exports.pow = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.numeric = exports.number = exports.nthRoot = exports.not = exports.norm = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.mode = exports.mod = exports.min = void 0;
+exports.mean = exports.max = exports.matrix = exports.mapSlices = exports.map = exports.mad = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.index = exports.hypot = exports.hasNumericValue = exports.gcd = exports.gamma = exports.format = exports.forEach = exports.floor = exports.fix = exports.filter = exports.factorial = exports.expm1 = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.e = exports.divideScalar = exports.divide = exports.deepEqual = exports.cumsum = exports.cube = exports.csch = exports.csc = exports.coth = exports.cot = exports.cosh = exports.cos = exports.corr = exports.composition = exports.compareText = exports.compareNatural = exports.compare = exports.combinationsWithRep = exports.combinations = exports.clone = exports.ceil = exports.cbrt = exports.catalan = exports.boolean = exports.bitXor = exports.bitOr = exports.bitNot = exports.bitAnd = exports.bigint = exports.bellNumbers = exports.atanh = exports.atan2 = exports.atan = exports.asinh = exports.asin = exports.asech = exports.asec = exports.apply = exports.and = exports.addScalar = exports.add = exports.acsch = exports.acsc = exports.acoth = exports.acot = exports.acosh = exports.acos = exports.abs = exports._true = exports._null = exports._false = exports._NaN = exports._Infinity = exports.SQRT2 = exports.SQRT1_2 = exports.ResultSet = exports.Range = exports.LOG2E = exports.LOG10E = exports.LN2 = exports.LN10 = void 0;
+exports.zeta = exports.xor = exports.xgcd = exports.version = exports.variance = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.tau = exports.tanh = exports.tan = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.std = exports.square = exports.sqrt = exports.smallerEq = exports.smaller = exports.size = exports.sinh = exports.sin = exports.sign = exports.sech = exports.sec = exports.round = exports.rightLogShift = exports.rightArithShift = exports.replacer = exports.range = exports.randomInt = exports.random = exports.quantileSeq = exports.prod = exports.print = exports.pow = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.numeric = exports.number = exports.nthRoot = exports.not = exports.norm = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.mode = exports.mod = exports.min = exports.median = void 0;
 var _configReadonly = require("./configReadonly.js");
 var _factoriesNumber = require("../factoriesNumber.js");
 /**
@@ -292,6 +292,11 @@ const lcm = exports.lcm = /* #__PURE__ */(0, _factoriesNumber.createLcm)({
 const log10 = exports.log10 = /* #__PURE__ */(0, _factoriesNumber.createLog10)({
   typed
 });
+const mapSlices = exports.mapSlices = /* #__PURE__ */(0, _factoriesNumber.createMapSlices)({
+  isInteger,
+  typed
+});
+const apply = exports.apply = mapSlices;
 const multiplyScalar = exports.multiplyScalar = /* #__PURE__ */(0, _factoriesNumber.createMultiplyScalar)({
   typed
 });
@@ -323,11 +328,7 @@ const unaryMinus = exports.unaryMinus = /* #__PURE__ */(0, _factoriesNumber.crea
 const acosh = exports.acosh = /* #__PURE__ */(0, _factoriesNumber.createAcosh)({
   typed
 });
-const apply = exports.apply = /* #__PURE__ */(0, _factoriesNumber.createApply)({
-  isInteger,
-  typed
-});
-const asinh = exports.asinh = /* #__PURE__ */(0, _factoriesNumber.createAsinh)({
+const atan2 = exports.atan2 = /* #__PURE__ */(0, _factoriesNumber.createAtan2)({
   typed
 });
 const bitOr = exports.bitOr = /* #__PURE__ */(0, _factoriesNumber.createBitOr)({
@@ -396,9 +397,9 @@ const subtractScalar = exports.subtractScalar = /* #__PURE__ */(0, _factoriesNum
 const _true = exports._true = /* #__PURE__ */(0, _factoriesNumber.createTrue)({});
 const variance = exports.variance = /* #__PURE__ */(0, _factoriesNumber.createVariance)({
   add,
-  apply,
   divide,
   isNaN,
+  mapSlices,
   multiply,
   subtract,
   typed
@@ -422,9 +423,6 @@ const zeta = exports.zeta = /* #__PURE__ */(0, _factoriesNumber.createZeta)({
 const acsch = exports.acsch = /* #__PURE__ */(0, _factoriesNumber.createAcsch)({
   typed
 });
-const atan2 = exports.atan2 = /* #__PURE__ */(0, _factoriesNumber.createAtan2)({
-  typed
-});
 const catalan = exports.catalan = /* #__PURE__ */(0, _factoriesNumber.createCatalan)({
   addScalar,
   combinations,
@@ -496,6 +494,7 @@ const quantileSeq = exports.quantileSeq = /* #__PURE__ */(0, _factoriesNumber.cr
   divide,
   isInteger,
   larger,
+  mapSlices,
   multiply,
   partitionSelect,
   smaller,
@@ -526,6 +525,9 @@ const stirlingS2 = exports.stirlingS2 = /* #__PURE__ */(0, _factoriesNumber.crea
   subtractScalar,
   typed
 });
+const asinh = exports.asinh = /* #__PURE__ */(0, _factoriesNumber.createAsinh)({
+  typed
+});
 const ceil = exports.ceil = /* #__PURE__ */(0, _factoriesNumber.createCeil)({
   config: _configReadonly.config,
   round,

package/lib/cjs/expression/embeddedDocs/embeddedDocs.js

@@ -117,6 +117,7 @@ var _and = require("./function/logical/and.js");
 var _not = require("./function/logical/not.js");
 var _or = require("./function/logical/or.js");
 var _xor = require("./function/logical/xor.js");
+var _mapSlices = require("./function/matrix/mapSlices.js");
 var _column = require("./function/matrix/column.js");
 var _concat = require("./function/matrix/concat.js");
 var _count = require("./function/matrix/count.js");
@@ -578,6 +579,7 @@ const embeddedDocs = exports.embeddedDocs = {
   or: _or.orDocs,
   xor: _xor.xorDocs,
   // functions - matrix
+  mapSlices: _mapSlices.mapSlicesDocs,
   concat: _concat.concatDocs,
   count: _count.countDocs,
   cross: _cross.crossDocs,

package/lib/cjs/expression/embeddedDocs/function/matrix/mapSlices.js

@@ -0,0 +1,17 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.mapSlicesDocs = void 0;
+const mapSlicesDocs = exports.mapSlicesDocs = {
+  name: 'mapSlices',
+  category: 'Matrix',
+  syntax: ['mapSlices(A, dim, callback)'],
+  description: 'Generate a matrix one dimension less than A by applying callback to ' + 'each slice of A along dimension dim.',
+  examples: ['A = [[1, 2], [3, 4]]', 'mapSlices(A, 1, sum)',
+  // returns [4, 6]
+  'mapSlices(A, 2, product)' // returns [2, 12]
+  ],
+  seealso: ['map', 'forEach']
+};

package/lib/cjs/expression/transform/index.transform.js

@@ -38,14 +38,14 @@ const createIndexTransform = exports.createIndexTransform = /* #__PURE__ */(0, _
             return v - 1;
           });
         }
-      } else if ((0, _is.isNumber)(arg)) {
+      } else if ((0, _is.isNumber)(arg) || (0, _is.isBigInt)(arg)) {
         arg--;
       } else if ((0, _is.isBigNumber)(arg)) {
         arg = arg.toNumber() - 1;
       } else if (typeof arg === 'string') {
         // leave as is
       } else {
-        throw new TypeError('Dimension must be an Array, Matrix, number, string, or Range');
+        throw new TypeError('Dimension must be an Array, Matrix, number, bigint, string, or Range');
       }
       args[i] = arg;
     }

package/lib/cjs/expression/transform/{apply.transform.js → mapSlices.transform.js}

@@ -3,33 +3,33 @@
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
-exports.createApplyTransform = void 0;
+exports.createMapSlicesTransform = void 0;
 var _errorTransform = require("./utils/errorTransform.js");
 var _factory = require("../../utils/factory.js");
-var _apply = require("../../function/matrix/apply.js");
+var _mapSlices = require("../../function/matrix/mapSlices.js");
 var _is = require("../../utils/is.js");
-const name = 'apply';
+const name = 'mapSlices';
 const dependencies = ['typed', 'isInteger'];
 
 /**
- * Attach a transform function to math.apply
+ * Attach a transform function to math.mapSlices
  * Adds a property transform containing the transform function.
  *
- * This transform changed the last `dim` parameter of function apply
+ * This transform changed the last `dim` parameter of function mapSlices
  * from one-based to zero based
  */
-const createApplyTransform = exports.createApplyTransform = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+const createMapSlicesTransform = exports.createMapSlicesTransform = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
     isInteger
   } = _ref;
-  const apply = (0, _apply.createApply)({
+  const mapSlices = (0, _mapSlices.createMapSlices)({
     typed,
     isInteger
   });
 
   // @see: comment of concat itself
-  return typed('apply', {
+  return typed('mapSlices', {
     '...any': function (args) {
       // change dim from one-based to zero-based
       const dim = args[1];
@@ -39,12 +39,13 @@ const createApplyTransform = exports.createApplyTransform = /* #__PURE__ */(0, _
         args[1] = dim.minus(1);
       }
       try {
-        return apply.apply(null, args);
+        return mapSlices.apply(null, args);
       } catch (err) {
         throw (0, _errorTransform.errorTransform)(err);
       }
     }
   });
 }, {
-  isTransformFunction: true
+  isTransformFunction: true,
+  ..._mapSlices.createMapSlices.meta
 });

package/lib/cjs/expression/transform/quantileSeq.transform.js

@@ -8,7 +8,7 @@ var _factory = require("../../utils/factory.js");
 var _quantileSeq = require("../../function/statistics/quantileSeq.js");
 var _lastDimToZeroBase = require("./utils/lastDimToZeroBase.js");
 const name = 'quantileSeq';
-const dependencies = ['typed', 'bignumber', 'add', 'subtract', 'divide', 'multiply', 'partitionSelect', 'compare', 'isInteger', 'smaller', 'smallerEq', 'larger'];
+const dependencies = ['typed', 'bignumber', 'add', 'subtract', 'divide', 'multiply', 'partitionSelect', 'compare', 'isInteger', 'smaller', 'smallerEq', 'larger', 'mapSlices'];
 
 /**
  * Attach a transform function to math.quantileSeq
@@ -30,7 +30,8 @@ const createQuantileSeqTransform = exports.createQuantileSeqTransform = /* #__PU
     isInteger,
     smaller,
     smallerEq,
-    larger
+    larger,
+    mapSlices
   } = _ref;
   const quantileSeq = (0, _quantileSeq.createQuantileSeq)({
     typed,
@@ -44,7 +45,8 @@ const createQuantileSeqTransform = exports.createQuantileSeqTransform = /* #__PU
     isInteger,
     smaller,
     smallerEq,
-    larger
+    larger,
+    mapSlices
   });
   return typed('quantileSeq', {
     'Array | Matrix, number | BigNumber': quantileSeq,

package/lib/cjs/expression/transform/variance.transform.js

@@ -9,7 +9,7 @@ var _errorTransform = require("./utils/errorTransform.js");
 var _variance = require("../../function/statistics/variance.js");
 var _lastDimToZeroBase = require("./utils/lastDimToZeroBase.js");
 const name = 'variance';
-const dependencies = ['typed', 'add', 'subtract', 'multiply', 'divide', 'apply', 'isNaN'];
+const dependencies = ['typed', 'add', 'subtract', 'multiply', 'divide', 'mapSlices', 'isNaN'];
 
 /**
  * Attach a transform function to math.var
@@ -25,7 +25,7 @@ const createVarianceTransform = exports.createVarianceTransform = /* #__PURE__ *
     subtract,
     multiply,
     divide,
-    apply,
+    mapSlices,
     isNaN
   } = _ref;
   const variance = (0, _variance.createVariance)({
@@ -34,7 +34,7 @@ const createVarianceTransform = exports.createVarianceTransform = /* #__PURE__ *
     subtract,
     multiply,
     divide,
-    apply,
+    mapSlices,
     isNaN
   });
   return typed(name, {

package/lib/cjs/factoriesAny.js

@@ -75,18 +75,6 @@ Object.defineProperty(exports, "createAndTransform", {
     return _andTransform.createAndTransform;
   }
 });
-Object.defineProperty(exports, "createApply", {
-  enumerable: true,
-  get: function () {
-    return _apply.createApply;
-  }
-});
-Object.defineProperty(exports, "createApplyTransform", {
-  enumerable: true,
-  get: function () {
-    return _applyTransform.createApplyTransform;
-  }
-});
 Object.defineProperty(exports, "createArg", {
   enumerable: true,
   get: function () {
@@ -1149,6 +1137,18 @@ Object.defineProperty(exports, "createMap", {
     return _map.createMap;
   }
 });
+Object.defineProperty(exports, "createMapSlices", {
+  enumerable: true,
+  get: function () {
+    return _mapSlices.createMapSlices;
+  }
+});
+Object.defineProperty(exports, "createMapSlicesTransform", {
+  enumerable: true,
+  get: function () {
+    return _mapSlicesTransform.createMapSlicesTransform;
+  }
+});
 Object.defineProperty(exports, "createMapTransform", {
   enumerable: true,
   get: function () {
@@ -2209,7 +2209,7 @@ var _splitUnit = require("./type/unit/function/splitUnit.js");
 var _unaryMinus = require("./function/arithmetic/unaryMinus.js");
 var _unaryPlus = require("./function/arithmetic/unaryPlus.js");
 var _abs = require("./function/arithmetic/abs.js");
-var _apply = require("./function/matrix/apply.js");
+var _mapSlices = require("./function/matrix/mapSlices.js");
 var _addScalar = require("./function/arithmetic/addScalar.js");
 var _subtractScalar = require("./function/arithmetic/subtractScalar.js");
 var _cbrt = require("./function/arithmetic/cbrt.js");
@@ -2446,7 +2446,7 @@ var _reviver = require("./json/reviver.js");
 var _replacer = require("./json/replacer.js");
 var _constants = require("./constants.js");
 var _physicalConstants = require("./type/unit/physicalConstants.js");
-var _applyTransform = require("./expression/transform/apply.transform.js");
+var _mapSlicesTransform = require("./expression/transform/mapSlices.transform.js");
 var _columnTransform = require("./expression/transform/column.transform.js");
 var _filterTransform = require("./expression/transform/filter.transform.js");
 var _forEachTransform = require("./expression/transform/forEach.transform.js");

package/lib/cjs/factoriesNumber.js

@@ -11,18 +11,6 @@ Object.defineProperty(exports, "createAccessorNode", {
   }
 });
 exports.createAnd = exports.createAddScalar = exports.createAdd = exports.createAcsch = exports.createAcsc = exports.createAcoth = exports.createAcot = exports.createAcosh = exports.createAcos = void 0;
-Object.defineProperty(exports, "createApply", {
-  enumerable: true,
-  get: function () {
-    return _apply.createApply;
-  }
-});
-Object.defineProperty(exports, "createApplyTransform", {
-  enumerable: true,
-  get: function () {
-    return _applyTransform.createApplyTransform;
-  }
-});
 Object.defineProperty(exports, "createArrayNode", {
   enumerable: true,
   get: function () {
@@ -366,7 +354,8 @@ Object.defineProperty(exports, "createLargerEq", {
     return _largerEq.createLargerEqNumber;
   }
 });
-exports.createLog2 = exports.createLog1p = exports.createLog10 = exports.createLog = exports.createLgamma = exports.createLeftShift = exports.createLcm = void 0;
+exports.createLeftShift = exports.createLcm = void 0;
+exports.createLog2 = exports.createLog1p = exports.createLog10 = exports.createLog = exports.createLgamma = void 0;
 Object.defineProperty(exports, "createMad", {
   enumerable: true,
   get: function () {
@@ -379,6 +368,18 @@ Object.defineProperty(exports, "createMap", {
     return _map.createMap;
   }
 });
+Object.defineProperty(exports, "createMapSlices", {
+  enumerable: true,
+  get: function () {
+    return _mapSlices.createMapSlices;
+  }
+});
+Object.defineProperty(exports, "createMapSlicesTransform", {
+  enumerable: true,
+  get: function () {
+    return _mapSlicesTransform.createMapSlicesTransform;
+  }
+});
 Object.defineProperty(exports, "createMapTransform", {
   enumerable: true,
   get: function () {
@@ -845,7 +846,7 @@ var _parse = require("./expression/parse.js");
 var _compile = require("./expression/function/compile.js");
 var _evaluate = require("./expression/function/evaluate.js");
 var _Parser = require("./expression/Parser.js");
-var _apply = require("./function/matrix/apply.js");
+var _mapSlices = require("./function/matrix/mapSlices.js");
 var _filter = require("./function/matrix/filter.js");
 var _forEach = require("./function/matrix/forEach.js");
 var _map = require("./function/matrix/map.js");
@@ -888,7 +889,7 @@ var _std = require("./function/statistics/std.js");
 var _corr = require("./function/statistics/corr.js");
 var _format = require("./function/string/format.js");
 var _print = require("./function/string/print.js");
-var _applyTransform = require("./expression/transform/apply.transform.js");
+var _mapSlicesTransform = require("./expression/transform/mapSlices.transform.js");
 var _filterTransform = require("./expression/transform/filter.transform.js");
 var _forEachTransform = require("./expression/transform/forEach.transform.js");
 var _mapTransform = require("./expression/transform/map.transform.js");

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

@@ -15,29 +15,34 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matAlgo14xDs = require("../../type/matrix/utils/matAlgo14xDs.js");
 const name = 'ceil';
 const dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar', 'zeros', 'DenseMatrix'];
+const bigTen = new _decimal.default(10);
 const createCeilNumber = exports.createCeilNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'config', 'round'], _ref => {
   let {
     typed,
     config,
     round
   } = _ref;
+  function _ceilNumber(x) {
+    // See ./floor.js _floorNumber for rationale here
+    const c = Math.ceil(x);
+    const r = round(x);
+    if (c === r) return c;
+    if ((0, _number.nearlyEqual)(x, r, config.relTol, config.absTol) && !(0, _number.nearlyEqual)(x, c, config.relTol, config.absTol)) {
+      return r;
+    }
+    return c;
+  }
   return typed(name, {
-    number: function (x) {
-      if ((0, _number.nearlyEqual)(x, round(x), config.relTol, config.absTol)) {
-        return round(x);
-      } else {
-        return Math.ceil(x);
-      }
-    },
+    number: _ceilNumber,
     'number, number': function (x, n) {
-      if ((0, _number.nearlyEqual)(x, round(x, n), config.relTol, config.absTol)) {
-        return round(x, n);
-      } else {
-        let [number, exponent] = `${x}e`.split('e');
-        const result = Math.ceil(Number(`${number}e${Number(exponent) + n}`));
-        [number, exponent] = `${result}e`.split('e');
-        return Number(`${number}e${Number(exponent) - n}`);
+      if (!(0, _number.isInteger)(n)) {
+        throw new RangeError('number of decimals in function ceil must be an integer');
+      }
+      if (n < 0 || n > 15) {
+        throw new RangeError('number of decimals in ceil number must be in range 0-15');
       }
+      const shift = 10 ** n;
+      return _ceilNumber(x * shift) / shift;
     }
   });
 });
@@ -67,6 +72,15 @@ const createCeil = exports.createCeil = /* #__PURE__ */(0, _factory.factory)(nam
     config,
     round
   });
+  function _bigCeil(x) {
+    // see ./floor.js _floorNumber for rationale
+    const bne = (a, b) => (0, _nearlyEqual.nearlyEqual)(a, b, config.relTol, config.absTol);
+    const c = x.ceil();
+    const r = round(x);
+    if (c.eq(r)) return c;
+    if (bne(x, r) && !bne(x, c)) return r;
+    return c;
+  }
   /**
    * Round a value towards plus infinity
    * If `x` is complex, both real and imaginary part are rounded towards plus infinity.
@@ -125,20 +139,14 @@ const createCeil = exports.createCeil = /* #__PURE__ */(0, _factory.factory)(nam
     'Complex, BigNumber': function (x, n) {
       return x.ceil(n.toNumber());
     },
-    BigNumber: function (x) {
-      if ((0, _nearlyEqual.nearlyEqual)(x, round(x), config.relTol, config.absTol)) {
-        return round(x);
-      } else {
-        return x.ceil();
-      }
-    },
+    BigNumber: _bigCeil,
     'BigNumber, BigNumber': function (x, n) {
-      if ((0, _nearlyEqual.nearlyEqual)(x, round(x, n), config.relTol, config.absTol)) {
-        return round(x, n);
-      } else {
-        return x.toDecimalPlaces(n.toNumber(), _decimal.default.ROUND_CEIL);
-      }
+      const shift = bigTen.pow(n);
+      return _bigCeil(x.mul(shift)).div(shift);
     },
+    bigint: b => b,
+    'bigint, number': (b, _dummy) => b,
+    'bigint, BigNumber': (b, _dummy) => b,
     Fraction: function (x) {
       return x.ceil();
     },

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

@@ -112,6 +112,9 @@ const createFix = exports.createFix = /* #__PURE__ */(0, _factory.factory)(name,
     'BigNumber, number | BigNumber': function (x, n) {
       return x.isNegative() ? ceil(x, n) : floor(x, n);
     },
+    bigint: b => b,
+    'bigint, number': (b, _dummy) => b,
+    'bigint, BigNumber': (b, _dummy) => b,
     Fraction: function (x) {
       return x.s < 0n ? x.ceil() : x.floor();
     },

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

@@ -15,29 +15,41 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matAlgo14xDs = require("../../type/matrix/utils/matAlgo14xDs.js");
 const name = 'floor';
 const dependencies = ['typed', 'config', 'round', 'matrix', 'equalScalar', 'zeros', 'DenseMatrix'];
+const bigTen = new _decimal.default(10);
 const createFloorNumber = exports.createFloorNumber = /* #__PURE__ */(0, _factory.factory)(name, ['typed', 'config', 'round'], _ref => {
   let {
     typed,
     config,
     round
   } = _ref;
+  function _floorNumber(x) {
+    // First, if the floor and the round are identical we can be
+    // quite comfortable that is the best answer:
+    const f = Math.floor(x);
+    const r = round(x);
+    if (f === r) return f;
+    // OK, they are different. If x is truly distinct from f but
+    // appears indistinguishable from r, presume it really is just
+    // the integer r with rounding/computation error, and return that
+    if ((0, _number.nearlyEqual)(x, r, config.relTol, config.absTol) && !(0, _number.nearlyEqual)(x, f, config.relTol, config.absTol)) {
+      return r;
+    }
+    // Otherwise (x distinct from both r and f, or indistinguishable from
+    // both r and f) may as well just return f, as that's the best
+    // candidate we can discern:
+    return f;
+  }
   return typed(name, {
-    number: function (x) {
-      if ((0, _number.nearlyEqual)(x, round(x), config.relTol, config.absTol)) {
-        return round(x);
-      } else {
-        return Math.floor(x);
-      }
-    },
+    number: _floorNumber,
     'number, number': function (x, n) {
-      if ((0, _number.nearlyEqual)(x, round(x, n), config.relTol, config.absTol)) {
-        return round(x, n);
-      } else {
-        let [number, exponent] = `${x}e`.split('e');
-        const result = Math.floor(Number(`${number}e${Number(exponent) + n}`));
-        [number, exponent] = `${result}e`.split('e');
-        return Number(`${number}e${Number(exponent) - n}`);
+      if (!(0, _number.isInteger)(n)) {
+        throw new RangeError('number of decimals in function floor must be an integer');
+      }
+      if (n < 0 || n > 15) {
+        throw new RangeError('number of decimals in floor number must be in range 0 - 15');
       }
+      const shift = 10 ** n;
+      return _floorNumber(x * shift) / shift;
     }
   });
 });
@@ -67,6 +79,15 @@ const createFloor = exports.createFloor = /* #__PURE__ */(0, _factory.factory)(n
     config,
     round
   });
+  function _bigFloor(x) {
+    // see _floorNumber above for rationale
+    const bne = (a, b) => (0, _nearlyEqual.nearlyEqual)(a, b, config.relTol, config.absTol);
+    const f = x.floor();
+    const r = round(x);
+    if (f.eq(r)) return f;
+    if (bne(x, r) && !bne(x, f)) return r;
+    return f;
+  }
   /**
    * Round a value towards minus infinity.
    * For matrices, the function is evaluated element wise.
@@ -128,20 +149,14 @@ const createFloor = exports.createFloor = /* #__PURE__ */(0, _factory.factory)(n
     'Complex, BigNumber': function (x, n) {
       return x.floor(n.toNumber());
     },
-    BigNumber: function (x) {
-      if ((0, _nearlyEqual.nearlyEqual)(x, round(x), config.relTol, config.absTol)) {
-        return round(x);
-      } else {
-        return x.floor();
-      }
-    },
+    BigNumber: _bigFloor,
     'BigNumber, BigNumber': function (x, n) {
-      if ((0, _nearlyEqual.nearlyEqual)(x, round(x, n), config.relTol, config.absTol)) {
-        return round(x, n);
-      } else {
-        return x.toDecimalPlaces(n.toNumber(), _decimal.default.ROUND_FLOOR);
-      }
+      const shift = bigTen.pow(n);
+      return _bigFloor(x.mul(shift)).div(shift);
     },
+    bigint: b => b,
+    'bigint, number': (b, _dummy) => b,
+    'bigint, BigNumber': (b, _dummy) => b,
     Fraction: function (x) {
       return x.floor();
     },

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

@@ -143,6 +143,10 @@ const createRound = exports.createRound = /* #__PURE__ */(0, _factory.factory)(n
       const xSelected = (0, _nearlyEqual.nearlyEqual)(x, xEpsilon, config.relTol, config.absTol) ? xEpsilon : x;
       return xSelected.toDecimalPlaces(n.toNumber());
     },
+    // bigints can't be rounded
+    bigint: b => b,
+    'bigint, number': (b, _dummy) => b,
+    'bigint, BigNumber': (b, _dummy) => b,
     Fraction: function (x) {
       return x.round();
     },

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

@@ -26,7 +26,7 @@ const createCtranspose = exports.createCtranspose = /* #__PURE__ */(0, _factory.
    * Examples:
    *
    *     const A = [[1, 2, 3], [4, 5, math.complex(6,7)]]
-   *     math.ctranspose(A)               // returns [[1, 4], [2, 5], [3, {re:6,im:7}]]
+   *     math.ctranspose(A)      // returns [[1, 4], [2, 5], [3, {re:6,im:-7}]]
    *
    * See also:
    *