Fixed type prediction issue in *(StaticMatrix, Vector)

Author: Andy Ferris

src/matrix_multiply.jl

@@ -176,7 +176,7 @@ end
     s = (sA[1],)
     T = promote_op(matprod, TA, Tb)
 
-    if T == Tb
+    if T == TA
         newtype = similar_type(A, s)
     else
         newtype = similar_type(A, T, s)
@@ -204,7 +204,7 @@ end
     s = (sA[1],)
     T = promote_op(matprod, TA, Tb)
 
-    if T == Tb
+    if T == TA
         newtype = similar_type(A, s)
     else
         newtype = similar_type(A, T, s)

test/matrix_multiply.jl

@@ -27,6 +27,11 @@
         m5 = @SMatrix [1.0 2.0; 3.0 4.0]
         v7 = [1.0, 2.0]
         @test (m5*v7)::SVector ≈ @SVector [5.0, 11.0]
+
+        m6 = @SMatrix Float32[1.0 2.0; 3.0 4.0]
+        v8 = Float64[1.0, 2.0]
+        @test (m6*v8)::SVector{2,Float64} ≈ @SVector [5.0, 11.0]
+
     end
 
     @testset "Vector-matrix" begin