import { factory } from '../../utils/factory.js'; import { createMatAlgo01xDSid } from '../../type/matrix/utils/matAlgo01xDSid.js'; import { createMatAlgo04xSidSid } from '../../type/matrix/utils/matAlgo04xSidSid.js'; import { createMatAlgo10xSids } from '../../type/matrix/utils/matAlgo10xSids.js'; import { createMatrixAlgorithmSuite } from '../../type/matrix/utils/matrixAlgorithmSuite.js'; import { gcdNumber } from '../../plain/number/index.js'; var name = 'gcd'; var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix']; export var createGcd = /* #__PURE__ */factory(name, dependencies, _ref => { var { typed, matrix, equalScalar, BigNumber, DenseMatrix } = _ref; var matAlgo01xDSid = createMatAlgo01xDSid({ typed }); var matAlgo04xSidSid = createMatAlgo04xSidSid({ typed, equalScalar }); var matAlgo10xSids = createMatAlgo10xSids({ typed, DenseMatrix }); var matrixAlgorithmSuite = createMatrixAlgorithmSuite({ typed, matrix }); var gcdTypes = 'number | BigNumber | Fraction | Matrix | Array'; var gcdManySignature = {}; gcdManySignature["".concat(gcdTypes, ", ").concat(gcdTypes, ", ...").concat(gcdTypes)] = typed.referToSelf(self => (a, b, args) => { var res = self(a, b); for (var i = 0; i < args.length; i++) { res = self(res, args[i]); } return res; }); /** * Calculate the greatest common divisor for two or more values or arrays. * * For matrices, the function is evaluated element wise. * * Syntax: * * math.gcd(a, b) * math.gcd(a, b, c, ...) * * Examples: * * math.gcd(8, 12) // returns 4 * math.gcd(-4, 6) // returns 2 * math.gcd(25, 15, -10) // returns 5 * * math.gcd([8, -4], [12, 6]) // returns [4, 2] * * See also: * * lcm, xgcd * * @param {... number | BigNumber | Fraction | Array | Matrix} args Two or more integer numbers * @return {number | BigNumber | Fraction | Array | Matrix} The greatest common divisor */ return typed(name, { 'number, number': gcdNumber, 'BigNumber, BigNumber': _gcdBigNumber, 'Fraction, Fraction': (x, y) => x.gcd(y) }, matrixAlgorithmSuite({ SS: matAlgo04xSidSid, DS: matAlgo01xDSid, Ss: matAlgo10xSids }), gcdManySignature); /** * Calculate gcd for BigNumbers * @param {BigNumber} a * @param {BigNumber} b * @returns {BigNumber} Returns greatest common denominator of a and b * @private */ function _gcdBigNumber(a, b) { if (!a.isInt() || !b.isInt()) { throw new Error('Parameters in function gcd must be integer numbers'); } // https://en.wikipedia.org/wiki/Euclidean_algorithm var zero = new BigNumber(0); while (!b.isZero()) { var r = a.mod(b); a = b; b = r; } return a.lt(zero) ? a.neg() : a; } });