Index: include/experimental/numeric =================================================================== --- include/experimental/numeric +++ include/experimental/numeric @@ -0,0 +1,79 @@ +// -*- C++ -*- +//===--------------------------- numeric ----------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +#ifndef _LIBCPP_EXPERIMENTAL_NUMERIC +#define _LIBCPP_EXPERIMENTAL_NUMERIC +/* + experimental/numeric synopsis + +// C++1z +namespace std { +namespace experimental { +inline namespace fundamentals_v2 { + + // 13.1.2, Greatest common divisor + template + constexpr common_type_t gcd(M m, N n); + + // 13.1.3, Least common multiple + template + constexpr common_type_t lcm(M m, N n); + +} // namespace fundamentals_v2 +} // namespace experimental +} // namespace std + + */ + +#include +#include +#include // is_integral +#include // abs + +#if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER) +#pragma GCC system_header +#endif + +_LIBCPP_BEGIN_NAMESPACE_EXPERIMENTAL + +template ::value> struct __abs; + +template +struct __abs<_Tp, true> { + _LIBCPP_INLINE_VISIBILITY _Tp operator()(_Tp __t) const noexcept { return _VSTD::abs(__t); } +}; + +template +struct __abs<_Tp, false> { + _LIBCPP_INLINE_VISIBILITY _Tp operator()(_Tp __t) const noexcept { return __t; } +}; + + +template +constexpr common_type_t<_Tp,_Up> _LIBCPP_INLINE_VISIBILITY +gcd(_Tp __m, _Up __n) +{ + static_assert((is_integral<_Tp>::value && is_integral<_Up>::value), "Arguments to gcd must be integer types"); + return __n == 0 ? __abs<_Tp>()(__m) : gcd(__n, __abs<_Tp>()(__m) % __abs<_Up>()(__n)); +} + +template +constexpr common_type_t<_Tp,_Up> _LIBCPP_INLINE_VISIBILITY +lcm(_Tp __m, _Up __n) +{ + static_assert((is_integral<_Tp>::value && is_integral<_Up>::value), "Arguments to lcm must be integer types"); + if (__m == 0 || __n == 0) + return 0; + return (__abs<_Tp>()(__m) / gcd(__m,__n)) * __abs<_Up>()(__n); +} + +_LIBCPP_END_NAMESPACE_EXPERIMENTAL + +#endif /* _LIBCPP_EXPERIMENTAL_MAP */ Index: test/std/experimental/numeric/numeric.ops.overview/nothing_to_do.pass.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops.overview/nothing_to_do.pass.cpp +++ test/std/experimental/numeric/numeric.ops.overview/nothing_to_do.pass.cpp @@ -0,0 +1,12 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +#include + +int main () {} Index: test/std/experimental/numeric/numeric.ops/nothing_to_do.pass.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/nothing_to_do.pass.cpp +++ test/std/experimental/numeric/numeric.ops/nothing_to_do.pass.cpp @@ -0,0 +1,12 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +#include + +int main () {} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral1.fail.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral1.fail.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral1.fail.cpp @@ -0,0 +1,24 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> gcd(_M __m, _N __n) + +// Remarks: If either M or N is not an integer type, the program is ill-formed. + +#include + + +int main() +{ + std::experimental::gcd(2.0, 4); +} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral2.fail.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral2.fail.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.not_integral2.fail.cpp @@ -0,0 +1,24 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> gcd(_M __m, _N __n) + +// Remarks: If either M or N is not an integer type, the program is ill-formed. + +#include + + +int main() +{ + std::experimental::gcd(4, 6.0); +} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.pass.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.pass.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.gcd/gcd.pass.cpp @@ -0,0 +1,251 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> gcd(_M __m, _N __n) + +#include +#include +#include + +template +void test0(Input1 in1, Input2 in2, Output out) +{ + static_assert((std::is_same::value), "" ); + assert(out == std::experimental::gcd(in1, in2)); +} + + +template +void test_same_signed() +{ + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + test0(-1, 0, 1); + test0( 0, -1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); + + test0( -1, 1, 1); + test0( -2, 3, 1); + test0( -2, 4, 2); + test0(-36, 17, 1); + test0(-36, 18, 18); + + test0( 1, -1, 1); + test0( 2, -3, 1); + test0( 2, -4, 2); + test0( 36, -17, 1); + test0( 36, -18, 18); + + test0( -1, -1, 1); + test0( -2, -3, 1); + test0( -2, -4, 2); + test0(-36, -17, 1); + test0(-36, -18, 18); +} + +template +void test_same_unsigned() +{ + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); +} + + +template +void test_different_signed() +{ + static_assert(( std::is_signed::value), "" ); + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + test0(-1, 0, 1); + test0( 0, -1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); + + test0( -1, 1, 1); + test0( -2, 3, 1); + test0( -2, 4, 2); + test0(-36, 17, 1); + test0(-36, 18, 18); + + test0( 1, -1, 1); + test0( 2, -3, 1); + test0( 2, -4, 2); + test0( 36, -17, 1); + test0( 36, -18, 18); + + test0( -1, -1, 1); + test0( -2, -3, 1); + test0( -2, -4, 2); + test0(-36, -17, 1); + test0(-36, -18, 18); +} + +template +void test_different_unsigned() +{ + static_assert((!std::is_signed::value), "" ); + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); +} + +template +void test_mixed_first_signed() +{ + static_assert(( std::is_signed::value), "" ); + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( -1, 1, 1); + test0( -2, 3, 1); // relatively prime + test0( -2, 4, 2); + test0(-36, 17, 1); // relatively prime + test0(-36, 18, 18); +} + +template +void test_mixed_second_signed() +{ + static_assert((!std::is_signed::value), "" ); + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when m and n are both zero. + test0( 0, 0, 0); + test0( 1, 0, 1); + test0( 0, 1, 1); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 1); // relatively prime + test0( 2, 4, 2); + test0( 36, 17, 1); // relatively prime + test0( 36, 18, 18); + +// Otherwise, returns the greatest common divisor of |m| and |n|. + test0( 1, -1, 1); + test0( 2, -3, 1); // relatively prime + test0( 2, -4, 2); + test0( 36, -17, 1); // relatively prime + test0( 36, -18, 18); +} + +int main() +{ + test_same_signed(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + + test_same_signed< int8_t>(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + + test_same_unsigned< uint8_t>(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); +} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral1.fail.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral1.fail.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral1.fail.cpp @@ -0,0 +1,24 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> lcm(_M __m, _N __n) + +// Remarks: If either M or N is not an integer type, the program is ill-formed. + +#include + + +int main() +{ + std::experimental::lcm(2.0, 4); +} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral2.fail.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral2.fail.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.not_integral2.fail.cpp @@ -0,0 +1,24 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> lcm(_M __m, _N __n) + +// Remarks: If either M or N is not an integer type, the program is ill-formed. + +#include + + +int main() +{ + std::experimental::lcm(4, 6.0); +} Index: test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.pass.cpp =================================================================== --- test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.pass.cpp +++ test/std/experimental/numeric/numeric.ops/numeric.ops.lcm/lcm.pass.cpp @@ -0,0 +1,252 @@ +//===----------------------------------------------------------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// + +// UNSUPPORTED: c++98, c++03, c++11 +// + +// template +// constexpr common_type_t<_M,_N> gcd(_M __m, _N __n) + +#include +#include +#include + +template +void test0(Input1 in1, Input2 in2, Output out) +{ + static_assert((std::is_same::value), "" ); + assert(out == std::experimental::lcm(in1, in2)); +} + + +template +void test_same_signed() +{ + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0( 0, 1, 0); + test0(-1, 0, 0); + test0( 0, -1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); + + test0( -1, 1, 1); + test0( -2, 3, 6); + test0( -2, 4, 4); + test0( -3, 17, 51); + test0(-36, 18, 36); + + test0( 1, -1, 1); + test0( 2, -3, 6); + test0( 2, -4, 4); + test0( 3, -17, 51); + test0( 36, -18, 36); + + test0( -1, -1, 1); + test0( -2, -3, 6); + test0( -2, -4, 4); + test0( -3, -17, 51); + test0(-36, -18, 36); +} + +template +void test_same_unsigned() +{ + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0( 0, 1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); +} + + +template +void test_different_signed() +{ + static_assert(( std::is_signed::value), "" ); + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0( 0, 1, 0); + test0(-1, 0, 0); + test0( 0, -1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); + + test0( -1, 1, 1); + test0( -2, 3, 6); + test0( -2, 4, 4); + test0( -3, 17, 51); + test0(-36, 18, 36); + + test0( 1, -1, 1); + test0( 2, -3, 6); + test0( 2, -4, 4); + test0( 3, -17, 51); + test0( 36, -18, 36); + + test0( -1, -1, 1); + test0( -2, -3, 6); + test0( -2, -4, 4); + test0( -3, -17, 51); + test0(-36, -18, 36); +} + +template +void test_different_unsigned() +{ + static_assert((!std::is_signed::value), "" ); + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0( 0, 1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); +} + +template +void test_mixed_first_signed() +{ + static_assert(( std::is_signed::value), "" ); + static_assert((!std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0(-1, 0, 0); + test0( 0, 1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); + + test0( -1, 1, 1); + test0( -2, 3, 6); + test0( -2, 4, 4); + test0( -3, 17, 51); + test0(-36, 18, 36); +} + +template +void test_mixed_second_signed() +{ + static_assert((!std::is_signed::value), "" ); + static_assert(( std::is_signed::value), "" ); + +// Returns: zero when either m or n is zero. + test0( 0, 0, 0); + test0( 1, 0, 0); + test0( 0, 1, 0); + test0( 0, -1, 0); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, 1, 1); + test0( 2, 3, 6); + test0( 2, 4, 4); + test0( 3, 17, 51); + test0( 36, 18, 36); + +// Otherwise, returns the least common multiple of |m| and |n|. + test0( 1, -1, 1); + test0( 2, -3, 6); + test0( 2, -4, 4); + test0( 3, -17, 51); + test0( 36, -18, 36); +} + +int main() +{ + test_same_signed(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + + test_same_signed< int8_t>(); + test_same_signed(); + test_same_signed(); + test_same_signed(); + + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + + test_same_unsigned< uint8_t>(); + test_same_unsigned(); + test_same_unsigned(); + test_same_unsigned(); + + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + test_different_signed(); + + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + test_different_unsigned(); + + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + test_mixed_first_signed(); + + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); + test_mixed_second_signed(); +}