3c21913415
The current std::gcd and std::chrono::duration::_S_gcd algorithms are both recursive. This is potentially expensive to evaluate in constant expressions, because each level of recursion makes a new copy of the function to evaluate. The maximum number of steps is bounded (proportional to the number of decimal digits in the smaller value) and so unlikely to exceed the limit for constexpr nesting, but the memory usage is still suboptimal. By using an iterative algorithm we avoid that compile-time cost. Because looping in constexpr functions is not allowed until C++14, we need to keep the recursive implementation in duration::_S_gcd for C++11 mode. For std::gcd we can also optimise runtime performance by using the binary GCD algorithm. libstdc++-v3/ChangeLog: * include/std/chrono (duration::_S_gcd): Use iterative algorithm for C++14 and later. * include/std/numeric (__detail::__gcd): Replace recursive Euclidean algorithm with iterative version of binary GCD algorithm. * testsuite/26_numerics/gcd/1.cc: Test additional inputs. * testsuite/26_numerics/gcd/gcd_neg.cc: Adjust dg-error lines. * testsuite/26_numerics/lcm/lcm_neg.cc: Likewise. * testsuite/experimental/numeric/gcd.cc: Test additional inputs. * testsuite/26_numerics/gcd/2.cc: New test.
54 lines
2.2 KiB
C++
54 lines
2.2 KiB
C++
// Copyright (C) 2015-2020 Free Software Foundation, Inc.
|
|
//
|
|
// This file is part of the GNU ISO C++ Library. This library is free
|
|
// software; you can redistribute it and/or modify it under the
|
|
// terms of the GNU General Public License as published by the
|
|
// Free Software Foundation; either version 3, or (at your option)
|
|
// any later version.
|
|
|
|
// This library is distributed in the hope that it will be useful,
|
|
// but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
// GNU General Public License for more details.
|
|
|
|
// You should have received a copy of the GNU General Public License along
|
|
// with this library; see the file COPYING3. If not see
|
|
// <http://www.gnu.org/licenses/>.
|
|
|
|
// { dg-do compile { target c++17 } }
|
|
|
|
#include <numeric>
|
|
|
|
#ifndef __cpp_lib_gcd_lcm
|
|
# error "Feature-test macro for gcd missing"
|
|
#elif __cpp_lib_gcd_lcm != 201606
|
|
# error "Feature-test macro for gcd has wrong value"
|
|
#endif
|
|
|
|
using std::gcd;
|
|
using std::is_same_v;
|
|
|
|
static_assert( gcd(1071, 462) == 21, "" );
|
|
static_assert( gcd(2000, 20) == 20, "" );
|
|
static_assert( gcd(2011, 17) == 1, "GCD of two primes is 1" );
|
|
static_assert( gcd(200, 200) == 200, "GCD of equal numbers is that number" );
|
|
static_assert( gcd(0, 13) == 13, "GCD of any number and 0 is that number" );
|
|
static_assert( gcd(29, 0) == 29, "GCD of any number and 0 is that number" );
|
|
static_assert( gcd(0, 0) == 0, "Zarro Boogs found" );
|
|
|
|
static_assert(gcd(1u, 2) == 1, "unsigned and signed");
|
|
static_assert(gcd(9, 6u) == 3, "unsigned and signed");
|
|
static_assert(gcd(3, 4u) == 1, "signed and unsigned");
|
|
static_assert(gcd(32u, 24) == 8, "signed and unsigned");
|
|
static_assert(gcd(1u, -2) == 1, "unsigned and negative");
|
|
static_assert(gcd(-21, 28u) == 7, "unsigned and negative");
|
|
static_assert(gcd(-3, 4u) == 1, "negative and unsigned");
|
|
static_assert(gcd(33u, -44) == 11, "negative and unsigned");
|
|
static_assert(gcd(5u, 6u) == 1, "unsigned and unsigned");
|
|
static_assert(gcd(54u, 36u) == 18, "unsigned and unsigned");
|
|
static_assert(gcd(-5, -6) == 1, "negative and negative");
|
|
static_assert(gcd(-50, -60) == 10, "negative and negative");
|
|
|
|
static_assert( is_same_v<decltype(gcd(1l, 1)), long> );
|
|
static_assert( is_same_v<decltype(gcd(1ul, 1ull)), unsigned long long> );
|