Почему _mm256_load_pd скомпилирован в MOVUPD вместо MOVAPD?

Почему следующий код приводит к невыровненным инструкциям AVX (MOVUPD вместо MOVAPD)? Я скомпилировал это на Visual Studio 2015. Как я могу сказать компилятору, что мои данные действительно выровнены?Почему _mm256_load_pd скомпилирован в MOVUPD вместо MOVAPD?

const size_t ALIGN_SIZE = 64; 
    const size_t ARRAY_SIZE = 1024; 

    double __declspec(align(ALIGN_SIZE)) a[ARRAY_SIZE]; 
    double __declspec(align(ALIGN_SIZE)) b[ARRAY_SIZE]; 

    //Calculate the dotproduct 
    __m256d ymm0 = _mm256_set1_pd(0.0); 
    for (int i = 0; i < ARRAY_SIZE; i += 8) 
    { 
     __m256d ymm1 = _mm256_load_pd(a + i); 
     __m256d ymm2 = _mm256_load_pd(b + i); 
     __m256d ymm3 = _mm256_mul_pd(ymm1, ymm2); 
     ymm0 = _mm256_add_pd(ymm3, ymm0); 

     __m256d ymm4 = _mm256_load_pd(a + i + 4); 
     __m256d ymm5 = _mm256_load_pd(b + i + 4); 
     __m256d ymm6 = _mm256_mul_pd(ymm4, ymm5); 
     ymm0 = _mm256_add_pd(ymm6, ymm0); 
    } 



Assembly of the loop: 
00007FF7AC7A1400 vmovupd  ymm1,ymmword ptr [rbp+rax*8+2020h] 
00007FF7AC7A1409 vmulpd  ymm3,ymm1,ymmword ptr [rbp+rax*8+20h] 
00007FF7AC7A140F vmovupd  ymm2,ymmword ptr [rbp+rax*8] 
00007FF7AC7A1415 vmulpd  ymm0,ymm2,ymmword ptr b[rax*8] 
00007FF7AC7A141E add   r8d,8 
00007FF7AC7A1422 movsxd  rax,r8d 
00007FF7AC7A1425 vaddpd  ymm1,ymm0,ymm4 
00007FF7AC7A1429 vaddpd  ymm4,ymm1,ymm3 
00007FF7AC7A142D cmp   rax,400h 
00007FF7AC7A1433 jb   main+70h (07FF7AC7A1400h)

источник

2016-04-19 Laci

Это на самом деле не важно, - нет практически никакого штрафа за использование невыровненных нагрузок с выровненными данными в современных процессорах - сценаристы компилятора, вероятно, просто решил всегда использовать невыровненные нагрузки вместо того, дополнительная логика, чтобы решить, когда использовать выровненные по сравнению с неуравновешенными нагрузками. –

FWIW gcc * et al * делать правильные вещи, поэтому похоже, что это просто специфический для Microsoft причуда. –

@PaulR, зачем использовать слово virtual? Вообще-то я не знаю, что мне известно. 'vmovapd' устарел. 'mvovapd' по-прежнему полезен для nehalem, потому что' movupd' не может сбрасываться с другими операциями, но я сомневаюсь, что это имеет большое значение на практике. Возможно, это то, что вы подразумеваете под виртуальным, но в этом случае это относится только к Nehalem, и этот ответ явно не скомпилирован для Nehalem. –

Существует способ решить эту проблему (это позволяет использовать команду VMOVDQA (аналог MOVAPD) вместо MOVUPD):

inline __m256d Load(const double * p) 
{ 
#ifdef _MSC_VER 
    return _mm256_castsi256_pd(_mm256_load_si256((__m256i*)p)); 
#else 
    return _mm256_load_pd(p); 
#endif 
}

Аналогичное решение для типа поплавка:

inline __m256 Load(const float * p) 
{ 
#ifdef _MSC_VER 
    return _mm256_castsi256_ps(_mm256_load_si256((__m256i*)p)); 
#else 
    return _mm256_load_ps(p); 
#endif 
}

Но для обмана компилятора Visual Studio вам нужно использовать динамически выделенные указатели. В противном случае компилятор не использует инструкцию VMOVDQA.

#include <immintrin.h> 

int main() 
{ 
    float * ps = (float*)_mm_malloc(40, 32); 
    double * pd = (double*)_mm_malloc(40, 32); 

    __m256 s = Load(ps); 
//00007FF79FF81325 vmovdqa  ymm1,ymmword ptr [rdi] 
    __m256d d = Load(pd); 
//00007FF79FF8132F vmovdqa  ymm0,ymmword ptr [rax] 

    _mm256_storeu_ps(ps, s); 
    _mm256_storeu_pd(pd, d); 

    _mm_free(ps); 
    _mm_free(pd); 
}

источник

2017-03-29 15:03:57 ErmIg

Почему _mm256_load_pd скомпилирован в MOVUPD вместо MOVAPD?

ответ

Смежные вопросы