dsq = 1. / N * np.ones((N, D)) * dvar