Sparse spectral-tau method for the three-dimensional helically reduced wave equation on two-center domains