Calculating the discrete spectrum of the transport operator with arbitrary order anisotropic scattering

We consider here the numerical determination of the discrete spectrum of the neutron transport operator in homogeneous subcritical infinite medium in the case of highly anisotropic scattering. Under mild assumptions on the scattering coefficients, one knows that the discrete spectrum is a finite set whose elements are real with absolute values larger than one. Under these assumptions, we describe a method to compute all eigenvalues for an arbitrary order of anisotropic scattering. The method has been implemented in a C++ library using IEEE double precision arithmetic, allowing for fast calculations. Numerical results are compared to values in literature and to high-precision calculations using MAPLE