docs: protected double computeDCovar(double[] dVarMatrixA, double[] dVarMatrixB, int n) {
      double result = 0.;
      for(int i = 0, c = 0; i < n; i++) {
        for(int j = 0; j < i; j++) {
          result += 2. * dVarMatrixA[c] * dVarMatrixB[c];
          c++;
        }
        // Diagonal entry.
        result += dVarMatrixA[c] * dVarMatrixB[c];
        c++;
      }
      return result / (n * n);
    }