Is there really no better way to write this? It really seems like there ought to be (like, in the sense of, if it isn't possible then this seems like a flaw in numpy, given that mathematically you're doing exactly the same thing in both branches). I guess the obvious thing would be `np.sqrt((lhs * lhs.conj()).sum(1))`, or `np.sqrt(np.square(np.abs(lhs)).sum(1))`, but I suppose that both of those require making two temporaries instead of just one.