improve \ override

Author: alyst

src/dual.jl

@@ -837,14 +837,15 @@ function _map_dual_components!(fvalue!, fpartial!, y::AbstractArray{DT}, x::Abst
     return y
 end
 
-function Base.:\(m::Union{LowerTriangular{<:LinearAlgebra.BlasFloat},
-                          UpperTriangular{<:LinearAlgebra.BlasFloat}},
-                 x::StridedVector{<:Dual})
-    T = valtype(eltype(x))
-    res = copy(x)
-    ldiv!(m, reinterpret(reshape, T, res)')
-    return res
-end
+LinearAlgebra.ldiv!(y::StridedVector{T},
+                    m::Union{LowerTriangular{<:LinearAlgebra.BlasFloat},
+                             UpperTriangular{<:LinearAlgebra.BlasFloat}},
+                    x::StridedVector{T}) where T <: Dual =
+    (ldiv!(reinterpret(reshape, valtype(T), y)', m, reinterpret(reshape, valtype(T), x)'); y)
+
+Base.:\(m::Union{LowerTriangular{<:LinearAlgebra.BlasFloat},
+                 UpperTriangular{<:LinearAlgebra.BlasFloat}},
+        x::StridedVector{<:Dual}) = ldiv!(similar(x), m, x)
 
 for MT in (StridedMatrix{<:LinearAlgebra.BlasFloat},
            LowerTriangular{<:LinearAlgebra.BlasFloat},