docs: public static double jensenShannonDivergence(double[] a, double[] b) {
      double[] average = pairwiseAdd(a, b);
      multiplyInPlace(average, .5);
      return .5 * klDivergence(a, average) + .5 * klDivergence(b, average);
    }