Eu preciso implementar um algoritmo de soma de prefixo e precisaria que fosse o mais rápido possível. Ex:

[3, 1,  7,  0,  4,  1,  6,  3]
should give
[3, 4, 11, 11, 15, 16, 22, 25]

Existe uma maneira de fazer isso usando SSE / mmx / SIMD cpu instrução?

Minha primeira idéia é somar cada par em paralelo recursivamente até que todas as quantias tenham sido computadas como abaixo!

       //in parallel do 
        for (int i = 0; i<z.length; i++){
            z[i] = x[i<<1] + x[(i<<1)+1];
        }

Para tornar o algoritmo um pouco mais claro, "z" não é a saída final

mas, em vez disso, usado para calcular a saída

        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];
        }