gpt4 book ai didi

c++ - GCC 5.3.1 C++ 在编译 Variadic 模板时停止

转载 作者:行者123 更新时间:2023-11-30 00:46:32 25 4
gpt4 key购买 nike

我写了下面的代码:

#include<array>
#include<type_traits>

namespace math{
namespace detail{

template<std::size_t... Is> struct seq{};

template<std::size_t N, std::size_t... Is>
struct gen_seq : gen_seq<N-1, N-1, Is...>{};

template<std::size_t... Is>
struct gen_seq<0, Is...> : seq<Is...>{};

template<class T,std::size_t N>
struct sin_coeffs{
using array_type = std::array<T,N>;
constexpr static inline T coeff(std::size_t n){
return power(-1, n-1) * inverse((T)factorial((2 * n)-1));
}
template<std::size_t...NS>
constexpr static array_type _coeffs(seq<NS...>){
return {{coeff(NS)...}};
}
constexpr static array_type coeffs=_coeffs(gen_seq<N>{});
};
}

template<class T,std::size_t N = max_factorial, class dcy = std::decay_t<T>>
constexpr std::enable_if_t<std::is_floating_point<dcy>::value,dcy> sin(T x) noexcept{
constexpr std::array<dcy,N>& coeffs = detail::sin_coeffs<dcy,N>::coeffs;
const dcy x_2 = x*x;
dcy pow = x;
dcy result = 0;
for(std::size_t i=0;i<N;++i){
result += coeffs[i] * pow;
pow*=x_2;
}
return result;
}
}

主要示例:

int main()
{
constexpr double d = math::sin(0.0);
}

代码被设计成一个 constexpr sin 函数,它使用一个包含系数的 constexpr 数组来进行所需的计算。

所有未列出的函数都存在于单独的头文件中,编译没有问题。

我正在尝试使用“indicies”技巧使用基于 this answer to another question 的 constexpr 函数填充数组.

我正在使用带有标志的 GCC 5.3.1 进行编译

--std=c++1z -pthread -g -O3 -MMD -MP -Wall -pedantic

编译器没有向我的代码发出错误,但在编译时停止了。

我已经让编译运行了几分钟,但它什么也没做。

我已经测试了代码中使用的所有函数,它们都可以独立于本节进行编译。

int main()
{
math::detail::sin_coeffs<double,20>::coeffs[0];
}

此代码片段还重现了使我相信它与 sin 函数本身无关但与 sin_coeffs 结构无关的问题。

编辑:以下是要求的其他功能:

#include <type_traits>
namespace math{

template<class T,class dcy = std::decay_t<T>>
constexpr inline std::enable_if_t<std::is_floating_point<T>::value,dcy> inverse(T value){
return (value == 0) ? 0.0 : 1.0 / value;
}
template <class T>
constexpr inline std::decay_t<T> sign(T value) {
return value < 0 ? -1 : 1;
}
template <typename T>
constexpr inline std::decay_t<T> abs(T value) {
return value * sign(value);
}
template<class T>
constexpr inline std::decay_t<T> power(T const& base, std::size_t const& pow){
if(pow==0){return 1;}
else if(pow == 1){return base;}
else{
T result = base;
for(std::size_t i=1;i<pow;++i){
result*=base;
}
return result;
}
}
constexpr std::intmax_t factorial(std::intmax_t const& n){
if(n==0){return 1;}
std::intmax_t result = n;
for(std::intmax_t i=n-1;i>0;--i){
result *=i;
}
return result;
}
constexpr static std::size_t max_factorial = 20;//replace with calculated version later
}

最佳答案

如何修复您的代码?

  1. 您缺少一个 const限定符在这里:
constexpr const std::array<dcy,N>& coeffs = /* ... */ ;
^^^^^
  1. 你的 gen_seq0 生成值至 N - 1 , 但关于你的 coeff函数,您需要来自 1 的值至 N .您可以通过以下任一方式解决此问题:
  • 更改基本模板以从 1 生成值至 N (详细解释见本答案底部):
