mathjs 12.3.1 → 12.4.0

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

@@ -601,6 +601,7 @@ var rightArithShift = exports.rightArithShift = /* #__PURE__ */(0, _factoriesAny
 var round = exports.round = /* #__PURE__ */(0, _factoriesAny.createRound)({
   BigNumber: BigNumber,
   DenseMatrix: DenseMatrix,
+  config: _configReadonly.config,
   equalScalar: equalScalar,
   matrix: matrix,
   typed: typed,

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

@@ -13,6 +13,7 @@ var _inherits2 = _interopRequireDefault(require("@babel/runtime/helpers/inherits
 var _defineProperty2 = _interopRequireDefault(require("@babel/runtime/helpers/defineProperty"));
 var _is = require("../../utils/is.js");
 var _array = require("../../utils/array.js");
+var _scope = require("../../utils/scope.js");
 var _string = require("../../utils/string.js");
 var _customs = require("../../utils/customs.js");
 var _operators = require("../operators.js");
@@ -311,7 +312,7 @@ var createOperatorNode = exports.createOperatorNode = /* #__PURE__ */(0, _factor
           // "raw" evaluation
           var rawArgs = this.args;
           return function evalOperatorNode(scope, args, context) {
-            return fn(rawArgs, math, scope);
+            return fn(rawArgs, math, (0, _scope.createSubScope)(scope, args));
           };
         } else if (evalArgs.length === 1) {
           var evalArg0 = evalArgs[0];

package/lib/cjs/expression/parse.js

@@ -1416,8 +1416,10 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
           while (state.token === ';') {
             // eslint-disable-line no-unmodified-loop-condition
             getToken(state);
-            params[rows] = parseRow(state);
-            rows++;
+            if (state.token !== ']') {
+              params[rows] = parseRow(state);
+              rows++;
+            }
           }
           if (state.token !== ']') {
             throw createSyntaxError(state, 'End of matrix ] expected');
@@ -1465,8 +1467,10 @@ var createParse = exports.createParse = /* #__PURE__ */(0, _factory.factory)(nam
       getToken(state);
 
       // parse expression
-      params[len] = parseAssignment(state);
-      len++;
+      if (state.token !== ']' && state.token !== ';') {
+        params[len] = parseAssignment(state);
+        len++;
+      }
     }
     return new ArrayNode(params);
   }

package/lib/cjs/expression/transform/utils/compileInlineExpression.js

@@ -5,12 +5,12 @@ Object.defineProperty(exports, "__esModule", {
 });
 exports.compileInlineExpression = compileInlineExpression;
 var _is = require("../../../utils/is.js");
-var _scope = require("../../../utils/scope.js");
+var _map = require("../../../utils/map.js");
 /**
  * Compile an inline expression like "x > 0"
  * @param {Node} expression
  * @param {Object} math
- * @param {Object} scope
+ * @param {Map} scope
  * @return {function} Returns a function with one argument which fills in the
  *                    undefined variable (like "x") and evaluates the expression
  */
@@ -25,10 +25,11 @@ function compileInlineExpression(expression, math, scope) {
 
   // create a test function for this equation
   var name = symbol.name; // variable name
-  var subScope = (0, _scope.createSubScope)(scope);
+  var argsScope = new Map();
+  var subScope = new _map.PartitionedMap(scope, argsScope, new Set([name]));
   var eq = expression.compile();
   return function inlineExpression(x) {
-    subScope.set(name, x);
+    argsScope.set(name, x);
     return eq.evaluate(subScope);
   };
 }

package/lib/cjs/function/arithmetic/multiply.js

@@ -119,11 +119,11 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     // a dense
     var adata = a._data;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a.getDataType();
     // b dense
     var bdata = b._data;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b.getDataType();
     // rows & columns
     var alength = asize[0];
     var bcolumns = bsize[1];
@@ -136,7 +136,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var mf = multiplyScalar;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -163,7 +163,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     return a.createDenseMatrix({
       data: c,
       size: [bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -207,10 +207,10 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     // a dense
     var adata = a._data;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a.getDataType();
     // b dense
     var bdata = b._data;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b.getDataType();
     // rows & columns
     var arows = asize[0];
     var acolumns = asize[1];
@@ -223,7 +223,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var mf = multiplyScalar;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -252,7 +252,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     return a.createDenseMatrix({
       data: c,
       size: [arows],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -265,14 +265,15 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
    * @return {Matrix}             DenseMatrix    (MxC)
    */
   function _multiplyDenseMatrixDenseMatrix(a, b) {
+    // getDataType()
     // a dense
     var adata = a._data;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a.getDataType();
     // b dense
     var bdata = b._data;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b.getDataType();
     // rows & columns
     var arows = asize[0];
     var acolumns = asize[1];
@@ -286,7 +287,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var mf = multiplyScalar;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -320,7 +321,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     return a.createDenseMatrix({
       data: c,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -336,13 +337,13 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     // a dense
     var adata = a._data;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a.getDataType();
     // b sparse
     var bvalues = b._values;
     var bindex = b._index;
     var bptr = b._ptr;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b._data === undefined ? b._datatype : b.getDataType();
     // validate b matrix
     if (!bvalues) {
       throw new Error('Cannot multiply Dense Matrix times Pattern only Matrix');
@@ -363,7 +364,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var zero = 0;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -384,7 +385,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // loop b columns
@@ -448,14 +449,14 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var avalues = a._values;
     var aindex = a._index;
     var aptr = a._ptr;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // validate a matrix
     if (!avalues) {
       throw new Error('Cannot multiply Pattern only Matrix times Dense Matrix');
     }
     // b dense
     var bdata = b._data;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b.getDataType();
     // rows & columns
     var arows = a._size[0];
     var brows = b._size[0];
@@ -476,7 +477,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var zero = 0;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -529,13 +530,13 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     // update ptr
     cptr[1] = cindex.length;
 
-    // return sparse matrix
+    // matrix to return
     return a.createSparseMatrix({
       values: cvalues,
       index: cindex,
       ptr: cptr,
       size: [arows, 1],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -552,14 +553,14 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var avalues = a._values;
     var aindex = a._index;
     var aptr = a._ptr;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // validate a matrix
     if (!avalues) {
       throw new Error('Cannot multiply Pattern only Matrix times Dense Matrix');
     }
     // b dense
     var bdata = b._data;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b.getDataType();
     // rows & columns
     var arows = a._size[0];
     var brows = b._size[0];
@@ -577,7 +578,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var zero = 0;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -598,7 +599,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // workspace
@@ -665,12 +666,12 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var avalues = a._values;
     var aindex = a._index;
     var aptr = a._ptr;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // b sparse
     var bvalues = b._values;
     var bindex = b._index;
     var bptr = b._ptr;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b._data === undefined ? b._datatype : b.getDataType();
 
     // rows & columns
     var arows = a._size[0];
@@ -686,7 +687,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
     var mf = multiplyScalar;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signatures that matches (dt, dt)
@@ -704,7 +705,7 @@ var createMultiply = exports.createMultiply = /* #__PURE__ */(0, _factory.factor
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // workspace

package/lib/cjs/function/arithmetic/round.js

@@ -6,19 +6,22 @@ Object.defineProperty(exports, "__esModule", {
 exports.createRound = void 0;
 var _factory = require("../../utils/factory.js");
 var _collection = require("../../utils/collection.js");
+var _number = require("../../utils/number.js");
+var _nearlyEqual = require("../../utils/bignumber/nearlyEqual.js");
 var _matAlgo11xS0s = require("../../type/matrix/utils/matAlgo11xS0s.js");
 var _matAlgo12xSfs = require("../../type/matrix/utils/matAlgo12xSfs.js");
 var _matAlgo14xDs = require("../../type/matrix/utils/matAlgo14xDs.js");
 var _index = require("../../plain/number/index.js");
 var NO_INT = 'Number of decimals in function round must be an integer';
 var name = 'round';
-var dependencies = ['typed', 'matrix', 'equalScalar', 'zeros', 'BigNumber', 'DenseMatrix'];
+var dependencies = ['typed', 'config', 'matrix', 'equalScalar', 'zeros', 'BigNumber', 'DenseMatrix'];
 var createRound = exports.createRound = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
   var typed = _ref.typed,
+    config = _ref.config,
     matrix = _ref.matrix,
     equalScalar = _ref.equalScalar,
     zeros = _ref.zeros,
-    BigNumber = _ref.BigNumber,
+    _BigNumber = _ref.BigNumber,
     DenseMatrix = _ref.DenseMatrix;
   var matAlgo11xS0s = (0, _matAlgo11xS0s.createMatAlgo11xS0s)({
     typed: typed,
@@ -31,6 +34,9 @@ var createRound = exports.createRound = /* #__PURE__ */(0, _factory.factory)(nam
   var matAlgo14xDs = (0, _matAlgo14xDs.createMatAlgo14xDs)({
     typed: typed
   });
+  function toExponent(epsilon) {
+    return Math.abs((0, _number.splitNumber)(epsilon).exponent);
+  }
 
   /**
    * Round a value towards the nearest rounded value.
@@ -77,13 +83,27 @@ var createRound = exports.createRound = /* #__PURE__ */(0, _factory.factory)(nam
    * @return {number | BigNumber | Fraction | Complex | Array | Matrix} Rounded value
    */
   return typed(name, {
-    number: _index.roundNumber,
-    'number, number': _index.roundNumber,
+    number: function number(x) {
+      // Handle round off errors by first rounding to epsilon precision
+      var xEpsilon = (0, _index.roundNumber)(x, toExponent(config.epsilon));
+      var xSelected = (0, _number.nearlyEqual)(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return (0, _index.roundNumber)(xSelected);
+    },
+    'number, number': function numberNumber(x, n) {
+      // Same as number: unless user specifies more decimals than epsilon
+      var epsilonExponent = toExponent(config.epsilon);
+      if (n >= epsilonExponent) {
+        return (0, _index.roundNumber)(x, n);
+      }
+      var xEpsilon = (0, _index.roundNumber)(x, epsilonExponent);
+      var xSelected = (0, _number.nearlyEqual)(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return (0, _index.roundNumber)(xSelected, n);
+    },
     'number, BigNumber': function numberBigNumber(x, n) {
       if (!n.isInteger()) {
         throw new TypeError(NO_INT);
       }
-      return new BigNumber(x).toDecimalPlaces(n.toNumber());
+      return new _BigNumber(x).toDecimalPlaces(n.toNumber());
     },
     Complex: function Complex(x) {
       return x.round();
@@ -102,13 +122,24 @@ var createRound = exports.createRound = /* #__PURE__ */(0, _factory.factory)(nam
       return x.round(_n);
     },
     BigNumber: function BigNumber(x) {
-      return x.toDecimalPlaces(0);
+      // Handle round off errors by first rounding to epsilon precision
+      var xEpsilon = new _BigNumber(x).toDecimalPlaces(toExponent(config.epsilon));
+      var xSelected = (0, _nearlyEqual.nearlyEqual)(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return xSelected.toDecimalPlaces(0);
     },
     'BigNumber, BigNumber': function BigNumberBigNumber(x, n) {
       if (!n.isInteger()) {
         throw new TypeError(NO_INT);
       }
-      return x.toDecimalPlaces(n.toNumber());
+
+      // Same as BigNumber: unless user specifies more decimals than epsilon
+      var epsilonExponent = toExponent(config.epsilon);
+      if (n >= epsilonExponent) {
+        return x.toDecimalPlaces(n.toNumber());
+      }
+      var xEpsilon = x.toDecimalPlaces(epsilonExponent);
+      var xSelected = (0, _nearlyEqual.nearlyEqual)(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return xSelected.toDecimalPlaces(n.toNumber());
     },
     Fraction: function Fraction(x) {
       return x.round();

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

@@ -66,9 +66,9 @@ var createDot = exports.createDot = /* #__PURE__ */(0, _factory.factory)(name, d
   function _denseDot(a, b) {
     var N = _validateDim(a, b);
     var adata = (0, _is.isMatrix)(a) ? a._data : a;
-    var adt = (0, _is.isMatrix)(a) ? a._datatype : undefined;
+    var adt = (0, _is.isMatrix)(a) ? a._datatype || a.getDataType() : undefined;
     var bdata = (0, _is.isMatrix)(b) ? b._data : b;
-    var bdt = (0, _is.isMatrix)(b) ? b._datatype : undefined;
+    var bdt = (0, _is.isMatrix)(b) ? b._datatype || b.getDataType() : undefined;
 
     // are these 2-dimensional column vectors? (as opposed to 1-dimensional vectors)
     var aIsColumn = _size(a).length === 2;
@@ -77,7 +77,7 @@ var createDot = exports.createDot = /* #__PURE__ */(0, _factory.factory)(name, d
     var mul = multiplyScalar;
 
     // process data types
-    if (adt && bdt && adt === bdt && typeof adt === 'string') {
+    if (adt && bdt && adt === bdt && typeof adt === 'string' && adt !== 'mixed') {
       var dt = adt;
       // find signatures that matches (dt, dt)
       add = typed.find(addScalar, [dt, dt]);

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 12.3.1
- * @date    2024-02-01
+ * @version 12.4.0
+ * @date    2024-02-22
  *
  * @license
  * Copyright (C) 2013-2024 Jos de Jong <wjosdejong@gmail.com>

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

@@ -33,13 +33,13 @@ var createMatAlgo01xDSid = exports.createMatAlgo01xDSid = /* #__PURE__ */(0, _fa
     // dense matrix arrays
     var adata = denseMatrix._data;
     var asize = denseMatrix._size;
-    var adt = denseMatrix._datatype;
+    var adt = denseMatrix._datatype || denseMatrix.getDataType();
     // sparse matrix arrays
     var bvalues = sparseMatrix._values;
     var bindex = sparseMatrix._index;
     var bptr = sparseMatrix._ptr;
     var bsize = sparseMatrix._size;
-    var bdt = sparseMatrix._datatype;
+    var bdt = sparseMatrix._datatype || sparseMatrix._data === undefined ? sparseMatrix._datatype : sparseMatrix.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -61,7 +61,7 @@ var createMatAlgo01xDSid = exports.createMatAlgo01xDSid = /* #__PURE__ */(0, _fa
     var columns = asize[1];
 
     // process data types
-    var dt = typeof adt === 'string' && adt === bdt ? adt : undefined;
+    var dt = typeof adt === 'string' && adt !== 'mixed' && adt === bdt ? adt : undefined;
     // callback function
     var cf = dt ? typed.find(callback, [dt, dt]) : callback;
 
@@ -110,7 +110,7 @@ var createMatAlgo01xDSid = exports.createMatAlgo01xDSid = /* #__PURE__ */(0, _fa
     return denseMatrix.createDenseMatrix({
       data: cdata,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -34,13 +34,13 @@ var createMatAlgo02xDS0 = exports.createMatAlgo02xDS0 = /* #__PURE__ */(0, _fact
     // dense matrix arrays
     var adata = denseMatrix._data;
     var asize = denseMatrix._size;
-    var adt = denseMatrix._datatype;
+    var adt = denseMatrix._datatype || denseMatrix.getDataType();
     // sparse matrix arrays
     var bvalues = sparseMatrix._values;
     var bindex = sparseMatrix._index;
     var bptr = sparseMatrix._ptr;
     var bsize = sparseMatrix._size;
-    var bdt = sparseMatrix._datatype;
+    var bdt = sparseMatrix._datatype || sparseMatrix._data === undefined ? sparseMatrix._datatype : sparseMatrix.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -71,7 +71,7 @@ var createMatAlgo02xDS0 = exports.createMatAlgo02xDS0 = /* #__PURE__ */(0, _fact
     var cf = callback;
 
     // process data types
-    if (typeof adt === 'string' && adt === bdt) {
+    if (typeof adt === 'string' && adt === bdt && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signature that matches (dt, dt)
@@ -114,7 +114,7 @@ var createMatAlgo02xDS0 = exports.createMatAlgo02xDS0 = /* #__PURE__ */(0, _fact
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -33,13 +33,13 @@ var createMatAlgo03xDSf = exports.createMatAlgo03xDSf = /* #__PURE__ */(0, _fact
     // dense matrix arrays
     var adata = denseMatrix._data;
     var asize = denseMatrix._size;
-    var adt = denseMatrix._datatype;
+    var adt = denseMatrix._datatype || denseMatrix.getDataType();
     // sparse matrix arrays
     var bvalues = sparseMatrix._values;
     var bindex = sparseMatrix._index;
     var bptr = sparseMatrix._ptr;
     var bsize = sparseMatrix._size;
-    var bdt = sparseMatrix._datatype;
+    var bdt = sparseMatrix._datatype || sparseMatrix._data === undefined ? sparseMatrix._datatype : sparseMatrix.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -68,7 +68,7 @@ var createMatAlgo03xDSf = exports.createMatAlgo03xDSf = /* #__PURE__ */(0, _fact
     var cf = callback;
 
     // process data types
-    if (typeof adt === 'string' && adt === bdt) {
+    if (typeof adt === 'string' && adt === bdt && adt !== 'mixed') {
       // datatype
       dt = adt;
       // convert 0 to the same datatype
@@ -120,7 +120,7 @@ var createMatAlgo03xDSf = exports.createMatAlgo03xDSf = /* #__PURE__ */(0, _fact
     return denseMatrix.createDenseMatrix({
       data: cdata,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -35,13 +35,13 @@ var createMatAlgo04xSidSid = exports.createMatAlgo04xSidSid = /* #__PURE__ */(0,
     var aindex = a._index;
     var aptr = a._ptr;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // sparse matrix arrays
     var bvalues = b._values;
     var bindex = b._index;
     var bptr = b._ptr;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b._data === undefined ? b._datatype : b.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -67,7 +67,7 @@ var createMatAlgo04xSidSid = exports.createMatAlgo04xSidSid = /* #__PURE__ */(0,
     var cf = callback;
 
     // process data types
-    if (typeof adt === 'string' && adt === bdt) {
+    if (typeof adt === 'string' && adt === bdt && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signature that matches (dt, dt)
@@ -177,7 +177,7 @@ var createMatAlgo04xSidSid = exports.createMatAlgo04xSidSid = /* #__PURE__ */(0,
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

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

@@ -35,13 +35,13 @@ var createMatAlgo05xSfSf = exports.createMatAlgo05xSfSf = /* #__PURE__ */(0, _fa
     var aindex = a._index;
     var aptr = a._ptr;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // sparse matrix arrays
     var bvalues = b._values;
     var bindex = b._index;
     var bptr = b._ptr;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b._data === undefined ? b._datatype : b.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -67,7 +67,7 @@ var createMatAlgo05xSfSf = exports.createMatAlgo05xSfSf = /* #__PURE__ */(0, _fa
     var cf = callback;
 
     // process data types
-    if (typeof adt === 'string' && adt === bdt) {
+    if (typeof adt === 'string' && adt === bdt && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signature that matches (dt, dt)
@@ -169,7 +169,7 @@ var createMatAlgo05xSfSf = exports.createMatAlgo05xSfSf = /* #__PURE__ */(0, _fa
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

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

@@ -34,11 +34,11 @@ var createMatAlgo06xS0S0 = exports.createMatAlgo06xS0S0 = /* #__PURE__ */(0, _fa
     // sparse matrix arrays
     var avalues = a._values;
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // sparse matrix arrays
     var bvalues = b._values;
     var bsize = b._size;
-    var bdt = b._datatype;
+    var bdt = b._datatype || b._data === undefined ? b._datatype : b.getDataType();
 
     // validate dimensions
     if (asize.length !== bsize.length) {
@@ -64,7 +64,7 @@ var createMatAlgo06xS0S0 = exports.createMatAlgo06xS0S0 = /* #__PURE__ */(0, _fa
     var cf = callback;
 
     // process data types
-    if (typeof adt === 'string' && adt === bdt) {
+    if (typeof adt === 'string' && adt === bdt && adt !== 'mixed') {
       // datatype
       dt = adt;
       // find signature that matches (dt, dt)
@@ -151,7 +151,7 @@ var createMatAlgo06xS0S0 = exports.createMatAlgo06xS0S0 = /* #__PURE__ */(0, _fa
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });