mathjs 12.3.1 → 12.4.0

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

@@ -115,11 +115,11 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     // 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];
@@ -132,7 +132,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -159,7 +159,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     return a.createDenseMatrix({
       data: c,
       size: [bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -203,10 +203,10 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     // 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];
@@ -219,7 +219,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -248,7 +248,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     return a.createDenseMatrix({
       data: c,
       size: [arows],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -261,14 +261,15 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
    * @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];
@@ -282,7 +283,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -316,7 +317,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     return a.createDenseMatrix({
       data: c,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   }
 
@@ -332,13 +333,13 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     // 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');
@@ -359,7 +360,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -380,7 +381,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // loop b columns
@@ -444,14 +445,14 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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];
@@ -472,7 +473,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -525,13 +526,13 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     // 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
     });
   }
 
@@ -548,14 +549,14 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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];
@@ -573,7 +574,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -594,7 +595,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // workspace
@@ -661,12 +662,12 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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];
@@ -682,7 +683,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
     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)
@@ -700,7 +701,7 @@ export var createMultiply = /* #__PURE__ */factory(name, dependencies, _ref => {
       index: cindex,
       ptr: cptr,
       size: [arows, bcolumns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
 
     // workspace

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

@@ -1,19 +1,22 @@
 import { factory } from '../../utils/factory.js';
 import { deepMap } from '../../utils/collection.js';
+import { nearlyEqual, splitNumber } from '../../utils/number.js';
+import { nearlyEqual as bigNearlyEqual } from '../../utils/bignumber/nearlyEqual.js';
 import { createMatAlgo11xS0s } from '../../type/matrix/utils/matAlgo11xS0s.js';
 import { createMatAlgo12xSfs } from '../../type/matrix/utils/matAlgo12xSfs.js';
 import { createMatAlgo14xDs } from '../../type/matrix/utils/matAlgo14xDs.js';
 import { roundNumber } from '../../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'];
 export var createRound = /* #__PURE__ */factory(name, dependencies, _ref => {
   var {
     typed,
+    config,
     matrix,
     equalScalar,
     zeros,
-    BigNumber,
+    BigNumber: _BigNumber,
     DenseMatrix
   } = _ref;
   var matAlgo11xS0s = createMatAlgo11xS0s({
@@ -27,6 +30,9 @@ export var createRound = /* #__PURE__ */factory(name, dependencies, _ref => {
   var matAlgo14xDs = createMatAlgo14xDs({
     typed
   });
+  function toExponent(epsilon) {
+    return Math.abs(splitNumber(epsilon).exponent);
+  }
 
   /**
    * Round a value towards the nearest rounded value.
@@ -73,13 +79,27 @@ export var createRound = /* #__PURE__ */factory(name, dependencies, _ref => {
    * @return {number | BigNumber | Fraction | Complex | Array | Matrix} Rounded value
    */
   return typed(name, {
-    number: roundNumber,
-    'number, number': roundNumber,
+    number: function number(x) {
+      // Handle round off errors by first rounding to epsilon precision
+      var xEpsilon = roundNumber(x, toExponent(config.epsilon));
+      var xSelected = nearlyEqual(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return 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 roundNumber(x, n);
+      }
+      var xEpsilon = roundNumber(x, epsilonExponent);
+      var xSelected = nearlyEqual(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return 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();
@@ -98,13 +118,24 @@ export var createRound = /* #__PURE__ */factory(name, dependencies, _ref => {
       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 = bigNearlyEqual(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 = bigNearlyEqual(x, xEpsilon, config.epsilon) ? xEpsilon : x;
+      return xSelected.toDecimalPlaces(n.toNumber());
     },
     Fraction: function Fraction(x) {
       return x.round();

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

@@ -62,9 +62,9 @@ export var createDot = /* #__PURE__ */factory(name, dependencies, _ref => {
   function _denseDot(a, b) {
     var N = _validateDim(a, b);
     var adata = isMatrix(a) ? a._data : a;
-    var adt = isMatrix(a) ? a._datatype : undefined;
+    var adt = isMatrix(a) ? a._datatype || a.getDataType() : undefined;
     var bdata = isMatrix(b) ? b._data : b;
-    var bdt = isMatrix(b) ? b._datatype : undefined;
+    var bdt = 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;
@@ -73,7 +73,7 @@ export var createDot = /* #__PURE__ */factory(name, dependencies, _ref => {
     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/esm/type/matrix/utils/matAlgo01xDSid.js

@@ -29,13 +29,13 @@ export var createMatAlgo01xDSid = /* #__PURE__ */factory(name, dependencies, _re
     // 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) {
@@ -57,7 +57,7 @@ export var createMatAlgo01xDSid = /* #__PURE__ */factory(name, dependencies, _re
     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;
 
@@ -106,7 +106,7 @@ export var createMatAlgo01xDSid = /* #__PURE__ */factory(name, dependencies, _re
     return denseMatrix.createDenseMatrix({
       data: cdata,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -30,13 +30,13 @@ export var createMatAlgo02xDS0 = /* #__PURE__ */factory(name, dependencies, _ref
     // 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) {
@@ -67,7 +67,7 @@ export var createMatAlgo02xDS0 = /* #__PURE__ */factory(name, dependencies, _ref
     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)
@@ -110,7 +110,7 @@ export var createMatAlgo02xDS0 = /* #__PURE__ */factory(name, dependencies, _ref
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -29,13 +29,13 @@ export var createMatAlgo03xDSf = /* #__PURE__ */factory(name, dependencies, _ref
     // 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) {
@@ -64,7 +64,7 @@ export var createMatAlgo03xDSf = /* #__PURE__ */factory(name, dependencies, _ref
     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
@@ -116,7 +116,7 @@ export var createMatAlgo03xDSf = /* #__PURE__ */factory(name, dependencies, _ref
     return denseMatrix.createDenseMatrix({
       data: cdata,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === denseMatrix._datatype && bdt === sparseMatrix._datatype ? dt : undefined
     });
   };
 });

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

@@ -31,13 +31,13 @@ export var createMatAlgo04xSidSid = /* #__PURE__ */factory(name, dependencies, _
     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) {
@@ -63,7 +63,7 @@ export var createMatAlgo04xSidSid = /* #__PURE__ */factory(name, dependencies, _
     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)
@@ -173,7 +173,7 @@ export var createMatAlgo04xSidSid = /* #__PURE__ */factory(name, dependencies, _
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

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

@@ -31,13 +31,13 @@ export var createMatAlgo05xSfSf = /* #__PURE__ */factory(name, dependencies, _re
     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) {
@@ -63,7 +63,7 @@ export var createMatAlgo05xSfSf = /* #__PURE__ */factory(name, dependencies, _re
     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)
@@ -165,7 +165,7 @@ export var createMatAlgo05xSfSf = /* #__PURE__ */factory(name, dependencies, _re
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

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

@@ -30,11 +30,11 @@ export var createMatAlgo06xS0S0 = /* #__PURE__ */factory(name, dependencies, _re
     // 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) {
@@ -60,7 +60,7 @@ export var createMatAlgo06xS0S0 = /* #__PURE__ */factory(name, dependencies, _re
     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)
@@ -147,7 +147,7 @@ export var createMatAlgo06xS0S0 = /* #__PURE__ */factory(name, dependencies, _re
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

package/lib/esm/type/matrix/utils/matAlgo07xSSf.js

@@ -24,10 +24,10 @@ export var createMatAlgo07xSSf = /* #__PURE__ */factory(name, dependencies, _ref
   return function matAlgo07xSSf(a, b, callback) {
     // sparse matrix arrays
     var asize = a._size;
-    var adt = a._datatype;
+    var adt = a._datatype || a._data === undefined ? a._datatype : a.getDataType();
     // sparse matrix arrays
     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) {
@@ -51,7 +51,7 @@ export var createMatAlgo07xSSf = /* #__PURE__ */factory(name, dependencies, _ref
     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
@@ -99,7 +99,7 @@ export var createMatAlgo07xSSf = /* #__PURE__ */factory(name, dependencies, _ref
     return new DenseMatrix({
       data: cdata,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
   function _scatter(m, j, w, x, mark) {

package/lib/esm/type/matrix/utils/matAlgo08xS0Sid.js

@@ -31,13 +31,13 @@ export var createMatAlgo08xS0Sid = /* #__PURE__ */factory(name, dependencies, _r
     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) {
@@ -68,7 +68,7 @@ export var createMatAlgo08xS0Sid = /* #__PURE__ */factory(name, dependencies, _r
     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)
@@ -148,7 +148,7 @@ export var createMatAlgo08xS0Sid = /* #__PURE__ */factory(name, dependencies, _r
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });

package/lib/esm/type/matrix/utils/matAlgo09xS0Sf.js

@@ -31,13 +31,13 @@ export var createMatAlgo09xS0Sf = /* #__PURE__ */factory(name, dependencies, _re
     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) {
@@ -63,7 +63,7 @@ export var createMatAlgo09xS0Sf = /* #__PURE__ */factory(name, dependencies, _re
     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)
@@ -136,7 +136,7 @@ export var createMatAlgo09xS0Sf = /* #__PURE__ */factory(name, dependencies, _re
       index: cindex,
       ptr: cptr,
       size: [rows, columns],
-      datatype: dt
+      datatype: adt === a._datatype && bdt === b._datatype ? dt : undefined
     });
   };
 });