Hi Julians,
When trying this, I get (Julia 0.3.3):
*julia> **b = sub(a,1:2,1:2)*
*2x2
SubArray{Float64,2,SparseMatrixCSC{Float64,Int64},(UnitRange{Int64},UnitRange{Int64})}:*
* 1.0 0.0*
* 0.0 1.0*
*julia> **b*ones(2)*
*ERROR: `A_mul_B!` has no method matching A_mul_B!(::Array{Float64,1},
::SubArray{Float64,2,SparseMatrixCSC{Float64,Int64},(UnitRange{Int64},UnitRange{Int64})},
::Array{Float64,1})*
* in * at linalg/matmul.jl:72*
Am I doing something wrong? If it's unsupported, does multiplication of
sparse subarrays work in 0.4?
Best regards,
Robert
