mathjs 11.5.1 → 11.6.0

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

@@ -1,7 +1,7 @@
 import { arraySize } from '../../utils/array.js';
 import { factory } from '../../utils/factory.js';
 var name = 'fft';
-var dependencies = ['typed', 'matrix', 'addScalar', 'multiplyScalar', 'divideScalar', 'exp', 'tau', 'i'];
+var dependencies = ['typed', 'matrix', 'addScalar', 'multiplyScalar', 'divideScalar', 'exp', 'tau', 'i', 'dotDivide', 'conj', 'pow', 'ceil', 'log2'];
 export var createFft = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
@@ -11,7 +11,12 @@ export var createFft = /* #__PURE__ */factory(name, dependencies, _ref => {
     divideScalar,
     exp,
     tau,
-    i: I
+    i: I,
+    dotDivide,
+    conj,
+    pow,
+    ceil,
+    log2
   } = _ref;
   /**
    * Calculate N-dimensional fourier transform
@@ -70,7 +75,32 @@ export var createFft = /* #__PURE__ */factory(name, dependencies, _ref => {
     }
     return _transpose(_1dFft(_transpose(arr), 1));
   }
-
+  /**
+   * Perform an 1-dimensional non-power-of-2 Fourier transform using Chirp-Z Transform
+   *
+   * @param {Array} arr      The array
+   * @return {Array}         resulting array
+   */
+  function _czt(arr) {
+    var n = arr.length;
+    var w = exp(divideScalar(multiplyScalar(-1, multiplyScalar(I, tau)), n));
+    var chirp = [];
+    for (var i = 1 - n; i < n; i++) {
+      chirp.push(pow(w, divideScalar(pow(i, 2), 2)));
+    }
+    var N2 = pow(2, ceil(log2(n + n - 1)));
+    var xp = [...new Array(n).fill(0).map((_, i) => multiplyScalar(arr[i], chirp[n - 1 + i])), ...new Array(N2 - n).fill(0)];
+    var ichirp = [...new Array(n + n - 1).fill(0).map((_, i) => divideScalar(1, chirp[i])), ...new Array(N2 - (n + n - 1)).fill(0)];
+    var fftXp = _fft(xp);
+    var fftIchirp = _fft(ichirp);
+    var fftProduct = new Array(N2).fill(0).map((_, i) => multiplyScalar(fftXp[i], fftIchirp[i]));
+    var ifftProduct = dotDivide(conj(_ndFft(conj(fftProduct))), N2);
+    var ret = [];
+    for (var _i = n - 1; _i < n + n - 1; _i++) {
+      ret.push(multiplyScalar(ifftProduct[_i], chirp[_i]));
+    }
+    return ret;
+  }
   /**
    * Perform an 1-dimensional Fourier transform
    *
@@ -89,7 +119,10 @@ export var createFft = /* #__PURE__ */factory(name, dependencies, _ref => {
         ret[k + len / 2] = addScalar(p, multiplyScalar(-1, q));
       }
       return ret;
+    } else {
+      // use chirp-z transform for non-power-of-2 FFT
+      return _czt(arr);
     }
-    throw new Error('Can only calculate FFT of power-of-two size');
+    // throw new Error('Can only calculate FFT of power-of-two size')
   }
 });

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

