Suma prefiksów SIMD w procesorze Intel
Muszę zaimplementować algorytm sumy przedrostkowej i będzie on musiał być tak szybki, jak to możliwe. Dawny:
[3, 1, 7, 0, 4, 1, 6, 3]
should give
[3, 4, 11, 11, 15, 16, 22, 25]
Czy można to zrobić za pomocą instrukcji procesora SSE / mmx / SIMD?
Moim pierwszym pomysłem jest sumowanie każdej pary równolegle rekurencyjnie, aż cała suma zostanie obliczona jak poniżej!
//in parallel do
for (int i = 0; i<z.length; i++){
z[i] = x[i<<1] + x[(i<<1)+1];
}
Aby algorytm stał się bardziej przejrzysty, „z” nie jest końcowym wyjściem
ale zamiast tego służy do obliczania wydajności
int[] w = computePrefixSum(z);
for (int i = 1; i<ouput.length; i++){
ouput[i] = (i%2==0) ? (x[i] + ouput[i-1]) : w[(i-1)>>1];
}