mathjs 14.4.0 → 14.5.0

package/lib/cjs/constants.js

@@ -76,47 +76,47 @@ const createLOG2E = exports.createLOG2E = /* #__PURE__ */recreateFactory('LOG2E'
   } = _ref9;
   return config.number === 'BigNumber' ? new BigNumber(1).div(new BigNumber(2).ln()) : Math.LOG2E;
 });
-const createLOG10E = exports.createLOG10E = /* #__PURE__ */recreateFactory('LOG10E', ['config', '?BigNumber'], _ref10 => {
+const createLOG10E = exports.createLOG10E = /* #__PURE__ */recreateFactory('LOG10E', ['config', '?BigNumber'], _ref0 => {
   let {
     config,
     BigNumber
-  } = _ref10;
+  } = _ref0;
   return config.number === 'BigNumber' ? new BigNumber(1).div(new BigNumber(10).ln()) : Math.LOG10E;
 });
 const createSQRT1_2 = exports.createSQRT1_2 = /* #__PURE__ */recreateFactory(
 // eslint-disable-line camelcase
-'SQRT1_2', ['config', '?BigNumber'], _ref11 => {
+'SQRT1_2', ['config', '?BigNumber'], _ref1 => {
   let {
     config,
     BigNumber
-  } = _ref11;
+  } = _ref1;
   return config.number === 'BigNumber' ? new BigNumber('0.5').sqrt() : Math.SQRT1_2;
 });
-const createSQRT2 = exports.createSQRT2 = /* #__PURE__ */recreateFactory('SQRT2', ['config', '?BigNumber'], _ref12 => {
+const createSQRT2 = exports.createSQRT2 = /* #__PURE__ */recreateFactory('SQRT2', ['config', '?BigNumber'], _ref10 => {
   let {
     config,
     BigNumber
-  } = _ref12;
+  } = _ref10;
   return config.number === 'BigNumber' ? new BigNumber(2).sqrt() : Math.SQRT2;
 });
-const createI = exports.createI = /* #__PURE__ */recreateFactory('i', ['Complex'], _ref13 => {
+const createI = exports.createI = /* #__PURE__ */recreateFactory('i', ['Complex'], _ref11 => {
   let {
     Complex
-  } = _ref13;
+  } = _ref11;
   return Complex.I;
 });
 
 // for backward compatibility with v5
-const createUppercasePi = exports.createUppercasePi = /* #__PURE__ */(0, _factory.factory)('PI', ['pi'], _ref14 => {
+const createUppercasePi = exports.createUppercasePi = /* #__PURE__ */(0, _factory.factory)('PI', ['pi'], _ref12 => {
   let {
     pi
-  } = _ref14;
+  } = _ref12;
   return pi;
 });
-const createUppercaseE = exports.createUppercaseE = /* #__PURE__ */(0, _factory.factory)('E', ['e'], _ref15 => {
+const createUppercaseE = exports.createUppercaseE = /* #__PURE__ */(0, _factory.factory)('E', ['e'], _ref13 => {
   let {
     e
-  } = _ref15;
+  } = _ref13;
   return e;
 });
 const createVersion = exports.createVersion = /* #__PURE__ */(0, _factory.factory)('version', [], () => _version.version);

package/lib/cjs/core/create.js

@@ -17,8 +17,7 @@ var emitter = _interopRequireWildcard(require("./../utils/emitter.js"));
 var _config = require("./config.js");
 var _config2 = require("./function/config.js");
 var _import = require("./function/import.js");
-function _getRequireWildcardCache(e) { if ("function" != typeof WeakMap) return null; var r = new WeakMap(), t = new WeakMap(); return (_getRequireWildcardCache = function (e) { return e ? t : r; })(e); }
-function _interopRequireWildcard(e, r) { if (!r && e && e.__esModule) return e; if (null === e || "object" != typeof e && "function" != typeof e) return { default: e }; var t = _getRequireWildcardCache(r); if (t && t.has(e)) return t.get(e); var n = { __proto__: null }, a = Object.defineProperty && Object.getOwnPropertyDescriptor; for (var u in e) if ("default" !== u && {}.hasOwnProperty.call(e, u)) { var i = a ? Object.getOwnPropertyDescriptor(e, u) : null; i && (i.get || i.set) ? Object.defineProperty(n, u, i) : n[u] = e[u]; } return n.default = e, t && t.set(e, n), n; }
+function _interopRequireWildcard(e, t) { if ("function" == typeof WeakMap) var r = new WeakMap(), n = new WeakMap(); return (_interopRequireWildcard = function (e, t) { if (!t && e && e.__esModule) return e; var o, i, f = { __proto__: null, default: e }; if (null === e || "object" != typeof e && "function" != typeof e) return f; if (o = t ? n : r) { if (o.has(e)) return o.get(e); o.set(e, f); } for (const t in e) "default" !== t && {}.hasOwnProperty.call(e, t) && ((i = (o = Object.defineProperty) && Object.getOwnPropertyDescriptor(e, t)) && (i.get || i.set) ? o(f, t, i) : f[t] = e[t]); return f; })(e, t); }
 /**
  * Create a mathjs instance from given factory functions and optionally config
  *

package/lib/cjs/defaultInstance.js

@@ -6,6 +6,5 @@ Object.defineProperty(exports, "__esModule", {
 exports.default = void 0;
 var all = _interopRequireWildcard(require("./factoriesAny.js"));
 var _create = require("./core/create.js");
-function _getRequireWildcardCache(e) { if ("function" != typeof WeakMap) return null; var r = new WeakMap(), t = new WeakMap(); return (_getRequireWildcardCache = function (e) { return e ? t : r; })(e); }
-function _interopRequireWildcard(e, r) { if (!r && e && e.__esModule) return e; if (null === e || "object" != typeof e && "function" != typeof e) return { default: e }; var t = _getRequireWildcardCache(r); if (t && t.has(e)) return t.get(e); var n = { __proto__: null }, a = Object.defineProperty && Object.getOwnPropertyDescriptor; for (var u in e) if ("default" !== u && {}.hasOwnProperty.call(e, u)) { var i = a ? Object.getOwnPropertyDescriptor(e, u) : null; i && (i.get || i.set) ? Object.defineProperty(n, u, i) : n[u] = e[u]; } return n.default = e, t && t.set(e, n), n; }
+function _interopRequireWildcard(e, t) { if ("function" == typeof WeakMap) var r = new WeakMap(), n = new WeakMap(); return (_interopRequireWildcard = function (e, t) { if (!t && e && e.__esModule) return e; var o, i, f = { __proto__: null, default: e }; if (null === e || "object" != typeof e && "function" != typeof e) return f; if (o = t ? n : r) { if (o.has(e)) return o.get(e); o.set(e, f); } for (const t in e) "default" !== t && {}.hasOwnProperty.call(e, t) && ((i = (o = Object.defineProperty) && Object.getOwnPropertyDescriptor(e, t)) && (i.get || i.set) ? o(f, t, i) : f[t] = e[t]); return f; })(e, t); }
 var _default = exports.default = (0, _create.create)(all);

package/lib/cjs/entry/allFactoriesAny.js

@@ -5,8 +5,7 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.all = void 0;
 var allFactories = _interopRequireWildcard(require("../factoriesAny.js"));
-function _getRequireWildcardCache(e) { if ("function" != typeof WeakMap) return null; var r = new WeakMap(), t = new WeakMap(); return (_getRequireWildcardCache = function (e) { return e ? t : r; })(e); }
-function _interopRequireWildcard(e, r) { if (!r && e && e.__esModule) return e; if (null === e || "object" != typeof e && "function" != typeof e) return { default: e }; var t = _getRequireWildcardCache(r); if (t && t.has(e)) return t.get(e); var n = { __proto__: null }, a = Object.defineProperty && Object.getOwnPropertyDescriptor; for (var u in e) if ("default" !== u && {}.hasOwnProperty.call(e, u)) { var i = a ? Object.getOwnPropertyDescriptor(e, u) : null; i && (i.get || i.set) ? Object.defineProperty(n, u, i) : n[u] = e[u]; } return n.default = e, t && t.set(e, n), n; }
+function _interopRequireWildcard(e, t) { if ("function" == typeof WeakMap) var r = new WeakMap(), n = new WeakMap(); return (_interopRequireWildcard = function (e, t) { if (!t && e && e.__esModule) return e; var o, i, f = { __proto__: null, default: e }; if (null === e || "object" != typeof e && "function" != typeof e) return f; if (o = t ? n : r) { if (o.has(e)) return o.get(e); o.set(e, f); } for (const t in e) "default" !== t && {}.hasOwnProperty.call(e, t) && ((i = (o = Object.defineProperty) && Object.getOwnPropertyDescriptor(e, t)) && (i.get || i.set) ? o(f, t, i) : f[t] = e[t]); return f; })(e, t); }
 // creating all factories here in a separate file is needed to get tree-shaking working
 
 const all = exports.all = allFactories;

package/lib/cjs/entry/allFactoriesNumber.js

@@ -5,8 +5,7 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.all = void 0;
 var allFactories = _interopRequireWildcard(require("../factoriesNumber.js"));
-function _getRequireWildcardCache(e) { if ("function" != typeof WeakMap) return null; var r = new WeakMap(), t = new WeakMap(); return (_getRequireWildcardCache = function (e) { return e ? t : r; })(e); }
-function _interopRequireWildcard(e, r) { if (!r && e && e.__esModule) return e; if (null === e || "object" != typeof e && "function" != typeof e) return { default: e }; var t = _getRequireWildcardCache(r); if (t && t.has(e)) return t.get(e); var n = { __proto__: null }, a = Object.defineProperty && Object.getOwnPropertyDescriptor; for (var u in e) if ("default" !== u && {}.hasOwnProperty.call(e, u)) { var i = a ? Object.getOwnPropertyDescriptor(e, u) : null; i && (i.get || i.set) ? Object.defineProperty(n, u, i) : n[u] = e[u]; } return n.default = e, t && t.set(e, n), n; }
+function _interopRequireWildcard(e, t) { if ("function" == typeof WeakMap) var r = new WeakMap(), n = new WeakMap(); return (_interopRequireWildcard = function (e, t) { if (!t && e && e.__esModule) return e; var o, i, f = { __proto__: null, default: e }; if (null === e || "object" != typeof e && "function" != typeof e) return f; if (o = t ? n : r) { if (o.has(e)) return o.get(e); o.set(e, f); } for (const t in e) "default" !== t && {}.hasOwnProperty.call(e, t) && ((i = (o = Object.defineProperty) && Object.getOwnPropertyDescriptor(e, t)) && (i.get || i.set) ? o(f, t, i) : f[t] = e[t]); return f; })(e, t); }
 // creating all factories here in a separate file is needed to get tree-shaking working
 
 const all = exports.all = allFactories;

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

@@ -10,5 +10,5 @@ const indexDocs = exports.indexDocs = {
   syntax: ['[start]', '[start:end]', '[start:step:end]', '[start1, start 2, ...]', '[start1:end1, start2:end2, ...]', '[start1:step1:end1, start2:step2:end2, ...]'],
   description: 'Create an index to get or replace a subset of a matrix',
   examples: ['A = [1, 2, 3; 4, 5, 6]', 'A[1, :]', 'A[1, 2] = 50', 'A[1:2, 1:2] = 1', 'B = [1, 2, 3]', 'B[B>1 and B<3]'],
-  seealso: ['bignumber', 'boolean', 'complex', 'matrix,', 'number', 'range', 'string', 'unit']
+  seealso: ['bignumber', 'boolean', 'complex', 'matrix', 'number', 'range', 'string', 'unit']
 };

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

@@ -110,6 +110,9 @@ var _conj = require("./function/complex/conj.js");
 var _im = require("./function/complex/im.js");
 var _re = require("./function/complex/re.js");
 var _evaluate = require("./function/expression/evaluate.js");
+var _parser = require("./function/expression/parser.js");
+var _parse = require("./function/expression/parse.js");
+var _compile = require("./function/expression/compile.js");
 var _help = require("./function/expression/help.js");
 var _distance = require("./function/geometry/distance.js");
 var _intersect = require("./function/geometry/intersect.js");
@@ -570,6 +573,9 @@ const embeddedDocs = exports.embeddedDocs = {
   // functions - expression
   evaluate: _evaluate.evaluateDocs,
   help: _help.helpDocs,
+  parse: _parse.parseDocs,
+  parser: _parser.parserDocs,
+  compile: _compile.compileDocs,
   // functions - geometry
   distance: _distance.distanceDocs,
   intersect: _intersect.intersectDocs,

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

@@ -0,0 +1,14 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.compileDocs = void 0;
+const compileDocs = exports.compileDocs = {
+  name: 'compile',
+  category: 'Expression',
+  syntax: ['compile(expr) ', 'compile([expr1, expr2, expr3, ...])'],
+  description: 'Parse and compile an expression. Returns a an object with a function evaluate([scope]) to evaluate the compiled expression.',
+  examples: ['code1 = compile("sqrt(3^2 + 4^2)")', 'code1.evaluate() ', 'code2 = compile("a * b")', 'code2.evaluate({a: 3, b: 4})'],
+  seealso: ['parser', 'parse', 'evaluate']
+};

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

@@ -10,5 +10,5 @@ const evaluateDocs = exports.evaluateDocs = {
   syntax: ['evaluate(expression)', 'evaluate(expression, scope)', 'evaluate([expr1, expr2, expr3, ...])', 'evaluate([expr1, expr2, expr3, ...], scope)'],
   description: 'Evaluate an expression or an array with expressions.',
   examples: ['evaluate("2 + 3")', 'evaluate("sqrt(16)")', 'evaluate("2 inch to cm")', 'evaluate("sin(x * pi)", { "x": 1/2 })', 'evaluate(["width=2", "height=4","width*height"])'],
-  seealso: []
+  seealso: ['parser', 'parse', 'compile']
 };

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

@@ -0,0 +1,14 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.parseDocs = void 0;
+const parseDocs = exports.parseDocs = {
+  name: 'parse',
+  category: 'Expression',
+  syntax: ['parse(expr)', 'parse(expr, options)', 'parse([expr1, expr2, expr3, ...])', 'parse([expr1, expr2, expr3, ...], options)'],
+  description: 'Parse an expression. Returns a node tree, which can be evaluated by invoking node.evaluate() or transformed into a functional object via node.compile().',
+  examples: ['node1 = parse("sqrt(3^2 + 4^2)")', 'node1.evaluate()', 'code1 = node1.compile()', 'code1.evaluate()', 'scope = {a: 3, b: 4}', 'node2 = parse("a * b")', 'node2.evaluate(scope)', 'code2 = node2.compile()', 'code2.evaluate(scope)'],
+  seealso: ['parser', 'evaluate', 'compile']
+};

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

@@ -0,0 +1,14 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.parserDocs = void 0;
+const parserDocs = exports.parserDocs = {
+  name: 'parser',
+  category: 'Expression',
+  syntax: ['parser()'],
+  description: 'Create a parser object that keeps a context of variables and their values, allowing the evaluation of expressions in that context.',
+  examples: ['myParser = parser()', 'myParser.evaluate("sqrt(3^2 + 4^2)")', 'myParser.set("x", 3)', 'myParser.evaluate("y = x + 3")', 'myParser.evaluate(["y = x + 3", "y = y + 1"])', 'myParser.get("y")'],
+  seealso: ['evaluate', 'parse', 'compile']
+};

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

@@ -10,5 +10,5 @@ const sumDocs = exports.sumDocs = {
   syntax: ['sum(a, b, c, ...)', 'sum(A)', 'sum(A, dimension)'],
   description: 'Compute the sum of all values.',
   examples: ['sum(2, 3, 4, 1)', 'sum([2, 3, 4, 1])', 'sum([2, 5; 4, 3])'],
-  seealso: ['max', 'mean', 'median', 'min', 'prod', 'std', 'sum', 'variance']
+  seealso: ['max', 'mean', 'median', 'min', 'prod', 'std', 'variance']
 };

package/lib/cjs/expression/function/parser.js

@@ -13,7 +13,7 @@ const createParser = exports.createParser = /* #__PURE__ */(0, _factory.factory)
     Parser
   } = _ref;
   /**
-   * Create a parser. The function creates a new `math.Parser` object.
+   * Create a `math.Parser` object that keeps a context of variables and their values, allowing the evaluation of expressions in that context.
    *
    * Syntax:
    *

package/lib/cjs/expression/parse.js

@@ -36,7 +36,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
   } = _ref;
   /**
    * Parse an expression. Returns a node tree, which can be evaluated by
-   * invoking node.evaluate().
+   * invoking node.evaluate() or transformed into a functional object via node.compile().
    *
    * Note the evaluating arbitrary expressions may involve security risks,
    * see [https://mathjs.org/docs/expressions/security.html](https://mathjs.org/docs/expressions/security.html) for more information.
@@ -55,6 +55,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
    *
    *     let scope = {a:3, b:4}
    *     const node2 = math.parse('a * b') // 12
+   *     node2.evaluate(scope) // 12
    *     const code2 = node2.compile()
    *     code2.evaluate(scope) // 12
    *     scope.a = 5
@@ -877,7 +878,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
       name = state.token;
       fn = operators[name];
       getTokenSkipNewline(state);
-      if (name === 'in' && state.token === '') {
+      if (name === 'in' && '])},;'.includes(state.token)) {
         // end of expression -> this is the unit 'in' ('inch')
         node = new OperatorNode('*', 'multiply', [node, new SymbolNode('in')], true);
       } else {

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

@@ -38,39 +38,39 @@ const createAdd = exports.createAdd = /* #__PURE__ */(0, _factory.factory)(name,
     concat
   });
   /**
-  * Add two or more values, `x + y`.
-  * For matrices, the function is evaluated element wise.
-  *
-  * Syntax:
-  *
-  *    math.add(x, y)
-  *    math.add(x, y, z, ...)
-  *
-  * Examples:
-  *
-  *    math.add(2, 3)               // returns number 5
-  *    math.add(2, 3, 4)            // returns number 9
-  *
-  *    const a = math.complex(2, 3)
-  *    const b = math.complex(-4, 1)
-  *    math.add(a, b)               // returns Complex -2 + 4i
-  *
-  *    math.add([1, 2, 3], 4)       // returns Array [5, 6, 7]
-  *
-  *    const c = math.unit('5 cm')
-  *    const d = math.unit('2.1 mm')
-  *    math.add(c, d)               // returns Unit 52.1 mm
-  *
-  *    math.add("2.3", "4")         // returns number 6.3
-  *
-  * See also:
-  *
-  *    subtract, sum
-  *
-  * @param  {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} x First value to add
-  * @param  {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} y Second value to add
-  * @return {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} Sum of `x` and `y`
-  */
+   * Add two or more values, `x + y`.
+   * For matrices, the function is evaluated element wise.
+   *
+   * Syntax:
+   *
+   *    math.add(x, y)
+   *    math.add(x, y, z, ...)
+   *
+   * Examples:
+   *
+   *    math.add(2, 3)               // returns number 5
+   *    math.add(2, 3, 4)            // returns number 9
+   *
+   *    const a = math.complex(2, 3)
+   *    const b = math.complex(-4, 1)
+   *    math.add(a, b)               // returns Complex -2 + 4i
+   *
+   *    math.add([1, 2, 3], 4)       // returns Array [5, 6, 7]
+   *
+   *    const c = math.unit('5 cm')
+   *    const d = math.unit('2.1 mm')
+   *    math.add(c, d)               // returns Unit 52.1 mm
+   *
+   *    math.add("2.3", "4")         // returns number 6.3
+   *
+   * See also:
+   *
+   *    subtract, sum
+   *
+   * @param  {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} x First value to add
+   * @param  {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} y Second value to add
+   * @return {number | BigNumber | bigint | Fraction | Complex | Unit | Array | Matrix} Sum of `x` and `y`
+   */
   return typed(name, {
     'any, any': addScalar,
     'any, any, ...any': typed.referToSelf(self => (x, y, rest) => {

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.4.0
- * @date    2025-03-28
+ * @version 14.5.0
+ * @date    2025-05-21
  *
  * @license
  * Copyright (C) 2013-2025 Jos de Jong <wjosdejong@gmail.com>

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

@@ -248,36 +248,46 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
       }
 
       // retrieve submatrix
-      // TODO: more efficient when creating an empty matrix and setting _data and _size manually
-      return new DenseMatrix(_getSubmatrix(matrix._data, index, size.length, 0), matrix._datatype);
+      const returnMatrix = new DenseMatrix([]);
+      const submatrix = _getSubmatrix(matrix._data, index);
+      returnMatrix._size = submatrix.size;
+      returnMatrix._datatype = matrix._datatype;
+      returnMatrix._data = submatrix.data;
+      return returnMatrix;
     }
   }
 
   /**
-   * Recursively get a submatrix of a multi dimensional matrix.
+   * Get a submatrix of a multi dimensional matrix.
    * Index is not checked for correct number or length of dimensions.
    * @memberof DenseMatrix
    * @param {Array} data
    * @param {Index} index
-   * @param {number} dims   Total number of dimensions
-   * @param {number} dim    Current dimension
    * @return {Array} submatrix
    * @private
    */
-  function _getSubmatrix(data, index, dims, dim) {
-    const last = dim === dims - 1;
-    const range = index.dimension(dim);
-    if (last) {
-      return range.map(function (i) {
-        (0, _array.validateIndex)(i, data.length);
-        return data[i];
-      }).valueOf();
-    } else {
-      return range.map(function (i) {
-        (0, _array.validateIndex)(i, data.length);
-        const child = data[i];
-        return _getSubmatrix(child, index, dims, dim + 1);
-      }).valueOf();
+  function _getSubmatrix(data, index) {
+    const maxDepth = index.size().length - 1;
+    const size = Array(maxDepth);
+    return {
+      data: getSubmatrixRecursive(data),
+      size
+    };
+    function getSubmatrixRecursive(data) {
+      let depth = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+      const ranges = index.dimension(depth);
+      size[depth] = ranges.size()[0];
+      if (depth < maxDepth) {
+        return ranges.map(rangeIndex => {
+          (0, _array.validateIndex)(rangeIndex, data.length);
+          return getSubmatrixRecursive(data[rangeIndex], depth + 1);
+        }).valueOf();
+      } else {
+        return ranges.map(rangeIndex => {
+          (0, _array.validateIndex)(rangeIndex, data.length);
+          return data[rangeIndex];
+        }).valueOf();
+      }
     }
   }
 
@@ -365,9 +375,7 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
       _fit(matrix, size, defaultValue);
 
       // insert the sub matrix
-      const dims = iSize.length;
-      const dim = 0;
-      _setSubmatrix(matrix._data, index, submatrix, dims, dim);
+      _setSubmatrix(matrix._data, index, submatrix);
     }
     return matrix;
   }
@@ -378,23 +386,25 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
    * @param {Array} data
    * @param {Index} index
    * @param {Array} submatrix
-   * @param {number} dims   Total number of dimensions
-   * @param {number} dim
    * @private
    */
-  function _setSubmatrix(data, index, submatrix, dims, dim) {
-    const last = dim === dims - 1;
-    const range = index.dimension(dim);
-    if (last) {
-      range.forEach(function (dataIndex, subIndex) {
-        (0, _array.validateIndex)(dataIndex);
-        data[dataIndex] = submatrix[subIndex[0]];
-      });
-    } else {
-      range.forEach(function (dataIndex, subIndex) {
-        (0, _array.validateIndex)(dataIndex);
-        _setSubmatrix(data[dataIndex], index, submatrix[subIndex[0]], dims, dim + 1);
-      });
+  function _setSubmatrix(data, index, submatrix) {
+    const maxDepth = index.size().length - 1;
+    setSubmatrixRecursive(data, submatrix);
+    function setSubmatrixRecursive(data, submatrix) {
+      let depth = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : 0;
+      const range = index.dimension(depth);
+      if (depth < maxDepth) {
+        range.forEach((rangeIndex, i) => {
+          (0, _array.validateIndex)(rangeIndex, data.length);
+          setSubmatrixRecursive(data[rangeIndex], submatrix[i[0]], depth + 1);
+        });
+      } else {
+        range.forEach((rangeIndex, i) => {
+          (0, _array.validateIndex)(rangeIndex, data.length);
+          data[rangeIndex] = submatrix[i[0]];
+        });
+      }
     }
   }
 
@@ -525,84 +535,75 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
   };
 
   /**
-   * Applies a callback function to a reference to each element of the matrix
+   * Create a new matrix with the results of the callback function executed on
+   * each entry of the matrix.
    * @memberof DenseMatrix
    * @param {Function} callback   The callback function is invoked with three
-   *                              parameters: the array containing the element,
-   *                              the index of the element within that array (as an integer),
-   *                              and for non unarry callbacks copy of the current index (as an array of integers).
+   *                              parameters: the value of the element, the index
+   *                              of the element, and the Matrix being traversed.
+   * @param {boolean} skipZeros   If true, the callback function is invoked only for non-zero entries
+   * @param {boolean} isUnary     If true, the callback function is invoked with one parameter
+   *
+   * @return {DenseMatrix} matrix
    */
-  DenseMatrix.prototype._forEach = function (callback) {
-    const isUnary = callback.length === 2; // callback has 2 parameters: value, index
-    const maxDepth = this._size.length - 1;
-    if (maxDepth < 0) return;
-    if (isUnary) {
-      iterateUnary(this._data);
-      return;
+  DenseMatrix.prototype.map = function (callback) {
+    let skipZeros = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : false;
+    let isUnary = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : false;
+    const me = this;
+    const maxDepth = me._size.length - 1;
+    if (maxDepth < 0) return me.clone();
+    const fastCallback = (0, _optimizeCallback.optimizeCallback)(callback, me, 'map', isUnary);
+    const fastCallbackFn = fastCallback.fn;
+    const result = me.create(undefined, me._datatype);
+    result._size = me._size;
+    if (isUnary || fastCallback.isUnary) {
+      result._data = iterateUnary(me._data);
+      return result;
     }
     if (maxDepth === 0) {
-      for (let i = 0; i < this._data.length; i++) {
-        callback(this._data, i, [i]);
+      const inputData = me.valueOf();
+      const data = Array(inputData.length);
+      for (let i = 0; i < inputData.length; i++) {
+        data[i] = fastCallbackFn(inputData[i], [i], me);
       }
-      return;
+      result._data = data;
+      return result;
     }
-    const index = new Array(maxDepth + 1);
-    iterate(this._data);
+    const index = [];
+    result._data = iterate(me._data);
+    return result;
     function iterate(data) {
       let depth = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+      const result = Array(data.length);
       if (depth < maxDepth) {
         for (let i = 0; i < data.length; i++) {
           index[depth] = i;
-          iterate(data[i], depth + 1);
+          result[i] = iterate(data[i], depth + 1);
         }
       } else {
         for (let i = 0; i < data.length; i++) {
           index[depth] = i;
-          callback(data, i, index.slice());
+          result[i] = fastCallbackFn(data[i], index.slice(), me);
         }
       }
+      return result;
     }
     function iterateUnary(data) {
       let depth = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+      const result = Array(data.length);
       if (depth < maxDepth) {
         for (let i = 0; i < data.length; i++) {
-          iterateUnary(data[i], depth + 1);
+          result[i] = iterateUnary(data[i], depth + 1);
         }
       } else {
         for (let i = 0; i < data.length; i++) {
-          callback(data, i);
+          result[i] = fastCallbackFn(data[i]);
         }
       }
+      return result;
     }
   };
 
-  /**
-   * Create a new matrix with the results of the callback function executed on
-   * each entry of the matrix.
-   * @memberof DenseMatrix
-   * @param {Function} callback   The callback function is invoked with three
-   *                              parameters: the value of the element, the index
-   *                              of the element, and the Matrix being traversed.
-   * @param {boolean} skipZeros   If true, the callback function is invoked only for non-zero entries
-   * @param {boolean} isUnary     If true, the callback function is invoked with one parameter
-   *
-   * @return {DenseMatrix} matrix
-   */
-  DenseMatrix.prototype.map = function (callback) {
-    let skipZeros = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : false;
-    let isUnary = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : false;
-    const me = this;
-    const result = new DenseMatrix(me);
-    const fastCallback = (0, _optimizeCallback.optimizeCallback)(callback, me._data, 'map', isUnary);
-    const applyCallback = isUnary || fastCallback.isUnary ? (arr, i) => {
-      arr[i] = fastCallback.fn(arr[i]);
-    } : (arr, i, index) => {
-      arr[i] = fastCallback.fn(arr[i], index, me);
-    };
-    result._forEach(applyCallback);
-    return result;
-  };
-
   /**
    * Execute a callback function on each entry of the matrix.
    * @memberof DenseMatrix
@@ -616,13 +617,48 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
     let skipZeros = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : false;
     let isUnary = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : false;
     const me = this;
-    const fastCallback = (0, _optimizeCallback.optimizeCallback)(callback, me._data, 'map', isUnary);
-    const applyCallback = isUnary || fastCallback.isUnary ? (arr, i) => {
-      fastCallback.fn(arr[i]);
-    } : (arr, i, index) => {
-      fastCallback.fn(arr[i], index, me);
-    };
-    me._forEach(applyCallback);
+    const maxDepth = me._size.length - 1;
+    if (maxDepth < 0) return;
+    const fastCallback = (0, _optimizeCallback.optimizeCallback)(callback, me, 'map', isUnary);
+    const fastCallbackFn = fastCallback.fn;
+    if (isUnary || fastCallback.isUnary) {
+      iterateUnary(me._data);
+      return;
+    }
+    if (maxDepth === 0) {
+      for (let i = 0; i < me._data.length; i++) {
+        fastCallbackFn(me._data[i], [i], me);
+      }
+      return;
+    }
+    const index = [];
+    iterate(me._data);
+    function iterate(data) {
+      let depth = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+      if (depth < maxDepth) {
+        for (let i = 0; i < data.length; i++) {
+          index[depth] = i;
+          iterate(data[i], depth + 1);
+        }
+      } else {
+        for (let i = 0; i < data.length; i++) {
+          index[depth] = i;
+          fastCallbackFn(data[i], index.slice(), me);
+        }
+      }
+    }
+    function iterateUnary(data) {
+      let depth = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0;
+      if (depth < maxDepth) {
+        for (let i = 0; i < data.length; i++) {
+          iterateUnary(data[i], depth + 1);
+        }
+      } else {
+        for (let i = 0; i < data.length; i++) {
+          fastCallbackFn(data[i]);
+        }
+      }
+    }
   };
 
   /**