SparseArrays.jl

SparseArrays.jl copied to clipboard

On Julia 1.6 and master, broadcasting the scaled sum of 5-6 complex sparse arrays sporadically produces eltypes of type Any. Examples:

julia> using SparseArrays

julia> A, B, C, D, E, F = ntuple(_ -> sprand(ComplexF64, 100, 100, 0.1), 6);

julia> (A .+ B .+ C .+ D .+ E .+ F) |> eltype
ComplexF64 (alias for Complex{Float64})

julia> (1 .* A .+ B .+ C .+ D .+ E .+ F) |> eltype
ComplexF64 (alias for Complex{Float64})

julia> (1 .* A .+ 2 .* B .+ C .+ D .+ E .+ F) |> eltype
Any

julia> (1 .* A .+ 2 .* B .+ C .+ D .+ E) |> eltype
Any

julia> (1 .* A .+ 2 .* B .+ C .+ D) |> eltype
ComplexF64 (alias for Complex{Float64})

This does not happen if the matrices have an eltype of Float64. This same issue happened with Julia 1.0-1.2. For Julia 1.3-1.5, the cases with eltype Any had eltype Complex (not ComplexF64).

Duplicate of JuliaLang/SparseArrays.jl#59.