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)); });