  There are several physically motivated density matrix functionals in the
literature, built from the knowledge of the natural orbitals and the occupation
numbers of the one-body reduced density matrix. With the help of the equivalent
phase-space formalism, we thoroughly test some of the most popular of those
functionals on a completely solvable model.
