mathjs 14.9.1 → 15.1.0

package/lib/esm/type/matrix/MatrixIndex.js

@@ -1,4 +1,4 @@
-import { isArray, isMatrix, isRange } from '../../utils/is.js';
+import { isArray, isMatrix, isRange, isNumber, isString } from '../../utils/is.js';
 import { clone } from '../../utils/object.js';
 import { isInteger } from '../../utils/number.js';
 import { factory } from '../../utils/factory.js';
@@ -31,7 +31,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
    * @Constructor Index
    * @param {...*} ranges
    */
-  function Index(ranges) {
+  function Index() {
     if (!(this instanceof Index)) {
       throw new SyntaxError('Constructor must be called with the new operator');
     }
@@ -39,7 +39,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     this._sourceSize = [];
     this._isScalar = true;
     for (var i = 0, ii = arguments.length; i < ii; i++) {
-      var arg = arguments[i];
+      var arg = i < 0 || arguments.length <= i ? undefined : arguments[i];
       var argIsArray = isArray(arg);
       var argIsMatrix = isMatrix(arg);
       var argType = typeof arg;
@@ -50,6 +50,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
       } else if (argIsArray || argIsMatrix) {
         // create matrix
         var m = void 0;
+        this._isScalar = false;
         if (getMatrixDataType(arg) === 'boolean') {
           if (argIsArray) m = _createImmutableMatrix(_booleansArrayToNumbersForIndex(arg).valueOf());
           if (argIsMatrix) m = _createImmutableMatrix(_booleansArrayToNumbersForIndex(arg._data).valueOf());
@@ -58,16 +59,10 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
           m = _createImmutableMatrix(arg.valueOf());
         }
         this._dimensions.push(m);
-        // size
-        var size = m.size();
-        // scalar
-        if (size.length !== 1 || size[0] !== 1 || sourceSize !== null) {
-          this._isScalar = false;
-        }
       } else if (argType === 'number') {
-        this._dimensions.push(_createImmutableMatrix([arg]));
+        this._dimensions.push(arg);
       } else if (argType === 'bigint') {
-        this._dimensions.push(_createImmutableMatrix([Number(arg)]));
+        this._dimensions.push(Number(arg));
       } else if (argType === 'string') {
         // object property (arguments.count should be 1)
         this._dimensions.push(arg);
@@ -87,12 +82,15 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
   function _createImmutableMatrix(arg) {
     // loop array elements
     for (var i = 0, l = arg.length; i < l; i++) {
-      if (typeof arg[i] !== 'number' || !isInteger(arg[i])) {
+      if (!isNumber(arg[i]) || !isInteger(arg[i])) {
         throw new TypeError('Index parameters must be positive integer numbers');
       }
     }
     // create matrix
-    return new ImmutableDenseMatrix(arg);
+    var matrix = new ImmutableDenseMatrix();
+    matrix._data = arg;
+    matrix._size = [arg.length];
+    return matrix;
   }
 
   /**
@@ -130,7 +128,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     var size = [];
     for (var i = 0, ii = this._dimensions.length; i < ii; i++) {
       var d = this._dimensions[i];
-      size[i] = typeof d === 'string' ? 1 : d.size()[0];
+      size[i] = isString(d) || isNumber(d) ? 1 : d.size()[0];
     }
     return size;
   };
@@ -144,7 +142,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     var values = [];
     for (var i = 0, ii = this._dimensions.length; i < ii; i++) {
       var range = this._dimensions[i];
-      values[i] = typeof range === 'string' ? range : range.max();
+      values[i] = isString(range) || isNumber(range) ? range : range.max();
     }
     return values;
   };
@@ -158,7 +156,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     var values = [];
     for (var i = 0, ii = this._dimensions.length; i < ii; i++) {
       var range = this._dimensions[i];
-      values[i] = typeof range === 'string' ? range : range.min();
+      values[i] = isString(range) || isNumber(range) ? range : range.min();
     }
     return values;
   };
@@ -183,10 +181,11 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
    * @returns {Range | null} range
    */
   Index.prototype.dimension = function (dim) {
-    if (typeof dim !== 'number') {
+    var _this$_dimensions$dim;
+    if (!isNumber(dim)) {
       return null;
     }
-    return this._dimensions[dim] || null;
+    return (_this$_dimensions$dim = this._dimensions[dim]) !== null && _this$_dimensions$dim !== void 0 ? _this$_dimensions$dim : null;
   };
 
   /**
@@ -194,7 +193,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
    * @returns {boolean} Returns true if the index is an object property
    */
   Index.prototype.isObjectProperty = function () {
-    return this._dimensions.length === 1 && typeof this._dimensions[0] === 'string';
+    return this._dimensions.length === 1 && isString(this._dimensions[0]);
   };
 
   /**
@@ -228,7 +227,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     var array = [];
     for (var i = 0, ii = this._dimensions.length; i < ii; i++) {
       var dimension = this._dimensions[i];
-      array.push(typeof dimension === 'string' ? dimension : dimension.toArray());
+      array.push(isString(dimension) || isNumber(dimension) ? dimension : dimension.toArray());
     }
     return array;
   };
@@ -250,7 +249,7 @@ export var createIndexClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     var strings = [];
     for (var i = 0, ii = this._dimensions.length; i < ii; i++) {
       var dimension = this._dimensions[i];
-      if (typeof dimension === 'string') {
+      if (isString(dimension)) {
         strings.push(JSON.stringify(dimension));
       } else {
         strings.push(dimension.toString());

package/lib/esm/type/matrix/SparseMatrix.js

@@ -15,9 +15,27 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     Matrix
   } = _ref;
   /**
-   * Sparse Matrix implementation. This type implements
-   * a [Compressed Column Storage](https://en.wikipedia.org/wiki/Sparse_matrix#Compressed_sparse_column_(CSC_or_CCS))
-   * format for two-dimensional sparse matrices.
+   * Sparse Matrix implementation. This type (currently) implements 2D
+   * matrices only via the format known as
+   * [Compressed Column Storage](https://en.wikipedia.org/wiki/Sparse_matrix#Compressed_sparse_column_(CSC_or_CCS)).
+   *
+   * The structure/invariants of the internal data should be:
+   * 1. _values is an array of the nonzero values in order from top to bottom
+   *    (of each column), left to right.
+   * 2. _index is an array of row numbers, of the same length as and
+   *    corresponding positionally to _values.
+   * 3. _ptr is an array of length one more than the number of columns. For j
+   *    less than the number of columns, the "half-open" span of indices
+   *    _ptr[j] to _ptr[j+1] (i.e. including _ptr[j] if it is less than
+   *    _ptr[j+1], but never including _ptr[j+1]) are the indices in _values
+   *    of the nonzero elements in column j. Note there are no nonzero elements
+   *    in column j exactly when _ptr[j] === _ptr[j+1], and that the final
+   *    entry in _ptr is always exactly the number of nonzero entries in the
+   *    matrix.
+   * 4. _size is a length-2 array consisting of the number of rows followed by
+   *    the number of columns.
+   * 5. _datatype, if set, is the mathjs typeOf value of all entries of the
+   *    SparseMatrix.
    * @class SparseMatrix
    */
   function SparseMatrix(data, datatype) {
@@ -135,7 +153,9 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
       } while (j < columns);
     }
     // store number of values in ptr
-    matrix._ptr.push(matrix._index.length);
+    while (matrix._ptr.length <= columns) {
+      matrix._ptr.push(matrix._index.length);
+    }
     // size
     matrix._size = [rows, columns];
   }
@@ -298,12 +318,13 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     var pv = [];
 
     // loop rows in resulting matrix
-    rows.forEach(function (i, r) {
+    function rowsCallback(i, r) {
       // update permutation vector
       pv[i] = r[0];
       // mark i in workspace
       w[i] = true;
-    });
+    }
+    if (Number.isInteger(rows)) rowsCallback(rows, [0]);else rows.forEach(rowsCallback);
 
     // result matrix arrays
     var values = mvalues ? [] : undefined;
@@ -311,7 +332,7 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
     var ptr = [];
 
     // loop columns in result matrix
-    columns.forEach(function (j) {
+    function columnsCallback(j) {
       // update ptr
       ptr.push(index.length);
       // loop values in column j
@@ -328,7 +349,8 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
           }
         }
       }
-    });
+    }
+    if (Number.isInteger(columns)) columnsCallback(columns);else columns.forEach(columnsCallback);
     // update ptr
     ptr.push(index.length);
 
