[WIP] Reference density state

This is a first idea of how a reference density implementation for different approximations could work.