> weirdly, it seems ``(a**2+b**2)**0.5`` is the winner in terms of speed (although close to `hypot`).

For large arrays, I get that `np.hypot` is a factor two slower than the
direct calculation, which is rather sad!  Clearly, it is not optimized
properly.  May want to change back (eventually).