Dissipation Support for sparse matrices Implement more efficient Liouvillian evolution with action of matrix exponential on a state, instead of explicit matrix exponentiation. General Unit testing