#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