Generalize fast colored iteration on banded matrix to array output

The function already uses 2D indexing so it's fast and safe for array output too, it just needs to make sure to get the sparsity information from sparsity (which it should have done in the first place, it just assumed the two were the same).