template<std::size_t N, std::size_t... Is>
struct gen_seq: gen_seq<N-1, N, Is...> {};
// ^--- N instead of N - 1
  • 改变调用 coeff 的方式(感谢@T.C.,查看评论):
template<std::size_t...NS>
constexpr static array_type _coeffs(seq<NS...>){
return {{coeff(NS + 1)...}};
// ^^^^
}
  1. 你不应该使用 power计算 -1 ** n ,使用三元条件:
constexpr static inline T coeff(std::size_t n){
return (n%2 ? 1 : -1) * inverse((T)factorial((2 * n)-1));
}

有了这个,我可以计算系数:

auto arr = math::detail::sin_coeffs<double, 10>::coeffs;
for (auto x: arr) {
std::cout << x << " ";
}

输出:

1 -0.166667 0.00833333 -0.000198413 2.75573e-06 -2.50521e-08 1.6059e-10 -7.64716e-13 ...

据我所知,这些是正确的系数( 1-1/3!1/5! ……)。请注意,我必须使用 N = 10 ,否则我会溢出 std::intmax_t (在我的架构上)- 如果溢出 std::intmax_t,Clang 会在编译时警告您与 factorial (多么好的编译器!)。

通过以上两个修改,您的代码可以正常工作(max_factorial 值除外,但您可以根据需要调整它)。

参见 rextester 上的代码: http://rextester.com/CRR35028

您的代码中的主要问题是什么?

你的 gen_seq<N>0 生成一个序列至 N - 1 , 所以你调用 coeff(0) ,正在调用 power(-1, static_cast<size_t>(0) - 1) ,实际上是 power(-1, 18446744073709551615) (在我的架构上),无法编译。在 power 中添加一个简单的案例“修复”编译(表明这是一个问题,但没有解决真正的问题):

template<class T>
constexpr inline std::decay_t<T> power(T const& base, std::size_t const& pow) {
if (pow == static_cast<size_t>(0) - 1) { // Stupid test
return 1;
}
/* ... */
}

此外,您的 max_factorial值(value)也可能太大了。更正后power , 我无法为 max_factorial > 11 编译(我可能有 32 位 std::intmax_t,所以你可以超过这个,但我认为 20 在所有情况下都太大了)。

此外,对于 future 的问题,clang似乎提供了关于为什么它不编译的更好信息:

  • 它对迭代次数有限制,所以我能够找出 power正在做一个(几乎)无限循环。
  • factorial的返回值给出警告溢出 intmax_t .

如何gen_seq把戏?

你的 seq struct 基本上是 std::size_t... 的结构持有者(因为您不能将其直接存储在变量中)。

gen_seq是构建gen_seq<N, ...> 的“递归”结构使用 gen_seq<N - 1, ...> .它是如何工作的:

gen_seq<3>: gen_seq<2, 2>
gen_seq<2, 2>: gen_seq<1, 1, 2>
gen_seq<1, 1, 2>: gen_seq<0, 0, 1, 2>
gen_seq<0, 0, 1, 2>: seq<0, 1, 2> // Specialized case

如你所见,既然你继承了gen_seq<N - 1, N - 1, ...> , 最后继承自 seq0值(value),你不想要的。因为什么是“发送”到seq是可变参数 std::size_t... , 你想避免 0 , 所以你改成 gen_seq<N - 1, N, ...> :

gen_seq<3>: gen_seq<2, 3>
gen_seq<2, 3>: gen_seq<1, 2, 3>
gen_seq<1, 2, 3>: gen_seq<0, 1, 2, 3>
gen_seq<0, 1, 2, 3>: seq<1, 2, 3> // Specialized case

现在,当您调用 _coeffs(gen_seq<N>{}) ,您允许编译器推断模板参数 NS_coeffs .由此,您可以使用 NS_coeffspack expansion :

{{coeff(NS)...}};

关于c++ - GCC 5.3.1 C++ 在编译 Variadic 模板时停止,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/38069092/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com