Mean squared jumping distance.

def msjd(theta):
    """Mean squared jumping distance.
    """
    s = 0. 
    for p in theta.dtype.names:
        s += np.sum(np.diff(theta[p], axis=0) ** 2)
    return s