I'm using healpy.synfast to create maps, but it seems that healpy does not have the "iseed" function (as in here: http://healpix.jpl.nasa.gov/html/facilitiesnode14.htm) which let me define the random seed to be used for the generation of alms from the power spectrum.
Could anyone tell me how to achieve the "iseed" function in healpy? Thanks!
healpy
internally usesnp.random.standard_normal
to generate the real and imaginary components of thealms
, seesphtfunc.py
.Therefore you can use the
numpy.random.seed
function to set the seed, as:before running
synfast
.