Qui ci sono funzioni libere che fanno lo stesso, ma nel primo caso il ciclo non è vettorizzato ma negli altri casi lo è. Perché?Perché la vettorizzazione si comporta in modo diverso per quasi lo stesso codice?
#include <vector>
typedef std::vector<double> Vec;
void update(Vec& a, const Vec& b, double gamma) {
const size_t K = a.size();
for (size_t i = 0; i < K; ++i) { // not vectorized
a[i] = b[i] * gamma - a[i];
}
}
void update2(Vec& a, const Vec& b, double gamma) {
for (size_t i = 0; i < a.size(); ++i) { // vectorized
a[i] = b[i] * gamma - a[i];
}
}
void update3(Vec& a, size_t K, const Vec& b, double gamma) {
for (size_t i = 0; i < K; ++i) { // vectorized
a[i] = b[i] * gamma - a[i];
}
}
int main(int argc, const char* argv[]) {
Vec a(argc), b;
update(a, b, 0.5);
update2(a, b, 0.5);
update3(a, a.size(), b, 0.5);
return 0;
}
messaggi rilevanti dal compilatore (VS2013):
1> c:\home\dima\trws\trw_s-v1.3\trws\test\vector.cpp(7) : info C5002: loop not vectorized due to reason '1200'
1> c:\home\dima\trws\trw_s-v1.3\trws\test\vector.cpp(13) : info C5001: loop vectorized
1> c:\home\dima\trws\trw_s-v1.3\trws\test\vector.cpp(19) : info C5001: loop vectorized
Dal commento di @tony
Motivo 1200: "Loop contiene dipendenze di dati di loop-effettuata che impediscono vettorializzazione. Diverse iterazioni del loop interferiscono con ciascuna in modo tale che la vettorizzazione del loop produce risposte errate e l'auto-vettorizzatore c annot dimostra a se stesso che non ci sono tali dipendenze ". source
Provare un compilatore diverso? Altri (gcc e clang) rendono tutte le 3 funzioni. –
Qual è la "ragione 1200" documentata come? –
Motivo 1200: "Il loop contiene dipendenze di dati trasportate da loop che impediscono la vettorizzazione.Eterni iterazioni del loop interferiscono tra loro in modo tale che la vettorizzazione del loop produce risposte errate e l'auto-vettore non può dimostrare a se stesso che non ci sono tali dati dipendenze ". [fonte] (https://msdn.microsoft.com/en-us/library/jj658585.aspx#BKMK_ReasonCode120x) – Tony