mathjs 11.3.3 → 11.4.0

package/lib/cjs/function/algebra/sylvester.js

@@ -0,0 +1,127 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createSylvester = void 0;
+var _factory = require("../../utils/factory.js");
+var name = 'sylvester';
+var dependencies = ['typed', 'schur', 'matrixFromColumns', 'matrix', 'multiply', 'range', 'concat', 'transpose', 'index', 'subset', 'add', 'subtract', 'identity', 'lusolve', 'abs'];
+var createSylvester = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
+  var typed = _ref.typed,
+    schur = _ref.schur,
+    matrixFromColumns = _ref.matrixFromColumns,
+    matrix = _ref.matrix,
+    multiply = _ref.multiply,
+    range = _ref.range,
+    concat = _ref.concat,
+    transpose = _ref.transpose,
+    index = _ref.index,
+    subset = _ref.subset,
+    add = _ref.add,
+    subtract = _ref.subtract,
+    identity = _ref.identity,
+    lusolve = _ref.lusolve,
+    abs = _ref.abs;
+  /**
+   *
+   * Solves the real-valued Sylvester equation AX+XB=C for X, where A, B and C are
+   * matrices of appropriate dimensions, being A and B squared. Notice that other
+   * equivalent definitions for the Sylvester equation exist and this function
+   * assumes the one presented in the original publication of the the Bartels-
+   * Stewart algorithm, which is implemented by this function.
+   * https://en.wikipedia.org/wiki/Sylvester_equation
+   *
+   * Syntax:
+   *
+   *     math.sylvester(A, B, C)
+   *
+   * Examples:
+   *
+   *     const A = [[-1, -2], [1, 1]]
+   *     const B = [[2, -1], [1, -2]]
+   *     const C = [[-3, 2], [3, 0]]
+   *     math.sylvester(A, B, C)      // returns DenseMatrix [[-0.25, 0.25], [1.5, -1.25]]
+   *
+   * See also:
+   *
+   *     schur, lyap
+   *
+   * @param {Matrix | Array} A  Matrix A
+   * @param {Matrix | Array} B  Matrix B
+   * @param {Matrix | Array} C  Matrix C
+   * @return {Matrix | Array}   Matrix X, solving the Sylvester equation
+   */
+  return typed(name, {
+    'Matrix, Matrix, Matrix': _sylvester,
+    'Array, Matrix, Matrix': function ArrayMatrixMatrix(A, B, C) {
+      return _sylvester(matrix(A), B, C);
+    },
+    'Array, Array, Matrix': function ArrayArrayMatrix(A, B, C) {
+      return _sylvester(matrix(A), matrix(B), C);
+    },
+    'Array, Matrix, Array': function ArrayMatrixArray(A, B, C) {
+      return _sylvester(matrix(A), B, matrix(C));
+    },
+    'Matrix, Array, Matrix': function MatrixArrayMatrix(A, B, C) {
+      return _sylvester(A, matrix(B), C);
+    },
+    'Matrix, Array, Array': function MatrixArrayArray(A, B, C) {
+      return _sylvester(A, matrix(B), matrix(C));
+    },
+    'Matrix, Matrix, Array': function MatrixMatrixArray(A, B, C) {
+      return _sylvester(A, B, matrix(C));
+    },
+    'Array, Array, Array': function ArrayArrayArray(A, B, C) {
+      return _sylvester(matrix(A), matrix(B), matrix(C)).toArray();
+    }
+  });
+  function _sylvester(A, B, C) {
+    var n = B.size()[0];
+    var m = A.size()[0];
+    var sA = schur(A);
+    var F = sA.T;
+    var U = sA.U;
+    var sB = schur(multiply(-1, B));
+    var G = sB.T;
+    var V = sB.U;
+    var D = multiply(multiply(transpose(U), C), V);
+    var all = range(0, m);
+    var y = [];
+    var hc = function hc(a, b) {
+      return concat(a, b, 1);
+    };
+    var vc = function vc(a, b) {
+      return concat(a, b, 0);
+    };
+    for (var k = 0; k < n; k++) {
+      if (k < n - 1 && abs(subset(G, index(k + 1, k))) > 1e-5) {
+        var RHS = vc(subset(D, index(all, k)), subset(D, index(all, k + 1)));
+        for (var j = 0; j < k; j++) {
+          RHS = add(RHS, vc(multiply(y[j], subset(G, index(j, k))), multiply(y[j], subset(G, index(j, k + 1)))));
+        }
+        var gkk = multiply(identity(m), multiply(-1, subset(G, index(k, k))));
+        var gmk = multiply(identity(m), multiply(-1, subset(G, index(k + 1, k))));
+        var gkm = multiply(identity(m), multiply(-1, subset(G, index(k, k + 1))));
+        var gmm = multiply(identity(m), multiply(-1, subset(G, index(k + 1, k + 1))));
+        var LHS = vc(hc(add(F, gkk), gmk), hc(gkm, add(F, gmm)));
+        var yAux = lusolve(LHS, RHS);
+        y[k] = yAux.subset(index(range(0, m), 0));
+        y[k + 1] = yAux.subset(index(range(m, 2 * m), 0));
+        k++;
+      } else {
+        var _RHS = subset(D, index(all, k));
+        for (var _j = 0; _j < k; _j++) {
+          _RHS = add(_RHS, multiply(y[_j], subset(G, index(_j, k))));
+        }
+        var _gkk = subset(G, index(k, k));
+        var _LHS = subtract(F, multiply(_gkk, identity(m)));
+        y[k] = lusolve(_LHS, _RHS);
+      }
+    }
+    var Y = matrix(matrixFromColumns.apply(void 0, y));
+    var X = multiply(U, multiply(Y, transpose(V)));
+    return X;
+  }
+});
+exports.createSylvester = createSylvester;

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

