mathjs 11.5.1 → 11.6.0

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

@@ -9,7 +9,7 @@ var _toConsumableArray2 = _interopRequireDefault(require("@babel/runtime/helpers
 var _array = require("../../utils/array.js");
 var _factory = require("../../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'];
 var createFft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     matrix = _ref.matrix,
@@ -18,7 +18,12 @@ var createFft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
     divideScalar = _ref.divideScalar,
     exp = _ref.exp,
     tau = _ref.tau,
-    I = _ref.i;
+    I = _ref.i,
+    dotDivide = _ref.dotDivide,
+    conj = _ref.conj,
+    pow = _ref.pow,
+    ceil = _ref.ceil,
+    log2 = _ref.log2;
   /**
    * Calculate N-dimensional fourier transform
    *
@@ -84,7 +89,38 @@ var createFft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
     }
     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 = [].concat((0, _toConsumableArray2["default"])(new Array(n).fill(0).map(function (_, i) {
+      return multiplyScalar(arr[i], chirp[n - 1 + i]);
+    })), (0, _toConsumableArray2["default"])(new Array(N2 - n).fill(0)));
+    var ichirp = [].concat((0, _toConsumableArray2["default"])(new Array(n + n - 1).fill(0).map(function (_, i) {
+      return divideScalar(1, chirp[i]);
+    })), (0, _toConsumableArray2["default"])(new Array(N2 - (n + n - 1)).fill(0)));
+    var fftXp = _fft(xp);
+    var fftIchirp = _fft(ichirp);
+    var fftProduct = new Array(N2).fill(0).map(function (_, i) {
+      return 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
    *
@@ -107,8 +143,11 @@ var createFft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
         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')
   }
 });
 exports.createFft = createFft;

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

@@ -25,7 +25,7 @@ var createGetMatrixDataType = /* #__PURE__ */(0, _factory.factory)(name, depende
    *
    *    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/cjs/function/matrix/row.js

@@ -5,6 +5,7 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.createRow = void 0;
 var _factory = require("../../utils/factory.js");
+var _is = require("../../utils/is.js");
 var _object = require("../../utils/object.js");
 var _array = require("../../utils/array.js");
 var name = 'row';
@@ -56,7 +57,8 @@ var createRow = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
     (0, _array.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 (0, _is.isMatrix)(result) ? result : matrix([[result]]);
   }
 });
 exports.createRow = createRow;

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

@@ -31,7 +31,7 @@ var createSubset = /* #__PURE__ */(0, _factory.factory)(name, dependencies, func
    *
    *     // 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/cjs/function/relational/compare.js

@@ -13,7 +13,7 @@ var _matAlgo05xSfSf = require("../../type/matrix/utils/matAlgo05xSfSf.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _compareUnits = require("./compareUnits.js");
 var name = 'compare';
-var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'BigNumber', 'Fraction', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'BigNumber', 'Fraction', 'DenseMatrix', 'concat'];
 var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
@@ -21,7 +21,8 @@ var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
     matrix = _ref.matrix,
     BigNumber = _ref.BigNumber,
     Fraction = _ref.Fraction,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -35,7 +36,8 @@ var createCompare = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
   var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed

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

@@ -8,14 +8,16 @@ var _string = require("../../utils/string.js");
 var _factory = require("../../utils/factory.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var name = 'compareText';
-var dependencies = ['typed', 'matrix'];
+var dependencies = ['typed', 'matrix', 'concat'];
 _string.compareText.signature = 'any, any';
 var createCompareText = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
-    matrix = _ref.matrix;
+    matrix = _ref.matrix,
+    concat = _ref.concat;
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
 
   /**

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

@@ -10,12 +10,13 @@ var _matAlgo07xSSf = require("../../type/matrix/utils/matAlgo07xSSf.js");
 var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var name = 'equal';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'equalScalar', 'DenseMatrix', 'concat'];
 var createEqual = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     matrix = _ref.matrix,
     equalScalar = _ref.equalScalar,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -29,7 +30,8 @@ var createEqual = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
 
   /**

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

@@ -13,12 +13,13 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _compareUnits = require("./compareUnits.js");
 var name = 'larger';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 var createLarger = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
     matrix = _ref.matrix,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -32,7 +33,8 @@ var createLarger = /* #__PURE__ */(0, _factory.factory)(name, dependencies, func
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
   var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed

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

@@ -13,12 +13,13 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _compareUnits = require("./compareUnits.js");
 var name = 'largerEq';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 var createLargerEq = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
     matrix = _ref.matrix,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -32,7 +33,8 @@ var createLargerEq = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fu
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
   var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed

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

@@ -13,12 +13,13 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _compareUnits = require("./compareUnits.js");
 var name = 'smaller';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 var createSmaller = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
     matrix = _ref.matrix,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -32,7 +33,8 @@ var createSmaller = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
   var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed

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

@@ -13,12 +13,13 @@ var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var _compareUnits = require("./compareUnits.js");
 var name = 'smallerEq';
-var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'DenseMatrix', 'concat'];
 var createSmallerEq = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
     matrix = _ref.matrix,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -32,7 +33,8 @@ var createSmallerEq = /* #__PURE__ */(0, _factory.factory)(name, dependencies, f
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
   var compareUnits = (0, _compareUnits.createCompareUnits)({
     typed: typed

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

@@ -10,13 +10,14 @@ var _matAlgo07xSSf = require("../../type/matrix/utils/matAlgo07xSSf.js");
 var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var name = 'unequal';
-var dependencies = ['typed', 'config', 'equalScalar', 'matrix', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'equalScalar', 'matrix', 'DenseMatrix', 'concat'];
 var createUnequal = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     config = _ref.config,
     equalScalar = _ref.equalScalar,
     matrix = _ref.matrix,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
     typed: typed
   });
@@ -30,7 +31,8 @@ var createUnequal = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
 
   /**

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

@@ -33,7 +33,7 @@ var createMax = /* #__PURE__ */(0, _factory.factory)(name, dependencies, functio
    *
    *     // 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/cjs/function/statistics/quantileSeq.js

@@ -110,11 +110,12 @@ var createQuantileSeq = /* #__PURE__ */(0, _factory.factory)(name, dependencies,
             return probArr;
           }
         }
-        if (Array.isArray(probOrN)) {
+        if ((0, _is.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 ((0, _is.isNumber)(currProb)) {
               if (currProb < 0 || currProb > 1) {
                 throw new Error('Probability must be between 0 and 1, inclusive');

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

@@ -12,13 +12,14 @@ var _matAlgo11xS0s = require("../../type/matrix/utils/matAlgo11xS0s.js");
 var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var name = 'atan2';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix'];
+var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix', 'concat'];
 var createAtan2 = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
     matrix = _ref.matrix,
     equalScalar = _ref.equalScalar,
     BigNumber = _ref.BigNumber,
-    DenseMatrix = _ref.DenseMatrix;
+    DenseMatrix = _ref.DenseMatrix,
+    concat = _ref.concat;
   var matAlgo02xDS0 = (0, _matAlgo02xDS.createMatAlgo02xDS0)({
     typed: typed,
     equalScalar: equalScalar
@@ -40,7 +41,8 @@ var createAtan2 = /* #__PURE__ */(0, _factory.factory)(name, dependencies, funct
   });
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
 
   /**

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

@@ -7,13 +7,15 @@ exports.createTo = void 0;
 var _factory = require("../../utils/factory.js");
 var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
 var name = 'to';
-var dependencies = ['typed', 'matrix'];
+var dependencies = ['typed', 'matrix', 'concat'];
 var createTo = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
-    matrix = _ref.matrix;
+    matrix = _ref.matrix,
+    concat = _ref.concat;
   var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
     typed: typed,
-    matrix: matrix
+    matrix: matrix,
+    concat: concat
   });
 
   /**
@@ -27,9 +29,9 @@ var createTo = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function
    *
    * 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/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 11.5.1
- * @date    2023-01-31
+ * @version 11.6.0
+ * @date    2023-02-24
  *
  * @license
  * Copyright (C) 2013-2023 Jos de Jong <wjosdejong@gmail.com>

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

@@ -0,0 +1,84 @@
+"use strict";
+
+var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createBroadcast = void 0;
+var _toConsumableArray2 = _interopRequireDefault(require("@babel/runtime/helpers/toConsumableArray"));
+var _factory = require("../../../utils/factory.js");
+var name = 'broadcast';
+var dependancies = ['concat'];
+var createBroadcast = /* #__PURE__ */(0, _factory.factory)(name, dependancies, function (_ref) {
+  var concat = _ref.concat;
+  /**
+  * 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(function (dim, i) {
+        return 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 [].concat((0, _toConsumableArray2["default"])(Array(N - shape.length).fill(filler)), (0, _toConsumableArray2["default"])(shape));
+  }
+  function _stretch(arrayToStretch, sizeToStretch, dimToStretch) {
+    // stretches a matrix up to a certain size in a certain dimension
+    return concat.apply(void 0, (0, _toConsumableArray2["default"])(Array(sizeToStretch).fill(arrayToStretch)).concat([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]));
+    }
+  }
+});
+exports.createBroadcast = createBroadcast;