def calc_kernel_old(Dl, Ds, dgr, dcr, rhog, fr):
gamma = pifrth
sigma = bva
alphal, betal = get_alpha_beta(Dl, dgr, dcr, rhog, fr)
alphas, betas = get_alpha_beta(Ds, dgr, dcr, rhog, fr)
vtl = tv.calc_vt(Dl, alphal, betal, gamma, sigma)
vts = tv.calc_vt(Ds, alphas, betas, gamma, sigma)
print Dl, Ds, vtl, vts, np.abs(vtl-vts)
return E_coll*pifrth*(Dl+Ds)**2*np.abs(vtl-vts)
def calc_vt_vec(D, dgr, dcr, rhog, fr):
alpha, beta = get_alpha_beta_vec(D, dgr, dcr, rhog, fr)
# gamma, sigma = get_gamma_sigma_vec(D, dgr, dcr, rhog, fr)
gamma = pifrth
sigma = bva
vt = tv.calc_vt(D, alpha, beta, gamma, sigma)
return vt