We devise fast and provably accurate algorithms to transform between an N × N × N Cartesian voxel representation of a three-dimensional function and its expansion into the ball harmonics, that is, the eigenbasis of the Dirichlet Laplacian on the unit... ...