[WIP] Reference density state

This is a first idea of how a reference density implementation for different approximations could work.
Double semicolon at the end of the line.