Add some specializations (#74)

* Add some specializations

* Add a test for different input eltypes

Author: kshyatt

src/fillalgebra.jl

@@ -99,7 +99,27 @@ end
 *(a::Diagonal, b::Zeros{<:Any,1}) = mult_zeros(a, b)
 *(a::Diagonal, b::Zeros{<:Any,2}) = mult_zeros(a, b)
 
+*(a::Adjoint{T, <:StridedMatrix{T}},   b::Fill{T, 1}) where T = reshape(sum(conj.(parent(a)); dims=1) .* b.value, size(parent(a), 2))
+*(a::Transpose{T, <:StridedMatrix{T}}, b::Fill{T, 1}) where T = reshape(sum(parent(a); dims=1) .* b.value, size(parent(a), 2))
+*(a::StridedMatrix{T}, b::Fill{T, 1}) where T         = reshape(sum(a; dims=2) .* b.value, size(a, 1))
 
+function *(a::Adjoint{T, <:StridedMatrix{T}}, b::Fill{T, 2}) where T
+    fB = similar(parent(a), size(b, 1), size(b, 2))
+    fill!(fB, b.value)
+    return a*fB
+end
+
+function *(a::Transpose{T, <:StridedMatrix{T}}, b::Fill{T, 2}) where T
+    fB = similar(parent(a), size(b, 1), size(b, 2))
+    fill!(fB, b.value)
+    return a*fB
+end
+
+function *(a::StridedMatrix{T}, b::Fill{T, 2}) where T
+    fB = similar(a, size(b, 1), size(b, 2))
+    fill!(fB, b.value)
+    return a*fB
+end
 function *(a::Adjoint{T, <:AbstractVector{T}}, b::Zeros{S, 1}) where {T, S}
     la, lb = length(a), length(b)
     if la ≠ lb

test/runtests.jl

@@ -846,4 +846,23 @@ end
         p in (-Inf, 0, 0.1, 1, 2, 3, Inf)
         @test norm(a,p) ≈ norm(Array(a),p)
     end
-end
+end
+
+@testset "multiplication" begin
+    for T in (Float64, ComplexF64)
+        fv = T == Float64 ? Float64(1.6) : ComplexF64(1.6, 1.3)
+        n  = 10
+        k  = 12
+        m  = 15
+        fillvec = Fill(fv, k)
+        fillmat = Fill(fv, k, m)
+        A  = rand(ComplexF64, n, k)
+        @test A*fillvec ≈ A*Array(fillvec)
+        @test A*fillmat ≈ A*Array(fillmat)
+        A  = rand(ComplexF64, k, n)
+        @test transpose(A)*fillvec ≈ transpose(A)*Array(fillvec)
+        @test transpose(A)*fillmat ≈ transpose(A)*Array(fillmat)
+        @test adjoint(A)*fillvec ≈ adjoint(A)*Array(fillvec)
+        @test adjoint(A)*fillmat ≈ adjoint(A)*Array(fillmat)
+    end
+end