def transitquick_circ_rprs(t,p):
""" Computes a transit lightcurve, normalized to unity, as a
function of time, t, usually in HJD or BJD.
Input parameters (x) are:
p[0] = tT = transit duration in days
p[1] = tG = ingress/egress duration in days
p[2] = RpRsSQ = (R_p/R_*)^2 ~ transit depth
p[3] = u1 = linear limb-darkening parameter
p[4] = u2 = quadratic limb-darkening coefficient
p[5] = F0 = uneclipsed flux
p[6] = Period = calculated in multi_transitquick
p[7] = T0.TN = Transit time
"""
RpRs = np.sqrt(p[2])
th = (2e0*np.pi/p[6])*(t-p[7])
vel = 2e0*np.sqrt(RpRs/(p[0]*p[1]))
aRs = p[6]*vel/(2e0*np.pi)
b = (np.sqrt(1e0-(p[0]/p[1])*RpRs))
inc = np.arccos(b/aRs)
z0 = aRs*np.sqrt(1e0 - (np.cos(th)*np.sin(inc))**2e0)
out = tqe.occultquad(z0,p[3],p[4],RpRs)
flux= out[:][0]*p[5]
#print RpRs, vel, aRs, b, inc, a5, a4, a0
return flux
def transitquick_circ_depth(t,p):
""" The python translation of transitquick.pro. Computes a transit lightcurve as a
function of time, t, usually in HJD or BJD.
Input parameters (x) are:
p[0] = tT = transit duration in days
p[1] = tG = ingress/egress duration in days
p[2] = D = (R_p/R_*)^2 * I(b)/I(0)= transit depth
p[3] = u1 = linear limb-darkening parameter
p[4] = u2 = quadratic limb-darkening coefficient
p[5] = F0 = uneclipsed flux
p[6] = Period = calculated in multi_transitquick
p[7] = T0.TN = Transit time
"""
a5 = -1.0e0*(p[3]+2.0e0*p[4])/(p[4]*np.sqrt(p[0]/p[1]))
a4 = -1.0e0*(1.0e0-p[3]-p[4])/(p[4]*p[0]/p[1])
a0 = p[2]*(1e0-p[3]/3e0-p[4]/6e0)/(p[4]*p[0]/p[1])
RpRs = (tqe.newtraph6(0.3,a5,a4,a0))**2e0
th = (2e0*np.pi/p[6])*(t-p[7])
vel = 2e0*np.sqrt(RpRs/(p[0]*p[1]))
aRs = p[6]*vel/(2e0*np.pi)
b = (np.sqrt(1e0-(p[0]/p[1])*RpRs))
inc = np.arccos(b/aRs)
z0 = aRs*np.sqrt(1e0 - (np.cos(th)*np.sin(inc))**2e0)
out = tqe.occultquad(z0,p[3],p[4],RpRs)
flux= out[:][0]*p[5]
#print RpRs, vel, aRs, b, inc, a5, a4, a0
return flux