1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
#include <math.h>
#if __FMA__
double fma(double x, double y, double z)
{
__asm__ ("vfmadd132sd %1, %2, %0" : "+x" (x) : "x" (y), "x" (z));
return x;
}
#elif __FMA4__
double fma(double x, double y, double z)
{
__asm__ ("vfmaddsd %3, %2, %1, %0" : "=x" (x) : "x" (x), "x" (y), "x" (z));
return x;
}
#else
#include "../fma.c"
#endif
|