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]