Pourquoi std :: inner_product est-il plus lent que l'implémentation naïve?

Question

Ceci est mon implémentation naïve du produit scalaire:

float simple_dot(int N, float *A, float *B) { 
    float dot = 0; 
    for(int i = 0; i < N; ++i) { 
    dot += A[i] * B[i]; 
    } 

    return dot; 
} 

Et cela est d'utiliser la bibliothèque C++:

float library_dot(int N, float *A, float *B) { 
    return std::inner_product(A, A+N, B, 0); 
} 

J'ai couru une référence (code est ici https://github.com/ijklr/sse), et la version de la bibliothèque est un beaucoup plus lent. Mon drapeau de compilateur est -Ofast -march=native

Answer 1

Vos deux fonctions ne font pas la même chose. L'algorithme utilise un accumulateur dont le type est déduit de la valeur initiale, qui dans votre cas (0) est int. Accumuler des flotteurs dans un int ne prend pas simplement plus de temps que de s'accumuler dans un flotteur, mais produit également un résultat différent. L'équivalent de votre code de boucle brute est d'utiliser la valeur initiale 0.0f, ou de manière équivalente float{}.

(Notez que std::accumulate est très similaire à cet égard.)