def METHOD_NAME(self):
    for fan in ["in", "out"]:
        expected = connect_test_base.get_expected_degrees_bernoulli(self.p, fan, self.N_s, self.N_t)
        pvalues = []
        for i in range(self.stat_dict["n_runs"]):
            connect_test_base.reset_seed(i + 1, self.nr_threads)
            self.setUpNetwork(conn_dict=self.conn_dict, N1=self.N_s, N2=self.N_t)
            degrees = connect_test_base.get_degrees(fan, self.pop1, self.pop2)
            degrees = connect_test_base.gather_data(degrees)
            # degrees = self.comm.gather(degrees, root=0)
            # if self.rank == 0:
            if degrees is not None:
                chi, p = connect_test_base.chi_squared_check(degrees, expected, self.rule)
                pvalues.append(p)
            connect_test_base.mpi_barrier()
        if degrees is not None:
            ks, p = scipy.stats.kstest(pvalues, "uniform")
            self.assertGreater(p, self.stat_dict["alpha2"])