The StatsCalculator interface defined a ``divergence_matrix`` function. Add this back in and test. It is possible to get all pairs easily enough using ``itertools.combinations()``, but this seems like a useful thing to have.