dev = np.sqrt(np.diag(Sc))  # std