Adds dispatch rules for column vector times transpose(row vector)

Fixes issue JuliaArrays#537

Author: zygmuntszpak

src/matrix_multiply.jl

@@ -10,16 +10,18 @@ import LinearAlgebra: BlasFloat, matprod, mul!
 @inline *(A::StaticVector, B::StaticMatrix) = *(reshape(A, Size(Size(A)[1], 1)), B)
 @inline *(A::StaticVector, B::Transpose{<:Any, <:StaticVector}) = _mul(Size(A), Size(B), A, B)
 @inline *(A::StaticVector, B::Adjoint{<:Any, <:StaticVector}) = _mul(Size(A), Size(B), A, B)
+@inline *(A::StaticMatrix, B::Adjoint{<:Any, <:StaticVector}) = *(reshape(A, Size(Size(A)[1],)), B)
 
 @inline mul!(dest::StaticVecOrMat, A::StaticMatrix, B::StaticVector) = _mul!(Size(dest), dest, Size(A), Size(B), A, B)
 @inline mul!(dest::StaticVecOrMat, A::StaticMatrix, B::StaticMatrix) = _mul!(Size(dest), dest, Size(A), Size(B), A, B)
 @inline mul!(dest::StaticVecOrMat, A::StaticVector, B::StaticMatrix) = mul!(dest, reshape(A, Size(Size(A)[1], 1)), B)
 @inline mul!(dest::StaticVecOrMat, A::StaticVector, B::Transpose{<:Any, <:StaticVector}) = _mul!(Size(dest), dest, Size(A), Size(B), A, B)
 @inline mul!(dest::StaticVecOrMat, A::StaticVector, B::Adjoint{<:Any, <:StaticVector}) = _mul!(Size(dest), dest, Size(A), Size(B), A, B)
-
+@inline mul!(dest::StaticVecOrMat, A::StaticMatrix, B::Adjoint{<:Any, <:StaticVector}) = mul!(dest, reshape(A, Size(Size(A)[1],)), B)
 #@inline *{TA<:LinearAlgebra.BlasFloat,Tb}(A::StaticMatrix{TA}, b::StaticVector{Tb})
 
 
+
 # Implementations
 
 @generated function _mul(::Size{sa}, a::StaticMatrix{<:Any, <:Any, Ta}, b::AbstractVector{Tb}) where {sa, Ta, Tb}

test/matrix_multiply.jl

@@ -77,6 +77,20 @@ using StaticArrays, Test, LinearAlgebra
         @test v4 * transpose(v5) === @SMatrix [3+0im 4+0im; 6+0im 8+0im]
     end
 
+    @testset "Column vector-vector" begin
+        cv_array = rand(4,1)
+        rv_array = rand(4)
+        a_array = cv_array * rv_array'
+
+        cv = SMatrix{4,1}(cv_array)
+        rv = SVector{4}(rv_array)
+        @test (cv*rv')::SMatrix ≈ a_array
+
+        cv = MMatrix{4,1}(cv_array)
+        rv = MVector{4}(rv_array)
+        @test (cv*rv')::SMatrix ≈ a_array
+    end
+
     @testset "Matrix-matrix" begin
         m = @SMatrix [1 2; 3 4]
         n = @SMatrix [2 3; 4 5]