WIP: Updates to steady state module

Added ARPACK eigensolver, updated power method to use RCM, start from
maximally mixed state, and not recalculate LU factors everytime solver
is called.
minor thin, but instead of qeye(...).full() we could directly use np.diag.