def evaluate(self, par, mjd, t0par=None, integrate=False):
"""Evaluate 'par' at given MJD(s), with zero point t0par
Parameters
----------
par : string
key into dictionary = parameter to evaluate. Takes into
account possible par+'DOT' (for 'F': F0, F1, F2, F3)
mjd : float or Time object
MJD at which to evaluate (or Time with .tdb.mjd attribute)
t0par : string or None
key into dictionary with zero point for par
default: None -> PEPOCH for 'F', 'TASC' for all others
integrate : bool
Whether to integrate the polynomial (e.g., to get mean
anomaly out of 'FB' or pulse phase out of 'F')
"""
if par == 'F':
t0par = t0par or 'PEPOCH'
parpol = Polynomial((self['F0'], self.get('F1', 0.),
self.get('F2', 0.), self.get('F3', 0.)))
else:
t0par = t0par or 'TASC'
parpol = Polynomial((self[par], self.get(par + 'DOT', 0.)))
if integrate:
parpol = parpol.integ()
# given time can be Time object
if hasattr(mjd, 'tdb'):
mjd = mjd.tdb.mjd
return parpol((mjd - self[t0par]) * 86400.)
def evaluate(self, par, mjd, t0par=None, integrate=False):
"""Evaluate 'par' at given MJD(s), with zero point t0par
Parameters
----------
par : string
key into dictionary = parameter to evaluate. Takes into
account possible par+'DOT' (for 'F': F0, F1, F2, F3)
mjd : float or Time object
MJD at which to evaluate (or Time with .tdb.mjd attribute)
t0par : string or None
key into dictionary with zero point for par
default: None -> PEPOCH for 'F', 'TASC' for all others
integrate : bool
Whether to integrate the polynomial (e.g., to get mean
anomaly out of 'FB' or pulse phase out of 'F')
"""
if par == 'F':
t0par = t0par or 'PEPOCH'
parpol = Polynomial((self['F0'], self.get('F1', 0.),
self.get('F2',0.), self.get('F3',0.)))
else:
t0par = t0par or 'TASC'
parpol = Polynomial((self[par], self.get(par+'DOT', 0.)))
if integrate:
parpol = parpol.integ()
# given time can be Time object
if hasattr(mjd, 'tdb'):
mjd = mjd.tdb.mjd
return parpol((mjd - self[t0par]) * 86400.)
def evaluate(self, par, mjd, t0par='TASC', integrate=False):
if par == 'F':
parpol = Polynomial((self['F'], self.get('F1', 0.),
self.get('F2',0.), self.get('F3',0.)))
else:
parpol = Polynomial((self[par], self.get(par+'DOT', 0.)))
if integrate:
parpol = parpol.integ()
dt = (mjd-self[t0par])*24.*3600.
return parpol(dt)
def evaluate(self, par, mjd, t0par='TASC', integrate=False):
parpol = Polynomial((self[par], self.get(par+'DOT', 0.)))
if integrate:
parpol = parpol.integ()
dt = (mjd-self[t0par])*24.*3600.
return parpol(dt)
def evaluate(self, par, mjd, t0par="TASC", integrate=False):
parpol = Polynomial((self[par], self.get(par + "DOT", 0.0)))
if integrate:
parpol = parpol.integ()
dt = (mjd - self[t0par]) * 24.0 * 3600.0
return parpol(dt)
