def _geq(self, X): numerator = self._log_proba(X) denominator = logsumexp(numerator, axis=1) denominator = denominator.reshape(len(denominator), 1) return numerator - denominator