Tabella di ricerca utilizzando SIMD

Question

Ho una grande funzione di elaborazione dei pixel che attualmente sto tentando di ottimizzare utilizzando le funzioni intrinseche.

Essendo un novizio SSE, non sono sicuro di come affrontare la parte del codice che coinvolge le tabelle di ricerca.

In sostanza, sto cercando di vettorizzare il seguente codice alla vaniglia C++:

//outside loop 
const float LUT_RATIO = 1000.0F; 

//in loop 
float v = ... //input value 
v = myLookupTable[static_cast<int>(v * LUT_RATIO)]; 

Quello che sto cercando:

//outside loop 
const __m128 LUT_RATIO = _mm_set1_ps(1000.0F); 

//in loop 
__m128 v = _mm_set_ps(v1, v2, v3, v4); //input values 
__m128i vI = _mm_cvtps_epi32(_mm_mul_ps(v, LUT_RATIO)); //multiply and convert to integers 
v = ??? // how to get vI indices of myLookupTable? 

edit: Ildjarn fa un punto che richiede un chiarimento da parte mia. Non sto cercando di ottenere l'accelerazione per il codice della tabella di ricerca, sto semplicemente cercando di evitare di archiviare i registri su float appositamente per fare la ricerca, poiché questa parte è inserita tra 2 altre parti che potrebbero teoricamente beneficiare di SSE.

Answer 1

Se è possibile attendere fino al prossimo anno, le CPU Haswell di Intel disporranno di AVX2 che include istruzioni per i carichi raccolti. Questo ti permette di fare ad es. 8 ricerche LUT parallele in un'unica istruzione (vedere ad esempio VGATHERDPS). A parte questo, sei sfortunato, a meno che le tue LUT siano piuttosto piccole (ad esempio 16 elementi), nel qual caso puoi usare PSHUFB.