better internal support for working with sparse matrices when calculating propagators

Work in progress. Direct superoperator support plan as well.