def setUpClass(self):
self.parfileB1855 = 'B1855+09_polycos.par'
self.timB1855 = 'B1855_polyco.tim'
self.toasB1855 = toa.get_TOAs(self.timB1855,
ephem="DE405",
planets=False)
self.modelB1855 = mb.get_model(self.parfileB1855)
# Read tempo style polycos.
self.plc = Polycos()
self.plc.read_polyco_file('B1855_polyco.dat', 'tempo')
def setUpClass(self):
self.parfileB1855 = 'B1855+09_polycos.par'
self.timB1855 = 'B1855_polyco.tim'
self.toasB1855 = toa.get_TOAs(self.timB1855, ephem="DE405", planets=False)
self.modelB1855 = mb.get_model(self.parfileB1855)
# Read tempo style polycos.
self.plc = Polycos()
self.plc.read_polyco_file('B1855_polyco.dat', 'tempo')
class TestD_phase_d_toa(unittest.TestCase):
@classmethod
def setUpClass(self):
self.parfileB1855 = 'B1855+09_polycos.par'
self.timB1855 = 'B1855_polyco.tim'
self.toasB1855 = toa.get_TOAs(self.timB1855, ephem="DE405", planets=False)
self.modelB1855 = mb.get_model(self.parfileB1855)
# Read tempo style polycos.
self.plc = Polycos()
self.plc.read_polyco_file('B1855_polyco.dat', 'tempo')
def TestD_phase_d_toa(self):
pint_d_phase_d_toa = self.modelB1855.d_phase_d_toa(self.toasB1855)
mjd = np.array([np.longdouble(t.jd1 - ut.DJM0)+np.longdouble(t.jd2) for t in self.toasB1855.table['mjd']])
tempo_d_phase_d_toa = self.plc.eval_spin_freq(mjd)
diff = pint_d_phase_d_toa - tempo_d_phase_d_toa
relative_diff = diff/tempo_d_phase_d_toa
assert np.all(relative_diff < 1e-8), 'd_phae_d_toa test filed.'
class TestD_phase_d_toa(unittest.TestCase):
@classmethod
def setUpClass(self):
self.parfileB1855 = 'B1855+09_polycos.par'
self.timB1855 = 'B1855_polyco.tim'
self.toasB1855 = toa.get_TOAs(self.timB1855,
ephem="DE405",
planets=False)
self.modelB1855 = mb.get_model(self.parfileB1855)
# Read tempo style polycos.
self.plc = Polycos()
self.plc.read_polyco_file('B1855_polyco.dat', 'tempo')
def TestD_phase_d_toa(self):
pint_d_phase_d_toa = self.modelB1855.d_phase_d_toa(self.toasB1855)
mjd = np.array([
np.longdouble(t.jd1 - ut.DJM0) + np.longdouble(t.jd2)
for t in self.toasB1855.table['mjd']
])
tempo_d_phase_d_toa = self.plc.eval_spin_freq(mjd)
diff = pint_d_phase_d_toa - tempo_d_phase_d_toa
relative_diff = diff / tempo_d_phase_d_toa
assert np.all(relative_diff < 1e-8), 'd_phae_d_toa test filed.'