mathjs 10.5.3 → 10.6.0

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

@@ -0,0 +1,128 @@
+"use strict";
+
+var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createFft = void 0;
+
+var _toConsumableArray2 = _interopRequireDefault(require("@babel/runtime/helpers/toConsumableArray"));
+
+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 createFft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+  var typed = _ref.typed,
+      matrix = _ref.matrix,
+      addScalar = _ref.addScalar,
+      multiplyScalar = _ref.multiplyScalar,
+      divideScalar = _ref.divideScalar,
+      exp = _ref.exp,
+      tau = _ref.tau,
+      I = _ref.i;
+
+  /**
+   * Calculate N-dimensional fourier transform
+   *
+   * Syntax:
+   *
+   *     math.fft(arr)
+   *
+   * Examples:
+   *
+   *    math.fft([[1, 0], [1, 0]]) // returns [[{re:2, im:0}, {re:2, im:0}], [{re:0, im:0}, {re:0, im:0}]]
+   *
+   *
+   * See Also:
+   *
+   *      ifft
+   *
+   * @param {Array | Matrix} arr    An array or matrix
+   * @return {Array | Matrix}       N-dimensional fourier transformation of the array
+   */
+  return typed(name, {
+    Array: _ndFft,
+    Matrix: function Matrix(matrix) {
+      return matrix.create(_ndFft(matrix.toArray()));
+    }
+  });
+  /**
+   * Perform an N-dimensional Fourier transform
+   *
+   * @param {Array} arr      The array
+   * @return {Array}         resulting array
+   */
+
+  function _ndFft(arr) {
+    var size = (0, _array.arraySize)(arr);
+    if (size.length === 1) return _fft(arr, size[0]); // ndFft along dimension 1,...,N-1 then 1dFft along dimension 0
+
+    return _1dFft(arr.map(function (slice) {
+      return _ndFft(slice, size.slice(1));
+    }), 0);
+  }
+  /**
+   * Perform an 1-dimensional Fourier transform
+   *
+   * @param {Array} arr      The array
+   * @param {number} dim     dimension of the array to perform on
+   * @return {Array}         resulting array
+   */
+
+
+  function _1dFft(arr, dim) {
+    var size = (0, _array.arraySize)(arr);
+    if (dim !== 0) return new Array(size[0]).fill(0).map(function (_, i) {
+      return _1dFft(arr[i], dim - 1);
+    });
+    if (size.length === 1) return _fft(arr);
+
+    function _transpose(arr) {
+      // Swap first 2 dimensions
+      var size = (0, _array.arraySize)(arr);
+      return new Array(size[1]).fill(0).map(function (_, j) {
+        return new Array(size[0]).fill(0).map(function (_, i) {
+          return arr[i][j];
+        });
+      });
+    }
+
+    return _transpose(_1dFft(_transpose(arr), 1));
+  }
+  /**
+   * Perform an 1-dimensional Fourier transform
+   *
+   * @param {Array} arr      The array
+   * @return {Array}         resulting array
+   */
+
+
+  function _fft(arr) {
+    var len = arr.length;
+    if (len === 1) return [arr[0]];
+
+    if (len % 2 === 0) {
+      var ret = [].concat((0, _toConsumableArray2.default)(_fft(arr.filter(function (_, i) {
+        return i % 2 === 0;
+      }), len / 2)), (0, _toConsumableArray2.default)(_fft(arr.filter(function (_, i) {
+        return i % 2 === 1;
+      }), len / 2)));
+
+      for (var k = 0; k < len / 2; k++) {
+        var p = ret[k];
+        var q = multiplyScalar(ret[k + len / 2], exp(multiplyScalar(multiplyScalar(tau, I), divideScalar(-k, len))));
+        ret[k] = addScalar(p, q);
+        ret[k + len / 2] = addScalar(p, multiplyScalar(-1, q));
+      }
+
+      return ret;
+    }
+
+    throw new Error('Can only calculate FFT of power-of-two size');
+  }
+});
+exports.createFft = createFft;

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

