dest='seed_diffpointing',
required=True,
type=int,
help='Seed to generate differential pointing.')
if __name__ == "__main__":
"""
Launch the pipeline!
"""
parser = argparse.ArgumentParser(description='MPI version of s4cmb')
addargs(parser)
args = parser.parse_args(None)
## Import parameters from the user parameter file
params = import_string_as_module(args.inifile)
## Overwrite ini file params with params pass to the App directly
for key in args.__dict__.keys():
new = getattr(args, key)
if key in params.__dict__.keys():
old = getattr(params, key)
if new in [None, 'None'] and old is not None:
continue
else:
print("Overwriting {} with new value: {} -> {}".format(
key, old, new))
setattr(params, key, new)
rank = MPI.COMM_WORLD.rank
dest='seed_diffpointing',
required=True,
type=int,
help='Seed to generate differential pointing.')
if __name__ == "__main__":
"""
Launch the pipeline!
"""
parser = argparse.ArgumentParser(description='MPI version of s4cmb')
addargs(parser)
args = parser.parse_args(None)
## Import parameters from the user parameter file
params = import_string_as_module(args.inifile)
## Overwrite ini file params with params pass to the App directly
for key in args.__dict__.keys():
new = getattr(args, key)
if key in params.__dict__.keys():
old = getattr(params, key)
if new in [None, 'None'] and old is not None:
continue
else:
print("Overwriting {} with new value: {} -> {}".format(
key, old, new))
setattr(params, key, new)
rank = MPI.COMM_WORLD.rank
## Other packages needed
import os
import healpy as hp
import numpy as np
from tqdm import *
if __name__ == "__main__":
"""
Launch the pipeline!
"""
## Import parameters from the user parameter file
#params = import_string_as_module(args.inifile)
#params = import_string_as_module("./examples/inifiles/simple_parameters_whwp.py")
params = import_string_as_module("./examples/inifiles/so_sac_parameters.py")
params.tag = 'test'
rank = MPI.COMM_WORLD.rank
size = MPI.COMM_WORLD.size
##################################################################
## START OF THE SIMULATION
## Flow is the following:
## sky -> instrument -> scanning strategy ->
## MAP2TOD -> (systematics) -> TOD2MAP
##################################################################
## Initialise our input maps
sky_in = HealpixFitsMap(params.input_filename,
fwhm_in=params.fwhm_in,