WIP, to be continued later, on a density matrix PDP solver