I’m playing with making a simple optimization package generic to any type of input. It is already generic to inputs in the form vectors, matrices, number types (precision), units, etc.
But I could not yet make it generic for sparse matrices as inputs, because iterators like
for i in eachindex(x)
iterate over all indexes of the matrices, stored or not.
Is there a generic iterator to iterate over stored-only elements of sparse matrices (specifically
SparseArrays). I have seen that some of that was discussed previously and led to the
ArrayIteration.jl package, but it seems that it was not released?