@@ -0,0 +1,49 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createIfft = void 0;
+
+var _array = require("../../utils/array.js");
+
+var _factory = require("../../utils/factory.js");
+
+var _is = require("../../utils/is.js");
+
+var name = 'ifft';
+var dependencies = ['typed', 'fft', 'dotDivide', 'conj'];
+var createIfft = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+  var typed = _ref.typed,
+      fft = _ref.fft,
+      dotDivide = _ref.dotDivide,
+      conj = _ref.conj;
+
+  /**
+   * Calculate N-dimensional inverse fourier transform
+   *
+   * Syntax:
+   *
+   *     math.ifft(arr)
+   *
+   * Examples:
+   *
+   *    math.ifft([[2, 2], [0, 0]]) // returns [[{re:1, im:0}, {re:0, im:0}], [{re:1, im:0}, {re:0, im:0}]]
+   *
+   * See Also:
+   *
+   *      fft
+   *
+   * @param {Array | Matrix} arr    An array or matrix
+   * @return {Array | Matrix}       N-dimensional fourier transformation of the array
+   */
+  return typed(name, {
+    'Array | Matrix': function ArrayMatrix(arr) {
+      var size = (0, _is.isMatrix)(arr) ? arr.size() : (0, _array.arraySize)(arr);
+      return dotDivide(conj(fft(conj(arr))), size.reduce(function (acc, curr) {
+        return acc * curr;
+      }, 1));
+    }
+  });
+});
+exports.createIfft = createIfft;

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 10.5.3
- * @date    2022-05-11
+ * @version 10.6.0
+ * @date    2022-05-24
  *
  * @license
  * Copyright (C) 2013-2022 Jos de Jong <wjosdejong@gmail.com>

package/lib/cjs/version.js

@@ -4,7 +4,7 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.version = void 0;
-var version = '10.5.3'; // Note: This file is automatically generated when building math.js.
+var version = '10.6.0'; // Note: This file is automatically generated when building math.js.
 // Changes made in this file will be overwritten.
 
 exports.version = version;

package/lib/esm/entry/dependenciesAny/dependenciesFft.generated.js

@@ -0,0 +1,24 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { addScalarDependencies } from './dependenciesAddScalar.generated.js';
+import { divideScalarDependencies } from './dependenciesDivideScalar.generated.js';
+import { expDependencies } from './dependenciesExp.generated.js';
+import { iDependencies } from './dependenciesI.generated.js';
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { multiplyScalarDependencies } from './dependenciesMultiplyScalar.generated.js';
+import { tauDependencies } from './dependenciesTau.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createFft } from '../../factoriesAny.js';
+export var fftDependencies = {
+  addScalarDependencies,
+  divideScalarDependencies,
+  expDependencies,
+  iDependencies,
+  matrixDependencies,
+  multiplyScalarDependencies,
+  tauDependencies,
+  typedDependencies,
+  createFft
+};

package/lib/esm/entry/dependenciesAny/dependenciesIfft.generated.js

@@ -0,0 +1,16 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { conjDependencies } from './dependenciesConj.generated.js';
+import { dotDivideDependencies } from './dependenciesDotDivide.generated.js';
+import { fftDependencies } from './dependenciesFft.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createIfft } from '../../factoriesAny.js';
+export var ifftDependencies = {
+  conjDependencies,
+  dotDivideDependencies,
+  fftDependencies,
+  typedDependencies,
+  createIfft
+};

package/lib/esm/entry/dependenciesAny.generated.js

@@ -111,6 +111,7 @@ export { factorialDependencies } from './dependenciesAny/dependenciesFactorial.g
 export { falseDependencies } from './dependenciesAny/dependenciesFalse.generated.js';
 export { faradayDependencies } from './dependenciesAny/dependenciesFaraday.generated.js';
 export { fermiCouplingDependencies } from './dependenciesAny/dependenciesFermiCoupling.generated.js';