@@ -21,7 +21,7 @@ export var createGetMatrixDataType = /* #__PURE__ */factory(name, dependencies,
    *
    *    const x = [ [1, 2, 3], [4, 5, 6] ]
    *    const mixedX = [ [1, true], [2, 3] ]
-   *    const fractionX = [ [math.fraction(1, 3)], [math.fraction(1, 3] ]
+   *    const fractionX = [ [math.fraction(1, 3)], [math.fraction(1, 3)] ]
    *    const unitX = [ [math.unit('5cm')], [math.unit('5cm')] ]
    *    const bigNumberX = [ [math.bignumber(1)], [math.bignumber(0)] ]
    *    const sparse = math.sparse(x)

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

@@ -1,4 +1,5 @@
 import { factory } from '../../utils/factory.js';
+import { isMatrix } from '../../utils/is.js';
 import { clone } from '../../utils/object.js';
 import { validateIndex } from '../../utils/array.js';
 var name = 'row';
@@ -52,6 +53,7 @@ export var createRow = /* #__PURE__ */factory(name, dependencies, _ref => {
     validateIndex(row, value.size()[0]);
     var columnRange = range(0, value.size()[1]);
     var index = new Index(row, columnRange);
-    return value.subset(index);
+    var result = value.subset(index);
+    return isMatrix(result) ? result : matrix([[result]]);
   }
 });

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

@@ -27,7 +27,7 @@ export var createSubset = /* #__PURE__ */factory(name, dependencies, _ref => {
    *
    *     // replace a subset
    *     const e = []
-   *     const f = math.subset(e, math.index(0, [0, 2]), [5, 6])  // f = [[5, 6]] and e = [[5, 0, 6]]
+   *     const f = math.subset(e, math.index(0, [0, 2]), [5, 6])  // f = [[5, 6]]
    *     const g = math.subset(f, math.index(1, 1), 7, 0)         // g = [[5, 6], [0, 7]]
    *
    *     // get submatrix using ranges

package/lib/esm/function/relational/compare.js

@@ -7,7 +7,7 @@ import { createMatAlgo05xSfSf } from '../../type/matrix/utils/matAlgo05xSfSf.js'
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 import { createCompareUnits } from './compareUnits.js';
 var name = 'compare';
-var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'BigNumber', 'Fraction', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'BigNumber', 'Fraction', 'DenseMatrix', 'concat'];
 export var createCompare = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
@@ -16,7 +16,8 @@ export var createCompare = /* #__PURE__ */factory(name, dependencies, _ref => {
     matrix,
     BigNumber,
     Fraction,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -31,7 +32,8 @@ export var createCompare = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
   var compareUnits = createCompareUnits({
     typed

package/lib/esm/function/relational/compareText.js

@@ -2,16 +2,18 @@ import { compareText as _compareText } from '../../utils/string.js';
 import { factory } from '../../utils/factory.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'compareText';
-var dependencies = ['typed', 'matrix'];
+var dependencies = ['typed', 'matrix', 'concat'];
 _compareText.signature = 'any, any';
 export var createCompareText = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
-    matrix
+    matrix,
+    concat
   } = _ref;
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
 
   /**

package/lib/esm/function/relational/equal.js

@@ -4,13 +4,14 @@ import { createMatAlgo07xSSf } from '../../type/matrix/utils/matAlgo07xSSf.js';
 import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'equal';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'equalScalar', 'DenseMatrix', 'concat'];
 export var createEqual = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     equalScalar,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -25,7 +26,8 @@ export var createEqual = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
 
   /**

package/lib/esm/function/relational/larger.js

@@ -7,13 +7,14 @@ import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 import { createCompareUnits } from './compareUnits.js';
 var name = 'larger';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 export var createLarger = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     matrix,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -28,7 +29,8 @@ export var createLarger = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
   var compareUnits = createCompareUnits({
     typed

package/lib/esm/function/relational/largerEq.js

@@ -7,13 +7,14 @@ import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 import { createCompareUnits } from './compareUnits.js';
 var name = 'largerEq';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 export var createLargerEq = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     matrix,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -28,7 +29,8 @@ export var createLargerEq = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
   var compareUnits = createCompareUnits({
     typed

package/lib/esm/function/relational/smaller.js

@@ -7,13 +7,14 @@ import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 import { createCompareUnits } from './compareUnits.js';
 var name = 'smaller';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 export var createSmaller = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     matrix,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -28,7 +29,8 @@ export var createSmaller = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
   var compareUnits = createCompareUnits({
     typed

package/lib/esm/function/relational/smallerEq.js

@@ -7,13 +7,14 @@ import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 import { createCompareUnits } from './compareUnits.js';
 var name = 'smallerEq';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 export var createSmallerEq = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     matrix,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -28,7 +29,8 @@ export var createSmallerEq = /* #__PURE__ */factory(name, dependencies, _ref =>
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
   var compareUnits = createCompareUnits({
     typed

package/lib/esm/function/relational/unequal.js

@@ -4,14 +4,15 @@ import { createMatAlgo07xSSf } from '../../type/matrix/utils/matAlgo07xSSf.js';
 import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'unequal';
-var dependencies = ['typed', 'config', 'equalScalar', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'equalScalar', 'matrix', 'DenseMatrix', 'concat'];
 export var createUnequal = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     config,
     equalScalar,
     matrix,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo03xDSf = createMatAlgo03xDSf({
     typed
@@ -26,7 +27,8 @@ export var createUnequal = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
 
   /**

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

@@ -29,7 +29,7 @@ export var createMax = /* #__PURE__ */factory(name, dependencies, _ref => {
    *
    *     // maximum over a specified dimension (zero-based)
    *     math.max([[2, 5], [4, 3], [1, 7]], 0) // returns [4, 7]
-   *     math.max([[2, 5], [4, 3]], [1, 7], 1) // returns [5, 4, 7]
+   *     math.max([[2, 5], [4, 3], [1, 7]], 1) // returns [5, 4, 7]
    *
    *     math.max(2.7, 7.1, -4.5, 2.0, 4.1)    // returns 7.1
    *     math.min(2.7, 7.1, -4.5, 2.0, 4.1)    // returns -4.5

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

@@ -106,11 +106,12 @@ export var createQuantileSeq = /* #__PURE__ */factory(name, dependencies, _ref =
             return probArr;
           }
         }
-        if (Array.isArray(probOrN)) {
+        if (isCollection(probOrN)) {
           // quantileSeq([a, b, c, d, ...], [prob1, prob2, ...][,sorted])
-          probArr = new Array(probOrN.length);
+          var probOrNArr = probOrN.valueOf();
+          probArr = new Array(probOrNArr.length);
           for (var _i2 = 0; _i2 < probArr.length; ++_i2) {
-            var currProb = probOrN[_i2];
+            var currProb = probOrNArr[_i2];
             if (isNumber(currProb)) {
               if (currProb < 0 || currProb > 1) {
                 throw new Error('Probability must be between 0 and 1, inclusive');

package/lib/esm/function/trigonometry/atan2.js

@@ -6,14 +6,15 @@ import { createMatAlgo11xS0s } from '../../type/matrix/utils/matAlgo11xS0s.js';
 import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'atan2';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix', 'concat'];
 export var createAtan2 = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
     matrix,
     equalScalar,
     BigNumber,
-    DenseMatrix
+    DenseMatrix,
+    concat
   } = _ref;
   var matAlgo02xDS0 = createMatAlgo02xDS0({
     typed,
@@ -36,7 +37,8 @@ export var createAtan2 = /* #__PURE__ */factory(name, dependencies, _ref => {
   });
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
 
   /**

package/lib/esm/function/unit/to.js

@@ -1,15 +1,17 @@
 import { factory } from '../../utils/factory.js';
 import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js';
 var name = 'to';
-var dependencies = ['typed', 'matrix'];
+var dependencies = ['typed', 'matrix', 'concat'];
 export var createTo = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
-    matrix
+    matrix,
+    concat
   } = _ref;
   var matrixAlgorithmSuite = createMatrixAlgorithmSuite({
     typed,
-    matrix
+    matrix,
+    concat
   });
 
   /**
@@ -23,9 +25,9 @@ export var createTo = /* #__PURE__ */factory(name, dependencies, _ref => {
    *
    * Examples:
    *
-   *    math.to(math.unit('2 inch'), 'cm')                   // returns Unit 5.08 cm
-   *    math.to(math.unit('2 inch'), math.unit(null, 'cm'))  // returns Unit 5.08 cm
-   *    math.to(math.unit(16, 'bytes'), 'bits')              // returns Unit 128 bits
+   *    math.to(math.unit('2 inch'), 'cm')             // returns Unit 5.08 cm
+   *    math.to(math.unit('2 inch'), math.unit('cm'))  // returns Unit 5.08 cm
+   *    math.to(math.unit(16, 'bytes'), 'bits')        // returns Unit 128 bits
    *
    * See also:
    *

package/lib/esm/type/matrix/utils/broadcast.js

@@ -0,0 +1,75 @@
+import { factory } from '../../../utils/factory.js';
+var name = 'broadcast';
+var dependancies = ['concat'];
+export var createBroadcast = /* #__PURE__ */factory(name, dependancies, _ref => {
+  var {
+    concat
+  } = _ref;
+  /**
+  * Broadcasts two matrices, and return both in an array
+  * It checks if it's possible with broadcasting rules
+  *
+  * @param {Matrix}   A      First Matrix
+  * @param {Matrix}   B      Second Matrix
+  *
+  * @return {Matrix[]}      [ broadcastedA, broadcastedB ]
+  */
+  return function (A, B) {
+    var N = Math.max(A._size.length, B._size.length); // max number of dims
+    if (A._size.length === B._size.length) {
+      if (A._size.every((dim, i) => dim === B._size[i])) {
+        // If matrices have the same size return them
+        return [A, B];
+      }
+    }
+    var sizeA = _padLeft(A._size, N, 0); // pad to the left to align dimensions to the right
+    var sizeB = _padLeft(B._size, N, 0); // pad to the left to align dimensions to the right
+
+    // calculate the max dimensions
+    var sizeMax = [];
+    for (var dim = 0; dim < N; dim++) {
+      sizeMax[dim] = Math.max(sizeA[dim], sizeB[dim]);
+    }
+
+    // check if the broadcasting rules applyes for both matrices
+    for (var _dim = 0; _dim < N; _dim++) {
+      _checkRules(sizeA, sizeMax, _dim);
+      _checkRules(sizeB, sizeMax, _dim);
+    }
+
+    // reshape A or B if needed to make them ready for concat
+    var AA = A.clone();
+    var BB = B.clone();
+    if (AA._size.length < N) {
+      AA.reshape(_padLeft(AA._size, N, 1));
+    } else if (BB._size.length < N) {
+      BB.reshape(_padLeft(BB._size, N, 1));
+    }
+
+    // stretches the matrices on each dimension to make them the same size
+    for (var _dim2 = 0; _dim2 < N; _dim2++) {
+      if (AA._size[_dim2] < sizeMax[_dim2]) {
+        AA = _stretch(AA, sizeMax[_dim2], _dim2);
+      }
+      if (BB._size[_dim2] < sizeMax[_dim2]) {
+        BB = _stretch(BB, sizeMax[_dim2], _dim2);
+      }
+    }
+
+    // return the array with the two broadcasted matrices
+    return [AA, BB];
+  };
+  function _padLeft(shape, N, filler) {
+    // pads an array of dimensions with numbers to the left, unitl the number of dimensions is N
+    return [...Array(N - shape.length).fill(filler), ...shape];
+  }
+  function _stretch(arrayToStretch, sizeToStretch, dimToStretch) {
+    // stretches a matrix up to a certain size in a certain dimension
+    return concat(...Array(sizeToStretch).fill(arrayToStretch), dimToStretch);
+  }
+  function _checkRules(shape, sizeMax, dim) {
+    if (shape[dim] < sizeMax[dim] & shape[dim] > 1) {
+      throw new Error("shape missmatch: missmatch is found in arg with shape (".concat(shape, ") not possible to broadcast dimension ").concat(dim, " with size ").concat(shape[dim], " to size ").concat(sizeMax[dim]));
+    }
+  }
+});

package/lib/esm/type/matrix/utils/matrixAlgorithmSuite.js

@@ -2,12 +2,14 @@ import { factory } from '../../../utils/factory.js';
 import { extend } from '../../../utils/object.js';
 import { createMatAlgo13xDD } from './matAlgo13xDD.js';
 import { createMatAlgo14xDs } from './matAlgo14xDs.js';
+import { createBroadcast } from './broadcast.js';
 var name = 'matrixAlgorithmSuite';
-var dependencies = ['typed', 'matrix'];
+var dependencies = ['typed', 'matrix', 'concat'];
 export var createMatrixAlgorithmSuite = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
-    matrix
+    matrix,
+    concat
   } = _ref;
   var matAlgo13xDD = createMatAlgo13xDD({
     typed
@@ -15,6 +17,9 @@ export var createMatrixAlgorithmSuite = /* #__PURE__ */factory(name, dependencie
   var matAlgo14xDs = createMatAlgo14xDs({
     typed
   });
+  var broadcast = createBroadcast({
+    concat
+  });
 
   /**
    * Return a signatures object with the usual boilerplate of
@@ -40,60 +45,60 @@ export var createMatrixAlgorithmSuite = /* #__PURE__ */factory(name, dependencie
     if (elop) {
       // First the dense ones
       matrixSignatures = {
-        'DenseMatrix, DenseMatrix': (x, y) => matAlgo13xDD(x, y, elop),
-        'Array, Array': (x, y) => matAlgo13xDD(matrix(x), matrix(y), elop).valueOf(),
-        'Array, DenseMatrix': (x, y) => matAlgo13xDD(matrix(x), y, elop),
-        'DenseMatrix, Array': (x, y) => matAlgo13xDD(x, matrix(y), elop)
+        'DenseMatrix, DenseMatrix': (x, y) => matAlgo13xDD(...broadcast(x, y), elop),
+        'Array, Array': (x, y) => matAlgo13xDD(...broadcast(matrix(x), matrix(y)), elop).valueOf(),
+        'Array, DenseMatrix': (x, y) => matAlgo13xDD(...broadcast(matrix(x), y), elop),
+        'DenseMatrix, Array': (x, y) => matAlgo13xDD(...broadcast(x, matrix(y)), elop)
       };
       // Now incorporate sparse matrices
       if (options.SS) {
-        matrixSignatures['SparseMatrix, SparseMatrix'] = (x, y) => options.SS(x, y, elop, false);
+        matrixSignatures['SparseMatrix, SparseMatrix'] = (x, y) => options.SS(...broadcast(x, y), elop, false);
       }
       if (options.DS) {
-        matrixSignatures['DenseMatrix, SparseMatrix'] = (x, y) => options.DS(x, y, elop, false);
-        matrixSignatures['Array, SparseMatrix'] = (x, y) => options.DS(matrix(x), y, elop, false);
+        matrixSignatures['DenseMatrix, SparseMatrix'] = (x, y) => options.DS(...broadcast(x, y), elop, false);
+        matrixSignatures['Array, SparseMatrix'] = (x, y) => options.DS(...broadcast(matrix(x), y), elop, false);
       }
       if (SD) {
-        matrixSignatures['SparseMatrix, DenseMatrix'] = (x, y) => SD(y, x, elop, true);
-        matrixSignatures['SparseMatrix, Array'] = (x, y) => SD(matrix(y), x, elop, true);
+        matrixSignatures['SparseMatrix, DenseMatrix'] = (x, y) => SD(...broadcast(y, x), elop, true);
+        matrixSignatures['SparseMatrix, Array'] = (x, y) => SD(...broadcast(matrix(y), x), elop, true);
       }
     } else {
       // No elop, use this
       // First the dense ones
       matrixSignatures = {
         'DenseMatrix, DenseMatrix': typed.referToSelf(self => (x, y) => {
-          return matAlgo13xDD(x, y, self);
+          return matAlgo13xDD(...broadcast(x, y), self);
         }),
         'Array, Array': typed.referToSelf(self => (x, y) => {
-          return matAlgo13xDD(matrix(x), matrix(y), self).valueOf();
+          return matAlgo13xDD(...broadcast(matrix(x), matrix(y)), self).valueOf();
         }),
         'Array, DenseMatrix': typed.referToSelf(self => (x, y) => {
-          return matAlgo13xDD(matrix(x), y, self);
+          return matAlgo13xDD(...broadcast(matrix(x), y), self);
         }),
         'DenseMatrix, Array': typed.referToSelf(self => (x, y) => {
-          return matAlgo13xDD(x, matrix(y), self);
+          return matAlgo13xDD(...broadcast(x, matrix(y)), self);
         })
       };
       // Now incorporate sparse matrices
       if (options.SS) {
         matrixSignatures['SparseMatrix, SparseMatrix'] = typed.referToSelf(self => (x, y) => {
-          return options.SS(x, y, self, false);
+          return options.SS(...broadcast(x, y), self, false);
         });
       }
       if (options.DS) {
         matrixSignatures['DenseMatrix, SparseMatrix'] = typed.referToSelf(self => (x, y) => {
-          return options.DS(x, y, self, false);
+          return options.DS(...broadcast(x, y), self, false);
         });
         matrixSignatures['Array, SparseMatrix'] = typed.referToSelf(self => (x, y) => {
-          return options.DS(matrix(x), y, self, false);
+          return options.DS(...broadcast(matrix(x), y), self, false);
         });
       }
       if (SD) {
         matrixSignatures['SparseMatrix, DenseMatrix'] = typed.referToSelf(self => (x, y) => {
-          return SD(y, x, self, true);
+          return SD(...broadcast(y, x), self, true);
         });
         matrixSignatures['SparseMatrix, Array'] = typed.referToSelf(self => (x, y) => {
-          return SD(matrix(y), x, self, true);
+          return SD(...broadcast(matrix(y), x), self, true);
         });
       }
     }

package/lib/esm/version.js

@@ -1,3 +1,3 @@
-export var version = '11.5.1';
+export var version = '11.6.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": "11.5.1",
+  "version": "11.6.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",
@@ -25,7 +25,7 @@
     "unit"
   ],
   "dependencies": {
-    "@babel/runtime": "^7.20.13",
+    "@babel/runtime": "^7.21.0",
     "complex.js": "^2.1.1",
     "decimal.js": "^10.4.3",
     "escape-latex": "^1.2.0",
@@ -36,24 +36,24 @@
     "typed-function": "^4.1.0"
   },
   "devDependencies": {
-    "@babel/core": "7.20.12",
+    "@babel/core": "7.21.0",
     "@babel/plugin-transform-object-assign": "7.18.6",
-    "@babel/plugin-transform-runtime": "7.19.6",
+    "@babel/plugin-transform-runtime": "7.21.0",
     "@babel/preset-env": "7.20.2",
-    "@babel/register": "7.18.9",
+    "@babel/register": "7.21.0",
     "@types/assert": "1.5.6",
     "@types/mocha": "10.0.1",
-    "@typescript-eslint/eslint-plugin": "5.50.0",
-    "@typescript-eslint/parser": "5.50.0",
+    "@typescript-eslint/eslint-plugin": "5.53.0",
+    "@typescript-eslint/parser": "5.53.0",
     "assert": "2.0.0",
     "babel-loader": "9.1.2",
     "benchmark": "2.1.4",
-    "c8": "7.12.0",
+    "c8": "7.13.0",
     "codecov": "3.8.3",
-    "core-js": "3.27.2",
+    "core-js": "3.28.0",
     "del": "6.1.1",
     "dtslint": "4.2.1",
-    "eslint": "8.33.0",
+    "eslint": "8.34.0",
     "eslint-config-prettier": "8.6.0",
     "eslint-config-standard": "17.0.0",
     "eslint-plugin-import": "2.27.5",
@@ -85,7 +85,7 @@
     "ndarray-pack": "1.2.1",
     "numericjs": "1.2.6",
     "pad-right": "0.2.2",
-    "prettier": "2.8.3",
+    "prettier": "2.8.4",
     "process": "0.11.10",
     "sylvester": "0.0.21",
     "ts-node": "10.9.1",