Implements pressure and stress correction for multiphase model with density contrast. Pressure is regularized and hence different from the form used by Safari and Fakahri. Surface tension is not implemented.