+export { fftDependencies } from './dependenciesAny/dependenciesFft.generated.js';
 export { FibonacciHeapDependencies } from './dependenciesAny/dependenciesFibonacciHeapClass.generated.js';
 export { filterDependencies } from './dependenciesAny/dependenciesFilter.generated.js';
 export { filterTransformDependencies } from './dependenciesAny/dependenciesFilterTransform.generated.js';
@@ -140,6 +141,7 @@ export { hexDependencies } from './dependenciesAny/dependenciesHex.generated.js'
 export { hypotDependencies } from './dependenciesAny/dependenciesHypot.generated.js';
 export { iDependencies } from './dependenciesAny/dependenciesI.generated.js';
 export { identityDependencies } from './dependenciesAny/dependenciesIdentity.generated.js';
+export { ifftDependencies } from './dependenciesAny/dependenciesIfft.generated.js';
 export { imDependencies } from './dependenciesAny/dependenciesIm.generated.js';
 export { ImmutableDenseMatrixDependencies } from './dependenciesAny/dependenciesImmutableDenseMatrixClass.generated.js';
 export { indexDependencies } from './dependenciesAny/dependenciesIndex.generated.js';

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

@@ -5,9 +5,9 @@ import _extends from "@babel/runtime/helpers/extends";
  * DON'T MAKE CHANGES HERE
  */
 import { config } from './configReadonly.js';
