def running_varsum(v_k1, x_k, m_k1, m_k):
    """ Computing the running variance numerator.

    Ref: https://www.johndcook.com/blog/standard_deviation/
    """
    return v_k1 + (x_k - m_k1) * (x_k - m_k)