[Numpy-discussion] avoiding loops when downsampling arrays

[Sturla Molden]

>           for i in range(m):
>               for j in range(n):
>                   found[i//4,j//4] = cond(x[i,j])
>

Blah, that should be

     found[i//4,j//4] |= cond(x[i,j])


Sturla