使用模板类型减少OpenMP

Question

template <typename T, std::size_t N>
static T sum(const std::array<T, N>& a)
{
    T result;

    // type of result (T) is not determined when pre-process?
    #pragma omp parallel for reduction(+: result)
    for(int i = 0; i < static_cast<int>(N); i++)
    {
        result += a[i];
    }
    return result;
}

我可以用MSVC和gcc编译并运行上面的代码。 是的，这太棒了！

但我的问题是在代码评论中; “因为在预处理'#pragma'时没有确定结果类型（T），编译器如何验证结果类型是否适合OpenMP减少？” 。

我确定如果T = double，那就没关系，如果T = std :: string则是NG，但是预处理器如何知道T的类型？

我记得很久以前我用一些小的c ++编译器无法编译上面的代码。

让我问一下在C ++ / OpenMP规范的上下文中哪种行为（可编译或不可编译）是正确的。

Answer 1

它未指定 （对于OpenMP 3.0或更高版本）或未定义 （对于OpenMP 2.5）

reduction是数据共享属性子句之一， OpenMP应用程序接口版本2.5说：

2.8.3数据共享属性子句
---- C / C ++ ----
如果数据共享属性子句中引用的变量具有从模板派生的类型，并且程序中没有对该变量的其他引用，则与该变量相关的任何行为都是未定义的。
---- C / C ++ ----

OpenMP应用程序接口版本3.0说：

2.9.3数据共享属性子句
---- C / C ++ ----
如果数据共享属性子句中引用的变量具有从模板派生的类型，并且程序中没有对该变量的其他引用，则未指定与该变量相关的任何行为。
---- C / C ++ ----