Use matprod_dest to obtain destination

Author: jishnub

src/diagonal.jl

@@ -337,16 +337,18 @@ function _diag_adj_mul(A::AdjOrTransAbsMat, D::Diagonal)
     copy(adj(adj(D) * adj(A)))
 end
 function _diag_adj_mul(A::AdjOrTransAbsMat{<:Number, <:StridedMatrix}, D::Diagonal{<:Number})
-    Ac = copy_similar(A, promote_op(*, eltype(A), eltype(D)))
-    rmul!(Ac, D)
+    TS = promote_op(matprod, eltype(A), eltype(D))
+    C = matprod_dest(A, D, TS)
+    mul!(C, A, D)
 end
 function _diag_adj_mul(D::Diagonal, A::AdjOrTransAbsMat)
     adj = wrapperop(A)
     copy(adj(adj(A) * adj(D)))
 end
 function _diag_adj_mul(D::Diagonal{<:Number}, A::AdjOrTransAbsMat{<:Number, <:StridedMatrix})
-    Ac = copy_similar(A, promote_op(*, eltype(A), eltype(D)))
-    lmul!(D, Ac)
+    TS = promote_op(matprod, eltype(A), eltype(D))
+    C = matprod_dest(D, A, TS)
+    mul!(C, D, A)
 end
 
 function (*)(A::AdjOrTransAbsMat, D::Diagonal)