def calc_anoms(params, MJD): """ Calculate the mean, eccentric, and true anomalies at the barycentric epoch MJD (written by Scott Ransom) Input ----- params - a list of parameter values MJD - the barycentric epoch MJD Output ------ mean_anom - mean anomaly in radians ecc_anom - eccentric enomaly in radians true_anom - the true anomaly in radians """ MJD = myasarray(MJD) difft = (MJD - params[3]) * SECPERDAY sec_since_peri = fmod(difft, params[1]) sec_since_peri[sec_since_peri < 0.0] += params[1] mean_anom = sec_since_peri / params[1] * TWOPI ecc_anom = eccentric_anomaly( [params[0], params[1], params[2], params[3], params[4], params[5]], mean_anom) true_anom = psr_utils.true_anomaly(ecc_anom, params[4]) return (mean_anom, ecc_anom, true_anom)
def calc_anoms(params, MJD): """ Calculate the mean, eccentric, and true anomalies at the barycentric epoch MJD (written by Scott Ransom) Input ----- params - a list of parameter values MJD - the barycentric epoch MJD Output ------ mean_anom - mean anomaly in radians ecc_anom - eccentric enomaly in radians true_anom - the true anomaly in radians """ MJD = myasarray(MJD) difft = (MJD - params[3])*SECPERDAY sec_since_peri = fmod(difft, params[1]) sec_since_peri[sec_since_peri < 0.0] += params[1] mean_anom = sec_since_peri/params[1]*TWOPI ecc_anom = eccentric_anomaly([params[0],params[1], params[2],params[3],params[4], params[5]], mean_anom) true_anom = psr_utils.true_anomaly(ecc_anom, params[4]) return (mean_anom, ecc_anom, true_anom)
def calc_anoms(self, MJD): """ calc_anoms(MJD): Return a tuple of the mean, eccentric, and true anomalies (all in radians) at the barycentric epoch MJD(s). """ MJD = myasarray(MJD) difft = (MJD - self.T0) * SECPERDAY sec_since_peri = Num.fmod(difft, self.PBsec) sec_since_peri[sec_since_peri < 0.0] += self.PBsec mean_anom = sec_since_peri / self.PBsec * TWOPI ecc_anom = self.eccentric_anomaly(mean_anom) true_anom = psr_utils.true_anomaly(ecc_anom, self.par.E) return (mean_anom, ecc_anom, true_anom)
def calc_anoms(self, MJD): """ calc_anoms(MJD): Return a tuple of the mean, eccentric, and true anomalies (all in radians) at the barycentric epoch MJD(s). """ MJD = myasarray(MJD) difft = (MJD - self.T0)*SECPERDAY sec_since_peri = Num.fmod(difft, self.PBsec) sec_since_peri[sec_since_peri < 0.0] += self.PBsec mean_anom = sec_since_peri/self.PBsec*TWOPI ecc_anom = self.eccentric_anomaly(mean_anom) true_anom = psr_utils.true_anomaly(ecc_anom, self.par.E) return (mean_anom, ecc_anom, true_anom)