Я пытался векторизовать цикл, содержащий использование функции "pow" в математической библиотеке. Я знаю, что компилятор Intel поддерживает использование "pow" для инструкций sse, но я не могу заставить его работать с gcc (я думаю). Это тот случай, с которым я работаю:
int main(){
int i=0;
float a[256],
b[256];
float x= 2.3;
for (i =0 ; i<256; i++){
a[i]=1.5;
}
for (i=0; i<256; i++){
b[i]=pow(a[i],x);
}
for (i=0; i<256; i++){
b[i]=a[i]*a[i];
}
return 0;
}
Я компилирую со следующим:
gcc -O3 -Wall -ftree-vectorize -msse2 -ftree-vectorizer-verbose=5 code.c -o runthis
Это относится к os X 10.5.8 с использованием gcc версии 4.2 (я также использовал 4.5 и не мог сказать, что он что-то векторизовал - поскольку он вообще ничего не выводил). Похоже, что ни одна из векторов в векторе не существует - есть ли проблема выделения или какая-либо другая проблема, которую мне нужно использовать? Если я напишу одну из циклов как функцию, я получу немного более подробный вывод (код):
void pow2(float *a, float * b, int n) {
int i;
for (i=0; i<n; i++){
b[i]=a[i]*a[i];
}
}
(используя подробный вывод уровня 7):
note: not vectorized: can't determine dependence between *D.2878_13 and *D.2877_8
bad data dependence.
Я просмотрел страницу gcc auto-vectorization, но это не помогло. Если невозможно использовать pow в версии gcc, где бы я мог найти ресурс для выполнения функции, эквивалентной pow (в основном я имею дело с целыми степенями).
Отредактируйте, так что я просто копался в другом источнике - как он векторизовал это?!:
void array_op(double * d,int len,double value,void (*f)(double*,double*) ) {
for ( int i = 0; i < len; i++ ){
f(&d[i],&value);
}
};
Соответствующий вывод gcc:
note: Profitability threshold is 3 loop iterations.
note: LOOP VECTORIZED.
Ну, теперь я в недоумении - 'd' и 'value' изменяются функцией, которую gcc не знает - странно? Возможно, мне нужно проверить эту часть немного более тщательно, чтобы убедиться, что результаты верны для векторизованной части. Все еще ищете векторную математическую библиотеку - почему нет открытых исходных текстов?