fma is only available on recent x86_64 cpus and it is much faster than
a software fma, so this should be done with a runtime check, however
that requires more changes, this patch just adds the code so it can be
tested when musl is compiled with -mfma or -mfma4.
--- /dev/null
+#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
--- /dev/null
+#include <math.h>
+
+#if __FMA__
+
+float fmaf(float x, float y, float z)
+{
+ __asm__ ("vfmadd132ss %1, %2, %0" : "+x" (x) : "x" (y), "x" (z));
+ return x;
+}
+
+#elif __FMA4__
+
+float fmaf(float x, float y, float z)
+{
+ __asm__ ("vfmaddss %3, %2, %1, %0" : "=x" (x) : "x" (x), "x" (y), "x" (z));
+ return x;
+}
+
+#else
+
+#include "../fmaf.c"
+
+#endif
--- /dev/null
+#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
--- /dev/null
+#include <math.h>
+
+#if __FMA__
+
+float fmaf(float x, float y, float z)
+{
+ __asm__ ("vfmadd132ss %1, %2, %0" : "+x" (x) : "x" (y), "x" (z));
+ return x;
+}
+
+#elif __FMA4__
+
+float fmaf(float x, float y, float z)
+{
+ __asm__ ("vfmaddss %3, %2, %1, %0" : "=x" (x) : "x" (x), "x" (y), "x" (z));
+ return x;
+}
+
+#else
+
+#include "../fmaf.c"
+
+#endif