(size=100, sample_count=10000)
| 17 | |
| 18 | |
| 19 | def simulate(size=100, sample_count=10000): |
| 20 | truth = gen_prob_dist(size) |
| 21 | accept, alias = create_alias_table(truth) |
| 22 | |
| 23 | sampled = np.zeros(size) |
| 24 | for _ in range(sample_count): |
| 25 | sampled[alias_sample(accept, alias)] += 1 |
| 26 | return sampled / np.sum(sampled), truth |
| 27 | |
| 28 | |
| 29 | def main(smoke=False, show=True): |
no test coverage detected