mathjs 14.6.0 → 14.8.0

package/lib/cjs/core/function/config.js

@@ -23,7 +23,6 @@ function configFactory(config, emit) {
    *
    * Examples:
    *
-   *
    *     import { create, all } from 'mathjs'
    *
    *     // create a mathjs instance

package/lib/cjs/entry/dependenciesAny/dependenciesNullish.generated.js

@@ -0,0 +1,25 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.nullishDependencies = void 0;
+var _dependenciesDeepEqualGenerated = require("./dependenciesDeepEqual.generated.js");
+var _dependenciesFlattenGenerated = require("./dependenciesFlatten.generated.js");
+var _dependenciesMatrixGenerated = require("./dependenciesMatrix.generated.js");
+var _dependenciesSizeGenerated = require("./dependenciesSize.generated.js");
+var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
+var _factoriesAny = require("../../factoriesAny.js");
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+
+const nullishDependencies = exports.nullishDependencies = {
+  deepEqualDependencies: _dependenciesDeepEqualGenerated.deepEqualDependencies,
+  flattenDependencies: _dependenciesFlattenGenerated.flattenDependencies,
+  matrixDependencies: _dependenciesMatrixGenerated.matrixDependencies,
+  sizeDependencies: _dependenciesSizeGenerated.sizeDependencies,
+  typedDependencies: _dependenciesTypedGenerated.typedDependencies,
+  createNullish: _factoriesAny.createNullish
+};

package/lib/cjs/entry/dependenciesAny/dependenciesNullishTransform.generated.js

@@ -0,0 +1,25 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.nullishTransformDependencies = void 0;
+var _dependenciesDeepEqualGenerated = require("./dependenciesDeepEqual.generated.js");
+var _dependenciesFlattenGenerated = require("./dependenciesFlatten.generated.js");
+var _dependenciesMatrixGenerated = require("./dependenciesMatrix.generated.js");
+var _dependenciesSizeGenerated = require("./dependenciesSize.generated.js");
+var _dependenciesTypedGenerated = require("./dependenciesTyped.generated.js");
+var _factoriesAny = require("../../factoriesAny.js");
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+
+const nullishTransformDependencies = exports.nullishTransformDependencies = {
+  deepEqualDependencies: _dependenciesDeepEqualGenerated.deepEqualDependencies,
+  flattenDependencies: _dependenciesFlattenGenerated.flattenDependencies,
+  matrixDependencies: _dependenciesMatrixGenerated.matrixDependencies,
+  sizeDependencies: _dependenciesSizeGenerated.sizeDependencies,
+  typedDependencies: _dependenciesTypedGenerated.typedDependencies,
+  createNullishTransform: _factoriesAny.createNullishTransform
+};

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

@@ -1443,6 +1443,18 @@ Object.defineProperty(exports, "nullDependencies", {
     return _dependenciesNullGenerated.nullDependencies;
   }
 });
