非定常输运问题适应于消息传递并行编程环境的香农熵计算方法?

在多计算步的非定常输运问题的蒙特卡罗模拟中，为自动调整每一步的样本数以获得较高的计算效率，可以有多种准则。一种可选的方法是在每一步每隔若干样本监测一次系统中未死亡粒子属性分布对应的香农熵的收敛情况以决定何时停止追加样本，此种方法需要在每一步频繁计算香农熵值。由于在MPI消息传递并行编程环境下香农熵的经典计算方法必须广播大量的数据，导致每一步的计算时间随香农熵计算频率的提高而快速增大，这显然是不能满足实际需求的。本文提出了一种适应于消息传递并行编程环境的香农熵计算新方法，该方法计算得到的香农熵值并不等价于经典方法，但二者之间的差别会随着样本数的增加而趋于零。新方法的最大优势是高频计算香农熵值的时间代价大为降低，为最终实现基于香农熵收敛判断的每步样本数的自动调整奠定了必要的基础。