[Numpy-discussion] Can this be done efficiently with array functions?