@@ -37,7 +37,7 @@ var createForEach = /* #__PURE__ */(0, _factory.factory)(name, dependencies, fun
   return typed(name, {
     'Array, function': _forEach,
     'Matrix, function': function MatrixFunction(x, callback) {
-      return x.forEach(callback);
+      x.forEach(callback);
     }
   });
 });

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.3.3
- * @date    2022-11-07
+ * @version 11.4.0
+ * @date    2022-11-18
  *
  * @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 = '11.3.3';
+var version = '11.4.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/dependenciesLyap.generated.js

@@ -0,0 +1,18 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { multiplyDependencies } from './dependenciesMultiply.generated.js';
+import { sylvesterDependencies } from './dependenciesSylvester.generated.js';
+import { transposeDependencies } from './dependenciesTranspose.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createLyap } from '../../factoriesAny.js';
+export var lyapDependencies = {
+  matrixDependencies,
+  multiplyDependencies,
+  sylvesterDependencies,
+  transposeDependencies,
+  typedDependencies,
+  createLyap
+};

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

@@ -0,0 +1,34 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { addDependencies } from './dependenciesAdd.generated.js';
+import { cbrtDependencies } from './dependenciesCbrt.generated.js';
+import { divideDependencies } from './dependenciesDivide.generated.js';
+import { equalScalarDependencies } from './dependenciesEqualScalar.generated.js';
+import { imDependencies } from './dependenciesIm.generated.js';
+import { isZeroDependencies } from './dependenciesIsZero.generated.js';
+import { multiplyDependencies } from './dependenciesMultiply.generated.js';
+import { reDependencies } from './dependenciesRe.generated.js';
+import { sqrtDependencies } from './dependenciesSqrt.generated.js';
+import { subtractDependencies } from './dependenciesSubtract.generated.js';
+import { typeOfDependencies } from './dependenciesTypeOf.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { unaryMinusDependencies } from './dependenciesUnaryMinus.generated.js';
+import { createPolynomialRoot } from '../../factoriesAny.js';
+export var polynomialRootDependencies = {
+  addDependencies,
+  cbrtDependencies,
+  divideDependencies,
+  equalScalarDependencies,
+  imDependencies,
+  isZeroDependencies,
+  multiplyDependencies,
+  reDependencies,
+  sqrtDependencies,
+  subtractDependencies,
+  typeOfDependencies,
+  typedDependencies,
+  unaryMinusDependencies,
+  createPolynomialRoot
+};

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

@@ -0,0 +1,22 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { identityDependencies } from './dependenciesIdentity.generated.js';
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { multiplyDependencies } from './dependenciesMultiply.generated.js';
+import { normDependencies } from './dependenciesNorm.generated.js';
+import { qrDependencies } from './dependenciesQr.generated.js';
+import { subtractDependencies } from './dependenciesSubtract.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createSchur } from '../../factoriesAny.js';
+export var schurDependencies = {
+  identityDependencies,
+  matrixDependencies,
+  multiplyDependencies,
+  normDependencies,
+  qrDependencies,
+  subtractDependencies,
+  typedDependencies,
+  createSchur
+};

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

