diff options
author | Cedric Nugteren <web@cedricnugteren.nl> | 2018-10-31 10:37:21 +0100 |
---|---|---|
committer | GitHub <noreply@github.com> | 2018-10-31 10:37:21 +0100 |
commit | 4215bbe62a03a84c09d3abb6c24b9e6eee0b5007 (patch) | |
tree | fdc57080f0ea3da59699907129cbf13f1374cea0 /src/utilities/utilities.cpp | |
parent | d45911b61dedafcbd74f65df263b4197697d6a81 (diff) | |
parent | 0b3d04f70902e00f86c572a5e3c379f9335b216f (diff) |
Merge pull request #330 from vbkaisetsu/CLBlast-270-col2im
Add col2im function
Diffstat (limited to 'src/utilities/utilities.cpp')
-rw-r--r-- | src/utilities/utilities.cpp | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/src/utilities/utilities.cpp b/src/utilities/utilities.cpp index a8fdaa19..a6cd82e7 100644 --- a/src/utilities/utilities.cpp +++ b/src/utilities/utilities.cpp @@ -489,4 +489,30 @@ std::string GetDeviceName(const Device& device) { } // ================================================================================================= + +// Solve Bezout's identity +// a * p + b * q = r = GCD(a, b) +void EuclidGCD(int a, int b, int &p, int &q, int &r) { + p = 0; + q = 1; + int p_1 = 1; + int q_1 = 0; + for (;;) { + const int c = a % b; + if (c == 0) { + break; + } + const int p_2 = p_1; + const int q_2 = q_1; + p_1 = p; + q_1 = q; + p = p_2 - p_1 * (a / b); + q = q_2 - q_1 * (a / b); + a = b; + b = c; + } + r = b; +} + +// ================================================================================================= } // namespace clblast |