我正在使用healpy.synfast创建地图,但似乎healpy没有“iseed”功能(如此处:http://healpix.jpl.nasa.gov/html/facilitiesnode14.htm),这让我可以定义用于从功率谱中产生施舍。
有谁能告诉我如何在healpy中实现“iseed”功能?谢谢!
答案 0 :(得分:2)
healpy
在内部使用np.random.standard_normal
生成alms
的实部和虚部,请参阅sphtfunc.py
。
因此,您可以使用numpy.random.seed
函数设置种子,如:
numpy.random.seed(1234)
在运行synfast
之前。