fixing grad() and div() calls in pressure sovler; replacing inhomo_solver::psi() with solver_common::state()