Return covariance between X and X2.

def K(self, X, X2, target):
        """Return covariance between X and X2."""
        if (X2 is None) or (X2 is X):
            target[np.diag_indices_from(target)] += self._Kdiag(X)