@@ -402,7 +424,7 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
       if (iSize.length === 1) {
         // if the replacement index only has 1 dimension, go trough each one and set its value
         var range = index.dimension(0);
-        range.forEach(function (dataIndex, subIndex) {
+        _forEachIndex(range, (dataIndex, subIndex) => {
           validateIndex(dataIndex);
           matrix.set([dataIndex, 0], submatrix[subIndex[0]], defaultValue);
         });
@@ -410,9 +432,9 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
         // if the replacement index has 2 dimensions, go through each one and set the value in the correct index
         var firstDimensionRange = index.dimension(0);
         var secondDimensionRange = index.dimension(1);
-        firstDimensionRange.forEach(function (firstDataIndex, firstSubIndex) {
+        _forEachIndex(firstDimensionRange, (firstDataIndex, firstSubIndex) => {
           validateIndex(firstDataIndex);
-          secondDimensionRange.forEach(function (secondDataIndex, secondSubIndex) {
+          _forEachIndex(secondDimensionRange, (secondDataIndex, secondSubIndex) => {
             validateIndex(secondDataIndex);
             matrix.set([firstDataIndex, secondDataIndex], submatrix[firstSubIndex[0]][secondSubIndex[0]], defaultValue);
           });
@@ -420,6 +442,10 @@ export var createSparseMatrixClass = /* #__PURE__ */factory(name, dependencies,
       }
     }
     return matrix;
+    function _forEachIndex(index, callback) {
+      // iterate cases where index is a Matrix or a Number
+      if (isNumber(index)) callback(index, [0]);else index.forEach(callback);
+    }
   }
 
   /**

package/lib/esm/type/unit/Unit.js

@@ -9,7 +9,7 @@ import { endsWith } from '../../utils/string.js';
 import { clone, hasOwnProperty } from '../../utils/object.js';
 import { createBigNumberPi as createPi } from '../../utils/bignumber/constants.js';
 var name = 'Unit';
-var dependencies = ['?on', 'config', 'addScalar', 'subtractScalar', 'multiplyScalar', 'divideScalar', 'pow', 'abs', 'fix', 'round', 'equal', 'isNumeric', 'format', 'toBest', 'number', 'Complex', 'BigNumber', 'Fraction'];
+var dependencies = ['?on', 'config', 'addScalar', 'subtractScalar', 'multiplyScalar', 'divideScalar', 'pow', 'abs', 'fix', 'round', 'equal', 'isNumeric', 'format', 'number', 'Complex', 'BigNumber', 'Fraction'];
 export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     on,
@@ -25,13 +25,15 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     equal,
     isNumeric,
     format,
-    toBest,
     number: _number,
     Complex,
     BigNumber: _BigNumber,
     Fraction: _Fraction
   } = _ref;
   var toNumber = _number;
+  var fixPrefixDefault = false;
+  var skipAutomaticSimplificationDefault = true;
+
   /**
    * A unit can be constructed in the following ways:
    *
@@ -58,13 +60,13 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
     if (!(value === null || value === undefined || isNumeric(value) || isComplex(value))) {
       throw new TypeError('First parameter in Unit constructor must be number, BigNumber, Fraction, Complex, or undefined');
     }
-    this.fixPrefix = false; // if true, function format will not search for the
+    this.fixPrefix = fixPrefixDefault; // if true, function format will not search for the
     // best prefix but leave it as initially provided.
     // fixPrefix is set true by the method Unit.to
 
     // The justification behind this is that if the constructor is explicitly called,
     // the caller wishes the units to be returned exactly as supplied.
-    this.skipAutomaticSimplification = true;
+    this.skipAutomaticSimplification = skipAutomaticSimplificationDefault;
     if (valuelessUnit === undefined) {
       this.units = [];
       this.dimensions = BASE_DIMENSIONS.map(x => 0);
@@ -855,14 +857,15 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
    * Get a JSON representation of the unit
    * @memberof Unit
    * @returns {Object} Returns a JSON object structured as:
-   *                   `{"mathjs": "Unit", "value": 2, "unit": "cm", "fixPrefix": false}`
+   *                   `{"mathjs": "Unit", "value": 2, "unit": "cm", "fixPrefix": false, "skipSimp": true}`
    */
   Unit.prototype.toJSON = function () {
     return {
       mathjs: 'Unit',
       value: this._denormalize(this.value),
       unit: this.units.length > 0 ? this.formatUnits() : null,
-      fixPrefix: this.fixPrefix
+      fixPrefix: this.fixPrefix,
+      skipSimp: this.skipAutomaticSimplification
     };
   };
 
@@ -874,9 +877,10 @@ export var createUnitClass = /* #__PURE__ */factory(name, dependencies, _ref =>
    * @return {Unit}
    */
   Unit.fromJSON = function (json) {
-    var _json$unit;
+    var _json$unit, _json$fixPrefix, _json$skipSimp;
     var unit = new Unit(json.value, (_json$unit = json.unit) !== null && _json$unit !== void 0 ? _json$unit : undefined);
-    unit.fixPrefix = json.fixPrefix || false;
+    unit.fixPrefix = (_json$fixPrefix = json.fixPrefix) !== null && _json$fixPrefix !== void 0 ? _json$fixPrefix : fixPrefixDefault;
+    unit.skipAutomaticSimplification = (_json$skipSimp = json.skipSimp) !== null && _json$skipSimp !== void 0 ? _json$skipSimp : skipAutomaticSimplificationDefault;
     return unit;
   };
 

package/lib/esm/utils/number.js

@@ -13,7 +13,7 @@ export function isInteger(value) {
   if (typeof value === 'boolean') {
     return true;
   }
-  return isFinite(value) ? value === Math.round(value) : false;
+  return Number.isFinite(value) ? value === Math.round(value) : false;
 }
 
 /**
@@ -105,7 +105,7 @@ export var cbrt = Math.cbrt || function cbrt(x) {
   if (negate) {
     x = -x;
   }
-  if (isFinite(x)) {
+  if (Number.isFinite(x)) {
     result = Math.exp(Math.log(x) / 3);
     // from https://en.wikipedia.org/wiki/Cube_root#Numerical_methods
     result = (x / (result * result) + 2 * result) / 3;
@@ -381,7 +381,7 @@ export function splitNumber(value) {
  * @param {number} [precision]        Optional number of significant figures to return.
  */
 export function toEngineering(value, precision) {
-  if (isNaN(value) || !isFinite(value)) {
+  if (isNaN(value) || !Number.isFinite(value)) {
     return String(value);
   }
   var split = splitNumber(value);
@@ -430,7 +430,7 @@ export function toEngineering(value, precision) {
  *                                        decimal point. null by default.
  */
 export function toFixed(value, precision) {
-  if (isNaN(value) || !isFinite(value)) {
+  if (isNaN(value) || !Number.isFinite(value)) {
     return String(value);
   }
   var splitValue = splitNumber(value);
@@ -465,7 +465,7 @@ export function toFixed(value, precision) {
  *                              is used.
  */
 export function toExponential(value, precision) {
-  if (isNaN(value) || !isFinite(value)) {
+  if (isNaN(value) || !Number.isFinite(Number(value))) {
     return String(value);
   }
 
@@ -496,7 +496,7 @@ export function toExponential(value, precision) {
  * @return {string}
  */
 export function toPrecision(value, precision, options) {
-  if (isNaN(value) || !isFinite(value)) {
+  if (isNaN(value) || !Number.isFinite(value)) {
     return String(value);
   }
 
@@ -635,7 +635,7 @@ export function nearlyEqual(a, b) {
   if (isNaN(a) || isNaN(b)) {
     return false;
   }
-  if (!isFinite(a) || !isFinite(b)) {
+  if (!Number.isFinite(a) || !Number.isFinite(b)) {
     return a === b;
   }
   if (a === b) {

package/lib/esm/utils/optimizeCallback.js

@@ -17,7 +17,19 @@ export function optimizeCallback(callback, array, name, isUnary) {
     if (isUnary) {
       numberOfArguments = 1;
     } else {
-      var firstIndex = (array.isMatrix ? array.size() : arraySize(array)).map(() => 0);
+      var size = array.isMatrix ? array.size() : arraySize(array);
+
+      // Check the size of the last dimension to see if the array/matrix is empty
+      var isEmpty = size.length ? size[size.length - 1] === 0 : true;
+      if (isEmpty) {
+        // don't optimize callbacks for empty arrays/matrix, as they will never be called
+        // and in fact will throw an exception when we try to access the first element below
+        return {
+          isUnary,
+          fn: callback
+        };
+      }
+      var firstIndex = size.map(() => 0);
       var firstValue = array.isMatrix ? array.get(firstIndex) : get(array, firstIndex);
       numberOfArguments = _findNumberOfArgumentsTyped(callback, firstValue, firstIndex, array);
     }

package/lib/esm/version.js

@@ -1,3 +1,3 @@
-export var version = '14.9.1';
+export var version = '15.1.0';
 // Note: This file is automatically generated when building math.js.
 // Changes made in this file will be overwritten.

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "mathjs",
-  "version": "14.9.1",
+  "version": "15.1.0",
   "description": "Math.js is an extensive math library for JavaScript and Node.js. It features a flexible expression parser with support for symbolic computation, comes with a large set of built-in functions and constants, and offers an integrated solution to work with different data types like numbers, big numbers, complex numbers, fractions, units, and matrices.",
   "author": "Jos de Jong <wjosdejong@gmail.com> (https://github.com/josdejong)",
   "homepage": "https://mathjs.org",
@@ -36,16 +36,16 @@
     "typed-function": "^4.2.1"
   },
   "devDependencies": {
-    "@babel/core": "7.28.4",
+    "@babel/core": "7.28.5",
     "@babel/plugin-transform-object-assign": "7.27.1",
     "@babel/plugin-transform-optional-catch-binding": "7.27.1",
-    "@babel/plugin-transform-runtime": "7.28.3",
-    "@babel/preset-env": "7.28.3",
+    "@babel/plugin-transform-runtime": "7.28.5",
+    "@babel/preset-env": "7.28.5",
     "@babel/register": "7.28.3",
     "@types/assert": "1.5.11",
     "@types/mocha": "10.0.10",
-    "@typescript-eslint/eslint-plugin": "8.46.0",
-    "@typescript-eslint/parser": "8.46.0",
+    "@typescript-eslint/eslint-plugin": "8.46.3",
+    "@typescript-eslint/parser": "8.46.3",
     "assert": "2.1.0",
     "babel-loader": "10.0.0",
     "c8": "10.1.3",
@@ -77,7 +77,7 @@
     "karma-webdriver-launcher": "1.0.8",
     "karma-webpack": "5.0.1",
     "mkdirp": "3.0.1",
-    "mocha": "11.7.4",
+    "mocha": "11.7.5",
     "mocha-junit-reporter": "2.2.1",
     "ndarray": "1.0.19",
     "ndarray-determinant": "1.0.0",
@@ -89,7 +89,7 @@
     "process": "0.11.10",
     "sinon": "21.0.0",
     "sylvester": "0.0.21",
-    "tinybench": "5.0.1",
+    "tinybench": "5.1.0",
     "ts-node": "10.9.2",
     "typescript": "5.9.3",
     "webpack": "5.102.1",