def resample_omega(self, x):
        """
        Resample omega from its conditional Polya-gamma distribution
        :return:
        """
        assert x.ndim != 2
        N = N_vec(x)

        #  Sum the N's (i.e. the b's in the denominator)
        NN = N.sum(0).astype(np.float)
        ppg.pgdrawvpar(self.ppgs, NN, self.psi, self.omega)
dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py/PGMultinomial.resample_omega/operator `==` -> `!=`