Todo List

Member Boundaries::getSatDensWater_T (const Temperature &T)

Limit getSatDensWater_T from being called with T <= TB_LOW

Member Boundaries::getSatDensSteam_T (const Temperature &T)

Limit getSatDensSteam_T from being called with T <= TB_LOW

Member Region3::set_pT (SteamCalculator &c, const Pressure &p, const Temperature &T, Num x)

: remove redundent 'x' from Region3::set_pT function?

Member Region4::set_pT (SteamCalculator &c, const Pressure &p, const Temperature &T, Num x)

Work out a way to do this without using a redundant parameter 'p'

Member SteamCalculator::setRegion3_rhoT (const Density &rho, const Temperature &T)

Ensure rho, T are valid for region 3

Member SteamCalculator::quality () const

Check this, make sure only legal values are possible in region 3. There is a bit of confusion in the code about 'region 4 inside region 3'.

Class Units< M, L, T, K, I >

Get the units stuff working with a specific units-tracking flag. Should be possible to say that all pressures should be output in MPa, all temperatures be output in °C, etc etc.

Member LOOP_REG2I

Make gampitau0 not require a steamcalculator argument

Member LOOP_REG2I

Convert LOOP_REG2I etc etc to expression template format a la Blitz++ etc