accelerate rho calculation using PSSS idea
modify rho0 calculation
comments added