p.storePrevIter();
    rho.storePrevIter();