The search logic has been updated to avoid unnecessary additional calls to the equation of state. Additionally, the neutral position is now calculated analyticaly for the PLM and PPM cases where a linearization of the coefficients are used,