worked a bit on handling of density matrices
latest work