lmul!(0, A) should replace NaN with 0.0 #560
Description
Standard arrays have the feature that lmul!(0.0, A) replaces NaNs with 0.0, consistent with BLAS:
julia> lmul!(0.0,fill(NaN,3,3))
3×3 Array{Float64,2}:
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0This is gone when generic_lmul! is called:
julia> lmul!(0,fill(NaN,3,3))
3×3 Array{Float64,2}:
NaN NaN NaN
NaN NaN NaN
NaN NaN NaNI'd propose fixing this inconsistency by making generic_lmul! follow the BLAS convention.