@@ -0,0 +1,38 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { absDependencies } from './dependenciesAbs.generated.js';
+import { addDependencies } from './dependenciesAdd.generated.js';
+import { concatDependencies } from './dependenciesConcat.generated.js';
+import { identityDependencies } from './dependenciesIdentity.generated.js';
+import { indexDependencies } from './dependenciesIndex.generated.js';
+import { lusolveDependencies } from './dependenciesLusolve.generated.js';
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { matrixFromColumnsDependencies } from './dependenciesMatrixFromColumns.generated.js';
+import { multiplyDependencies } from './dependenciesMultiply.generated.js';
+import { rangeDependencies } from './dependenciesRange.generated.js';
+import { schurDependencies } from './dependenciesSchur.generated.js';
+import { subsetDependencies } from './dependenciesSubset.generated.js';
+import { subtractDependencies } from './dependenciesSubtract.generated.js';
+import { transposeDependencies } from './dependenciesTranspose.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createSylvester } from '../../factoriesAny.js';
+export var sylvesterDependencies = {
+  absDependencies,
+  addDependencies,
+  concatDependencies,
+  identityDependencies,
+  indexDependencies,
+  lusolveDependencies,
+  matrixDependencies,
+  matrixFromColumnsDependencies,
+  multiplyDependencies,
+  rangeDependencies,
+  schurDependencies,
+  subsetDependencies,
+  subtractDependencies,
+  transposeDependencies,
+  typedDependencies,
+  createSylvester
+};

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

@@ -182,6 +182,7 @@ export { lsolveDependencies } from './dependenciesAny/dependenciesLsolve.generat
 export { lsolveAllDependencies } from './dependenciesAny/dependenciesLsolveAll.generated.js';
 export { lupDependencies } from './dependenciesAny/dependenciesLup.generated.js';
 export { lusolveDependencies } from './dependenciesAny/dependenciesLusolve.generated.js';
+export { lyapDependencies } from './dependenciesAny/dependenciesLyap.generated.js';
 export { madDependencies } from './dependenciesAny/dependenciesMad.generated.js';
 export { magneticConstantDependencies } from './dependenciesAny/dependenciesMagneticConstant.generated.js';
 export { magneticFluxQuantumDependencies } from './dependenciesAny/dependenciesMagneticFluxQuantum.generated.js';
@@ -240,6 +241,7 @@ export { planckLengthDependencies } from './dependenciesAny/dependenciesPlanckLe
 export { planckMassDependencies } from './dependenciesAny/dependenciesPlanckMass.generated.js';
 export { planckTemperatureDependencies } from './dependenciesAny/dependenciesPlanckTemperature.generated.js';
 export { planckTimeDependencies } from './dependenciesAny/dependenciesPlanckTime.generated.js';
+export { polynomialRootDependencies } from './dependenciesAny/dependenciesPolynomialRoot.generated.js';
 export { powDependencies } from './dependenciesAny/dependenciesPow.generated.js';
 export { printDependencies } from './dependenciesAny/dependenciesPrint.generated.js';
 export { prodDependencies } from './dependenciesAny/dependenciesProd.generated.js';
@@ -274,6 +276,7 @@ export { rydbergDependencies } from './dependenciesAny/dependenciesRydberg.gener
 export { SQRT1_2Dependencies } from './dependenciesAny/dependenciesSQRT1_2.generated.js'; // eslint-disable-line camelcase
 export { SQRT2Dependencies } from './dependenciesAny/dependenciesSQRT2.generated.js';
 export { sackurTetrodeDependencies } from './dependenciesAny/dependenciesSackurTetrode.generated.js';
+export { schurDependencies } from './dependenciesAny/dependenciesSchur.generated.js';
 export { secDependencies } from './dependenciesAny/dependenciesSec.generated.js';
 export { sechDependencies } from './dependenciesAny/dependenciesSech.generated.js';
 export { secondRadiationDependencies } from './dependenciesAny/dependenciesSecondRadiation.generated.js';
@@ -317,6 +320,7 @@ export { subsetTransformDependencies } from './dependenciesAny/dependenciesSubse
 export { subtractDependencies } from './dependenciesAny/dependenciesSubtract.generated.js';
 export { sumDependencies } from './dependenciesAny/dependenciesSum.generated.js';
 export { sumTransformDependencies } from './dependenciesAny/dependenciesSumTransform.generated.js';
