mathjs 14.2.1 → 14.3.0

package/lib/esm/expression/embeddedDocs/function/statistics/max.js

@@ -2,7 +2,7 @@ export var maxDocs = {
   name: 'max',
   category: 'Statistics',
   syntax: ['max(a, b, c, ...)', 'max(A)', 'max(A, dimension)'],
-  description: 'Compute the maximum value of a list of values.',
+  description: 'Compute the maximum value of a list of values. If any NaN values are found, the function yields the last NaN in the input.',
   examples: ['max(2, 3, 4, 1)', 'max([2, 3, 4, 1])', 'max([2, 5; 4, 3])', 'max([2, 5; 4, 3], 1)', 'max([2, 5; 4, 3], 2)', 'max(2.7, 7.1, -4.5, 2.0, 4.1)', 'min(2.7, 7.1, -4.5, 2.0, 4.1)'],
   seealso: ['mean', 'median', 'min', 'prod', 'std', 'sum', 'variance']
 };

package/lib/esm/expression/embeddedDocs/function/statistics/min.js

@@ -2,7 +2,7 @@ export var minDocs = {
   name: 'min',
   category: 'Statistics',
   syntax: ['min(a, b, c, ...)', 'min(A)', 'min(A, dimension)'],
-  description: 'Compute the minimum value of a list of values.',
+  description: 'Compute the minimum value of a list of values. If any NaN values are found, the function yields the last NaN in the input.',
   examples: ['min(2, 3, 4, 1)', 'min([2, 3, 4, 1])', 'min([2, 5; 4, 3])', 'min([2, 5; 4, 3], 1)', 'min([2, 5; 4, 3], 2)', 'min(2.7, 7.1, -4.5, 2.0, 4.1)', 'max(2.7, 7.1, -4.5, 2.0, 4.1)'],
   seealso: ['max', 'mean', 'median', 'prod', 'std', 'sum', 'variance']
 };

package/lib/esm/expression/transform/max.transform.js

@@ -3,19 +3,21 @@ import { errorTransform } from './utils/errorTransform.js';
 import { createMax } from '../../function/statistics/max.js';
 import { lastDimToZeroBase } from './utils/lastDimToZeroBase.js';
 var name = 'max';
