gcd.js 4.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132
  1. "use strict";
  2. var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
  3. Object.defineProperty(exports, "__esModule", {
  4. value: true
  5. });
  6. exports.createGcd = void 0;
  7. var _defineProperty2 = _interopRequireDefault(require("@babel/runtime/helpers/defineProperty"));
  8. var _toConsumableArray2 = _interopRequireDefault(require("@babel/runtime/helpers/toConsumableArray"));
  9. var _factory = require("../../utils/factory.js");
  10. var _matAlgo01xDSid = require("../../type/matrix/utils/matAlgo01xDSid.js");
  11. var _matAlgo04xSidSid = require("../../type/matrix/utils/matAlgo04xSidSid.js");
  12. var _matAlgo10xSids = require("../../type/matrix/utils/matAlgo10xSids.js");
  13. var _matrixAlgorithmSuite = require("../../type/matrix/utils/matrixAlgorithmSuite.js");
  14. var _index = require("../../plain/number/index.js");
  15. var _ArgumentsError = require("../../error/ArgumentsError.js");
  16. var name = 'gcd';
  17. var dependencies = ['typed', 'matrix', 'equalScalar', 'BigNumber', 'DenseMatrix', 'concat'];
  18. var gcdTypes = 'number | BigNumber | Fraction | Matrix | Array';
  19. var gcdManyTypesSignature = "".concat(gcdTypes, ", ").concat(gcdTypes, ", ...").concat(gcdTypes);
  20. function is1d(array) {
  21. return !array.some(function (element) {
  22. return Array.isArray(element);
  23. });
  24. }
  25. var createGcd = /* #__PURE__ */(0, _factory.factory)(name, dependencies, function (_ref) {
  26. var _typed;
  27. var typed = _ref.typed,
  28. matrix = _ref.matrix,
  29. equalScalar = _ref.equalScalar,
  30. BigNumber = _ref.BigNumber,
  31. DenseMatrix = _ref.DenseMatrix,
  32. concat = _ref.concat;
  33. var matAlgo01xDSid = (0, _matAlgo01xDSid.createMatAlgo01xDSid)({
  34. typed: typed
  35. });
  36. var matAlgo04xSidSid = (0, _matAlgo04xSidSid.createMatAlgo04xSidSid)({
  37. typed: typed,
  38. equalScalar: equalScalar
  39. });
  40. var matAlgo10xSids = (0, _matAlgo10xSids.createMatAlgo10xSids)({
  41. typed: typed,
  42. DenseMatrix: DenseMatrix
  43. });
  44. var matrixAlgorithmSuite = (0, _matrixAlgorithmSuite.createMatrixAlgorithmSuite)({
  45. typed: typed,
  46. matrix: matrix,
  47. concat: concat
  48. });
  49. /**
  50. * Calculate the greatest common divisor for two or more values or arrays.
  51. *
  52. * For matrices, the function is evaluated element wise.
  53. *
  54. * Syntax:
  55. *
  56. * math.gcd(a, b)
  57. * math.gcd(a, b, c, ...)
  58. *
  59. * Examples:
  60. *
  61. * math.gcd(8, 12) // returns 4
  62. * math.gcd(-4, 6) // returns 2
  63. * math.gcd(25, 15, -10) // returns 5
  64. *
  65. * math.gcd([8, -4], [12, 6]) // returns [4, 2]
  66. *
  67. * See also:
  68. *
  69. * lcm, xgcd
  70. *
  71. * @param {... number | BigNumber | Fraction | Array | Matrix} args Two or more integer numbers
  72. * @return {number | BigNumber | Fraction | Array | Matrix} The greatest common divisor
  73. */
  74. return typed(name, {
  75. 'number, number': _index.gcdNumber,
  76. 'BigNumber, BigNumber': _gcdBigNumber,
  77. 'Fraction, Fraction': function FractionFraction(x, y) {
  78. return x.gcd(y);
  79. }
  80. }, matrixAlgorithmSuite({
  81. SS: matAlgo04xSidSid,
  82. DS: matAlgo01xDSid,
  83. Ss: matAlgo10xSids
  84. }), (_typed = {}, (0, _defineProperty2["default"])(_typed, gcdManyTypesSignature, typed.referToSelf(function (self) {
  85. return function (a, b, args) {
  86. var res = self(a, b);
  87. for (var i = 0; i < args.length; i++) {
  88. res = self(res, args[i]);
  89. }
  90. return res;
  91. };
  92. })), (0, _defineProperty2["default"])(_typed, "Array", typed.referToSelf(function (self) {
  93. return function (array) {
  94. if (array.length === 1 && Array.isArray(array[0]) && is1d(array[0])) {
  95. return self.apply(void 0, (0, _toConsumableArray2["default"])(array[0]));
  96. }
  97. if (is1d(array)) {
  98. return self.apply(void 0, (0, _toConsumableArray2["default"])(array));
  99. }
  100. throw new _ArgumentsError.ArgumentsError('gcd() supports only 1d matrices!');
  101. };
  102. })), (0, _defineProperty2["default"])(_typed, "Matrix", typed.referToSelf(function (self) {
  103. return function (matrix) {
  104. return self(matrix.toArray());
  105. };
  106. })), _typed));
  107. /**
  108. * Calculate gcd for BigNumbers
  109. * @param {BigNumber} a
  110. * @param {BigNumber} b
  111. * @returns {BigNumber} Returns greatest common denominator of a and b
  112. * @private
  113. */
  114. function _gcdBigNumber(a, b) {
  115. if (!a.isInt() || !b.isInt()) {
  116. throw new Error('Parameters in function gcd must be integer numbers');
  117. }
  118. // https://en.wikipedia.org/wiki/Euclidean_algorithm
  119. var zero = new BigNumber(0);
  120. while (!b.isZero()) {
  121. var r = a.mod(b);
  122. a = b;
  123. b = r;
  124. }
  125. return a.lt(zero) ? a.neg() : a;
  126. }
  127. });
  128. exports.createGcd = createGcd;