[julia-users] sparse and full matrices give different answer

[Tamas Papp]

Is this a bug, or am I missing something obvious?

julia> versioninfo()
Julia Version 0.3.5
Platform Info:
  System: Linux (x86_64-linux-gnu)
  CPU: Intel(R) Core(TM) i7-3520M CPU @ 2.90GHz
  WORD_SIZE: 64
  BLAS: libblas
  LAPACK: liblapack
  LIBM: libopenlibm
  LLVM: libLLVM-3.5
  
julia> A
3x5 sparse matrix with 5 Float64 entries:
        [3, 1]  =  1.0
        [3, 2]  =  1.0
        [2, 3]  =  1.0
        [1, 4]  =  1.0
        [1, 5]  =  1.0

julia> p
5-element Array{Float64,1}:
 0.0      
 0.138946 
 0.0203463
 0.648933 
 0.191775 

julia> full(A * spdiagm(p))
3x5 Array{Float64,2}:
 0.0  0.0       0.0        0.648933  0.0     
 0.0  0.0       0.0203463  0.0       0.191775
 0.0  0.138946  0.0        0.0       0.0     

julia> full(A) * spdiagm(p)
3x5 Array{Float64,2}:
 0.0  0.0       0.0        0.648933  0.191775
 0.0  0.0       0.0203463  0.0       0.0     
 0.0  0.138946  0.0        0.0       0.0     

Best,

Tamas