Residual based solver added for P,u, and v. Working for just P