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.