1#include <math.h> 2 3double fma(double x, double y, double z) 4{ 5 __asm__ ("fmadd %0, %1, %2, %3" : "=d"(x) : "d"(x), "d"(y), "d"(z)); 6 return x; 7} 8