The density can be constructed from the (occupied block of) the C matrices, which are accessed by the Ca_subset function and put onto C_left in exactly the lines of code you quoted. It is more efficient not to use the density directly, but to decompose the density into the product of C matrices.
If that doesn’t answer the question, I’ll need you to explain more precisely.