+export { sylvesterDependencies } from './dependenciesAny/dependenciesSylvester.generated.js';
 export { SymbolNodeDependencies } from './dependenciesAny/dependenciesSymbolNode.generated.js';
 export { symbolicEqualDependencies } from './dependenciesAny/dependenciesSymbolicEqual.generated.js';
 export { tanDependencies } from './dependenciesAny/dependenciesTan.generated.js';

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

@@ -4,10 +4,10 @@ import _extends from "@babel/runtime/helpers/extends";
  * DON'T MAKE CHANGES HERE
  */
 import { config } from './configReadonly.js';
-import { createNode, createObjectNode, createOperatorNode, createParenthesisNode, createRelationalNode, createArrayNode, createBlockNode, createConditionalNode, createConstantNode, createRangeNode, createReviver, createChainClass, createFunctionAssignmentNode, createChain, createAccessorNode, createAssignmentNode, createIndexNode, createSymbolNode, createFunctionNode, createParse, createResolve, createSimplifyConstant, createCompile, createHelpClass, createLeafCount, createSimplifyCore, createEvaluate, createHelp, createParserClass, createParser, createSimplify, createSymbolicEqual, createDerivative, createRationalize, createFilterTransform, createForEachTransform, createMapTransform, createApplyTransform, createDiffTransform, createSubsetTransform, createConcatTransform, createMaxTransform, createMinTransform, createRangeTransform, createSumTransform, createCumSumTransform, createRowTransform, createColumnTransform, createIndexTransform, createMeanTransform, createVarianceTransform, createStdTransform } from '../factoriesAny.js';
+import { createNode, createObjectNode, createOperatorNode, createParenthesisNode, createRelationalNode, createArrayNode, createBlockNode, createConditionalNode, createConstantNode, createRangeNode, createReviver, createChainClass, createFunctionAssignmentNode, createChain, createAccessorNode, createAssignmentNode, createIndexNode, createSymbolNode, createFunctionNode, createParse, createResolve, createSimplifyConstant, createCompile, createHelpClass, createLeafCount, createSimplifyCore, createEvaluate, createHelp, createParserClass, createSimplify, createSymbolicEqual, createDerivative, createParser, createRationalize, 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, 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, variance, classicalElectronRadius, lusolve, median, setUnion, std, gasConstant, mad, 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, secondRadiation, setSymDifference, speedOfLight, stefanBoltzmann, thomsonCrossSection, avogadro, bohrRadius, coulomb, divide, electronMass, factorial, gravity, inverseConductanceQuantum, lup, magneticFluxQuantum, molarMassC12, multinomial, permutations, planckMass, polynomialRoot, setIntersect, stirlingS2, unit, bellNumbers, eigs, fermiCoupling, mean, molarVolume, planckConstant, quantumOfCirculation, setUnion, variance, classicalElectronRadius, lusolve, median, kldivergence, mad, std, gasConstant, norm, rotationMatrix, planckTime, schur, rotate, sylvester, lyap } from './pureFunctionsAny.generated.js';
 var math = {}; // NOT pure!
 var mathWithTransform = {}; // NOT pure!
 var classes = {}; // NOT pure!
