Routines# self_consistency_loop(hamiltonian, kgrid, ...) Self-consistently solves the gap equation for a given Hamiltonian.