docs: private void expandPsiMinusPsiSum(Vec input, double sum, Vec output)
      {
          double psiSum = FastMath.digamma(sum);
          for(int i = 0; i < input.length(); i++)
              output.set(i, FastMath.digamma(input.get(i))-psiSum);
      }