overload mul!(y::DualMtx, m::FloatMtx, x::DualMtx)

Author: alyst

src/dual.jl

@@ -840,9 +840,10 @@ for MT in (StridedMatrix{<:LinearAlgebra.BlasFloat},
 
     @eval Base.:*(m::$MT, x::StridedVector{<:Dual}) = mul!(similar(x, (size(m, 1),)), m, x)
 
-    @eval Base.:*(m::$MT, x::StridedMatrix{<:Dual}) =
-        _map_dual_components!((y, x) -> mul!(y, m, x), (y, x, _) -> mul!(y, m, x),
-                              similar(x, (size(m, 1), size(x, 2))), x)
+    @eval LinearAlgebra.mul!(y::StridedMatrix{T}, m::$MT, x::StridedMatrix{T}) where T <: Dual =
+        _map_dual_components!((y, x) -> mul!(y, m, x), (y, x, _) -> mul!(y, m, x), y, x)
+
+    @eval Base.:*(m::$MT, x::StridedMatrix{<:Dual}) = mul!(similar(x, (size(m, 1), size(x, 2))), m, x)
 end
 
 ###################