“…The problem has been solved in two different ways, either by direct construction of valid density matrices from wavefunctions 3 or by an algebraic method 4 . The complexity is large in both methods, as presently implemented, and at least the algebraic method is so far computationally very demanding 3,5 . The wavefunction method is very direct, but gives no information about whether density matrices and densities may or may not be added.…”