fma

Author: serge-sans-paille

include/xsimd/arch/xsimd_altivec.hpp

@@ -301,6 +301,19 @@ namespace xsimd
             }
         }
 
+        // fma
+        template <class A>
+        XSIMD_INLINE batch<float, A> fma(batch<float, A> const& x, batch<float, A> const& y, batch<float, A> const& z, requires_arch<altivec>) noexcept
+        {
+            return vec_madd(x, y, z);
+        }
+
+        template <class A>
+        XSIMD_INLINE batch<double, A> fma(batch<double, A> const& x, batch<double, A> const& y, batch<double, A> const& z, requires_arch<altivec>) noexcept
+        {
+            return vec_madd(x, y, z);
+        }
+
         // eq
         template <class A, class T, class = typename std::enable_if<std::is_scalar<T>::value, void>::type>
         XSIMD_INLINE batch_bool<T, A> eq(batch<T, A> const& self, batch<T, A> const& other, requires_arch<altivec>) noexcept
@@ -459,7 +472,6 @@ namespace xsimd
         XSIMD_INLINE batch<T, A> mul(batch<T, A> const& self, batch<T, A> const& other, requires_arch<altivec>) noexcept
         {
             return self.data * other.data;
-            // return vec_mul(self.data, other.data);
         }
 
         // neg