Despite its deceptively simple functional form, finding highly efficient methods for evaluating interactions involving the Coulomb operator has proven to be a formidable challenge. This difficulty stems from a combination of a slowly decaying tail and a singularity at the coalescence point. However, there have been some successful approaches to this problem that proceed by partitioning the operator into a short-ranged (SR) part that encapsulates the singularity and is readily treated by direct summation, with a complementary long-ranged (LR) term that is amenable to factorization: