Calculate matrix profile in debug

Returns profile reduction ratio in debug mode
Fixed solution time for power method
All solvers return residual norm
residual to return_info only
Use logging module
Merge pull request #250 from nonhermitian/steady_updates

Updates to Steady state module