+Object.defineProperty(exports, "nullishDependencies", {
+  enumerable: true,
+  get: function () {
+    return _dependenciesNullishGenerated.nullishDependencies;
+  }
+});
+Object.defineProperty(exports, "nullishTransformDependencies", {
+  enumerable: true,
+  get: function () {
+    return _dependenciesNullishTransformGenerated.nullishTransformDependencies;
+  }
+});
 Object.defineProperty(exports, "numberDependencies", {
   enumerable: true,
   get: function () {
@@ -2416,6 +2428,8 @@ var _dependenciesNthRootGenerated = require("./dependenciesAny/dependenciesNthRo
 var _dependenciesNthRootsGenerated = require("./dependenciesAny/dependenciesNthRoots.generated.js");
 var _dependenciesNuclearMagnetonGenerated = require("./dependenciesAny/dependenciesNuclearMagneton.generated.js");
 var _dependenciesNullGenerated = require("./dependenciesAny/dependenciesNull.generated.js");
+var _dependenciesNullishGenerated = require("./dependenciesAny/dependenciesNullish.generated.js");
+var _dependenciesNullishTransformGenerated = require("./dependenciesAny/dependenciesNullishTransform.generated.js");
 var _dependenciesNumberGenerated = require("./dependenciesAny/dependenciesNumber.generated.js");
 var _dependenciesNumericGenerated = require("./dependenciesAny/dependenciesNumeric.generated.js");
 var _dependenciesObjectNodeGenerated = require("./dependenciesAny/dependenciesObjectNode.generated.js");

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

@@ -416,6 +416,7 @@ const help = exports.help = (0, _factoriesAny.createHelp)({
   min: _pureFunctionsAnyGenerated.min,
   mod: _pureFunctionsAnyGenerated.mod,
   nthRoots: _pureFunctionsAnyGenerated.nthRoots,
+  nullish: _pureFunctionsAnyGenerated.nullish,
   or: _pureFunctionsAnyGenerated.or,
   partitionSelect: _pureFunctionsAnyGenerated.partitionSelect,
   qr: _pureFunctionsAnyGenerated.qr,
@@ -631,6 +632,13 @@ const help = exports.help = (0, _factoriesAny.createHelp)({
     smaller: _pureFunctionsAnyGenerated.smaller,
     typed: _pureFunctionsAnyGenerated.typed
   }),
+  nullish: (0, _factoriesAny.createNullishTransform)({
+    deepEqual: _pureFunctionsAnyGenerated.deepEqual,
+    flatten: _pureFunctionsAnyGenerated.flatten,
+    matrix: _pureFunctionsAnyGenerated.matrix,
+    size: _pureFunctionsAnyGenerated.size,
+    typed: _pureFunctionsAnyGenerated.typed
+  }),
   subset: (0, _factoriesAny.createSubsetTransform)({
     add: _pureFunctionsAnyGenerated.add,
     matrix: _pureFunctionsAnyGenerated.matrix,

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

@@ -4,9 +4,9 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.dotPow = exports.dotMultiply = exports.dotDivide = exports.dot = exports.divideScalar = exports.divide = exports.distance = exports.diff = exports.diag = exports.deuteronMass = exports.det = exports.deepEqual = exports.cumsum = exports.cube = exports.ctranspose = exports.csch = exports.csc = exports.cross = exports.createUnit = exports.count = exports.coulombConstant = exports.coulomb = exports.coth = exports.cot = exports.cosh = exports.cos = exports.corr = exports.conj = exports.conductanceQuantum = exports.concat = exports.composition = exports.complex = exports.compareText = exports.compareNatural = exports.compare = exports.combinationsWithRep = exports.combinations = exports.column = exports.clone = exports.classicalElectronRadius = exports.ceil = exports.cbrt = exports.catalan = exports.boolean = exports.boltzmann = exports.bohrRadius = exports.bohrMagneton = exports.bitXor = exports.bitOr = exports.bitNot = exports.bitAnd = exports.bin = exports.bignumber = exports.bigint = exports.bellNumbers = exports.avogadro = exports.atomicMass = exports.atanh = exports.atan2 = exports.atan = exports.asinh = exports.asin = exports.asech = exports.asec = exports.arg = exports.apply = exports.and = exports.addScalar = exports.add = exports.acsch = exports.acsc = exports.acoth = exports.acot = exports.acosh = exports.acos = exports.abs = exports._true = exports._null = exports._false = exports._NaN = exports._Infinity = exports.Unit = exports.SparseMatrix = exports.Spa = exports.SQRT2 = exports.SQRT1_2 = exports.ResultSet = exports.Range = exports.Matrix = exports.LOG2E = exports.LOG10E = exports.LN2 = exports.LN10 = exports.Index = exports.ImmutableDenseMatrix = exports.Fraction = exports.FibonacciHeap = exports.DenseMatrix = exports.Complex = exports.BigNumber = void 0;
-exports.number = exports.nuclearMagneton = exports.nthRoots = exports.nthRoot = exports.not = exports.norm = exports.neutronMass = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.molarVolume = exports.molarPlanckConstant = exports.molarMassC12 = exports.molarMass = exports.mode = exports.mod = exports.min = exports.median = exports.mean = exports.max = exports.matrixFromRows = exports.matrixFromFunction = exports.matrixFromColumns = exports.matrix = exports.mapSlices = exports.map = exports.magneticFluxQuantum = exports.magneticConstant = exports.mad = exports.lyap = exports.lusolve = exports.lup = exports.lsolveAll = exports.lsolve = exports.loschmidt = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.kron = exports.klitzing = exports.kldivergence = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.invmod = exports.inverseConductanceQuantum = exports.inv = exports.intersect = exports.index = exports.im = exports.ifft = exports.identity = exports.i = exports.hypot = exports.hex = exports.hasNumericValue = exports.hartreeEnergy = exports.gravity = exports.gravitationConstant = exports.getMatrixDataType = exports.gcd = exports.gasConstant = exports.gamma = exports.freqz = exports.fraction = exports.format = exports.forEach = exports.floor = exports.flatten = exports.fix = exports.firstRadiation = exports.fineStructure = exports.filter = exports.fft = exports.fermiCoupling = exports.faraday = exports.factorial = exports.expm1 = exports.expm = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.elementaryCharge = exports.electronMass = exports.electricConstant = exports.eigs = exports.efimovFactor = exports.e = void 0;
-exports.wienDisplacement = exports.weakMixingAngle = exports.version = exports.variance = exports.vacuumImpedance = exports.usolveAll = exports.usolve = exports.unit = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.transpose = exports.trace = exports.toBest = exports.to = exports.thomsonCrossSection = exports.tau = exports.tanh = exports.tan = exports.sylvester = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.stefanBoltzmann = exports.std = exports.squeeze = exports.square = exports.sqrtm = exports.sqrt = exports.splitUnit = exports.speedOfLight = exports.sparse = exports.sort = exports.solveODE = exports.smallerEq = exports.smaller = exports.slu = exports.size = exports.sinh = exports.sin = exports.sign = exports.setUnion = exports.setSymDifference = exports.setSize = exports.setPowerset = exports.setMultiplicity = exports.setIsSubset = exports.setIntersect = exports.setDistinct = exports.setDifference = exports.setCartesian = exports.secondRadiation = exports.sech = exports.sec = exports.schur = exports.sackurTetrode = exports.rydberg = exports.row = exports.round = exports.rotationMatrix = exports.rotate = exports.rightLogShift = exports.rightArithShift = exports.resize = exports.reshape = exports.replacer = exports.reducedPlanckConstant = exports.re = exports.range = exports.randomInt = exports.random = exports.quantumOfCirculation = exports.quantileSeq = exports.qr = exports.protonMass = exports.prod = exports.print = exports.pow = exports.polynomialRoot = exports.planckTime = exports.planckTemperature = exports.planckMass = exports.planckLength = exports.planckConstant = exports.planckCharge = exports.pinv = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.ones = exports.oct = exports.numeric = void 0;
-exports.zpk2tf = exports.zeta = exports.zeros = exports.xor = exports.xgcd = void 0;
+exports.nullish = exports.nuclearMagneton = exports.nthRoots = exports.nthRoot = exports.not = exports.norm = exports.neutronMass = exports.multiplyScalar = exports.multiply = exports.multinomial = exports.molarVolume = exports.molarPlanckConstant = exports.molarMassC12 = exports.molarMass = exports.mode = exports.mod = exports.min = exports.median = exports.mean = exports.max = exports.matrixFromRows = exports.matrixFromFunction = exports.matrixFromColumns = exports.matrix = exports.mapSlices = exports.map = exports.magneticFluxQuantum = exports.magneticConstant = exports.mad = exports.lyap = exports.lusolve = exports.lup = exports.lsolveAll = exports.lsolve = exports.loschmidt = exports.log2 = exports.log1p = exports.log10 = exports.log = exports.lgamma = exports.leftShift = exports.lcm = exports.largerEq = exports.larger = exports.kron = exports.klitzing = exports.kldivergence = exports.isZero = exports.isPrime = exports.isPositive = exports.isNumeric = exports.isNegative = exports.isNaN = exports.isInteger = exports.invmod = exports.inverseConductanceQuantum = exports.inv = exports.intersect = exports.index = exports.im = exports.ifft = exports.identity = exports.i = exports.hypot = exports.hex = exports.hasNumericValue = exports.hartreeEnergy = exports.gravity = exports.gravitationConstant = exports.getMatrixDataType = exports.gcd = exports.gasConstant = exports.gamma = exports.freqz = exports.fraction = exports.format = exports.forEach = exports.floor = exports.flatten = exports.fix = exports.firstRadiation = exports.fineStructure = exports.filter = exports.fft = exports.fermiCoupling = exports.faraday = exports.factorial = exports.expm1 = exports.expm = exports.exp = exports.erf = exports.equalText = exports.equalScalar = exports.equal = exports.elementaryCharge = exports.electronMass = exports.electricConstant = exports.eigs = exports.efimovFactor = exports.e = void 0;
+exports.weakMixingAngle = exports.version = exports.variance = exports.vacuumImpedance = exports.usolveAll = exports.usolve = exports.unit = exports.unequal = exports.unaryPlus = exports.unaryMinus = exports.typed = exports.typeOf = exports.transpose = exports.trace = exports.toBest = exports.to = exports.thomsonCrossSection = exports.tau = exports.tanh = exports.tan = exports.sylvester = exports.sum = exports.subtractScalar = exports.subtract = exports.subset = exports.string = exports.stirlingS2 = exports.stefanBoltzmann = exports.std = exports.squeeze = exports.square = exports.sqrtm = exports.sqrt = exports.splitUnit = exports.speedOfLight = exports.sparse = exports.sort = exports.solveODE = exports.smallerEq = exports.smaller = exports.slu = exports.size = exports.sinh = exports.sin = exports.sign = exports.setUnion = exports.setSymDifference = exports.setSize = exports.setPowerset = exports.setMultiplicity = exports.setIsSubset = exports.setIntersect = exports.setDistinct = exports.setDifference = exports.setCartesian = exports.secondRadiation = exports.sech = exports.sec = exports.schur = exports.sackurTetrode = exports.rydberg = exports.row = exports.round = exports.rotationMatrix = exports.rotate = exports.rightLogShift = exports.rightArithShift = exports.resize = exports.reshape = exports.replacer = exports.reducedPlanckConstant = exports.re = exports.range = exports.randomInt = exports.random = exports.quantumOfCirculation = exports.quantileSeq = exports.qr = exports.protonMass = exports.prod = exports.print = exports.pow = exports.polynomialRoot = exports.planckTime = exports.planckTemperature = exports.planckMass = exports.planckLength = exports.planckConstant = exports.planckCharge = exports.pinv = exports.pickRandom = exports.pi = exports.phi = exports.permutations = exports.partitionSelect = exports.or = exports.ones = exports.oct = exports.numeric = exports.number = void 0;
+exports.zpk2tf = exports.zeta = exports.zeros = exports.xor = exports.xgcd = exports.wienDisplacement = void 0;
 var _configReadonly = require("./configReadonly.js");
 var _factoriesAny = require("../factoriesAny.js");
 /**
@@ -779,6 +779,13 @@ const nthRoots = exports.nthRoots = /* #__PURE__ */(0, _factoriesAny.createNthRo
   divideScalar,
   typed
 });
+const nullish = exports.nullish = /* #__PURE__ */(0, _factoriesAny.createNullish)({
+  deepEqual,
+  flatten,
+  matrix,
+  size,
+  typed
+});
 const or = exports.or = /* #__PURE__ */(0, _factoriesAny.createOr)({
   DenseMatrix,
   concat,

package/lib/cjs/expression/Parser.js

@@ -5,6 +5,7 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.createParserClass = void 0;
 var _factory = require("../utils/factory.js");
+var _is = require("../utils/is.js");
 var _map = require("../utils/map.js");
 const name = 'Parser';
 const dependencies = ['evaluate', 'parse'];
@@ -160,7 +161,40 @@ const createParserClass = exports.createParserClass = /* #__PURE__ */(0, _factor
   Parser.prototype.clear = function () {
     this.scope.clear();
   };
+  Parser.prototype.toJSON = function () {
+    const json = {
+      mathjs: 'Parser',
+      variables: {},
+      functions: {}
+    };
+    for (const [name, value] of this.scope) {
+      if ((0, _is.isFunction)(value)) {
+        if (!isExpressionFunction(value)) {
+          throw new Error(`Cannot serialize external function ${name}`);
+        }
+        json.functions[name] = `${value.syntax} = ${value.expr}`;
+      } else {
+        json.variables[name] = value;
+      }
+    }
+    return json;
+  };
+  Parser.fromJSON = function (json) {
+    const parser = new Parser();
+    Object.entries(json.variables).forEach(_ref2 => {
+      let [name, value] = _ref2;
+      return parser.set(name, value);
+    });
+    Object.entries(json.functions).forEach(_ref3 => {
+      let [_name, fn] = _ref3;
+      return parser.evaluate(fn);
+    });
+    return parser;
+  };
   return Parser;
 }, {
   isClass: true
-});
+});
+function isExpressionFunction(value) {
+  return typeof value === 'function' && typeof value.syntax === 'string' && typeof value.expr === 'string';
+}

package/lib/cjs/expression/embeddedDocs/embeddedDocs.js

@@ -118,6 +118,7 @@ var _distance = require("./function/geometry/distance.js");
 var _intersect = require("./function/geometry/intersect.js");
 var _and = require("./function/logical/and.js");
 var _not = require("./function/logical/not.js");
+var _nullish = require("./function/logical/nullish.js");
 var _or = require("./function/logical/or.js");
 var _xor = require("./function/logical/xor.js");
 var _mapSlices = require("./function/matrix/mapSlices.js");
@@ -587,6 +588,7 @@ const embeddedDocs = exports.embeddedDocs = {
   // functions - logical
   and: _and.andDocs,
   not: _not.notDocs,
+  nullish: _nullish.nullishDocs,
   or: _or.orDocs,
   xor: _xor.xorDocs,
   // functions - matrix

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

@@ -0,0 +1,14 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.nullishDocs = void 0;
+const nullishDocs = exports.nullishDocs = {
+  name: 'nullish',
+  category: 'Logical',
+  syntax: ['x ?? y', 'nullish(x, y)'],
+  description: 'Nullish coalescing operator. Returns the right-hand operand when the left-hand operand is null or undefined, and otherwise returns the left-hand operand.',
+  examples: ['null ?? 42', 'undefined ?? 42', '0 ?? 42', 'false ?? 42', 'null ?? undefined ?? 42'],
+  seealso: ['and', 'or', 'not']
+};

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

@@ -105,7 +105,8 @@ const createFunctionAssignmentNode = exports.createFunctionAssignmentNode = /* #
       });
 
       // compile the function expression with the child args
-      const evalExpr = this.expr._compile(math, childArgNames);
+      const expr = this.expr;
+      const evalExpr = expr._compile(math, childArgNames);
       const name = this.name;
       const params = this.params;
       const signature = (0, _array.join)(this.types, ',');
@@ -121,6 +122,7 @@ const createFunctionAssignmentNode = exports.createFunctionAssignmentNode = /* #
         };
         const fn = typed(name, signatures);
         fn.syntax = syntax;
+        fn.expr = expr.toString();
         scope.set(name, fn);
         return fn;
       };

package/lib/cjs/expression/operators.js

@@ -232,6 +232,13 @@ const properties = exports.properties = [{
     associativity: 'right',
     associativeWith: []
   }
+}, {
+  // nullish coalescing
+  'OperatorNode:nullish': {
+    op: '??',
+    associativity: 'left',
+    associativeWith: []
+  }
 }, {
   // factorial
   'OperatorNode:factorial': {

package/lib/cjs/expression/parse.js

@@ -54,7 +54,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
    *     node1.compile().evaluate() // 5
    *
    *     let scope = {a:3, b:4}
-   *     const node2 = math.parse('a * b') // 12
+   *     const node2 = math.parse('a * b')
    *     node2.evaluate(scope) // 12
    *     const code2 = node2.compile()
    *     code2.evaluate(scope) // 12
@@ -136,6 +136,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
     '=': true,
     ':': true,
     '?': true,
+    '??': true,
     '==': true,
     '!=': true,
     '<': true,
@@ -1115,7 +1116,7 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
    */
   function parsePow(state) {
     let node, name, fn, params;
-    node = parseLeftHandOperators(state);
+    node = parseNullishCoalescing(state);
     if (state.token === '^' || state.token === '.^') {
       name = state.token;
       fn = name === '^' ? 'pow' : 'dotPow';
@@ -1126,6 +1127,21 @@ const createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(n
     return node;
   }
 
+  /**
+   * nullish coalescing operator
+   * @return {Node} node
+   * @private
+   */
+  function parseNullishCoalescing(state) {
+    let node = parseLeftHandOperators(state);
+    while (state.token === '??') {
+      // eslint-disable-line no-unmodified-loop-condition
+      getTokenSkipNewline(state);
+      node = new OperatorNode('??', 'nullish', [node, parseLeftHandOperators(state)]);
+    }
+    return node;
+  }
+
   /**
    * Left hand operators: factorial x!, ctranspose x'
    * @return {Node} node

package/lib/cjs/expression/transform/nullish.transform.js

@@ -0,0 +1,43 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createNullishTransform = void 0;
+var _nullish = require("../../function/logical/nullish.js");
+var _factory = require("../../utils/factory.js");
+var _is = require("../../utils/is.js");
+const name = 'nullish';
+const dependencies = ['typed', 'matrix', 'size', 'flatten', 'deepEqual'];
+const createNullishTransform = exports.createNullishTransform = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed,
+    matrix,
+    size,
+    flatten,
+    deepEqual
+  } = _ref;
+  const nullish = (0, _nullish.createNullish)({
+    typed,
+    matrix,
+    size,
+    flatten,
+    deepEqual
+  });
+  function nullishTransform(args, math, scope) {
+    const left = args[0].compile().evaluate(scope);
+
+    // If left is not a collection and not nullish, short-circuit and return it
+    if (!(0, _is.isCollection)(left) && left != null && left !== undefined) {
+      return left;
+    }
+
+    // Otherwise evaluate right and apply full nullish semantics (incl. element-wise)
+    const right = args[1].compile().evaluate(scope);
+    return nullish(left, right);
+  }
+  nullishTransform.rawArgs = true;
+  return nullishTransform;
+}, {
+  isTransformFunction: true
+});

package/lib/cjs/factoriesAny.js

@@ -1341,6 +1341,18 @@ Object.defineProperty(exports, "createNull", {
     return _constants.createNull;
   }
 });
+Object.defineProperty(exports, "createNullish", {
+  enumerable: true,
+  get: function () {
+    return _nullish.createNullish;
+  }
+});
+Object.defineProperty(exports, "createNullishTransform", {
+  enumerable: true,
+  get: function () {
+    return _nullishTransform.createNullishTransform;
+  }
+});
 Object.defineProperty(exports, "createNumber", {
   enumerable: true,
   get: function () {
@@ -2255,6 +2267,7 @@ var _conj = require("./function/complex/conj.js");
 var _im = require("./function/complex/im.js");
 var _re = require("./function/complex/re.js");
 var _not = require("./function/logical/not.js");
+var _nullish = require("./function/logical/nullish.js");
 var _or = require("./function/logical/or.js");
 var _xor = require("./function/logical/xor.js");
 var _concat = require("./function/matrix/concat.js");
@@ -2481,5 +2494,6 @@ var _varianceTransform = require("./expression/transform/variance.transform.js")
 var _printTransform = require("./expression/transform/print.transform.js");
 var _andTransform = require("./expression/transform/and.transform.js");
 var _orTransform = require("./expression/transform/or.transform.js");
+var _nullishTransform = require("./expression/transform/nullish.transform.js");
 var _bitAndTransform = require("./expression/transform/bitAnd.transform.js");
 var _bitOrTransform = require("./expression/transform/bitOr.transform.js");

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

@@ -0,0 +1,87 @@
+"use strict";
+
+Object.defineProperty(exports, "__esModule", {
+  value: true
+});
+exports.createNullish = void 0;
+var _factory = require("../../utils/factory.js");
+var _matAlgo03xDSf = require("../../type/matrix/utils/matAlgo03xDSf.js");
+var _matAlgo14xDs = require("../../type/matrix/utils/matAlgo14xDs.js");
+var _matAlgo13xDD = require("../../type/matrix/utils/matAlgo13xDD.js");
+var _DimensionError = require("../../error/DimensionError.js");
+const name = 'nullish';
+const dependencies = ['typed', 'matrix', 'size', 'flatten', 'deepEqual'];
+const createNullish = exports.createNullish = /* #__PURE__ */(0, _factory.factory)(name, dependencies, _ref => {
+  let {
+    typed,
+    matrix,
+    size,
+    flatten,
+    deepEqual
+  } = _ref;
+  const matAlgo03xDSf = (0, _matAlgo03xDSf.createMatAlgo03xDSf)({
+    typed
+  });
+  const matAlgo14xDs = (0, _matAlgo14xDs.createMatAlgo14xDs)({
+    typed
+  });
+  const matAlgo13xDD = (0, _matAlgo13xDD.createMatAlgo13xDD)({
+    typed
+  });
+
+  /**
+   * Nullish coalescing operator (??). Returns the right-hand side operand
+   * when the left-hand side operand is null or undefined, and otherwise
+   * returns the left-hand side operand.
+   *
+   * For matrices, the function is evaluated element wise.
+   *
+   * Syntax:
+   *
+   *    math.nullish(x, y)
+   *
+   * Examples:
+   *
+   *    math.nullish(null, 42)        // returns 42
+   *    math.nullish(undefined, 42)   // returns 42
+   *    math.nullish(0, 42)           // returns 0
+   *    math.nullish(false, 42)       // returns false
+   *    math.nullish('', 42)          // returns ''
+   *
+   *    // Object property access with fallback
+   *    const obj = {foo: 7, bar: 3}
+   *    math.nullish(obj.baz, 0)      // returns 0
+   *
+   * See also:
+   *
+   *    and, or, not
+   *
+   * @param  {*} x First value to check
+   * @param  {*} y Fallback value
+   * @return {*} Returns y when x is null or undefined, otherwise returns x
+   */
+
+  return typed(name, {
+    // Scalar and SparseMatrix-first short-circuit handlers
+    'number|bigint|Complex|BigNumber|Fraction|Unit|string|boolean|SparseMatrix, any': (x, _y) => x,
+    'null, any': (_x, y) => y,
+    '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());
+      if (deepEqual(sx, sy)) return x;
+      throw new _DimensionError.DimensionError(sx, sy);
+    },
+    // DenseMatrix-first handlers (no broadcasting between collections)
+    'DenseMatrix, DenseMatrix': typed.referToSelf(self => (x, y) => matAlgo13xDD(x, y, self)),
+    'DenseMatrix, SparseMatrix': typed.referToSelf(self => (x, y) => matAlgo03xDSf(x, y, self, false)),
+    'DenseMatrix, Array': typed.referToSelf(self => (x, y) => matAlgo13xDD(x, matrix(y), self)),
+    'DenseMatrix, any': typed.referToSelf(self => (x, y) => matAlgo14xDs(x, y, self, false)),
+    // Array-first handlers (bridge via matrix() where needed)
+    'Array, Array': typed.referToSelf(self => (x, y) => matAlgo13xDD(matrix(x), matrix(y), self).valueOf()),
+    'Array, DenseMatrix': typed.referToSelf(self => (x, y) => matAlgo13xDD(matrix(x), y, self)),
+    'Array, SparseMatrix': typed.referToSelf(self => (x, y) => matAlgo03xDSf(matrix(x), y, self, false)),
+    'Array, any': typed.referToSelf(self => (x, y) => matAlgo14xDs(matrix(x), y, self, false).valueOf())
+  });
+});

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

@@ -22,7 +22,7 @@ const createRange = exports.createRange = /* #__PURE__ */(0, _factory.factory)(n
     isPositive
   } = _ref;
   /**
-   * Create an array from a range.
+   * Create a matrix or array containing a range of values.
    * By default, the range end is excluded. This can be customized by providing
    * an extra parameter `includeEnd`.
    *
@@ -50,10 +50,12 @@ const createRange = exports.createRange = /* #__PURE__ */(0, _factory.factory)(n
    * - `includeEnd: boolean`
    *   Option to specify whether to include the end or not. False by default.
    *
-   * Note that the return type of the range is taken from the type of
-   * the start/end. If only one these is a built-in `number` type, it will
-   * be promoted to the type of the other endpoint. However, in the case of
-   * Unit values, both endpoints must have compatible units, and the return
+   * The function returns a `DenseMatrix` when the library is configured with
+   * `config = { matrix: 'Matrix' }, and returns an Array otherwise.
+   * Note that the type of the returned values is taken from the type of the
+   * provided start/end value. If only one of these is a built-in `number` type,
+   * it will be promoted to the type of the other endpoint. However, in the case
+   * of Unit values, both endpoints must have compatible units, and the return
    * value will have compatible units as well.
    *
    * Examples:

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 14.6.0
- * @date    2025-07-25
+ * @version 14.8.0
+ * @date    2025-09-24
  *
  * @license
  * Copyright (C) 2013-2025 Jos de Jong <wjosdejong@gmail.com>

package/lib/cjs/type/matrix/DenseMatrix.js

@@ -679,24 +679,28 @@ const createDenseMatrixClass = exports.createDenseMatrixClass = /* #__PURE__ */(
       }
       return;
     }
-    const index = [];
-    const recurse = function* (value, depth) {
-      if (depth < maxDepth) {
-        for (let i = 0; i < value.length; i++) {
-          index[depth] = i;
-          yield* recurse(value[i], depth + 1);
-        }
-      } else {
-        for (let i = 0; i < value.length; i++) {
-          index[depth] = i;
-          yield {
-            value: value[i],
-            index: index.slice()
-          };
-        }
+
+    // Multi-dimensional matrix: iterate over all elements
+    const index = Array(maxDepth + 1).fill(0);
+    const totalElements = this._size.reduce((a, b) => a * b, 1);
+    for (let count = 0; count < totalElements; count++) {
+      // Traverse to the current element using indices
+      let current = this._data;
+      for (let d = 0; d < maxDepth; d++) {
+        current = current[index[d]];
       }
-    };
-    yield* recurse(this._data, 0);
+      yield {
+        value: current[index[maxDepth]],
+        index: index.slice()
+      };
+
+      // Increment indices for next element
+      for (let d = maxDepth; d >= 0; d--) {
+        index[d]++;
+        if (index[d] < this._size[d]) break;
+        index[d] = 0;
+      }
+    }
   };
 
   /**

package/lib/cjs/type/string.js

@@ -23,13 +23,13 @@ const createString = exports.createString = /* #__PURE__ */(0, _factory.factory)
    *
    * Examples:
    *
-   *    math.string(4.2)               // returns string '4.2'
-   *    math.string(math.complex(3, 2) // returns string '3 + 2i'
+   *    math.string(4.2)                // returns string '4.2'
+   *    math.string(math.complex(3, 2)) // returns string '3 + 2i'
    *
    *    const u = math.unit(5, 'km')
-   *    math.string(u.to('m'))         // returns string '5000 m'
+   *    math.string(u.to('m'))          // returns string '5000 m'
    *
-   *    math.string([true, false])     // returns ['true', 'false']
+   *    math.string([true, false])      // returns ['true', 'false']
    *
    * See also:
    *

package/lib/cjs/version.js

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

package/lib/esm/core/function/config.js

@@ -16,7 +16,6 @@ export function configFactory(config, emit) {
    *
    * Examples:
    *
-   *
    *     import { create, all } from 'mathjs'
    *
    *     // create a mathjs instance

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

@@ -0,0 +1,18 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { deepEqualDependencies } from './dependenciesDeepEqual.generated.js';
+import { flattenDependencies } from './dependenciesFlatten.generated.js';
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { sizeDependencies } from './dependenciesSize.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createNullish } from '../../factoriesAny.js';
+export var nullishDependencies = {
+  deepEqualDependencies,
+  flattenDependencies,
+  matrixDependencies,
+  sizeDependencies,
+  typedDependencies,
+  createNullish
+};

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

@@ -0,0 +1,18 @@
+/**
+ * THIS FILE IS AUTO-GENERATED
+ * DON'T MAKE CHANGES HERE
+ */
+import { deepEqualDependencies } from './dependenciesDeepEqual.generated.js';
+import { flattenDependencies } from './dependenciesFlatten.generated.js';
+import { matrixDependencies } from './dependenciesMatrix.generated.js';
+import { sizeDependencies } from './dependenciesSize.generated.js';
+import { typedDependencies } from './dependenciesTyped.generated.js';
+import { createNullishTransform } from '../../factoriesAny.js';
+export var nullishTransformDependencies = {
+  deepEqualDependencies,
+  flattenDependencies,
+  matrixDependencies,
+  sizeDependencies,
+  typedDependencies,
+  createNullishTransform
+};