A Scalable Parallel Poisson Solver in Three Dimensions with Infinite-Domain Boundary Conditions