Re: [Numpy-discussion] Slicing slower than matrix multiplication?

Charles R Harris wrote:

Sort of, it's actually (Xi.T*S).T, now that I think of it... I'll see if that is any faster. And if there is a neater way of doing it I'd love to hear about it.

Xi*S[:,newaxis]

Thanks! (Obviously doesn't matter much in terms of performance, as it's only the initialization, but it definitely does make for nicer code.)