*julia> **A = sprand(5,5,2/5)*
*5x5 sparse matrix with 9 Float64 entries:* * [3, 1] = 0.640396* * [3, 2] = 0.873544* * [3, 3] = 0.278582* * [4, 3] = 0.751313* * [1, 4] = 0.644883* * [2, 4] = 0.0778264* * [3, 4] = 0.848185* * [2, 5] = 0.0856352* * [3, 5] = 0.553206* *julia> **B = sprand(5,5,2/5)* *5x5 sparse matrix with 9 Float64 entries:* * [2, 1] = 0.404673* * [3, 1] = 0.736787* * [4, 1] = 0.953803* * [1, 2] = 0.0951856* * [4, 2] = 0.519675* * [1, 3] = 0.0135403* * [4, 3] = 0.303399* * [1, 5] = 0.702557* * [4, 5] = 0.596537* *julia> **C = max(A,B)* *5x5 Array{Float64,2}:* * 0.0 0.0951856 0.0135403 0.644883 0.702557* * 0.404673 0.0 0.0 0.0778264 0.0856352* * 0.736787 0.873544 0.278582 0.848185 0.553206* * 0.953803 0.519675 0.751313 0.0 0.596537* * 0.0 0.0 0.0 0.0 0.0 * Isn't this an undesired behavior of max() when called on 2 sparse matrices? Shouldn't the output be sparse as well?