def henshin(self, thetas, kplanets):
for i in range(kplanets):
Ask = thetas[:, i*5]
Pk = thetas[:, i*5 + 1]
Ack = thetas[:, i*5 + 2]
Sk = thetas[:, i*5 + 3]
Ck = thetas[:, i*5 + 4]
ecck = Sk ** 2 + Ck ** 2
Ak = Ask ** 2 + Ack ** 2
wk = sp.arccos(Ck / (ecck ** 0.5))
Phasek = sp.arccos(Ack / (Ak ** 0.5))
for j in range(len(Sk)):
if Sk[j] < 0:
wk[j] = 2 * sp.pi - sp.arccos(Ck[j] / (ecck[j] ** 0.5))
if Ask[j] < 0:
Phasek[j] = 2 * sp.pi - sp.arccos(Ack[j] / (Ak[j] ** 0.5))
thetas[:, i*5] = Ak
thetas[:, i*5 + 1] = sp.exp(Pk)
thetas[:, i*5 + 2] = Phasek
thetas[:, i*5 + 3] = wk
thetas[:, i*5 + 4] = ecck
return thetas
########################################
# los IC sirven para los mod_lims? NO!
# mod_lims sólo acotan el espacio donde buscar, excepto para periodo
评论列表
文章目录