_mm_fmadd_pd Программа получила сигнал SIGILL, недопустимая инструкция

Я получаю странную ошибку для следующего кода:

#include <assert.h>
#include <stdio.h>
#include <immintrin.h>

inline static double myfma(double x,double y, double z) {
double r; // result
__m128d xx, yy, zz,rr;

xx = _mm_set_sd(x);// xx[0]=x, xx[1]=undefined
yy = _mm_set_sd(y);// yy[0]=y, yy[1]=undefined
zz = _mm_set_sd(z);// zz[0]=z, zz[1]=undefined
r = _mm_cvtsd_f64(_mm_fmadd_pd(xx,yy,zz));

return r;
}

void testfma() {
double x, y, z, res;
x = 1.0;
y = 2.0;
z = 3.0;

res =  myfma(x,y,z);
printf("test: res = x*y + z \n");
printf("    x: %g\n", x);
printf("    y: %g\n", y);
printf("    z: %g\n", z);
assert(res == 5.0);
}int main() {
testfma();
return 0;
}

компилировать код как:

g++ test.cpp -o a.out -std=c++11 -mavx2 -mfma  -march=native -g

когда я запускаю исполняемый файл, я получаю сообщение:

Illegal instruction (core dumped)

Использование GDB для получения дополнительной информации:

gdb ./a.out
(gdb) r
(gdb) r
Starting program: ....

Program received signal SIGILL, Illegal instruction.
0x000000000040067d in _mm_fmadd_pd(double __vector(2), double __vector(2), double __vector(2)) (__C=..., __B=..., __A=...)
at /usr/lib/gcc/x86_64-linux-gnu/5/include/fmaintrin.h:42
42                                                 (__v2df)__C);

Однако при использовании valgrind следующим образом:

valgrind ./a.out
==9825== Memcheck, a memory error detector
==9825== Copyright (C) 2002-2015, and GNU GPL'd, by Julian Seward et al.
==9825== Using Valgrind-3.11.0 and LibVEX; rerun with -h for copyright

info
==9825== Command: ./helios.x
==9825==
test: res = x*y + z
x: 1
y: 2
z: 3
res: 5
==9825==
==9825== HEAP SUMMARY:
==9825==     in use at exit: 0 bytes in 0 blocks
==9825==   total heap usage: 1 allocs, 1 frees, 1,024 bytes allocated
==9825==
==9825== All heap blocks were freed -- no leaks are possible
==9825==
==9825== For counts of detected and suppressed errors, rerun with: -v
==9825== ERROR SUMMARY: 0 errors from 0 contexts (suppressed: 0 from 0)

Программа, кажется, работает. Что мне здесь не хватает? Как я могу использовать _mm_fmadd_pd надежным способом? Можно ли заставить пример работать независимо от того, работает ли он на процессоре Intel или AMD? Можно ли сделать его скомпилировать независимо от G + или ICPC?

0

Решение

Я предполагаю, что ваш процессор не поддерживает инструкции FMA. Причина, по которой он не терпит неудачу при использовании valgrind, заключается в том, что valgrind может эмулировать определенные инструкции.

Вы можете рассмотреть возможность использования std::fma если вы хотите только SISD. С gcc генерирует встроенную инструкцию FMA, но если вы скомпилируете для не-FMA цели, то это будет отступить к реализации без FMA.

0

Другие решения

Других решений пока нет …

По вопросам рекламы [email protected]