def residuals(self, ts, rvs, p):
"""Return the residuals for the rv model with parameters ``p``
and the observations of radial velocitys ``rv`` at times
``t``"""
if p.npl == 0:
return rvs
else:
rvmodel = np.sum(rv.rv_model(ts,p), axis=0)
return rvs - rvmodel
def residuals(self, ts, rvs, p):
"""Return the residuals for the rv model with parameters ``p``
and the observations of radial velocitys ``rv`` at times
``t``"""
if p.npl == 0:
return rvs
else:
rvmodel = np.sum(rv.rv_model(ts, p), axis=0)
return rvs - rvmodel
def generate(self):
"""Returns (ts, rvs), where ts is a list of arrays of
observation times (one array for each observatory), and rvs is
a corresponding list of total radial velocity measurements."""
ts=self.generate_tobs()
noise=self.generate_noise(ts)
rvs=[]
for t,n in zip(ts, noise):
rvs.append(n + np.sum(rv.rv_model(t, self.params), axis=0))
return ts,rvs
def generate(self):
"""Returns (ts, rvs), where ts is a list of arrays of
observation times (one array for each observatory), and rvs is
a corresponding list of total radial velocity measurements."""
ts = self.generate_tobs()
noise = self.generate_noise(ts)
rvs = []
for t, n in zip(ts, noise):
rvs.append(n + np.sum(rv.rv_model(t, self.params), axis=0))
return ts, rvs
def rvs(self, p):
p = self.to_params(p)
return rv.rv_model(self.ts, p['K'], p['e'], p['omega'], p['chi'],
p['P'])
def rvs(self, p):
p = self.to_params(p)
return rv.rv_model(self.ts, p['K'], p['e'], p['omega'], p['chi'], p['P'])
def rvs(self, p):
p = self.to_params(p)
return rv.rv_model(self.ts, p["K"], p["e"], p["omega"], p["chi"],
p["P"])
