Use of linear operators on a GPU

[skleinbo]

@shakedregev

Dvec = sparse(rows, cols, vec)

This needs to be a CUSPARSE array, i.e

Dvec = CuSparseMatrixCSC(sparse(rows, cols, vec))

By the way, vec is a symbol in Base. Trying to reassign it should throw an error.