[Numpy-discussion] array indexing in numarray

[Joost v. E. <joo...@gm...>]

Hi list,

I have two questions on numarray indexing:
1) does anyone know an elegant way to assign values to a non-contiguous
slice. Like:

        In [1]:from numarray import *
        
        In [2]:a = arange(25).resize((5,5))
        
        In [3]:a
        Out[3]:
        array([[ 0,  1,  2,  3,  4],
               [ 5,  6,  7,  8,  9],
               [10, 11, 12, 13, 14],
               [15, 16, 17, 18, 19],
               [20, 21, 22, 23, 24]])
        
        In [4]:i = [1,3]
        
        In [5]:transpose(a[i])[i]=1

	----------------------------------------------------------------
        exceptions.ValueError                                Traceback
        (most recent call last)
        
        ...
        
        ValueError: Invalid destination array: partial indices require
        contiguous non-byteswapped destination
        
2) And second: why didn't you choose to return all combinations of
indices:

        In [8]:a[i,i]
        Out[8]:array([ 6, 18])

Where I, in my humble opinion, would prefer to see all possible
combinations. e.g. a[[1,3],[1,3]] = a[[1,1,3,3],[1,3,1,3]]

Anyone has an idea how to easyly obtain this behaviour?

Greets+thanks,
Joost