Fix product of Complex with Hermitian (#265)

src/dispatch.jl

@@ -484,13 +484,6 @@ function Base.:*(α::Number, A::LinearAlgebra.Symmetric{<:AbstractMutable})
     return LinearAlgebra.Symmetric(B, c)
 end
 
-function Base.:*(α::Number, A::LinearAlgebra.Hermitian{<:AbstractMutable})
-    c = LinearAlgebra.sym_uplo(A.uplo)
-    B = c == :U ? _mult_upper(α, A) : _mult_lower(α, A)
-    return LinearAlgebra.Hermitian(B, c)
-end
-
-# Fix ambiguity identified by Aqua.jl.
 function Base.:*(α::Real, A::LinearAlgebra.Hermitian{<:AbstractMutable})
     c = LinearAlgebra.sym_uplo(A.uplo)
     B = c == :U ? _mult_upper(α, A) : _mult_lower(α, A)

test/dispatch.jl

@@ -55,3 +55,11 @@ end
     @test all(MA.isequal_canonical.(2 * LinearAlgebra.Hermitian(A, :L), D))
     @test all(MA.isequal_canonical.(2 * LinearAlgebra.Hermitian(A, :U), D))
 end
+
+@testset "*(::Complex, ::Hermitian)" begin
+    A = BigInt[1 2; 2 3]
+    B = LinearAlgebra.Hermitian(DummyBigInt.(A))
+    C = 2im * A
+    @test 2im * B == C
+    @test C isa Matrix{Complex{BigInt}}
+end

test/dummy.jl

@@ -171,3 +171,5 @@ end
 function Base.:^(x::DummyBigInt, y::Union{Integer,UniformScaling{<:Integer}})
     return DummyBigInt(x.data^y)
 end
+
+Base.:*(x::Complex, y::DummyBigInt) = x * y.data