StepMatrixFinder or fast_evaluate() should have the option to generate sparse matrices

Something that quickly becomes desirable when you have large state vectors