rho as sqrt(eigenMin*eigenMax) for sparse, sqrt(eigenMax) for other formulations