标签： c distribution poisson
在我目前的代码中，我从Poisson分布生成两个数字（x和y），然后我从x中减去y。这两个操作几乎占我模拟总计算时间的50％！
我刚刚在wiki上发现两个不相符的Poisson分布的减法是Skellam分布。我希望通过让C从Skellam分布中只抽取一个随机数而不是从两个不同的Poissons分布中抽取两个数来获得计算时间。