78 lines
2.7 KiB
JavaScript
78 lines
2.7 KiB
JavaScript
|
import { factory } from '../../utils/factory.js';
|
||
|
|
import { extend } from '../../utils/object.js';
|
||
|
|
import { createAlgorithm11 } from '../../type/matrix/utils/algorithm11.js';
|
||
|
|
import { createAlgorithm14 } from '../../type/matrix/utils/algorithm14.js';
|
||
|
|
var name = 'divide';
|
||
|
|
var dependencies = ['typed', 'matrix', 'multiply', 'equalScalar', 'divideScalar', 'inv'];
|
||
|
|
export var createDivide = /* #__PURE__ */factory(name, dependencies, _ref => {
|
||
|
|
var {
|
||
|
|
typed,
|
||
|
|
matrix,
|
||
|
|
multiply,
|
||
|
|
equalScalar,
|
||
|
|
divideScalar,
|
||
|
|
inv
|
||
|
|
} = _ref;
|
||
|
|
var algorithm11 = createAlgorithm11({
|
||
|
|
typed,
|
||
|
|
equalScalar
|
||
|
|
});
|
||
|
|
var algorithm14 = createAlgorithm14({
|
||
|
|
typed
|
||
|
|
});
|
||
|
|
/**
|
||
|
|
* Divide two values, `x / y`.
|
||
|
|
* To divide matrices, `x` is multiplied with the inverse of `y`: `x * inv(y)`.
|
||
|
|
*
|
||
|
|
* Syntax:
|
||
|
|
*
|
||
|
|
* math.divide(x, y)
|
||
|
|
*
|
||
|
|
* Examples:
|
||
|
|
*
|
||
|
|
* math.divide(2, 3) // returns number 0.6666666666666666
|
||
|
|
*
|
||
|
|
* const a = math.complex(5, 14)
|
||
|
|
* const b = math.complex(4, 1)
|
||
|
|
* math.divide(a, b) // returns Complex 2 + 3i
|
||
|
|
*
|
||
|
|
* const c = [[7, -6], [13, -4]]
|
||
|
|
* const d = [[1, 2], [4, 3]]
|
||
|
|
* math.divide(c, d) // returns Array [[-9, 4], [-11, 6]]
|
||
|
|
*
|
||
|
|
* const e = math.unit('18 km')
|
||
|
|
* math.divide(e, 4.5) // returns Unit 4 km
|
||
|
|
*
|
||
|
|
* See also:
|
||
|
|
*
|
||
|
|
* multiply
|
||
|
|
*
|
||
|
|
* @param {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} x Numerator
|
||
|
|
* @param {number | BigNumber | Fraction | Complex | Array | Matrix} y Denominator
|
||
|
|
* @return {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} Quotient, `x / y`
|
||
|
|
*/
|
||
|
|
|
||
|
|
return typed('divide', extend({
|
||
|
|
// we extend the signatures of divideScalar with signatures dealing with matrices
|
||
|
|
'Array | Matrix, Array | Matrix': function ArrayMatrixArrayMatrix(x, y) {
|
||
|
|
// TODO: implement matrix right division using pseudo inverse
|
||
|
|
// https://www.mathworks.nl/help/matlab/ref/mrdivide.html
|
||
|
|
// https://www.gnu.org/software/octave/doc/interpreter/Arithmetic-Ops.html
|
||
|
|
// https://stackoverflow.com/questions/12263932/how-does-gnu-octave-matrix-division-work-getting-unexpected-behaviour
|
||
|
|
return multiply(x, inv(y));
|
||
|
|
},
|
||
|
|
'DenseMatrix, any': function DenseMatrixAny(x, y) {
|
||
|
|
return algorithm14(x, y, divideScalar, false);
|
||
|
|
},
|
||
|
|
'SparseMatrix, any': function SparseMatrixAny(x, y) {
|
||
|
|
return algorithm11(x, y, divideScalar, false);
|
||
|
|
},
|
||
|
|
'Array, any': function ArrayAny(x, y) {
|
||
|
|
// use matrix implementation
|
||
|
|
return algorithm14(matrix(x), y, divideScalar, false).valueOf();
|
||
|
|
},
|
||
|
|
'any, Array | Matrix': function anyArrayMatrix(x, y) {
|
||
|
|
return multiply(x, inv(y));
|
||
|
|
}
|
||
|
|
}, divideScalar.signatures));
|
||
|
|
});
|