rge.py 文件源码

def smeft_evolve(C_in, scale_high, scale_in, scale_out, **kwargs):
    def func(y, t0):
        return beta.beta_array(C=beta.C_array2dict(y.view(complex)),
                               HIGHSCALE=scale_high).view(float)/(16*pi**2)
    y0 = beta.C_dict2array(C_in).view(float)
    sol = odeint(func=func, y0=y0, t=(log(scale_in), log(scale_out)), **kwargs)
    return beta.C_array2dict(sol[1].view(complex))