astroemperor.py 文件源码

python
阅读 31 收藏 0 点赞 0 评论 0

项目:astroEMPEROR 作者: ReddTea 项目源码 文件源码
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
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号