mathjs 14.9.1 → 15.0.0

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

@@ -70,7 +70,8 @@ const tau = exports.tau = /* #__PURE__ */(0, _factoriesAny.createTau)({
 const _true = exports._true = /* #__PURE__ */(0, _factoriesAny.createTrue)({});
 const version = exports.version = /* #__PURE__ */(0, _factoriesAny.createVersion)({});
 const DenseMatrix = exports.DenseMatrix = /* #__PURE__ */(0, _factoriesAny.createDenseMatrixClass)({
-  Matrix
+  Matrix,
+  config: _configReadonly.config
 });
 const efimovFactor = exports.efimovFactor = /* #__PURE__ */(0, _factoriesAny.createEfimovFactor)({
   BigNumber,
@@ -288,6 +289,9 @@ const sign = exports.sign = /* #__PURE__ */(0, _factoriesAny.createSign)({
 const sin = exports.sin = /* #__PURE__ */(0, _factoriesAny.createSin)({
   typed
 });
+const size = exports.size = /* #__PURE__ */(0, _factoriesAny.createSize)({
+  typed
+});
 const SparseMatrix = exports.SparseMatrix = /* #__PURE__ */(0, _factoriesAny.createSparseMatrixClass)({
   Matrix,
   equalScalar,
@@ -343,6 +347,13 @@ const csch = exports.csch = /* #__PURE__ */(0, _factoriesAny.createCsch)({
   BigNumber,
   typed
 });
+const dot = exports.dot = /* #__PURE__ */(0, _factoriesAny.createDot)({
+  addScalar,
+  conj,
+  multiplyScalar,
+  size,
+  typed
+});
 const isNaN = exports.isNaN = /* #__PURE__ */(0, _factoriesAny.createIsNaN)({
   typed
 });
@@ -365,6 +376,14 @@ const matrixFromFunction = exports.matrixFromFunction = /* #__PURE__ */(0, _fact
   matrix,
   typed
 });
+const multiply = exports.multiply = /* #__PURE__ */(0, _factoriesAny.createMultiply)({
+  addScalar,
+  dot,
+  equalScalar,
+  matrix,
+  multiplyScalar,
+  typed
+});
 const ones = exports.ones = /* #__PURE__ */(0, _factoriesAny.createOnes)({
   BigNumber,
   config: _configReadonly.config,
@@ -505,6 +524,12 @@ const leftShift = exports.leftShift = /* #__PURE__ */(0, _factoriesAny.createLef
   typed,
   zeros
 });
+const matrixFromRows = exports.matrixFromRows = /* #__PURE__ */(0, _factoriesAny.createMatrixFromRows)({
+  flatten,
+  matrix,
+  size,
+  typed
+});
 const mode = exports.mode = /* #__PURE__ */(0, _factoriesAny.createMode)({
   isNaN,
   isNumeric,
@@ -549,11 +574,6 @@ const round = exports.round = /* #__PURE__ */(0, _factoriesAny.createRound)({
   typed,
   zeros
 });
-const size = exports.size = /* #__PURE__ */(0, _factoriesAny.createSize)({
-  matrix,
-  config: _configReadonly.config,
-  typed
-});
 const smaller = exports.smaller = /* #__PURE__ */(0, _factoriesAny.createSmaller)({
   DenseMatrix,
   SparseMatrix,
@@ -662,15 +682,6 @@ const divideScalar = exports.divideScalar = /* #__PURE__ */(0, _factoriesAny.cre
   numeric,
   typed
 });
-const dotDivide = exports.dotDivide = /* #__PURE__ */(0, _factoriesAny.createDotDivide)({
-  DenseMatrix,
-  SparseMatrix,
-  concat,
-  divideScalar,
-  equalScalar,
-  matrix,
-  typed
-});
 const equalText = exports.equalText = /* #__PURE__ */(0, _factoriesAny.createEqualText)({
   compareText,
   isZero,
@@ -749,13 +760,6 @@ const matrixFromColumns = exports.matrixFromColumns = /* #__PURE__ */(0, _factor
   size,
   typed
 });
-const max = exports.max = /* #__PURE__ */(0, _factoriesAny.createMax)({
-  config: _configReadonly.config,
-  isNaN,
-  larger,
-  numeric,
-  typed
-});
 const min = exports.min = /* #__PURE__ */(0, _factoriesAny.createMin)({
   config: _configReadonly.config,
   isNaN,
@@ -861,6 +865,13 @@ const usolve = exports.usolve = /* #__PURE__ */(0, _factoriesAny.createUsolve)({
   subtractScalar,
   typed
 });
+const zpk2tf = exports.zpk2tf = /* #__PURE__ */(0, _factoriesAny.createZpk2tf)({
+  Complex,
+  add,
+  multiply,
+  number,
+  typed
+});
 const catalan = exports.catalan = /* #__PURE__ */(0, _factoriesAny.createCatalan)({
   addScalar,
   combinations,
@@ -883,6 +894,21 @@ const composition = exports.composition = /* #__PURE__ */(0, _factoriesAny.creat
   larger,
   typed
 });
+const cross = exports.cross = /* #__PURE__ */(0, _factoriesAny.createCross)({
+  matrix,
+  multiply,
+  subtract,
+  typed
+});
+const det = exports.det = /* #__PURE__ */(0, _factoriesAny.createDet)({
+  divideScalar,
+  isZero,
+  matrix,
+  multiply,
+  subtractScalar,
+  typed,
+  unaryMinus
+});
 const diff = exports.diff = /* #__PURE__ */(0, _factoriesAny.createDiff)({
   matrix,
   number,
@@ -899,11 +925,13 @@ const distance = exports.distance = /* #__PURE__ */(0, _factoriesAny.createDista
   subtractScalar,
   typed
 });
-const dot = exports.dot = /* #__PURE__ */(0, _factoriesAny.createDot)({
-  addScalar,
-  conj,
-  multiplyScalar,
-  size,
+const dotDivide = exports.dotDivide = /* #__PURE__ */(0, _factoriesAny.createDotDivide)({
+  DenseMatrix,
+  SparseMatrix,
+  concat,
+  divideScalar,
+  equalScalar,
+  matrix,
   typed
 });
 const FibonacciHeap = exports.FibonacciHeap = /* #__PURE__ */(0, _factoriesAny.createFibonacciHeapClass)({
@@ -914,6 +942,23 @@ const index = exports.index = /* #__PURE__ */(0, _factoriesAny.createIndex)({
   Index,
   typed
 });
+const intersect = exports.intersect = /* #__PURE__ */(0, _factoriesAny.createIntersect)({
+  abs,
+  add,
+  addScalar,
+  config: _configReadonly.config,
+  divideScalar,
+  equalScalar,
+  flatten,
+  isNumeric,
+  isZero,
+  matrix,
+  multiply,
+  multiplyScalar,
+  smaller,
+  subtract,
+  typed
+});
 const invmod = exports.invmod = /* #__PURE__ */(0, _factoriesAny.createInvmod)({
   BigNumber,
   add,
@@ -947,18 +992,11 @@ const lsolveAll = exports.lsolveAll = /* #__PURE__ */(0, _factoriesAny.createLso
   subtractScalar,
   typed
 });
-const matrixFromRows = exports.matrixFromRows = /* #__PURE__ */(0, _factoriesAny.createMatrixFromRows)({
-  flatten,
-  matrix,
-  size,
-  typed
-});
-const multiply = exports.multiply = /* #__PURE__ */(0, _factoriesAny.createMultiply)({
-  addScalar,
-  dot,
-  equalScalar,
-  matrix,
-  multiplyScalar,
+const max = exports.max = /* #__PURE__ */(0, _factoriesAny.createMax)({
+  config: _configReadonly.config,
+  isNaN,
+  larger,
+  numeric,
   typed
 });
 const range = exports.range = /* #__PURE__ */(0, _factoriesAny.createRange)({
@@ -1043,13 +1081,6 @@ const usolveAll = exports.usolveAll = /* #__PURE__ */(0, _factoriesAny.createUso
   subtractScalar,
   typed
 });
-const zpk2tf = exports.zpk2tf = /* #__PURE__ */(0, _factoriesAny.createZpk2tf)({
-  Complex,
-  add,
-  multiply,
-  number,
-  typed
-});
 const and = exports.and = /* #__PURE__ */(0, _factoriesAny.createAnd)({
   concat,
   equalScalar,
@@ -1073,20 +1104,10 @@ const column = exports.column = /* #__PURE__ */(0, _factoriesAny.createColumn)({
   range,
   typed
 });
-const cross = exports.cross = /* #__PURE__ */(0, _factoriesAny.createCross)({
-  matrix,
-  multiply,
-  subtract,
-  typed
-});
-const det = exports.det = /* #__PURE__ */(0, _factoriesAny.createDet)({
-  divideScalar,
-  isZero,
-  matrix,
-  multiply,
-  subtractScalar,
+const cumsum = exports.cumsum = /* #__PURE__ */(0, _factoriesAny.createCumSum)({
+  add,
   typed,
-  unaryMinus
+  unaryPlus
 });
 const fix = exports.fix = /* #__PURE__ */(0, _factoriesAny.createFix)({
   Complex,
@@ -1098,6 +1119,40 @@ const fix = exports.fix = /* #__PURE__ */(0, _factoriesAny.createFix)({
   typed,
   zeros
 });
+const setDifference = exports.setDifference = /* #__PURE__ */(0, _factoriesAny.createSetDifference)({
+  DenseMatrix,
+  Index,
+  compareNatural,
+  size,
+  subset,
+  typed
+});
+const setMultiplicity = exports.setMultiplicity = /* #__PURE__ */(0, _factoriesAny.createSetMultiplicity)({
+  Index,
+  compareNatural,
+  size,
+  subset,
+  typed
+});
+const setSymDifference = exports.setSymDifference = /* #__PURE__ */(0, _factoriesAny.createSetSymDifference)({
+  Index,
+  concat,
+  setDifference,
+  size,
+  subset,
+  typed
+});
+const Spa = exports.Spa = /* #__PURE__ */(0, _factoriesAny.createSpaClass)({
+  FibonacciHeap,
+  addScalar,
+  equalScalar
+});
+const sum = exports.sum = /* #__PURE__ */(0, _factoriesAny.createSum)({
+  add,
+  config: _configReadonly.config,
+  numeric,
+  typed
+});
 const inv = exports.inv = /* #__PURE__ */(0, _factoriesAny.createInv)({
   abs,
   addScalar,
@@ -1109,6 +1164,21 @@ const inv = exports.inv = /* #__PURE__ */(0, _factoriesAny.createInv)({
   typed,
   unaryMinus
 });
+const lup = exports.lup = /* #__PURE__ */(0, _factoriesAny.createLup)({
+  DenseMatrix,
+  Spa,
+  SparseMatrix,
+  abs,
+  addScalar,
+  divideScalar,
+  equalScalar,
+  larger,
+  matrix,
+  multiplyScalar,
+  subtractScalar,
+  typed,
+  unaryMinus
+});
 const pinv = exports.pinv = /* #__PURE__ */(0, _factoriesAny.createPinv)({
   Complex,
   add,
@@ -1134,7 +1204,7 @@ const pow = exports.pow = /* #__PURE__ */(0, _factoriesAny.createPow)({
   number,
   typed
 });
-const setDifference = exports.setDifference = /* #__PURE__ */(0, _factoriesAny.createSetDifference)({
+const setIntersect = exports.setIntersect = /* #__PURE__ */(0, _factoriesAny.createSetIntersect)({
   DenseMatrix,
   Index,
   compareNatural,
@@ -1142,26 +1212,15 @@ const setDifference = exports.setDifference = /* #__PURE__ */(0, _factoriesAny.c
   subset,
   typed
 });
-const setMultiplicity = exports.setMultiplicity = /* #__PURE__ */(0, _factoriesAny.createSetMultiplicity)({
-  Index,
-  compareNatural,
-  size,
-  subset,
-  typed
-});
-const setSymDifference = exports.setSymDifference = /* #__PURE__ */(0, _factoriesAny.createSetSymDifference)({
+const setUnion = exports.setUnion = /* #__PURE__ */(0, _factoriesAny.createSetUnion)({
   Index,
   concat,
-  setDifference,
+  setIntersect,
+  setSymDifference,
   size,
   subset,
   typed
 });
-const Spa = exports.Spa = /* #__PURE__ */(0, _factoriesAny.createSpaClass)({
-  FibonacciHeap,
-  addScalar,
-  equalScalar
-});
 const sqrtm = exports.sqrtm = /* #__PURE__ */(0, _factoriesAny.createSqrtm)({
   abs,
   add,
@@ -1175,12 +1234,6 @@ const sqrtm = exports.sqrtm = /* #__PURE__ */(0, _factoriesAny.createSqrtm)({
   subtract,
   typed
 });
-const sum = exports.sum = /* #__PURE__ */(0, _factoriesAny.createSum)({
-  add,
-  config: _configReadonly.config,
-  numeric,
-  typed
-});
 const Unit = exports.Unit = /* #__PURE__ */(0, _factoriesAny.createUnitClass)({
   BigNumber,
   Complex,
@@ -1315,11 +1368,6 @@ const ifft = exports.ifft = /* #__PURE__ */(0, _factoriesAny.createIfft)({
   fft,
   typed
 });
-const inverseConductanceQuantum = exports.inverseConductanceQuantum = /* #__PURE__ */(0, _factoriesAny.createInverseConductanceQuantum)({
-  BigNumber,
-  Unit,
-  config: _configReadonly.config
-});
 const klitzing = exports.klitzing = /* #__PURE__ */(0, _factoriesAny.createKlitzing)({
   BigNumber,
   Unit,
@@ -1395,10 +1443,6 @@ const secondRadiation = exports.secondRadiation = /* #__PURE__ */(0, _factoriesA
   Unit,
   config: _configReadonly.config
 });
-const setSize = exports.setSize = /* #__PURE__ */(0, _factoriesAny.createSetSize)({
-  compareNatural,
-  typed
-});
 const speedOfLight = exports.speedOfLight = /* #__PURE__ */(0, _factoriesAny.createSpeedOfLight)({
   BigNumber,
   Unit,
@@ -1456,37 +1500,19 @@ const gravity = exports.gravity = /* #__PURE__ */(0, _factoriesAny.createGravity
   Unit,
   config: _configReadonly.config
 });
-const intersect = exports.intersect = /* #__PURE__ */(0, _factoriesAny.createIntersect)({
-  abs,
-  add,
-  addScalar,
-  config: _configReadonly.config,
-  divideScalar,
-  equalScalar,
-  flatten,
-  isNumeric,
-  isZero,
-  matrix,
-  multiply,
-  multiplyScalar,
-  smaller,
-  subtract,
-  typed
+const inverseConductanceQuantum = exports.inverseConductanceQuantum = /* #__PURE__ */(0, _factoriesAny.createInverseConductanceQuantum)({
+  BigNumber,
+  Unit,
+  config: _configReadonly.config
 });
-const lup = exports.lup = /* #__PURE__ */(0, _factoriesAny.createLup)({
+const lusolve = exports.lusolve = /* #__PURE__ */(0, _factoriesAny.createLusolve)({
   DenseMatrix,
-  Spa,
-  SparseMatrix,
-  abs,
-  addScalar,
-  divideScalar,
-  equalScalar,
-  larger,
+  lsolve,
+  lup,
   matrix,
-  multiplyScalar,
-  subtractScalar,
+  slu,
   typed,
-  unaryMinus
+  usolve
 });
 const magneticFluxQuantum = exports.magneticFluxQuantum = /* #__PURE__ */(0, _factoriesAny.createMagneticFluxQuantum)({
   BigNumber,
@@ -1531,12 +1557,8 @@ const polynomialRoot = exports.polynomialRoot = /* #__PURE__ */(0, _factoriesAny
   typed,
   unaryMinus
 });
-const setIntersect = exports.setIntersect = /* #__PURE__ */(0, _factoriesAny.createSetIntersect)({
-  DenseMatrix,
-  Index,
+const setSize = exports.setSize = /* #__PURE__ */(0, _factoriesAny.createSetSize)({
   compareNatural,
-  size,
-  subset,
   typed
 });
 const solveODE = exports.solveODE = /* #__PURE__ */(0, _factoriesAny.createSolveODE)({
@@ -1582,11 +1604,6 @@ const bellNumbers = exports.bellNumbers = /* #__PURE__ */(0, _factoriesAny.creat
   stirlingS2,
   typed
 });
-const cumsum = exports.cumsum = /* #__PURE__ */(0, _factoriesAny.createCumSum)({
-  add,
-  typed,
-  unaryPlus
-});
 const eigs = exports.eigs = /* #__PURE__ */(0, _factoriesAny.createEigs)({
   abs,
   add,
@@ -1643,15 +1660,6 @@ const kldivergence = exports.kldivergence = /* #__PURE__ */(0, _factoriesAny.cre
   sum,
   typed
 });
-const lusolve = exports.lusolve = /* #__PURE__ */(0, _factoriesAny.createLusolve)({
-  DenseMatrix,
-  lsolve,
-  lup,
-  matrix,
-  slu,
-  typed,
-  usolve
-});
 const mean = exports.mean = /* #__PURE__ */(0, _factoriesAny.createMean)({
   add,
   divide,
@@ -1682,15 +1690,6 @@ const quantileSeq = exports.quantileSeq = /* #__PURE__ */(0, _factoriesAny.creat
   subtract,
   typed
 });
-const setUnion = exports.setUnion = /* #__PURE__ */(0, _factoriesAny.createSetUnion)({
-  Index,
-  concat,
-  setIntersect,
-  setSymDifference,
-  size,
-  subset,
-  typed
-});
 const variance = exports.variance = /* #__PURE__ */(0, _factoriesAny.createVariance)({
   add,
   divide,

package/lib/cjs/entry/pureFunctionsNumber.generated.js

@@ -218,8 +218,6 @@ const sin = exports.sin = /* #__PURE__ */(0, _factoriesNumber.createSin)({
   typed
 });
 const size = exports.size = /* #__PURE__ */(0, _factoriesNumber.createSize)({
-  matrix,
-  config: _configReadonly.config,
   typed
 });
 const smallerEq = exports.smallerEq = /* #__PURE__ */(0, _factoriesNumber.createSmallerEq)({

package/lib/cjs/expression/node/IndexNode.js

@@ -86,7 +86,7 @@ const createIndexNode = exports.createIndexNode = /* #__PURE__ */(0, _factory.fa
             if (!(0, _is.isMatrix)(context) && !(0, _is.isArray)(context) && !(0, _is.isString)(context)) {
               throw new TypeError('Cannot resolve "end": ' + 'context must be a Matrix, Array, or string but is ' + (0, _is.typeOf)(context));
             }
-            const s = size(context).valueOf();
+            const s = size(context);
             const childArgs = Object.create(args);
             childArgs.end = s[i];
             return _evalDimension(scope, childArgs, context);

package/lib/cjs/expression/parse.js

@@ -329,7 +329,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
         next(state);
         state.token += currentCharacter(state);
         next(state);
-        while (parse.isHexDigit(currentCharacter(state))) {
+        while (parse.isAlpha(currentCharacter(state), prevCharacter(state), nextCharacter(state)) || parse.isDigit(currentCharacter(state))) {
           state.token += currentCharacter(state);
           next(state);
         }
@@ -338,7 +338,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
           state.token += '.';
           next(state);
           // get the digits after the radix
-          while (parse.isHexDigit(currentCharacter(state))) {
+          while (parse.isAlpha(currentCharacter(state), prevCharacter(state), nextCharacter(state)) || parse.isDigit(currentCharacter(state))) {
             state.token += currentCharacter(state);
             next(state);
           }
@@ -544,15 +544,6 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
     return c >= '0' && c <= '9';
   };
 
-  /**
-   * checks if the given char c is a hex digit
-   * @param {string} c   a string with one character
-   * @return {boolean}
-   */
-  parse.isHexDigit = function isHexDigit(c) {
-    return c >= '0' && c <= '9' || c >= 'a' && c <= 'f' || c >= 'A' && c <= 'F';
-  };
-
   /**
    * Start of the parse levels below, in order of precedence
    * @return {Node} node
@@ -941,7 +932,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
    */
   function parseAddSubtract(state) {
     let node, name, fn, params;
-    node = parseMultiplyDivideModulusPercentage(state);
+    node = parseMultiplyDivideModulus(state);
     const operators = {
       '+': 'add',
       '-': 'subtract'
@@ -950,7 +941,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
       name = state.token;
       fn = operators[name];
       getTokenSkipNewline(state);
-      const rightNode = parseMultiplyDivideModulusPercentage(state);
+      const rightNode = parseMultiplyDivideModulus(state);
       if (rightNode.isPercentage) {
         params = [node, new OperatorNode('*', 'multiply', [node, rightNode])];
       } else {
@@ -962,11 +953,11 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
   }
 
   /**
-   * multiply, divide, modulus, percentage
+   * multiply, divide, modulus
    * @return {Node} node
    * @private
    */
-  function parseMultiplyDivideModulusPercentage(state) {
+  function parseMultiplyDivideModulus(state) {
     let node, last, name, fn;
     node = parseImplicitMultiplication(state);
     last = node;
@@ -984,23 +975,8 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
         name = state.token;
         fn = operators[name];
         getTokenSkipNewline(state);
-        if (name === '%' && state.tokenType === TOKENTYPE.DELIMITER && state.token !== '(') {
-          // If the expression contains only %, then treat that as /100
-          if (state.token !== '' && operators[state.token]) {
-            const left = new OperatorNode('/', 'divide', [node, new ConstantNode(100)], false, true);
-            name = state.token;
-            fn = operators[name];
-            getTokenSkipNewline(state);
-            last = parseImplicitMultiplication(state);
-            node = new OperatorNode(name, fn, [left, last]);
-          } else {
-            node = new OperatorNode('/', 'divide', [node, new ConstantNode(100)], false, true);
-          }
-          // return node
-        } else {
-          last = parseImplicitMultiplication(state);
-          node = new OperatorNode(name, fn, [node, last]);
-        }
+        last = parseImplicitMultiplication(state);
+        node = new OperatorNode(name, fn, [node, last]);
       } else {
         break;
       }
@@ -1043,7 +1019,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
    * @private
    */
   function parseRule2(state) {
-    let node = parseUnary(state);
+    let node = parseUnaryPercentage(state);
     let last = node;
     const tokenStates = [];
     while (true) {
@@ -1065,7 +1041,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
             // Rewind once and build the "number / number" node; the symbol will be consumed later
             (0, _extends2.default)(state, tokenStates.pop());
             tokenStates.pop();
-            last = parseUnary(state);
+            last = parseUnaryPercentage(state);
             node = new OperatorNode('/', 'divide', [node, last]);
           } else {
             // Not a match, so rewind
@@ -1085,6 +1061,34 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
     return node;
   }
 
+  /**
+   * Unary percentage operator (treated as `value / 100`)
+   * @return {Node} node
+   * @private
+   */
+  function parseUnaryPercentage(state) {
+    let node = parseUnary(state);
+    if (state.token === '%') {
+      const previousState = (0, _extends2.default)({}, state);
+      getTokenSkipNewline(state);
+      // We need to decide if this is a unary percentage % or binary modulo %
+      // So we attempt to parse a unary expression at this point.
+      // If it fails, then the only possibility is that this is a unary percentage.
+      // If it succeeds, then we presume that this must be binary modulo, since the
+      // only things that parseUnary can handle are _higher_ precedence than unary %.
+      try {
+        parseUnary(state);
+        // Not sure if we could somehow use the result of that parseUnary? Without
+        // further analysis/testing, safer just to discard and let the parse proceed
+        (0, _extends2.default)(state, previousState);
+      } catch (_unused) {
+        // Not seeing a term at this point, so was a unary %
+        node = new OperatorNode('/', 'divide', [node, new ConstantNode(100)], false, true);
+      }
+    }
+    return node;
+  }
+
   /**
    * Unary plus and minus, and logical and bitwise not
    * @return {Node} node

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

@@ -23,7 +23,8 @@ const createSylvester = exports.createSylvester = /* #__PURE__ */(0, _factory.fa
     subtract,
     identity,
     lusolve,
-    abs
+    abs,
+    config
   } = _ref;
   /**
    *
@@ -94,7 +95,7 @@ const createSylvester = exports.createSylvester = /* #__PURE__ */(0, _factory.fa
     const vc = (a, b) => concat(a, b, 0);
     for (let k = 0; k < n; k++) {
       if (k < n - 1 && abs(subset(G, index(k + 1, k))) > 1e-5) {
-        let RHS = vc(subset(D, index(all, k)), subset(D, index(all, k + 1)));
+        let RHS = vc(subset(D, index(all, [k])), subset(D, index(all, [k + 1])));
         for (let 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)))));
         }
@@ -104,11 +105,11 @@ const createSylvester = exports.createSylvester = /* #__PURE__ */(0, _factory.fa
         const gmm = multiply(identity(m), multiply(-1, subset(G, index(k + 1, k + 1))));
         const LHS = vc(hc(add(F, gkk), gmk), hc(gkm, add(F, gmm)));
         const yAux = lusolve(LHS, RHS);
-        y[k] = yAux.subset(index(range(0, m), 0));
-        y[k + 1] = yAux.subset(index(range(m, 2 * m), 0));
+        y[k] = yAux.subset(index(range(0, m), [0]));
+        y[k + 1] = yAux.subset(index(range(m, 2 * m), [0]));
         k++;
       } else {
-        let RHS = subset(D, index(all, k));
+        let RHS = subset(D, index(all, [k]));
         for (let j = 0; j < k; j++) {
           RHS = add(RHS, multiply(y[j], subset(G, index(j, k))));
         }

package/lib/cjs/function/logical/nullish.js

@@ -68,8 +68,8 @@ const createNullish = exports.createNullish = /* #__PURE__ */(0, _factory.factor
     'undefined, any': (_x, y) => y,
     // SparseMatrix-first with collection RHS: enforce exact shape match
     'SparseMatrix, Array | Matrix': (x, y) => {
-      const sx = flatten(size(x).valueOf()); // work around #3529/#3530
-      const sy = flatten(size(y).valueOf());
+      const sx = size(x);
+      const sy = size(y);
       if (deepEqual(sx, sy)) return x;
       throw new _DimensionError.DimensionError(sx, sy);
     },

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

@@ -58,8 +58,9 @@ const createColumn = exports.createColumn = /* #__PURE__ */(0, _factory.factory)
     }
     (0, _array.validateIndex)(column, value.size()[1]);
     const rowRange = range(0, value.size()[0]);
-    const index = new Index(rowRange, column);
+    const index = new Index(rowRange, [column]);
     const result = value.subset(index);
+    // once config.legacySubset just return result
     return (0, _is.isMatrix)(result) ? result : matrix([[result]]);
   }
 });