Performance issues with deepMap

src/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'
 
 const name = 'forEach'
 const dependencies = ['typed']
@@ -45,5 +45,5 @@ export const createForEach = /* #__PURE__ */ factory(name, dependencies, ({ type
  * @private
  */
 function _forEach (array, callback) {
-  recurse(array, [], array, optimizeCallback(callback, array, name))
+  deepForEach(array, array, optimizeCallback(callback, array, name))
 }

src/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'
 
 const name = 'map'
@@ -151,6 +151,6 @@ export const createMap = /* #__PURE__ */ factory(name, dependencies, ({ typed })
  * @private
  */
   function _mapArray (array, callback) {
-    return recurse(array, [], array, optimizeCallback(callback, array, name))
+    return deepMap(array, array, optimizeCallback(callback, array, name))
   }
 })

src/type/matrix/DenseMatrix.js

@@ -1,5 +1,5 @@
 import { isArray, isBigNumber, isCollection, isIndex, isMatrix, isNumber, isString, typeOf } from '../../utils/is.js'
-import { arraySize, getArrayDataType, processSizesWildcard, reshape, resize, unsqueeze, validate, validateIndex, broadcastTo, get, recurse } from '../../utils/array.js'
+import { arraySize, getArrayDataType, processSizesWildcard, reshape, resize, unsqueeze, validate, validateIndex, broadcastTo, get, deepMap, deepForEach } from '../../utils/array.js'
 import { format } from '../../utils/string.js'
 import { isInteger } from '../../utils/number.js'
 import { clone, deepStrictEqual } from '../../utils/object.js'
@@ -541,7 +541,7 @@ export const createDenseMatrixClass = /* #__PURE__ */ factory(name, dependencies
 
     // determine the new datatype when the original matrix has datatype defined
     // TODO: should be done in matrix constructor instead
-    const data = recurse(me._data, [], me, fastCallback)
+    const data = deepMap(me._data, me, fastCallback)
     const datatype = me._datatype !== undefined
       ? getArrayDataType(data, typeOf)
       : undefined
@@ -559,7 +559,7 @@ export const createDenseMatrixClass = /* #__PURE__ */ factory(name, dependencies
     // matrix instance
     const me = this
     const fastCallback = optimizeCallback(callback, me._data, 'forEach')
-    recurse(this._data, [], me, fastCallback)
+    deepForEach(this._data, me, fastCallback)
   }
 
   /**

src/utils/array.js

@@ -4,6 +4,7 @@ import { format } from './string.js'
 import { DimensionError } from '../error/DimensionError.js'
 import { IndexError } from '../error/IndexError.js'
 import { deepStrictEqual } from './object.js'
+import { findNumberOfArguments } from './optimizeCallback.js'
 
 /**
  * Calculate the size of a multi dimensional array.
@@ -834,15 +835,124 @@ export function get (array, index) {
  * @param {Function} callback - Function that produces the element of the new Array, taking three arguments: the value of the element, the index of the element, and the Array being processed.
  * @returns {*} The new array with each element being the result of the callback function.
  */
-export function recurse (value, index, array, callback) {
-  if (Array.isArray(value)) {
-    return value.map(function (child, i) {
-      // we create a copy of the index array and append the new index value
-      return recurse(child, index.concat(i), array, callback)
-    })
-  } else {
-    // invoke the callback function with the right number of arguments
-    return callback(value, index, array)
+export function deepMap (value, array, callback) {
+  const numberOfArguments = findNumberOfArguments(callback, array)
+  switch (numberOfArguments) {
+    case 1:
+      return recurse1(value)
+    case 2:
+      return recurse2(value, [])
+    case 3:
+      return recurse3(value, [])
+    default:
+      return recurse3(value, [])
+  }
+
+  function recurse1 (value) {
+    if (Array.isArray(value)) {
+      return value.map(function (child) {
+        // we create a copy of the index array and append the new index value
+        const results = recurse1(child)
+        return results
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      return callback(value)
+    }
+  }
+
+  function recurse2 (value, index) {
+    if (Array.isArray(value)) {
+      return value.map(function (child, i) {
+        // we create a copy of the index array and append the new index value
+        index.push(i)
+        const results = recurse2(child, index)
+        index.pop()
+        return results
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      return callback(value, [...index])
+    }
+  }
+
+  function recurse3 (value, index) {
+    if (Array.isArray(value)) {
+      return value.map(function (child, i) {
+        // we create a copy of the index array and append the new index value
+        index.push(i)
+        const results = recurse3(child, index)
+        index.pop()
+        return results
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      return callback(value, [...index], array)
+    }
+  }
+}
+
+/**
+ * Recursive function to map a multi-dimensional array.
+ *
+ * @param {*} value - The current value being processed in the array.
+ * @param {Array} index - The index of the current value being processed in the array.
+ * @param {Array} array - The array being processed.
+ * @param {Function} callback - Function that produces the element of the new Array, taking three arguments: the value of the element, the index of the element, and the Array being processed.
+ * @returns {*} The new array with each element being the result of the callback function.
+ */
+export function deepForEach (value, array, callback) {
+  const numberOfArguments = findNumberOfArguments(callback, array)
+  switch (numberOfArguments) {
+    case 1:
+      recurse1(value)
+      break
+    case 2:
+      recurse2(value, [])
+      break
+    case 3:
+      recurse3(value, [])
+      break
+    default:
+      recurse3(value, [])
+      break
+  }
+
+  function recurse1 (value) {
+    if (Array.isArray(value)) {
+      value.forEach(function (child) {
+        recurse1(child)
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      callback(value)
+    }
+  }
+
+  function recurse2 (value, index) {
+    if (Array.isArray(value)) {
+      value.forEach(function (child, i) {
+        index.push(i)
+        recurse2(child, index)
+        index.pop()
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      callback(value, [...index])
+    }
+  }
+
+  function recurse3 (value, index) {
+    if (Array.isArray(value)) {
+      value.forEach(function (child, i) {
+        index.push(i)
+        recurse3(child, index)
+        index.pop()
+      })
+    } else {
+      // invoke the callback function with the right number of arguments
+      callback(value, [...index], array)
+    }
   }
 }
 

src/utils/collection.js

@@ -27,16 +27,15 @@ export function containsCollections (array) {
  */
 export function deepForEach (array, callback) {
   if (isMatrix(array)) {
-    array = array.valueOf()
+    recurse(array.valueOf())
+  } else {
+    recurse(array)
   }
-
-  for (let i = 0, ii = array.length; i < ii; i++) {
-    const value = array[i]
-
-    if (Array.isArray(value)) {
-      deepForEach(value, callback)
+  function recurse (array) {
+    if (Array.isArray(array)) {
+      array.forEach(value => recurse(value))
     } else {
-      callback(value)
+      callback(array)
     }
   }
 }
@@ -54,13 +53,28 @@ export function deepForEach (array, callback) {
  * @return {Array | Matrix} res
  */
 export function deepMap (array, callback, skipZeros) {
-  if (array && (typeof array.map === 'function')) {
-    // TODO: replace array.map with a for loop to improve performance
-    return array.map(function (x) {
-      return deepMap(x, callback, skipZeros)
-    })
+  if (skipZeros) {
+    const callbackSkip = (x) => x === 0 ? 0 : callback(x)
+    if (isMatrix(array)) {
+      return array.create(recurse(array.valueOf(), callbackSkip), array.datatype())
+    } else {
+      return recurse(array, callbackSkip)
+    }
   } else {
-    return callback(array)
+    if (isMatrix(array)) {
+      return array.create(recurse(array.valueOf(), callback), array.datatype())
+    } else {
+      return recurse(array, callback)
+    }
+  }
+  function recurse (array) {
+    if (Array.isArray(array)) {
+      return array.map(function (x) {
+        return recurse(x)
+      })
+    } else {
+      return callback(array)
+    }
   }
 }
 

src/utils/optimizeCallback.js

@@ -15,17 +15,27 @@ export function optimizeCallback (callback, array, name) {
     const firstIndex = (array.isMatrix ? array.size() : arraySize(array)).map(() => 0)
     const firstValue = array.isMatrix ? array.get(firstIndex) : get(array, firstIndex)
     const hasSingleSignature = Object.keys(callback.signatures).length === 1
-    const numberOfArguments = _findNumberOfArguments(callback, firstValue, firstIndex, array)
+    const numberOfArguments = _typedFindNumberOfArguments(callback, firstValue, firstIndex, array)
     const fastCallback = hasSingleSignature ? Object.values(callback.signatures)[0] : callback
     if (numberOfArguments >= 1 && numberOfArguments <= 3) {
-      return (...args) => _tryFunctionWithArgs(fastCallback, args.slice(0, numberOfArguments), name, callback.name)
+      return (...args) => tryFunctionWithArgs(fastCallback, args.slice(0, numberOfArguments), name, callback.name)
     }
-    return (...args) => _tryFunctionWithArgs(fastCallback, args, name, callback.name)
+    return (...args) => tryFunctionWithArgs(fastCallback, args, name, callback.name)
   }
   return callback
 }
 
-function _findNumberOfArguments (callback, value, index, array) {
+export function findNumberOfArguments (callback, array) {
+  if (typed.isTypedFunction(callback)) {
+    const firstIndex = (array.isMatrix ? array.size() : arraySize(array)).map(() => 0)
+    const firstValue = array.isMatrix ? array.get(firstIndex) : get(array, firstIndex)
+    return _typedFindNumberOfArguments(callback, firstValue, firstIndex, array)
+  } else {
+    return callback.length
+  }
+}
+
+function _typedFindNumberOfArguments (callback, value, index, array) {
   const testArgs = [value, index, array]
   for (let i = 3; i > 0; i--) {
     const args = testArgs.slice(0, i)
@@ -43,7 +53,7 @@ function _findNumberOfArguments (callback, value, index, array) {
    * @returns {*} Returns the return value of the invoked signature
    * @throws {TypeError} Throws an error when no matching signature was found
    */
-function _tryFunctionWithArgs (func, args, mappingFnName, callbackName) {
+export function tryFunctionWithArgs (func, args, mappingFnName, callbackName) {
   try {
     return func(...args)
   } catch (err) {