We report the solution of extremely large integral-equation problems involving electromagnetic scattering from conducting bodies. By orchestrating diverse activities, such as the multilevel fast multipole algorithm, iterative methods, preconditioning techniques, and parallelization, we are able to solve scattering problems that are discretized with tens of millions of unknowns. Specifically, we report the solution of a closed geometry containing 42 million unknowns and an open geometry containing 20 million unknowns, which are the largest problems of their classes, to the best of our knowledge.