I think the 1-D way of assignment is ok, because the most obvious way to get lot's of indices is nonzero(), which is also only 1-D. It would be useful to have a fast mapping function from 1-D to N-D indices and vice versa. I want to agree, that a standard form of doing it 1-D is better than nothing. And it's great, that something along these lines is coming. __Janko