healpy synfast：如何定义随机种子

Question

我正在使用healpy.synfast创建地图，但似乎healpy没有“iseed”功能（如此处：http://healpix.jpl.nasa.gov/html/facilitiesnode14.htm），这让我可以定义用于从功率谱中产生施舍。

有谁能告诉我如何在healpy中实现“iseed”功能？谢谢！

Answer 1

healpy在内部使用np.random.standard_normal生成alms的实部和虚部，请参阅sphtfunc.py。

因此，您可以使用numpy.random.seed函数设置种子，如：

numpy.random.seed(1234)

在运行synfast之前。