Algorithms for the fast and exact computation of Wigner matrices are described and their application to a fast and massively parallel 4 pi convolution code between a beam and a sky is also presented.