def RV_timeseries(self,ts,recalc=False):
if type(ts) != Quantity:<EOL><INDENT>ts *= u.day<EOL><DEDENT>if not recalc and hasattr(self,'<STR_LIT>'):<EOL><INDENT>if (ts == self.ts).all():<EOL><INDENT>return self.RV_measurements<EOL><DEDENT>else:<EOL><INDENT>pass<EOL><DEDENT><DEDENT>RVs = Quantity(np.zeros((len(ts),self.N)),unit='<STR_LIT>')<EOL>for i,t in enumerate(ts):<EOL><INDENT>RVs[i,:] = self.dRV(t,com=True)<EOL><DEDENT>self.RV_measurements = RVs<EOL>self.ts = ts<EOL>return RVs<EOL>
Radial Velocity time series for star 1 at given times ts.