Bugzilla – Bug 1271
Add view to nonZeros of SparseMatrix
Last modified: 2016-08-29 10:14:21 UTC
The ability to call a method that zeros the coefficients in a sparse matrix.
Currently the setZero() method zeros all allocated memory but it would be nice to set the Aij coefficients to zero, while keeping the row and column arrays allocated and populated.
Exposing a dense "view" on the non-zeros would solve this issue, e.g.:
while enabling many more efficient manipulations on the non zero coefficients. I've suggested this solution several times, e.g. in bug 1104.
Implementing it for compressed matrices is just a matter of returning a 1D dense Eigen::Map. For un-compressed ones, this is much more tricky because we would need to return a list of dense Map for which random access would be extremely inefficient.
I like the idea of providing a nonZerosView().
If it is to complicated for un-compressed matrices, I think simply restricting this feature to compressed matrices would be ok.
Minor question: Should the return value be a vector or an array? I'm tending towards array, since nonZerosView() would behave more like independent scalars.
I agree with the two suggestions (compressed-only + array).
Done for compressed mode: https://bitbucket.org/eigen/eigen/commits/614aeef
Supporting non compressed mode would require a more general expression to concatenate an arbitrary number of vectors/matrices.
Fantastic, suits my purposes. Many thanks!