-var dependencies = ['typed', 'config', 'numeric', 'larger'];
+var dependencies = ['typed', 'config', 'numeric', 'larger', 'isNaN'];
 export var createMaxTransform = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   } = _ref;
   var max = createMax({
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   });
 
   /**

package/lib/esm/expression/transform/min.transform.js

@@ -3,19 +3,21 @@ import { errorTransform } from './utils/errorTransform.js';
 import { createMin } from '../../function/statistics/min.js';
 import { lastDimToZeroBase } from './utils/lastDimToZeroBase.js';
 var name = 'min';
-var dependencies = ['typed', 'config', 'numeric', 'smaller'];
+var dependencies = ['typed', 'config', 'numeric', 'smaller', 'isNaN'];
 export var createMinTransform = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   } = _ref;
   var min = createMin({
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   });
 
   /**

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

@@ -20,7 +20,7 @@ export var createVarianceTransform = /* #__PURE__ */factory(name, dependencies,
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   } = _ref;
   var variance = createVariance({
     typed,
@@ -29,7 +29,7 @@ export var createVarianceTransform = /* #__PURE__ */factory(name, dependencies,
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   });
   return typed(name, {
     '...any': function any(args) {

package/lib/esm/function/matrix/flatten.js

@@ -32,7 +32,7 @@ export var createFlatten = /* #__PURE__ */factory(name, dependencies, _ref => {
     Matrix: function Matrix(x) {
       // Return the same matrix type as x (Dense or Sparse Matrix)
       // Return the same data type as x
-      return x.create(flattenArray(x.toArray()), x.datatype());
+      return x.create(flattenArray(x.valueOf()), x.datatype());
     }
   });
 });

package/lib/esm/function/matrix/forEach.js

@@ -1,6 +1,6 @@
 import { optimizeCallback } from '../../utils/optimizeCallback.js';
 import { factory } from '../../utils/factory.js';
-import { recurse } from '../../utils/array.js';
+import { deepForEach } from '../../utils/array.js';
 var name = 'forEach';
 var dependencies = ['typed'];
 export var createForEach = /* #__PURE__ */factory(name, dependencies, _ref => {
@@ -52,5 +52,5 @@ export var createForEach = /* #__PURE__ */factory(name, dependencies, _ref => {
  * @private
  */
 function _forEach(array, callback) {
-  recurse(array, [], array, optimizeCallback(callback, array, name));
+  deepForEach(array, optimizeCallback(callback, array, name));
 }

package/lib/esm/function/matrix/map.js

@@ -1,5 +1,5 @@
 import { optimizeCallback } from '../../utils/optimizeCallback.js';
-import { arraySize, broadcastSizes, broadcastTo, get, recurse } from '../../utils/array.js';
+import { arraySize, broadcastSizes, broadcastTo, get, deepMap } from '../../utils/array.js';
 import { factory } from '../../utils/factory.js';
 var name = 'map';
 var dependencies = ['typed'];
@@ -137,6 +137,6 @@ export var createMap = /* #__PURE__ */factory(name, dependencies, _ref => {
   * @private
   */
   function _mapArray(array, callback) {
-    return recurse(array, [], array, optimizeCallback(callback, array, name));
+    return deepMap(array, optimizeCallback(callback, array, name));
   }
 });

package/lib/esm/function/matrix/partitionSelect.js

@@ -7,7 +7,7 @@ export var createPartitionSelect = /* #__PURE__ */factory(name, dependencies, _r
   var {
     typed,
     isNumeric,
-    isNaN,
+    isNaN: mathIsNaN,
     compare
   } = _ref;
   var asc = compare;
@@ -98,7 +98,7 @@ export var createPartitionSelect = /* #__PURE__ */factory(name, dependencies, _r
 
     // check for NaN values since these can cause an infinite while loop
     for (var i = 0; i < arr.length; i++) {
-      if (isNumeric(arr[i]) && isNaN(arr[i])) {
+      if (isNumeric(arr[i]) && mathIsNaN(arr[i])) {
         return arr[i]; // return NaN
       }
     }

package/lib/esm/function/statistics/max.js

@@ -3,13 +3,14 @@ import { factory } from '../../utils/factory.js';
 import { safeNumberType } from '../../utils/number.js';
 import { improveErrorMessage } from './utils/improveErrorMessage.js';
 var name = 'max';
-var dependencies = ['typed', 'config', 'numeric', 'larger'];
+var dependencies = ['typed', 'config', 'numeric', 'larger', 'isNaN'];
 export var createMax = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the maximum value of a matrix or a  list with values.
@@ -83,8 +84,8 @@ export var createMax = /* #__PURE__ */factory(name, dependencies, _ref => {
     var res;
     deepForEach(array, function (value) {
       try {
-        if (typeof value === 'number' && isNaN(value)) {
-          res = NaN;
+        if (mathIsNaN(value)) {
+          res = value;
         } else if (res === undefined || larger(value, res)) {
           res = value;
         }

package/lib/esm/function/statistics/min.js

@@ -3,13 +3,14 @@ import { factory } from '../../utils/factory.js';
 import { safeNumberType } from '../../utils/number.js';
 import { improveErrorMessage } from './utils/improveErrorMessage.js';
 var name = 'min';
-var dependencies = ['typed', 'config', 'numeric', 'smaller'];
+var dependencies = ['typed', 'config', 'numeric', 'smaller', 'isNaN'];
 export var createMin = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the minimum value of a matrix or a  list of values.
@@ -83,8 +84,8 @@ export var createMin = /* #__PURE__ */factory(name, dependencies, _ref => {
     var min;
     deepForEach(array, function (value) {
       try {
-        if (typeof value === 'number' && isNaN(value)) {
-          min = NaN;
+        if (mathIsNaN(value)) {
+          min = value;
         } else if (min === undefined || smaller(value, min)) {
           min = value;
         }

package/lib/esm/function/statistics/mode.js

@@ -5,7 +5,7 @@ var dependencies = ['typed', 'isNaN', 'isNumeric'];
 export var createMode = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
-    isNaN,
+    isNaN: mathIsNaN,
     isNumeric
   } = _ref;
   /**
@@ -57,7 +57,7 @@ export var createMode = /* #__PURE__ */factory(name, dependencies, _ref => {
     var max = 0;
     for (var i = 0; i < values.length; i++) {
       var value = values[i];
-      if (isNumeric(value) && isNaN(value)) {
+      if (isNumeric(value) && mathIsNaN(value)) {
         throw new Error('Cannot calculate mode of an array containing NaN values');
       }
       if (!(value in count)) {

package/lib/esm/function/statistics/variance.js

@@ -13,7 +13,7 @@ export var createVariance = /* #__PURE__ */factory(name, dependencies, _ref => {
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the variance of a matrix or a  list with values.
@@ -123,7 +123,7 @@ export var createVariance = /* #__PURE__ */factory(name, dependencies, _ref => {
       var diff = subtract(value, mean);
       sum = sum === undefined ? multiply(diff, diff) : add(sum, multiply(diff, diff));
     });
-    if (isNaN(sum)) {
+    if (mathIsNaN(sum)) {
       return sum;
     }
     switch (normalization) {

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

@@ -525,56 +525,35 @@ export var createDenseMatrixClass = /* #__PURE__ */factory(name, dependencies, _
    *                              array, and the Matrix being traversed.
    */
   DenseMatrix.prototype._forEach = function (callback) {
-    // matrix instance
     var me = this;
     var s = me.size();
-
-    // if there is only one dimension, just loop through it
-    if (s.length === 1) {
-      for (var i = 0; i < s[0]; i++) {
+    var maxDepth = s.length - 1;
+    if (maxDepth < 0) {
+      return;
+    }
+    if (maxDepth === 0) {
+      var thisSize = s[0];
+      for (var i = 0; i < thisSize; i++) {
         callback(me._data, i, [i]);
       }
       return;
     }
-
-    // keep track of the current index permutation
-    var index = Array(s.length).fill(0);
-
-    // store a reference of each dimension of the matrix for faster access
-    var data = Array(s.length - 1);
-    var last = data.length - 1;
-    data[0] = me._data[0];
-    for (var _i = 0; _i < last; _i++) {
-      data[_i + 1] = data[_i][0];
-    }
-    index[last] = -1;
-    while (true) {
-      var _i2 = void 0;
-      for (_i2 = last; _i2 >= 0; _i2--) {
-        // march index to the next permutation
-        index[_i2]++;
-        if (index[_i2] === s[_i2]) {
-          index[_i2] = 0;
-          continue;
-        }
-
-        // update references to matrix dimensions
-        data[_i2] = _i2 === 0 ? me._data[index[_i2]] : data[_i2 - 1][index[_i2]];
-        for (var j = _i2; j < last; j++) {
-          data[j + 1] = data[j][0];
+    var index = Array(s.length);
+    function recurse(data, depth) {
+      var thisSize = s[depth];
+      if (depth < maxDepth) {
+        for (var _i = 0; _i < thisSize; _i++) {
+          index[depth] = _i;
+          recurse(data[_i], depth + 1);
         }
-
-        // loop through the last dimension and map each value
-        for (var _j = 0; _j < s[data.length]; _j++) {
-          index[data.length] = _j;
-          callback(data[last], _j, index.slice(0));
+      } else {
+        for (var _i2 = 0; _i2 < thisSize; _i2++) {
+          index[depth] = _i2;
+          callback(data, _i2, index.slice());
         }
-        break;
-      }
-      if (_i2 === -1) {
-        break;
       }
     }
+    recurse(me._data, 0);
   };
 
   /**
@@ -617,19 +596,37 @@ export var createDenseMatrixClass = /* #__PURE__ */factory(name, dependencies, _
    * @return {Iterable<{ value, index: number[] }>}
    */
   DenseMatrix.prototype[Symbol.iterator] = function* () {
-    var _recurse = function* recurse(value, index) {
-      if (isArray(value)) {
-        for (var i = 0; i < value.length; i++) {
-          yield* _recurse(value[i], index.concat(i));
-        }
-      } else {
+    var maxDepth = this._size.length - 1;
+    if (maxDepth < 0) {
+      return;
+    }
+    if (maxDepth === 0) {
+      for (var i = 0; i < this._data.length; i++) {
         yield {
-          value,
-          index
+          value: this._data[i],
+          index: [i]
         };
       }
+      return;
+    }
+    var index = [];
+    var _recurse = function* recurse(value, depth) {
+      if (depth < maxDepth) {
+        for (var _i3 = 0; _i3 < value.length; _i3++) {
+          index[depth] = _i3;
+          yield* _recurse(value[_i3], depth + 1);
+        }
+      } else {
+        for (var _i4 = 0; _i4 < value.length; _i4++) {
+          index[depth] = _i4;
+          yield {
+            value: value[_i4],
+            index: index.slice()
+          };
+        }
+      }
     };
-    yield* _recurse(this._data, []);
+    yield* _recurse(this._data, 0);
   };
 
   /**