We discuss iterative methods for computing criticality in nuclear reactors. In general this requires the solution of a generalised eigenvalue problem for an unsymmetric integro-differential operator in 6 independent variables, modelling transport, scattering and fission, where the dependent variable is the neutron angular flux. In engineering practice this problem is often solved iteratively, u...