WIP: Updates to mcsolver

- Attempt to optimize the mcsolve algorithms before an eventual move to cython.
Is this really correct? the cy_expect_rho_vec_csr is for the operation trace(matrix * matrix) and is intended for the case when calculating expectation values for density matrices in vector for. I don't think this will work on a ket input, but i might be missing something here..