def test_b_circular(self): """Impact parameter mapping test for circular orbits. Test impact parameter mapping to ensure that b == z(t0) for circular orbits. """ Is = o.i_from_ba(self.bs, self.a) zs = np.concatenate( [of.z_circular(self.t, self.t0, self.p, self.a, i, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([ of.z_eccentric_newton(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is ]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([ of.z_eccentric_iter(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is ]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([ of.z_eccentric_ps3(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is ]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([ of.z_eccentric_ps5(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is ]) npt.assert_array_almost_equal(self.bs, zs)
def compute_lc_model(self, pv): _tc, _p, _b = pv[0], pv[1], pv[3] _ld = pv[self.ldc_slice] _a = as_from_rhop(pv[2], pv[1]) _i = m.acos(_b/_a) k = np.sqrt(pv[self.k2_slice]) _k = k.mean() kf = (k/_k)**2 z = of.z_circular(self.time, _tc, _p, _a, _i, 4) flux = self.tmodel(z, _k, pv[self.ldc_slice]) return (kf*(flux-1.)+1.).reshape((self.npt,self.npb))
def compute_lc_model(self, pv): _tc, _p, _b = pv[0], pv[1], pv[3] _ld = pv[self.ldc_slice] _a = as_from_rhop(pv[2], pv[1]) _i = m.acos(_b/_a) self.ld[0::2] = 2.*np.sqrt(_ld[0::2])*_ld[1::2] self.ld[1::2] = np.sqrt(_ld[0::2])*(1.-2.*_ld[1::2]) k = np.sqrt(pv[self.k2_slice]) _k = k.mean() kf = (k/_k)**2 z = of.z_circular(self.time, _tc, _p, _a, _i, 4) flux = self.tmodel(z, _k, self.ld) return kf*(flux-1.)+1.
def test_b_circular(self): """Impact parameter mapping test for circular orbits. Test impact parameter mapping to ensure that b == z(t0) for circular orbits. """ Is = o.i_from_ba(self.bs, self.a) zs = np.concatenate([of.z_circular(self.t, self.t0, self.p, self.a, i, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([of.z_eccentric_newton(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([of.z_eccentric_iter(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([of.z_eccentric_ps3(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs) zs = np.concatenate([of.z_eccentric_ps5(self.t, self.t0, self.p, self.a, i, 0, 0, 1) for i in Is]) npt.assert_array_almost_equal(self.bs, zs)
def test_circular(self): for (i, t0), w in product(enumerate(self.t0s), self.ws): z = of.z_circular(self.t, t0, self.p, self.a, self.i, nth=1) npt.assert_array_almost_equal(z, self.ref_zs[i])
def test_circular(self): for (i,t0), w in product(enumerate(self.t0s), self.ws): z = of.z_circular(self.t, t0, self.p, self.a, self.i, nth=1) npt.assert_array_almost_equal(z, self.ref_zs[i])