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

[Andreas Noack]

I'm unable to reproduce this on 0.3.5, 0.3.7pre or master. How did you
construct A?

2015-02-18 9:17 GMT-05:00 Tamas Papp <tkp...@gmail.com>:

> 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
>