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}