lambda = randn(d, k);
psi = diag(s);
psi = diag(psi .* (.5+rand(size(psi))));