def toKepler(u, which = 'Pueyo', mass = 1, referenceTime = None):
"""
"""
if which == 'Pueyo':
res = np.zeros(6)
res[1] = u[1]
res[5] = u[5]
res[0] = semimajoraxis(math.exp(u[0]), starMass = mass)
res[2] = math.degrees(math.acos(u[2]))
res[3] = np.mod((u[3]-u[4])*0.5,360)
res[4] = np.mod((u[3]+u[4])*0.5,360)
return res
elif which == 'alternative':
res = np.zeros(6)
res[1] = u[1]
res[5] = u[5]
res[0] = semimajoraxis(math.exp(u[0]), starMass = mass)
res[2] = math.degrees(math.acos(u[2]))
res[3] = u[3]
res[4] = u[4]
return res
elif which == 'Chauvin':
stat = StatisticsMCMC()
res = stat.xFROMu(u,referenceTime,mass)
return res
return None
def toSynthetic(theta,which = 'Pueyo', mass = 1, referenceTime = None):
"""
"""
if which == 'Pueyo':
res = np.zeros(6)
temp = theta
res[0] = math.log(period(temp[0],starMass = mass))
res[2] = math.cos(math.radians(temp[2]))
res[3] = np.mod(temp[4]+temp[3],360)
res[4] = np.mod(temp[4]-temp[3],360)
res[1] = temp[1]
res[5] = temp[5]
elif which == 'alternative':
res = np.zeros(6)
temp = theta
res[0] = math.log(period(temp[0],starMass = mass))
res[2] = math.cos(math.radians(temp[2]))
res[3] = temp[3]
res[4] = temp[4]
res[1] = temp[1]
res[5] = temp[5]
elif which == 'Chauvin':
stat = StatisticsMCMC()
res = stat.uFROMx(theta,referenceTime,mass)
return res
