Finish #3

I made the uncomment for `setSat_p` but forgot for `setSat_T`, this is now corrected.

Projects/Sources/coolpropsolver.cpp

@@ -467,13 +467,9 @@ void CoolPropSolver::setSat_T(double &T, ExternalSaturationProperties *const pro
 		  // Ancillary equations - fast but not 100% consistent with the rest
 		  //properties->dl = state->saturation_ancillary(CoolProp::iDmass,0,CoolProp::iT,T);
 		  //properties->dl = state->saturation_ancillary(CoolProp::iDmolar,0,CoolProp::iT,T)/state->molar_mass();
-
-		  /* TODO:                                                                                                  **
-		  ** Uncoment the `specify_phase` when issue 656 is solved: https://github.com/CoolProp/CoolProp/issues/656 **
-		  ** This will allow for a full state update and a speed enhancement                                        */
 
 		  // At bubble line:
-		  // state->specify_phase(CoolProp::iphase_liquid);
+		  state->specify_phase(CoolProp::iphase_liquid);
 		  state->update(CoolProp::QT_INPUTS,0,T);
 		  properties->Tsat = T;
 		  properties->psat = state->p(); // At bubble line! (mather for pseudo-pure fluids)
@@ -486,7 +482,7 @@ void CoolPropSolver::setSat_T(double &T, ExternalSaturationProperties *const pro
 		  properties->sl    = state->smass();    // Specific entropy at bubble line (for pressure ps)
 
 		  // At dew line:
-		  // state->specify_phase(CoolProp::iphase_gas);
+		  state->specify_phase(CoolProp::iphase_gas);
 		  state->update(CoolProp::QT_INPUTS,1,T);
 		  properties->dv = state->rhomass();
 		  properties->hv = state->hmass();