宇宙大尺度结构数值模拟中的短程力

本文提出了宇宙大尺度结构P³M法数值模拟中短程力计算的改进方法。其改进主要在于对质点分块重新编号,从而大大减少检索所需时间。为了检验该算法的有效性,我们和国际上著名的Efstathiou等人的算法进行了比较,结果表明,在相同的初始条件下,在同类型机器上运算,我们的计算短程力算法要比他们的快一倍到五倍。