-import { createChainClass, createNode, createObjectNode, createOperatorNode, createParenthesisNode, createRelationalNode, createArrayNode, createBlockNode, createConditionalNode, createConstantNode, createRangeNode, createReviver, createChain, createFunctionAssignmentNode, createAccessorNode, createAssignmentNode, createIndexNode, createSymbolNode, createFunctionNode, createParse, createResolve, createSimplifyCore, createCompile, createHelpClass, createLeafCount, createEvaluate, createHelp, createParserClass, createSimplify, createParser, createRationalize, createDerivative, createSymbolicEqual, createFilterTransform, createForEachTransform, createMapTransform, createApplyTransform, createDiffTransform, createIndexTransform, createSubsetTransform, createConcatTransform, createMaxTransform, createMinTransform, createRangeTransform, createRowTransform, createSumTransform, createColumnTransform, createMeanTransform, createVarianceTransform, createStdTransform, createCumSumTransform } from '../factoriesAny.js';
+import { createChainClass, createNode, createObjectNode, createOperatorNode, createParenthesisNode, createRelationalNode, createArrayNode, createBlockNode, createConditionalNode, createConstantNode, createRangeNode, createReviver, createChain, createFunctionAssignmentNode, createAccessorNode, createAssignmentNode, createIndexNode, createSymbolNode, createFunctionNode, createParse, createResolve, createSimplifyCore, createCompile, createHelpClass, createLeafCount, createSimplify, createSymbolicEqual, createEvaluate, createHelp, createParserClass, createRationalize, createDerivative, createParser, createFilterTransform, createForEachTransform, createMapTransform, createApplyTransform, createDiffTransform, createSubsetTransform, createConcatTransform, createMaxTransform, createMinTransform, createRangeTransform, createSumTransform, createCumSumTransform, createRowTransform, createColumnTransform, createIndexTransform, createMeanTransform, createVarianceTransform, createStdTransform } from '../factoriesAny.js';
 import { BigNumber, Complex, e, _false, fineStructure, Fraction, i, _Infinity, LN10, LOG10E, Matrix, _NaN, _null, phi, Range, ResultSet, SQRT1_2, // eslint-disable-line camelcase
-sackurTetrode, tau, _true, version, DenseMatrix, efimovFactor, LN2, pi, replacer, SQRT2, typed, unaryPlus, weakMixingAngle, abs, acos, acot, acsc, addScalar, arg, asech, asinh, atan, atanh, bignumber, bitNot, boolean, clone, combinations, complex, conj, cosh, coth, csc, cube, equalScalar, erf, exp, expm1, filter, forEach, format, getMatrixDataType, hex, im, isInteger, isNegative, isPositive, isZero, LOG2E, lgamma, log10, log2, map, multiplyScalar, not, number, oct, pickRandom, print, random, re, sec, sign, sin, SparseMatrix, splitUnit, square, string, tan, typeOf, acosh, acsch, apply, asec, bin, combinationsWithRep, cos, csch, isNaN, isPrime, randomInt, sech, sinh, sparse, sqrt, tanh, unaryMinus, acoth, cot, fraction, isNumeric, matrix, matrixFromFunction, mod, nthRoot, numeric, or, prod, reshape, size, smaller, squeeze, subset, subtract, to, transpose, xgcd, zeros, and, bitAnd, bitXor, cbrt, compare, compareText, concat, count, ctranspose, diag, divideScalar, dotDivide, equal, flatten, gcd, hasNumericValue, hypot, ImmutableDenseMatrix, Index, kron, largerEq, leftShift, lsolve, matrixFromColumns, min, mode, nthRoots, ones, partitionSelect, resize, rightArithShift, round, smallerEq, unequal, usolve, xor, add, atan2, bitOr, catalan, compareNatural, cumsum, deepEqual, diff, dot, equalText, floor, identity, invmod, larger, log, lsolveAll, matrixFromRows, multiply, qr, range, rightLogShift, row, setCartesian, setDistinct, setIsSubset, setPowerset, slu, sum, trace, usolveAll, asin, ceil, column, composition, cross, det, distance, dotMultiply, FibonacciHeap, fix, index, inv, lcm, log1p, max, pinv, pow, quantileSeq, setDifference, setMultiplicity, setSymDifference, sort, sqrtm, Unit, vacuumImpedance, wienDisplacement, atomicMass, bohrMagneton, boltzmann, conductanceQuantum, createUnit, deuteronMass, dotPow, electricConstant, elementaryCharge, expm, faraday, firstRadiation, gamma, gravitationConstant, hartreeEnergy, intersect, klitzing, loschmidt, magneticConstant, molarMass, molarPlanckConstant, neutronMass, nuclearMagneton, planckCharge, planckLength, planckTemperature, protonMass, reducedPlanckConstant, rydberg, setIntersect, setUnion, Spa, stefanBoltzmann, thomsonCrossSection, avogadro, bohrRadius, coulomb, divide, electronMass, factorial, gravity, inverseConductanceQuantum, lup, magneticFluxQuantum, molarMassC12, multinomial, permutations, planckMass, quantumOfCirculation, secondRadiation, stirlingS2, unit, bellNumbers, eigs, fermiCoupling, mean, molarVolume, planckConstant, setSize, speedOfLight, variance, classicalElectronRadius, lusolve, median, gasConstant, mad, std, norm, rotationMatrix, kldivergence, rotate, planckTime } from './pureFunctionsAny.generated.js';
+sackurTetrode, tau, _true, version, DenseMatrix, efimovFactor, LN2, pi, replacer, SQRT2, typed, unaryPlus, weakMixingAngle, abs, acos, acot, acsc, addScalar, arg, asech, asinh, atan, atanh, bignumber, bitNot, boolean, clone, combinations, complex, conj, cosh, coth, csc, cube, equalScalar, erf, exp, expm1, filter, forEach, format, getMatrixDataType, hex, im, isInteger, isNegative, isPositive, isZero, LOG2E, lgamma, log10, log2, map, multiplyScalar, not, number, oct, pickRandom, print, random, re, sec, sign, sin, SparseMatrix, splitUnit, square, string, tan, typeOf, acosh, acsch, apply, asec, bin, combinationsWithRep, cos, csch, isNaN, isPrime, randomInt, sech, sinh, sparse, sqrt, tanh, unaryMinus, acoth, cot, fraction, isNumeric, matrix, matrixFromFunction, mod, nthRoot, numeric, or, prod, reshape, size, smaller, squeeze, subset, subtract, to, transpose, xgcd, zeros, and, bitAnd, bitXor, cbrt, compare, compareText, concat, count, ctranspose, diag, divideScalar, dotDivide, equal, fft, flatten, gcd, hasNumericValue, hypot, ifft, kron, largerEq, leftShift, lsolve, matrixFromColumns, min, mode, nthRoots, ones, partitionSelect, resize, rightArithShift, round, smallerEq, unequal, usolve, xor, add, atan2, bitOr, catalan, compareNatural, cumsum, deepEqual, diff, dot, equalText, floor, identity, invmod, larger, log, lsolveAll, matrixFromRows, multiply, qr, range, rightLogShift, setSize, slu, sum, trace, usolveAll, asin, ceil, composition, cross, det, distance, dotMultiply, FibonacciHeap, fix, ImmutableDenseMatrix, Index, intersect, lcm, log1p, max, quantileSeq, row, setCartesian, setDistinct, setIsSubset, setPowerset, sort, column, index, inv, pinv, pow, setDifference, setMultiplicity, Spa, sqrtm, Unit, vacuumImpedance, wienDisplacement, atomicMass, bohrMagneton, boltzmann, conductanceQuantum, createUnit, deuteronMass, dotPow, electricConstant, elementaryCharge, expm, faraday, firstRadiation, gamma, gravitationConstant, hartreeEnergy, klitzing, loschmidt, magneticConstant, molarMass, molarPlanckConstant, neutronMass, nuclearMagneton, planckCharge, planckLength, planckTemperature, protonMass, reducedPlanckConstant, rydberg, setIntersect, speedOfLight, stefanBoltzmann, thomsonCrossSection, avogadro, bohrRadius, coulomb, divide, electronMass, factorial, gravity, inverseConductanceQuantum, lup, magneticFluxQuantum, molarMassC12, multinomial, permutations, planckMass, quantumOfCirculation, secondRadiation, stirlingS2, unit, bellNumbers, eigs, fermiCoupling, mean, molarVolume, planckConstant, setSymDifference, classicalElectronRadius, lusolve, median, setUnion, variance, kldivergence, norm, planckTime, rotationMatrix, gasConstant, std, mad, rotate } from './pureFunctionsAny.generated.js';
 var math = {}; // NOT pure!
 
 var mathWithTransform = {}; // NOT pure!
