UPDATED: iterative stead-state solver now uses spilu preconditioner

Use sparse incomplete LU preconditioner for steady state solver.
Method is now stablized bi-conjugate gradient so that we do not need to
calculate dagger of preconditioner.