def derivative(x, mu, nf):
alpha_s = x[-2]
alpha_e = x[-1]
c_real = x[:-2]
c = c_real.view(np.complex)
d_alpha = betafunctions.beta_qcd_qed([alpha_s, alpha_e], mu, nf)
d_c = np.dot(adm(nf, alpha_s, alpha_e).T, c) / mu
d_c_real = d_c.view(np.float)
return np.append(d_c_real, d_alpha)
def derivative(x, mu, nf):
alpha_s = x[-2]
alpha_e = x[-1]
c_real = x[:-2]
c = c_real.view(np.complex)
d_alpha = betafunctions.beta_qcd_qed([alpha_s, alpha_e], mu, nf)
d_c = np.dot(adm(nf, alpha_s, alpha_e).T, c)/mu
d_c_real = d_c.view(np.float)
return np.append(d_c_real, d_alpha)
def _derivative_mq(x, mu, nf):
d_alphas = betafunctions.beta_qcd_qed([x[0],0], mu, nf)[0] # only alpha_s
d_m = masses.gamma_qcd(x[1], x[0], mu, nf)
return [ d_alphas, d_m ]
def _derivative_mq(x, mu, nf):
d_alphas = betafunctions.beta_qcd_qed([x[0], 0], mu, nf)[0] # only alpha_s
d_m = masses.gamma_qcd(x[1], x[0], mu, nf)
return [d_alphas, d_m]
