Workshop Programming of I will describe results from statistical studies of binary black holes in the post-Newtonian approximation using CUDA. The problem is embarrassing parallel and arithmetically intensive enough to lend itself to the efficient use of multiple GPUs. The implementation is simple yet the scientific payoff turned out to be significant, with new insights into the binary black hole problem.