more fixes

Author: dkarrasch

src/linalg.jl

@@ -51,13 +51,18 @@ const DenseViewWrappers{T,S} = Union{AdjOrTrans{T,S}, HermOrSym{T,S}, UpperOrLow
 const QuasiSparseMatrix = Union{SparseMatrixCSCUnion2, DenseViewWrappers{<:Any,<:SparseMatrixCSCUnion2}}
 const QuasiStridedMatrix = Union{StridedMatrix, DenseViewWrappers{<:Any,<:StridedMatrix}}
 
-matop_dest(::typeof(*), A::QuasiStridedMatrix, b::AbstractSparseVector) = Vector{TS}(undef, size(A, 1))
-matop_dest(::typeof(*), A, B::QuasiSparseMatrix) = similar(A, TS, (size(A, 1), size(B, 2)))
+matop_dest(::typeof(*), A::QuasiStridedMatrix, b::AbstractSparseVector) =
+    Vector{promote_op(matprod, eltype(A), eltype(b))}(undef, size(A, 1))
+matop_dest(::typeof(*), A, B::QuasiSparseMatrix) =
+    similar(A, promote_op(matprod, eltype(A), eltype(B)), (size(A, 1), size(B, 2)))
 # sparse products with banded matrices should return sparse arrays (Diagonal is handled by fallback)
-matop_dest(::typeof(*), ::BiTriSym, B::QuasiSparseMatrix) = similar(B, TS, size(B))
-matop_dest(::typeof(*), ::Diagonal, B::QuasiSparseMatrix) = similar(B, TS, size(B)) # disambiguation with LinearAlgebra
+matop_dest(::typeof(*), A::BiTriSym, B::QuasiSparseMatrix) =
+    similar(B, promote_op(matprod, eltype(A), eltype(B)), size(B))
+# needed for disambiguation with LinearAlgebra
+matop_dest(::typeof(*), A::Diagonal, B::QuasiSparseMatrix) =
+    similar(B, promote_op(matprod, eltype(A), eltype(B)), size(B))
 matop_dest(::typeof(*), A::QuasiSparseMatrix, B::BiTriSym) =
-    similar(A, TS, (size(A, 1), size(B, 2)))
+    similar(A, promote_op(matprod, eltype(A), eltype(B)), (size(A, 1), size(B, 2)))
 
 for op ∈ (:+, :-), Wrapper ∈ (:Hermitian, :Symmetric)
     @eval begin