Config.set("general", "analysis", analysis_section)
Config.set("general", "sims", sim_section)
Config.set("general", "kappa", kappa_section)
Config.set("general", "cosmology", cosmology_section)
Config.set("general", "cosmology", cosmology_section)
Config.set("general", "Nsims", str(Nsims))
Config.set("general", "Nsets", str(Nsets))
Config.set("general", "pol", str(not (args.skip_pol)))
Config.write(open(fout_dir + "config.ini", 'w'))
lens_order = Config.get(sim_section, 'lens_order')
pol = args.do_pol
if pol: raise NotImplementedError
shape_sim, wcs_sim, shape_dat, wcs_dat = aio.enmaps_from_config(
Config, sim_section, analysis_section, pol=pol)
if rank == 0:
print(("Dat shape : ", shape_dat))
print(("Sim shape : ", shape_sim))
analysis_resolution = np.min(
enmap.extent(shape_dat, wcs_dat) / shape_dat[-2:]) * 60. * 180. / np.pi
pixratio = analysis_resolution / Config.getfloat(sim_section, "pixel_arcmin")
min_ell = fmaps.minimum_ell(shape_dat, wcs_dat)
parray_dat = aio.patch_array_from_config(Config,
expf_name,
shape_dat,
wcs_dat,
dimensionless=False)
parray_sim = aio.patch_array_from_config(Config,
expf_name,
pol = args.pol
# Read config
iniFile = "input/recon.ini"
Config = SafeConfigParser()
Config.optionxform = str
Config.read(iniFile)
lens_order = Config.getint("sims", "lens_order")
pixratio = Config.getfloat("analysis", "pixel_arcmin") / Config.getfloat(
"sims", "pixel_arcmin")
pol = Config.getboolean("reconstruction", "pol")
shape_sim, wcs_sim, shape_dat, wcs_dat = aio.enmaps_from_config(Config,
"sims",
"analysis",
pol=pol)
parray_sim = aio.patch_array_from_config(Config,
exp_name,
shape_sim,
wcs_sim,
dimensionless=True)
parray_dat = aio.patch_array_from_config(Config,
expf_name,
shape_dat,
wcs_dat,
dimensionless=True)
lmax, tellmin, tellmax, pellmin, pellmax, kellmin, kellmax = aio.ellbounds_from_config(
Config, "reconstruction")
if cluster:
