Revert "explicitly SIMD muladd with duals (#562)"

This reverts commit ab0e239.

Author: KristofferC

src/ForwardDiff.jl

@@ -16,13 +16,6 @@ import SpecialFunctions
 import LogExpFunctions
 import CommonSubexpressions
 
-const SIMDFloat = Union{Float64, Float32}
-const SIMDInt = Union{
-                       Int128, Int64, Int32, Int16, Int8,
-                       UInt128, UInt64, UInt32, UInt16, UInt8,
-                     }
-const SIMDType = Union{SIMDFloat, SIMDInt}
-
 include("prelude.jl")
 include("partials.jl")
 include("dual.jl")

src/dual.jl

@@ -541,16 +541,6 @@ end
 # fma #
 #-----#
 
-@inline function calc_fma_xyz(x::Dual{T,V,N},
-                              y::Dual{T,V,N},
-                              z::Dual{T,V,N}) where {T, V<:SIMDFloat,N}
-    xv, yv, zv = value(x), value(y), value(z)
-    rv = fma(xv, yv, zv)
-    N == 0 && return Dual{T}(rv)
-    xp, yp, zp = Vec(partials(x).values), Vec(partials(y).values), Vec(partials(z).values)
-    parts = Tuple(fma(xv, yp, fma(yv, xp, zp)))
-    Dual{T}(rv, parts)
-end
 @generated function calc_fma_xyz(x::Dual{T,<:Any,N},
                                  y::Dual{T,<:Any,N},
                                  z::Dual{T,<:Any,N}) where {T,N}
@@ -593,16 +583,6 @@ end
 # muladd #
 #--------#
 
-@inline function calc_muladd_xyz(x::Dual{T,V,N},
-                                 y::Dual{T,V,N},
-                                 z::Dual{T,V,N}) where {T, V<:SIMDType,N}
-    xv, yv, zv = value(x), value(y), value(z)
-    rv = muladd(xv, yv, zv)
-    N == 0 && return Dual{T}(rv)
-    xp, yp, zp = Vec(partials(x).values), Vec(partials(y).values), Vec(partials(z).values)
-    parts = Tuple(muladd(xv, yp, muladd(yv, xp, zp)))
-    Dual{T}(rv, parts)
-end
 @generated function calc_muladd_xyz(x::Dual{T,<:Any,N},
                                     y::Dual{T,<:Any,N},
                                     z::Dual{T,<:Any,N}) where {T,N}

src/partials.jl

@@ -205,6 +205,13 @@ end
     return tupexpr(i -> :(rand(V)), N)
 end
 
+
+const SIMDFloat = Union{Float64, Float32}
+const SIMDInt = Union{
+                       Int128, Int64, Int32, Int16, Int8,
+                       UInt128, UInt64, UInt32, UInt16, UInt8,
+                     }
+const SIMDType = Union{SIMDFloat, SIMDInt}
 const NT{N,T} = NTuple{N,T}
 
 # SIMD implementation