Hi Christopher, same problem here (0.3.0-prerelease+495 (2013-12-15 17:39 UTC))! The problem arises already in sparsevec:
julia> sv=sparsevec([3,5,7],[0.1,0.0,3.2],4) 4x1 sparse matrix with 3 Float64 nonzeros: [3, 1] = 0.1 [5, 1] = 0.0 [7, 1] = 3.2 which shouldn't be 4x1. Interestingly sparse doesn't allow to construct a matrix in this way julia> sv=sparse([3,5,7],[1,1,1],[0.1,0.0,3.2],4,1) ERROR: BoundsError() in sparse at sparse/csparse.jl:30 in sparse at sparse/sparsematrix.jl:266 Seems the matrices are constructed in different ways in the two cases ... I guess you should file an issue at https://github.com/JuliaLang/julia/issues Cheers, Alex. On Wednesday, 18 December 2013 21:22:53 UTC+1, Christopher Rinderspacher wrote: > > I'm using: > 0.3.0-prerelease+407 (2013-12-12 15:50 UTC) > Commit 48e350f* > > As per the subject I've found that: > > sv=sparsevec([3,5,7],[0.1,0.0,3.2],4) > sv' > > segfaults. There is no such error for: > > sv=sparse([3,3,5,7],[1,2,4,4],[0.1,0.0,3.2,7.7],8,5) > sv' > > I presume boundaries aren't checked properly. > > Sincerely, > > Christopher >
