`np.var(X, axis=0).sum()` would spare a memory copy.