OMP for 循环条件不是简单的关系比较

Question

我的程序由于错误而无法编译：OpenMP for 循环的条件必须是关系比较（'<'、'<='、'>'、'>=' 或循环变量'i'的'!=')，参考for (size_t i = 2; i * i <= n; i++)。如何修改和修复它而不影响性能？这是由于 OpenMP 版本旧而导致的问题吗？ （因为我记得以前在另一台旧版本的计算机上遇到过不同的问题，现在已经解决了。）

#include <iostream>
#include <cstdio>

int main(int argc, char **argv)
{
    if (argc != 2)
        return EXIT_FAILURE;
    size_t n;
    if (sscanf(argv[1], "%zu", &n) == 0)
        return EXIT_FAILURE;
    auto *prime = new size_t[n + 1];
    for (size_t i = 2; i <= n; i++)
        prime[i] = i;
    #pragma omp parallel for
    for (size_t i = 2; i * i <= n; i++)
        for (size_t j = i * i; j <= n; j += i)
            prime[j] = 0;
    size_t N = 0;
    for (size_t i = 2; i <= n; i++)
        if (prime[i] != 0)
            N++;
    std::cout << N << '\n';
}

Answer 1

#pragma omp parallel for 之后的循环必须是 规范形式待确认。在您的情况下，问题在于测试表达式，它应该是以下之一：

var relational-op b
b relational-op var

因此，您必须使用@Yakk建议的内容：计算n的sqrt > 并将其与 i 进行比较：

const size_t max_i=sqrt(n);
#pragma omp parallel for
  for (size_t i = 2; i <= max_i; i++)
  ....