@@ -138,18 +138,6 @@ export var leafCount = createLeafCount({
   parse,
   typed
 });
-export var evaluate = createEvaluate({
-  parse,
-  typed
-});
-export var help = createHelp({
-  Help,
-  mathWithTransform,
-  typed
-});
-export var Parser = createParserClass({
-  evaluate
-});
 export var simplify = createSimplify({
   bignumber,
   fraction,
@@ -177,10 +165,24 @@ export var simplify = createSimplify({
   subtract,
   typed
 });
-export var parser = createParser({
-  Parser,
+export var symbolicEqual = createSymbolicEqual({
+  OperatorNode,
+  parse,
+  simplify,
+  typed
+});
+export var evaluate = createEvaluate({
+  parse,
   typed
 });
+export var help = createHelp({
+  Help,
+  mathWithTransform,
+  typed
+});
+export var Parser = createParserClass({
+  evaluate
+});
 export var rationalize = createRationalize({
   bignumber,
   fraction,
@@ -222,10 +224,8 @@ export var derivative = createDerivative({
   simplify,
   typed
 });
-export var symbolicEqual = createSymbolicEqual({
-  OperatorNode,
-  parse,
-  simplify,
+export var parser = createParser({
+  Parser,
   typed
 });
 
@@ -364,10 +364,12 @@ _extends(math, {
   divideScalar,
   dotDivide,
   equal,
+  fft,
   flatten,
   gcd,
   hasNumericValue,
   hypot,
+  ifft,
   kron,
   largerEq,
   leftShift,
@@ -406,36 +408,37 @@ _extends(math, {
   qr,
   range,
   rightLogShift,
-  row,
-  setCartesian,
-  setDistinct,
-  setIsSubset,
-  setPowerset,
+  setSize,
   slu,
   sum,
   trace,
   usolveAll,
   asin,
   ceil,
-  column,
   composition,
   cross,
   det,
   distance,
   dotMultiply,
   fix,
-  index,
-  inv,
+  intersect,
   lcm,
   log1p,
   max,
+  quantileSeq,
+  row,
+  setCartesian,
+  setDistinct,
+  setIsSubset,
+  setPowerset,
+  sort,
+  column,
+  index,
+  inv,
   pinv,
   pow,
-  quantileSeq,
   setDifference,
   setMultiplicity,
-  setSymDifference,
-  sort,
   sqrtm,
   vacuumImpedance,
   wienDisplacement,
@@ -454,7 +457,6 @@ _extends(math, {
   gamma,
   gravitationConstant,
   hartreeEnergy,
-  intersect,
   klitzing,
   loschmidt,
   magneticConstant,
@@ -469,7 +471,7 @@ _extends(math, {
   reducedPlanckConstant,
   rydberg,
   setIntersect,
-  setUnion,
+  speedOfLight,
   stefanBoltzmann,
   thomsonCrossSection,
   avogadro,
@@ -501,27 +503,27 @@ _extends(math, {
   mean,
   molarVolume,
   planckConstant,
-  setSize,
-  speedOfLight,
-  variance,
+  setSymDifference,
+  simplify,
+  symbolicEqual,
   classicalElectronRadius,
   evaluate,
   help,
   lusolve,
   median,
-  simplify,
-  gasConstant,
-  mad,
-  parser,
   rationalize,
-  std,
+  setUnion,
+  variance,
   derivative,
+  kldivergence,
   norm,
+  planckTime,
   rotationMatrix,
-  kldivergence,
+  gasConstant,
+  parser,
+  std,
+  mad,
   rotate,
-  planckTime,
-  symbolicEqual,
   config
 });
 
@@ -546,9 +548,6 @@ _extends(mathWithTransform, math, {
     subtract,
     typed
   }),
-  index: createIndexTransform({
-    Index
-  }),
   subset: createSubsetTransform({
     matrix,
     typed
@@ -580,24 +579,32 @@ _extends(mathWithTransform, math, {
     smallerEq,
     typed
   }),
-  row: createRowTransform({
-    Index,
-    matrix,
-    range,
-    typed
-  }),
   sum: createSumTransform({
     add,
     config,
     numeric,
     typed
   }),
+  cumsum: createCumSumTransform({
+    add,
+    typed,
+    unaryPlus
+  }),
+  row: createRowTransform({
+    Index,
+    matrix,
+    range,
+    typed
+  }),
   column: createColumnTransform({
     Index,
     matrix,
     range,
     typed
   }),
+  index: createIndexTransform({
+    Index
+  }),
   mean: createMeanTransform({
     add,
     divide,
@@ -616,11 +623,6 @@ _extends(mathWithTransform, math, {
     sqrt,
     typed,
     variance
-  }),
-  cumsum: createCumSumTransform({
-    add,
-    typed,
-    unaryPlus
   })
 });
 
@@ -647,12 +649,12 @@ _extends(classes, {
   SparseMatrix,
   AccessorNode,
   AssignmentNode,
-  ImmutableDenseMatrix,
-  Index,
   IndexNode,
   FibonacciHeap,
-  Unit,
+  ImmutableDenseMatrix,
+  Index,
   Spa,
+  Unit,
   SymbolNode,
   FunctionNode,
   Help,