A three-dimensional time-dependent multigroup SPN method for unstructured tetrahedral meshes is developed. The SPN equations are expressed in a canonical form that allows them to be solved using standard diffusion solution techniques in conjunction with source iteration, diffusion-synthetic acceleration, and fission-source acceleration. A computational comparison of this SPN method with an even-parity SN method is given.