################################################################################################################################
psr = [NX01_psr.PsrObj(t2psr[ii]) for ii in range(len(t2psr))]
[psr[ii].grab_all_vars() for ii in range(len(psr))]
psr_positions = [
np.array([psr[ii].psr_locs[0], np.pi / 2. - psr[ii].psr_locs[1]])
for ii in range(len(psr))
]
positions = np.array(psr_positions).copy()
CorrCoeff = np.array(anis.CorrBasis(
positions,
args.LMAX)) # computing all the correlation basis-functions for the array
harm_sky_vals = utils.SetupPriorSkyGrid(
args.LMAX) # computing the values of the spherical-harmonics up to order
# LMAX on a pre-specified grid
if args.anis_modefile is None:
gwfreqs_per_win = int(1. * args.nmodes / (1. * args.num_gwfreq_wins)
) # getting the number of GW frequencies per window
anis_modefreqs = np.arange(1, args.nmodes + 1)
anis_modefreqs = np.reshape(anis_modefreqs,
(args.num_gwfreq_wins, gwfreqs_per_win))
tmp_num_gwfreq_wins = args.num_gwfreq_wins
else:
tmp_modefreqs = np.loadtxt(args.anis_modefile)
tmp_num_gwfreq_wins = tmp_modefreqs.shape[0]
anis_modefreqs = []
for ii in range(tmp_num_gwfreq_wins):
