new program_state and fixpointiterator with makeLinearSolver, FunctionalFactory, NonlinearSolver working as intended
robust implementation of eigenvalue computation, prevents complex eigenvalues due to round-off errors