PARALLEL.RU

Может кто подскажет в чем дело.
Имеется следующее.

inline static Type Vec_Vec(const mtl::dense1D< Type > &A,const mtl::dense1D< Type > &B){
Type n = A[0]*B[0];
#pragma omp parallel for reduction(+:n)
for(int_type i = 1; i < B.size(); i++)
n += A[i] * B[i];
return(n);
}

Type может быть: double либо complex

при компиляции (Visual Studio 2008) возникает ошибка:
C3031 'n': variable in reduction clause must have scalar arithmetic type

Видимо, #pragma omp проверяет типы в момент компиляции, а не в момент инстанцирования. А до инстанцирования тип неизвестен, вот он и ругается. Так что стоит попробовать отказаться от шаблона, и сделать пару перегруженных функций.

попробовал.не работает.на функцию с данными complex ругается( и опять именно на reduction

А, ну да, все верно. Он же пишет, что хочет scalar type. Так что complex - не пойдет. И вроде как по стандарту он даже почти прав:
The operator specified in a reduction clause cannot be overloaded with respect to the list items that appear in that clause.
Так что с complex видимо не выйдет.

cпасибо.
Странно, то что с помощью прагмы удается, складывать вектора, умножать матрицу на вектор и не важно комплексные они или нет (вектора, матрицы), а вот именно reduction нельзя использовать при комплексных данных.

Ничего странного. Если вы имели в виду, что можно распараллелить цикл, внутри которого производятся операции с комплексными числами. Эта операция происходит внутри цикла с полном соответствии с синтаксисом C++, а в случае reduction это уже операция внутри OpenMP и зависит от реализации.

Можно попробовать обойти это ограничение, если в конце цикла присваивать Re и Im в отдельные переменные и редуцировать их. Правда в этом случае придётся делать синхронизацию в конце цикла (не получится с nowait) и это подойдёт только для сложения, а не для умножения.

спасибо. видимо. вы абсолютно правы.


тоже про это думал, но похоже выеденного яйца стоить не будет.

Оператор, указанный в пункте сокращение не может быть перегружен по отношению к пунктам списка, которые появляются в этой статье.
Так что с complex видимо не выйдет.

_________________
We offer best quality pmp certification test papers and test-king ielts exam preparation dumps materials. You can get our 100% guaranteed questions coloradomesa.edu to help you in passing the real exam. flsouthern.edu