Euclidian norm normalization, using broadcasting.

Python source code: normalize.py

import numpy as np

def normalize(v):
    return v / np.sqrt(np.sum(v**2, axis=-1))[..., None]