# make sure pulsar names are in correct order
# TODO: is this a very round about way to do this?
index = []
for ct,p in enumerate(pp):
if p.name == psr[ct].name:
index.append(ct)
else:
for ii in range(npsr):
if pp[ii].name == psr[ct].name:
index.append(ii)
pp = [pp[ii] for ii in index]
M = [PALutils.createQSDdesignmatrix(p.toas) for p in psr]
RQ = [PALutils.createRmatrix(M[ct], p.err) for ct, p in enumerate(psr)]
# construct noise matrix for new noise realizations
print 'Constructing noise cholesky decompositions'
L = []
for ct, p in enumerate(psr):
Amp = p.Amp
gam = p.gam
efac = p.efac
equad = p.equad
cequad = p.cequad
avetoas, U = PALutils.exploderMatrix(p.toas)
# make sure pulsar names are in correct order
# TODO: is this a very round about way to do this?
index = []
for ct, p in enumerate(pp):
if p.name == psr[ct].name:
index.append(ct)
else:
for ii in range(npsr):
if pp[ii].name == psr[ct].name:
index.append(ii)
pp = [pp[ii] for ii in index]
M = [PALutils.createQSDdesignmatrix(p.toas) for p in psr]
RQ = [PALutils.createRmatrix(M[ct], p.err) for ct, p in enumerate(psr)]
# construct noise matrix for new noise realizations
print 'Constructing noise cholesky decompositions'
L = []
for ct, p in enumerate(psr):
Amp = p.Amp
gam = p.gam
efac = p.efac
equad = p.equad
cequad = p.cequad
avetoas, U = PALutils.exploderMatrix(p.toas)
