mathjs 14.2.1 → 14.3.0

package/lib/cjs/entry/dependenciesAny/dependenciesMax.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.maxDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesLargerGenerated = require("./dependenciesLarger.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesAny = require("../../factoriesAny.js");
  */
 
 const maxDependencies = exports.maxDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   largerDependencies: _dependenciesLargerGenerated.largerDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesAny/dependenciesMaxTransform.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.maxTransformDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesLargerGenerated = require("./dependenciesLarger.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesAny = require("../../factoriesAny.js");
  */
 
 const maxTransformDependencies = exports.maxTransformDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   largerDependencies: _dependenciesLargerGenerated.largerDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesAny/dependenciesMin.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.minDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesSmallerGenerated = require("./dependenciesSmaller.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesAny = require("../../factoriesAny.js");
  */
 
 const minDependencies = exports.minDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   smallerDependencies: _dependenciesSmallerGenerated.smallerDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesAny/dependenciesMinTransform.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.minTransformDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesSmallerGenerated = require("./dependenciesSmaller.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesAny = require("../../factoriesAny.js");
  */
 
 const minTransformDependencies = exports.minTransformDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   smallerDependencies: _dependenciesSmallerGenerated.smallerDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesNumber/dependenciesMax.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.maxDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesLargerGenerated = require("./dependenciesLarger.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesNumber = require("../../factoriesNumber.js");
  */
 
 const maxDependencies = exports.maxDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   largerDependencies: _dependenciesLargerGenerated.largerDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesNumber/dependenciesMaxTransform.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.maxTransformDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesLargerGenerated = require("./dependenciesLarger.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesNumber = require("../../factoriesNumber.js");
  */
 
 const maxTransformDependencies = exports.maxTransformDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   largerDependencies: _dependenciesLargerGenerated.largerDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesNumber/dependenciesMin.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.minDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesSmallerGenerated = require("./dependenciesSmaller.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesNumber = require("../../factoriesNumber.js");
  */
 
 const minDependencies = exports.minDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   smallerDependencies: _dependenciesSmallerGenerated.smallerDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

package/lib/cjs/entry/dependenciesNumber/dependenciesMinTransform.generated.js

@@ -4,6 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.minTransformDependencies = void 0;
+var _dependenciesIsNaNGenerated = require("./dependenciesIsNaN.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesNumeric.generated.js");
 var _dependenciesSmallerGenerated = require("./dependenciesSmaller.generated.js");
 var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
@@ -14,6 +15,7 @@ var _factoriesNumber = require("../../factoriesNumber.js");
  */
 
 const minTransformDependencies = exports.minTransformDependencies = {
+  isNaNDependencies: _dependenciesIsNaNGenerated.isNaNDependencies,
   numericDependencies: _dependenciesNumericGenerated.numericDependencies,
   smallerDependencies: _dependenciesSmallerGenerated.smallerDependencies,
   typedDependencies: _dependenciesTypedGenerated.typedDependencies,

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

@@ -624,6 +624,7 @@ const help = exports.help = (0, _factoriesAny.createHelp)({
   }),
   min: (0, _factoriesAny.createMinTransform)({
     config: _configReadonly.config,
+    isNaN: _pureFunctionsAnyGenerated.isNaN,
     numeric: _pureFunctionsAnyGenerated.numeric,
     smaller: _pureFunctionsAnyGenerated.smaller,
     typed: _pureFunctionsAnyGenerated.typed
@@ -655,6 +656,7 @@ const help = exports.help = (0, _factoriesAny.createHelp)({
   }),
   max: (0, _factoriesAny.createMaxTransform)({
     config: _configReadonly.config,
+    isNaN: _pureFunctionsAnyGenerated.isNaN,
     larger: _pureFunctionsAnyGenerated.larger,
     numeric: _pureFunctionsAnyGenerated.numeric,
     typed: _pureFunctionsAnyGenerated.typed

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

@@ -462,12 +462,14 @@ const parser = exports.parser = (0, _factoriesNumber.createParser)({
   }),
   max: (0, _factoriesNumber.createMaxTransform)({
     config: _configReadonly.config,
+    isNaN: _pureFunctionsNumberGenerated.isNaN,
     larger: _pureFunctionsNumberGenerated.larger,
     numeric: _pureFunctionsNumberGenerated.numeric,
     typed: _pureFunctionsNumberGenerated.typed
   }),
   min: (0, _factoriesNumber.createMinTransform)({
     config: _configReadonly.config,
+    isNaN: _pureFunctionsNumberGenerated.isNaN,
     numeric: _pureFunctionsNumberGenerated.numeric,
     smaller: _pureFunctionsNumberGenerated.smaller,
     typed: _pureFunctionsNumberGenerated.typed

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

@@ -748,12 +748,14 @@ const matrixFromColumns = exports.matrixFromColumns = /* #__PURE__ */(0, _factor
 });
 const max = exports.max = /* #__PURE__ */(0, _factoriesAny.createMax)({
   config: _configReadonly.config,
+  isNaN,
   larger,
   numeric,
   typed
 });
 const min = exports.min = /* #__PURE__ */(0, _factoriesAny.createMin)({
   config: _configReadonly.config,
+  isNaN,
   numeric,
   smaller,
   typed

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

@@ -602,12 +602,14 @@ const corr = exports.corr = /* #__PURE__ */(0, _factoriesNumber.createCorr)({
 });
 const max = exports.max = /* #__PURE__ */(0, _factoriesNumber.createMax)({
   config: _configReadonly.config,
+  isNaN,
   larger,
   numeric,
   typed
 });
 const min = exports.min = /* #__PURE__ */(0, _factoriesNumber.createMin)({
   config: _configReadonly.config,
+  isNaN,
   numeric,
   smaller,
   typed

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

@@ -11,7 +11,7 @@ const mapSlicesDocs = exports.mapSlicesDocs = {
   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]
+  'mapSlices(A, 2, prod)' // returns [2, 12]
   ],
   seealso: ['map', 'forEach']
 };

package/lib/cjs/expression/embeddedDocs/function/relational/larger.js

@@ -8,7 +8,7 @@ const largerDocs = exports.largerDocs = {
   name: 'larger',
   category: 'Relational',
   syntax: ['x > y', 'larger(x, y)'],
-  description: 'Check if value x is larger than y. Returns true if x is larger than y, and false if not.',
+  description: 'Check if value x is larger than y. Returns true if x is larger than y, and false if not. Comparing a value with NaN returns false.',
   examples: ['2 > 3', '5 > 2*2', 'a = 3.3', 'b = 6-2.8', '(a > b)', '(b < a)', '5 cm > 2 inch'],
   seealso: ['equal', 'unequal', 'smaller', 'smallerEq', 'largerEq', 'compare']
 };

package/lib/cjs/expression/embeddedDocs/function/relational/smaller.js

@@ -8,7 +8,7 @@ const smallerDocs = exports.smallerDocs = {
   name: 'smaller',
   category: 'Relational',
   syntax: ['x < y', 'smaller(x, y)'],
-  description: 'Check if value x is smaller than value y. Returns true if x is smaller than y, and false if not.',
+  description: 'Check if value x is smaller than value y. Returns true if x is smaller than y, and false if not. Comparing a value with NaN returns false.',
   examples: ['2 < 3', '5 < 2*2', 'a = 3.3', 'b = 6-2.8', '(a < b)', '5 cm < 2 inch'],
   seealso: ['equal', 'unequal', 'larger', 'smallerEq', 'largerEq', 'compare']
 };

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

@@ -8,7 +8,7 @@ const maxDocs = exports.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/cjs/expression/embeddedDocs/function/statistics/min.js

@@ -8,7 +8,7 @@ const minDocs = exports.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/cjs/expression/transform/max.transform.js

@@ -9,19 +9,21 @@ var _errorTransform = require("./utils/errorTransform.js");
 var _max = require("../../function/statistics/max.js");
 var _lastDimToZeroBase = require("./utils/lastDimToZeroBase.js");
 const name = 'max';
-const dependencies = ['typed', 'config', 'numeric', 'larger'];
+const dependencies = ['typed', 'config', 'numeric', 'larger', 'isNaN'];
 const createMaxTransform = exports.createMaxTransform = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   } = _ref;
   const max = (0, _max.createMax)({
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   });
 
   /**

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

@@ -9,19 +9,21 @@ var _errorTransform = require("./utils/errorTransform.js");
 var _min = require("../../function/statistics/min.js");
 var _lastDimToZeroBase = require("./utils/lastDimToZeroBase.js");
 const name = 'min';
-const dependencies = ['typed', 'config', 'numeric', 'smaller'];
+const dependencies = ['typed', 'config', 'numeric', 'smaller', 'isNaN'];
 const createMinTransform = exports.createMinTransform = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   } = _ref;
   const min = (0, _min.createMin)({
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   });
 
   /**

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

@@ -26,7 +26,7 @@ const createVarianceTransform = exports.createVarianceTransform = /* #__PURE__ *
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   } = _ref;
   const variance = (0, _variance.createVariance)({
     typed,
@@ -35,7 +35,7 @@ const createVarianceTransform = exports.createVarianceTransform = /* #__PURE__ *
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   });
   return typed(name, {
     '...any': function (args) {

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

@@ -38,7 +38,7 @@ const createFlatten = exports.createFlatten = /* #__PURE__ */(0, _factory.factor
     Matrix: function (x) {
       // Return the same matrix type as x (Dense or Sparse Matrix)
       // Return the same data type as x
-      return x.create((0, _array.flatten)(x.toArray()), x.datatype());
+      return x.create((0, _array.flatten)(x.valueOf()), x.datatype());
     }
   });
 });

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

@@ -58,5 +58,5 @@ const createForEach = exports.createForEach = /* #__PURE__ */(0, _factory.factor
  * @private
  */
 function _forEach(array, callback) {
-  (0, _array.recurse)(array, [], array, (0, _optimizeCallback.optimizeCallback)(callback, array, name));
+  (0, _array.deepForEach)(array, (0, _optimizeCallback.optimizeCallback)(callback, array, name));
 }

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

@@ -143,6 +143,6 @@ const createMap = exports.createMap = /* #__PURE__ */(0, _factory.factory)(name,
   * @private
   */
   function _mapArray(array, callback) {
-    return (0, _array.recurse)(array, [], array, (0, _optimizeCallback.optimizeCallback)(callback, array, name));
+    return (0, _array.deepMap)(array, (0, _optimizeCallback.optimizeCallback)(callback, array, name));
   }
 });

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

@@ -13,7 +13,7 @@ const createPartitionSelect = exports.createPartitionSelect = /* #__PURE__ */(0,
   let {
     typed,
     isNumeric,
-    isNaN,
+    isNaN: mathIsNaN,
     compare
   } = _ref;
   const asc = compare;
@@ -104,7 +104,7 @@ const createPartitionSelect = exports.createPartitionSelect = /* #__PURE__ */(0,
 
     // check for NaN values since these can cause an infinite while loop
     for (let 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/cjs/function/statistics/max.js

@@ -9,13 +9,14 @@ var _factory = require("../../utils/factory.js");
 var _number = require("../../utils/number.js");
 var _improveErrorMessage = require("./utils/improveErrorMessage.js");
 const name = 'max';
-const dependencies = ['typed', 'config', 'numeric', 'larger'];
+const dependencies = ['typed', 'config', 'numeric', 'larger', 'isNaN'];
 const createMax = exports.createMax = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
     config,
     numeric,
-    larger
+    larger,
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the maximum value of a matrix or a  list with values.
@@ -89,8 +90,8 @@ const createMax = exports.createMax = /* #__PURE__ */(0, _factory.factory)(name,
     let res;
     (0, _collection.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/cjs/function/statistics/min.js

@@ -9,13 +9,14 @@ var _factory = require("../../utils/factory.js");
 var _number = require("../../utils/number.js");
 var _improveErrorMessage = require("./utils/improveErrorMessage.js");
 const name = 'min';
-const dependencies = ['typed', 'config', 'numeric', 'smaller'];
+const dependencies = ['typed', 'config', 'numeric', 'smaller', 'isNaN'];
 const createMin = exports.createMin = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
     config,
     numeric,
-    smaller
+    smaller,
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the minimum value of a matrix or a  list of values.
@@ -89,8 +90,8 @@ const createMin = exports.createMin = /* #__PURE__ */(0, _factory.factory)(name,
     let min;
     (0, _collection.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/cjs/function/statistics/mode.js

@@ -11,7 +11,7 @@ const dependencies = ['typed', 'isNaN', 'isNumeric'];
 const createMode = exports.createMode = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
   let {
     typed,
-    isNaN,
+    isNaN: mathIsNaN,
     isNumeric
   } = _ref;
   /**
@@ -63,7 +63,7 @@ const createMode = exports.createMode = /* #__PURE__ */(0, _factory.factory)(nam
     let max = 0;
     for (let i = 0; i < values.length; i++) {
       const 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/cjs/function/statistics/variance.js

@@ -19,7 +19,7 @@ const createVariance = exports.createVariance = /* #__PURE__ */(0, _factory.fact
     multiply,
     divide,
     mapSlices,
-    isNaN
+    isNaN: mathIsNaN
   } = _ref;
   /**
    * Compute the variance of a matrix or a  list with values.
@@ -129,7 +129,7 @@ const createVariance = exports.createVariance = /* #__PURE__ */(0, _factory.fact
       const 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/cjs/header.js

@@ -6,8 +6,8 @@
  * It features real and complex numbers, units, matrices, a large set of
  * mathematical functions, and a flexible expression parser.
  *
- * @version 14.2.1
- * @date    2025-02-05
+ * @version 14.3.0
+ * @date    2025-02-28
  *
  * @license
  * Copyright (C) 2013-2025 Jos de Jong <wjosdejong@gmail.com>

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

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