Fix Vector-AdjVector product (#155)

Author: odow

src/implementations/LinearAlgebra.jl

@@ -198,6 +198,13 @@ function promote_array_mul(
     return Matrix{promote_sum_mul(S, T)}
 end
 
+function promote_array_mul(
+    ::Type{<:AbstractVector{S}},
+    ::Type{<:LinearAlgebra.Adjoint{T,<:AbstractVector{T}}},
+) where {S,T}
+    return Matrix{promote_sum_mul(S, T)}
+end
+
 function promote_array_mul(
     ::Type{<:AbstractMatrix{S}},
     ::Type{<:AbstractVector{T}},

test/matmul.jl

@@ -312,3 +312,10 @@ Base.convert(::Type{Base.OneTo}, x::Issue65OneTo) = Base.OneTo(x.N)
     @test MA.operate(*, A, x[:, 1]) == x * x[:, 1]
     @test MA.operate(*, A, x) == x * x
 end
+
+@testset "Issue 154" begin
+    X = big.([1 2; 3 4])
+    c = big.([5, 6])
+    MA.operate!!(MA.add_mul, X, c, c')
+    @test X == [26 32; 33 40]
+end