@@ -169,10 +169,6 @@ export var help = createHelp({
 export var Parser = createParserClass({
   evaluate
 });
-export var parser = createParser({
-  Parser,
-  typed
-});
 export var simplify = createSimplify({
   bignumber,
   fraction,
@@ -221,6 +217,10 @@ export var derivative = createDerivative({
   simplify,
   typed
 });
+export var parser = createParser({
+  Parser,
+  typed
+});
 export var rationalize = createRationalize({
   bignumber,
   fraction,
@@ -490,7 +490,8 @@ _extends(math, {
   protonMass,
   reducedPlanckConstant,
   rydberg,
-  setIntersect,
+  secondRadiation,
+  setSymDifference,
   speedOfLight,
   stefanBoltzmann,
   thomsonCrossSection,
@@ -509,9 +510,9 @@ _extends(math, {
   parse,
   permutations,
   planckMass,
-  quantumOfCirculation,
+  polynomialRoot,
   resolve,
-  secondRadiation,
+  setIntersect,
   simplifyConstant,
   stirlingS2,
   unit,
@@ -523,7 +524,8 @@ _extends(math, {
   mean,
   molarVolume,
   planckConstant,
-  setSymDifference,
+  quantumOfCirculation,
+  setUnion,
   simplifyCore,
   variance,
   classicalElectronRadius,
@@ -531,20 +533,22 @@ _extends(math, {
   help,
   lusolve,
   median,
-  setUnion,
-  std,
-  gasConstant,
-  mad,
-  parser,
   simplify,
   symbolicEqual,
   derivative,
-  norm,
+  kldivergence,
+  mad,
+  parser,
   rationalize,
+  std,
+  gasConstant,
+  norm,
   rotationMatrix,
-  kldivergence,
-  rotate,
   planckTime,
+  schur,
+  rotate,
+  sylvester,
+  lyap,
   config
 });
 _extends(mathWithTransform, math, {

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

@@ -5,7 +5,7 @@
 import { config } from './configReadonly.js';
 import { createBigNumberClass, createComplexClass, createE, createFalse, createFineStructure, createFractionClass, createI, createInfinity, createLN10, createLOG10E, createMatrixClass, createNaN, createNull, createPhi, createRangeClass, createResultSet, createSQRT1_2,
 // eslint-disable-line camelcase
-createSackurTetrode, createTau, createTrue, createVersion, createDenseMatrixClass, createEfimovFactor, createLN2, createPi, createReplacer, createSQRT2, createTyped, createUnaryPlus, createWeakMixingAngle, createAbs, createAcos, createAcot, createAcsc, createAddScalar, createArg, createAsech, createAsinh, createAtan, createAtanh, createBignumber, createBitNot, createBoolean, createClone, createCombinations, createComplex, createConj, createCosh, createCoth, createCsc, createCube, createEqualScalar, createErf, createExp, createExpm1, createFilter, createForEach, createFormat, createGetMatrixDataType, createHex, createIm, createIsInteger, createIsNegative, createIsPositive, createIsZero, createLOG2E, createLgamma, createLog10, createLog2, createMap, createMultiplyScalar, createNot, createNumber, createOct, createPickRandom, createPrint, createRandom, createRe, createSec, createSign, createSin, createSparseMatrixClass, createSplitUnit, createSquare, createString, createTan, createTypeOf, createAcosh, createAcsch, createApply, createAsec, createBin, createCombinationsWithRep, createCos, createCsch, createIsNaN, createIsPrime, createRandomInt, createSech, createSinh, createSparse, createSqrt, createTanh, createUnaryMinus, createAcoth, createCot, createFraction, createIsNumeric, createMatrix, createMatrixFromFunction, createMod, createNthRoot, createNumeric, createOr, createProd, createReshape, createSize, createSmaller, createSqueeze, createSubset, createSubtract, createTo, createTranspose, createXgcd, createZeros, createAnd, createBitAnd, createBitXor, createCbrt, createCompare, createCompareText, createConcat, createCount, createCtranspose, createDiag, createDivideScalar, createDotDivide, createEqual, createFft, createFlatten, createGcd, createHasNumericValue, createHypot, createIfft, createKron, createLargerEq, createLeftShift, createLsolve, createMatrixFromColumns, createMin, createMode, createNthRoots, createOnes, createPartitionSelect, createResize, createRightArithShift, createRound, createSmallerEq, createUnequal, createUsolve, createXor, createAdd, createAtan2, createBitOr, createCatalan, createCompareNatural, createCumSum, createDeepEqual, createDiff, createDot, createEqualText, createFloor, createIdentity, createInvmod, createLarger, createLog, createLsolveAll, createMatrixFromRows, createMultiply, createQr, createRange, createRightLogShift, createSetSize, createSlu, createSum, createTrace, createUsolveAll, createAsin, createCeil, createComposition, createCross, createDet, createDistance, createDotMultiply, createFibonacciHeapClass, createFix, createImmutableDenseMatrixClass, createIndexClass, createIntersect, createLcm, createLog1p, createMax, createQuantileSeq, createRow, createSetCartesian, createSetDistinct, createSetIsSubset, createSetPowerset, createSort, createColumn, createIndex, createInv, createPinv, createPow, createSetDifference, createSetMultiplicity, createSpaClass, createSqrtm, createUnitClass, createVacuumImpedance, createWienDisplacement, createAtomicMass, createBohrMagneton, createBoltzmann, createConductanceQuantum, createCreateUnit, createDeuteronMass, createDotPow, createElectricConstant, createElementaryCharge, createExpm, createFaraday, createFirstRadiation, createGamma, createGravitationConstant, createHartreeEnergy, createKlitzing, createLoschmidt, createMagneticConstant, createMolarMass, createMolarPlanckConstant, createNeutronMass, createNuclearMagneton, createPlanckCharge, createPlanckLength, createPlanckTemperature, createProtonMass, createReducedPlanckConstant, createRydberg, createSetIntersect, createSpeedOfLight, createStefanBoltzmann, createThomsonCrossSection, createAvogadro, createBohrRadius, createCoulomb, createDivide, createElectronMass, createFactorial, createGravity, createInverseConductanceQuantum, createLup, createMagneticFluxQuantum, createMolarMassC12, createMultinomial, createPermutations, createPlanckMass, createQuantumOfCirculation, createSecondRadiation, createStirlingS2, createUnitFunction, createBellNumbers, createEigs, createFermiCoupling, createMean, createMolarVolume, createPlanckConstant, createSetSymDifference, createVariance, createClassicalElectronRadius, createLusolve, createMedian, createSetUnion, createStd, createGasConstant, createMad, createNorm, createRotationMatrix, createKldivergence, createRotate, createPlanckTime } from '../factoriesAny.js';
+createSackurTetrode, createTau, createTrue, createVersion, createDenseMatrixClass, createEfimovFactor, createLN2, createPi, createReplacer, createSQRT2, createTyped, createUnaryPlus, createWeakMixingAngle, createAbs, createAcos, createAcot, createAcsc, createAddScalar, createArg, createAsech, createAsinh, createAtan, createAtanh, createBignumber, createBitNot, createBoolean, createClone, createCombinations, createComplex, createConj, createCosh, createCoth, createCsc, createCube, createEqualScalar, createErf, createExp, createExpm1, createFilter, createForEach, createFormat, createGetMatrixDataType, createHex, createIm, createIsInteger, createIsNegative, createIsPositive, createIsZero, createLOG2E, createLgamma, createLog10, createLog2, createMap, createMultiplyScalar, createNot, createNumber, createOct, createPickRandom, createPrint, createRandom, createRe, createSec, createSign, createSin, createSparseMatrixClass, createSplitUnit, createSquare, createString, createTan, createTypeOf, createAcosh, createAcsch, createApply, createAsec, createBin, createCombinationsWithRep, createCos, createCsch, createIsNaN, createIsPrime, createRandomInt, createSech, createSinh, createSparse, createSqrt, createTanh, createUnaryMinus, createAcoth, createCot, createFraction, createIsNumeric, createMatrix, createMatrixFromFunction, createMod, createNthRoot, createNumeric, createOr, createProd, createReshape, createSize, createSmaller, createSqueeze, createSubset, createSubtract, createTo, createTranspose, createXgcd, createZeros, createAnd, createBitAnd, createBitXor, createCbrt, createCompare, createCompareText, createConcat, createCount, createCtranspose, createDiag, createDivideScalar, createDotDivide, createEqual, createFft, createFlatten, createGcd, createHasNumericValue, createHypot, createIfft, createKron, createLargerEq, createLeftShift, createLsolve, createMatrixFromColumns, createMin, createMode, createNthRoots, createOnes, createPartitionSelect, createResize, createRightArithShift, createRound, createSmallerEq, createUnequal, createUsolve, createXor, createAdd, createAtan2, createBitOr, createCatalan, createCompareNatural, createCumSum, createDeepEqual, createDiff, createDot, createEqualText, createFloor, createIdentity, createInvmod, createLarger, createLog, createLsolveAll, createMatrixFromRows, createMultiply, createQr, createRange, createRightLogShift, createSetSize, createSlu, createSum, createTrace, createUsolveAll, createAsin, createCeil, createComposition, createCross, createDet, createDistance, createDotMultiply, createFibonacciHeapClass, createFix, createImmutableDenseMatrixClass, createIndexClass, createIntersect, createLcm, createLog1p, createMax, createQuantileSeq, createRow, createSetCartesian, createSetDistinct, createSetIsSubset, createSetPowerset, createSort, createColumn, createIndex, createInv, createPinv, createPow, createSetDifference, createSetMultiplicity, createSpaClass, createSqrtm, createUnitClass, createVacuumImpedance, createWienDisplacement, createAtomicMass, createBohrMagneton, createBoltzmann, createConductanceQuantum, createCreateUnit, createDeuteronMass, createDotPow, createElectricConstant, createElementaryCharge, createExpm, createFaraday, createFirstRadiation, createGamma, createGravitationConstant, createHartreeEnergy, createKlitzing, createLoschmidt, createMagneticConstant, createMolarMass, createMolarPlanckConstant, createNeutronMass, createNuclearMagneton, createPlanckCharge, createPlanckLength, createPlanckTemperature, createProtonMass, createReducedPlanckConstant, createRydberg, createSecondRadiation, createSetSymDifference, createSpeedOfLight, createStefanBoltzmann, createThomsonCrossSection, createAvogadro, createBohrRadius, createCoulomb, createDivide, createElectronMass, createFactorial, createGravity, createInverseConductanceQuantum, createLup, createMagneticFluxQuantum, createMolarMassC12, createMultinomial, createPermutations, createPlanckMass, createPolynomialRoot, createSetIntersect, createStirlingS2, createUnitFunction, createBellNumbers, createEigs, createFermiCoupling, createMean, createMolarVolume, createPlanckConstant, createQuantumOfCirculation, createSetUnion, createVariance, createClassicalElectronRadius, createLusolve, createMedian, createKldivergence, createMad, createStd, createGasConstant, createNorm, createRotationMatrix, createPlanckTime, createSchur, createRotate, createSylvester, createLyap } from '../factoriesAny.js';
 export var BigNumber = /* #__PURE__ */createBigNumberClass({
   config
 });
@@ -1308,10 +1308,15 @@ export var rydberg = /* #__PURE__ */createRydberg({
   Unit,
   config
 });
-export var setIntersect = /* #__PURE__ */createSetIntersect({
-  DenseMatrix,
+export var secondRadiation = /* #__PURE__ */createSecondRadiation({
+  BigNumber,
+  Unit,
+  config
+});
+export var setSymDifference = /* #__PURE__ */createSetSymDifference({
   Index,
-  compareNatural,
+  concat,
+  setDifference,
   size,
   subset,
   typed
@@ -1416,15 +1421,28 @@ export var planckMass = /* #__PURE__ */createPlanckMass({
   Unit,
   config
 });
-export var quantumOfCirculation = /* #__PURE__ */createQuantumOfCirculation({
-  BigNumber,
-  Unit,
-  config
+export var polynomialRoot = /* #__PURE__ */createPolynomialRoot({
+  add,
+  cbrt,
+  divide,
+  equalScalar,
+  im,
+  isZero,
+  multiply,
+  re,
+  sqrt,
+  subtract,
+  typeOf,
+  typed,
+  unaryMinus
 });
-export var secondRadiation = /* #__PURE__ */createSecondRadiation({
-  BigNumber,
-  Unit,
-  config
+export var setIntersect = /* #__PURE__ */createSetIntersect({
+  DenseMatrix,
+  Index,
+  compareNatural,
+  size,
+  subset,
+  typed
 });
 export var stirlingS2 = /* #__PURE__ */createStirlingS2({
   bignumber,
@@ -1505,10 +1523,16 @@ export var planckConstant = /* #__PURE__ */createPlanckConstant({
   Unit,
   config
 });
-export var setSymDifference = /* #__PURE__ */createSetSymDifference({
+export var quantumOfCirculation = /* #__PURE__ */createQuantumOfCirculation({
+  BigNumber,
+  Unit,
+  config
+});
+export var setUnion = /* #__PURE__ */createSetUnion({
   Index,
   concat,
-  setDifference,
+  setIntersect,
+  setSymDifference,
   size,
   subset,
   typed
@@ -1543,13 +1567,22 @@ export var median = /* #__PURE__ */createMedian({
   partitionSelect,
   typed
 });
-export var setUnion = /* #__PURE__ */createSetUnion({
-  Index,
-  concat,
-  setIntersect,
-  setSymDifference,
-  size,
-  subset,
+export var kldivergence = /* #__PURE__ */createKldivergence({
+  divide,
+  dotDivide,
+  isNumeric,
+  log,
+  map,
+  matrix,
+  multiply,
+  sum,
+  typed
+});
+export var mad = /* #__PURE__ */createMad({
+  abs,
+  map,
+  median,
+  subtract,
   typed
 });
 export var std = /* #__PURE__ */createStd({
@@ -1563,13 +1596,6 @@ export var gasConstant = /* #__PURE__ */createGasConstant({
   Unit,
   config
 });
-export var mad = /* #__PURE__ */createMad({
-  abs,
-  map,
-  median,
-  subtract,
-  typed
-});
 export var norm = /* #__PURE__ */createNorm({
   abs,
   add,
@@ -1599,15 +1625,18 @@ export var rotationMatrix = /* #__PURE__ */createRotationMatrix({
   typed,
   unaryMinus
 });
-export var kldivergence = /* #__PURE__ */createKldivergence({
-  divide,
-  dotDivide,
-  isNumeric,
-  log,
-  map,
+export var planckTime = /* #__PURE__ */createPlanckTime({
+  BigNumber,
+  Unit,
+  config
+});
+export var schur = /* #__PURE__ */createSchur({
+  identity,
   matrix,
   multiply,
-  sum,
+  norm,
+  qr,
+  subtract,
   typed
 });
 export var rotate = /* #__PURE__ */createRotate({
@@ -1615,8 +1644,27 @@ export var rotate = /* #__PURE__ */createRotate({
   rotationMatrix,
   typed
 });
-export var planckTime = /* #__PURE__ */createPlanckTime({
-  BigNumber,
-  Unit,
-  config
+export var sylvester = /* #__PURE__ */createSylvester({
+  abs,
+  add,
+  concat,
+  identity,
+  index,
+  lusolve,
+  matrix,
+  matrixFromColumns,
+  multiply,
+  range,
+  schur,
+  subset,
+  subtract,
+  transpose,
+  typed
+});
+export var lyap = /* #__PURE__ */createLyap({
+  matrix,
+  multiply,
+  sylvester,
+  transpose,
+  typed
 });

package/lib/esm/expression/embeddedDocs/embeddedDocs.js

@@ -36,6 +36,7 @@ import { lsolveDocs } from './function/algebra/lsolve.js';
 import { lsolveAllDocs } from './function/algebra/lsolveAll.js';
 import { lupDocs } from './function/algebra/lup.js';
 import { lusolveDocs } from './function/algebra/lusolve.js';
+import { polynomialRootDocs } from './function/algebra/polynomialRoot.js';
 import { qrDocs } from './function/algebra/qr.js';
 import { rationalizeDocs } from './function/algebra/rationalize.js';
 import { resolveDocs } from './function/algebra/resolve.js';
@@ -78,6 +79,9 @@ import { roundDocs } from './function/arithmetic/round.js';
 import { signDocs } from './function/arithmetic/sign.js';
 import { sqrtDocs } from './function/arithmetic/sqrt.js';
 import { sqrtmDocs } from './function/arithmetic/sqrtm.js';
+import { sylvesterDocs } from './function/matrix/sylvester.js';
+import { schurDocs } from './function/matrix/schur.js';
+import { lyapDocs } from './function/matrix/lyap.js';
 import { squareDocs } from './function/arithmetic/square.js';
 import { subtractDocs } from './function/arithmetic/subtract.js';
 import { unaryMinusDocs } from './function/arithmetic/unaryMinus.js';
@@ -478,6 +482,7 @@ export var embeddedDocs = {
   lup: lupDocs,
   lusolve: lusolveDocs,
   leafCount: leafCountDocs,
+  polynomialRoot: polynomialRootDocs,
   resolve: resolveDocs,
   simplify: simplifyDocs,
   simplifyConstant: simplifyConstantDocs,
@@ -599,6 +604,9 @@ export var embeddedDocs = {
   zeros: zerosDocs,
   fft: fftDocs,
   ifft: ifftDocs,
+  sylvester: sylvesterDocs,
+  schur: schurDocs,
+  lyap: lyapDocs,
   // functions - probability
   combinations: combinationsDocs,
   combinationsWithRep: combinationsWithRepDocs,

package/lib/esm/expression/embeddedDocs/function/algebra/polynomialRoot.js

@@ -0,0 +1,8 @@
+export var polynomialRootDocs = {
+  name: 'polynomialRoot',
+  category: 'Algebra',
+  syntax: ['x=polynomialRoot(-6, 3)', 'x=polynomialRoot(4, -4, 1)', 'x=polynomialRoot(-8, 12, -6, 1)'],
+  description: 'Finds the roots of a univariate polynomial given by its coefficients starting from constant, linear, and so on, increasing in degree.',
+  examples: ['a = polynomialRoot(-6, 11, -6 1)'],
+  seealso: ['cbrt', 'sqrt']
+};