range(rank, scan.nces, size))
## Noise seeds
state_for_noise = np.random.RandomState(params.array_noise_seed)
seeds_for_noise = state_for_noise.randint(0, 1e6, scan.nces)
for pos_CES, CESnumber in enumerate(range(rank, scan.nces, size)):
if params.verbose:
print("Proc [{}] with seeds ".format(rank),
seeds_for_noise[CESnumber], seeds_for_noise)
tod = TimeOrderedDataPairDiff(
inst,
scan,
sky_in,
CESnumber=CESnumber,
projection=params.projection,
nside_out=params.nside_out,
pixel_size=params.pixel_size,
width=params.width,
array_noise_level=params.array_noise_level,
array_noise_seed=seeds_for_noise[CESnumber],
mapping_perpair=params.mapping_perpair)
## Initialise map containers for each processor
if pos_CES == 0:
sky_out_tot = OutputSkyMap(projection=tod.projection,
nside=tod.nside_out,
obspix=tod.obspix,
npixsky=tod.npixsky,
pixel_size=tod.pixel_size)
## Scan input map to get TODs
squid_ids = inst.focal_plane.get_indices('Sq')
bolo_ids = inst.focal_plane.bolo_index_in_squid
## Noise seeds
state_for_noise = np.random.RandomState(params.array_noise_seed)
seeds_for_noise = np.random.randint(0, 1e6, scan.nces)
for pos_CES, CESnumber in enumerate(range(rank, scan.nces, size)):
if params.verbose:
print("Proc [{}] with seeds ".format(rank),
seeds_for_noise[CESnumber], seeds_for_noise)
tod = TimeOrderedDataPairDiff(
inst,
scan,
sky_in,
CESnumber=CESnumber,
projection=params.projection,
nside_out=params.nside_out,
pixel_size=params.pixel_size,
width=params.width,
array_noise_level=params.array_noise_level,
array_noise_seed=seeds_for_noise[CESnumber],
mapping_perpair=False)
## Initialise map containers for each processor
if pos_CES == 0:
sky_out_tot = OutputSkyMap(projection=tod.projection,
nside=tod.nside_out,
obspix=tod.obspix,
npixsky=tod.npixsky,
pixel_size=tod.pixel_size)
## Scan input map to get TODs