UNPKG

3.09 kBJavaScriptView Raw
1"use strict";
2
3Object.defineProperty(exports, "__esModule", {
4 value: true
5});
6exports.createDivide = void 0;
7
8var _factory = require("../../utils/factory");
9
10var _object = require("../../utils/object");
11
12var _algorithm = require("../../type/matrix/utils/algorithm11");
13
14var _algorithm2 = require("../../type/matrix/utils/algorithm14");
15
16var name = 'divide';
17var dependencies = ['typed', 'matrix', 'multiply', 'equalScalar', 'divideScalar', 'inv'];
18var createDivide = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
19 var typed = _ref.typed,
20 matrix = _ref.matrix,
21 multiply = _ref.multiply,
22 equalScalar = _ref.equalScalar,
23 divideScalar = _ref.divideScalar,
24 inv = _ref.inv;
25 var algorithm11 = (0, _algorithm.createAlgorithm11)({
26 typed: typed,
27 equalScalar: equalScalar
28 });
29 var algorithm14 = (0, _algorithm2.createAlgorithm14)({
30 typed: typed
31 });
32 /**
33 * Divide two values, `x / y`.
34 * To divide matrices, `x` is multiplied with the inverse of `y`: `x * inv(y)`.
35 *
36 * Syntax:
37 *
38 * math.divide(x, y)
39 *
40 * Examples:
41 *
42 * math.divide(2, 3) // returns number 0.6666666666666666
43 *
44 * const a = math.complex(5, 14)
45 * const b = math.complex(4, 1)
46 * math.divide(a, b) // returns Complex 2 + 3i
47 *
48 * const c = [[7, -6], [13, -4]]
49 * const d = [[1, 2], [4, 3]]
50 * math.divide(c, d) // returns Array [[-9, 4], [-11, 6]]
51 *
52 * const e = math.unit('18 km')
53 * math.divide(e, 4.5) // returns Unit 4 km
54 *
55 * See also:
56 *
57 * multiply
58 *
59 * @param {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} x Numerator
60 * @param {number | BigNumber | Fraction | Complex | Array | Matrix} y Denominator
61 * @return {number | BigNumber | Fraction | Complex | Unit | Array | Matrix} Quotient, `x / y`
62 */
63
64 return typed('divide', (0, _object.extend)({
65 // we extend the signatures of divideScalar with signatures dealing with matrices
66 'Array | Matrix, Array | Matrix': function ArrayMatrixArrayMatrix(x, y) {
67 // TODO: implement matrix right division using pseudo inverse
68 // https://www.mathworks.nl/help/matlab/ref/mrdivide.html
69 // https://www.gnu.org/software/octave/doc/interpreter/Arithmetic-Ops.html
70 // https://stackoverflow.com/questions/12263932/how-does-gnu-octave-matrix-division-work-getting-unexpected-behaviour
71 return multiply(x, inv(y));
72 },
73 'DenseMatrix, any': function DenseMatrixAny(x, y) {
74 return algorithm14(x, y, divideScalar, false);
75 },
76 'SparseMatrix, any': function SparseMatrixAny(x, y) {
77 return algorithm11(x, y, divideScalar, false);
78 },
79 'Array, any': function ArrayAny(x, y) {
80 // use matrix implementation
81 return algorithm14(matrix(x), y, divideScalar, false).valueOf();
82 },
83 'any, Array | Matrix': function anyArrayMatrix(x, y) {
84 return multiply(x, inv(y));
85 }
86 }, divideScalar.signatures));
87});
88exports.createDivide = createDivide;
\No newline at end of file