def myloglike(cube, ndim, nparams): efac = np.ones(npsr) equad = np.zeros(npsr) A = 10 ** cube[0] gam = cube[1] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, Ared, gred, efac, equad, ORF) # print efac, rho, loglike return loglike
def myloglike(cube, ndim, nparams): efac = np.ones(npsr) equad = np.zeros(npsr) A = 10**cube[0] gam = cube[1] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \ Ared, gred, efac, equad, ORF) #print efac, rho, loglike return loglike
def loglike(x): A = x[0] * 1e-14 gam = x[1] efac = x[2:(2 + npsr)] equad = 10**x[(2 + npsr):(2 + 2 * npsr)] Ared = 10**x[(2 + 2 * npsr):(2 + 3 * npsr)] gred = x[(2 + 3 * npsr):(2 + 4 * npsr)] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \ Ared, gred, efac, equad, ORF) return loglike
def loglike(x): A = x[0]*1e-14 gam = x[1] efac = x[2:(2+npsr)] equad = 10**x[(2+npsr):(2+2*npsr)] Ared = 10**x[(2+2*npsr):(2+3*npsr)] gred = x[(2+3*npsr):(2+4*npsr)] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \ Ared, gred, efac, equad, ORF) return loglike
def myloglike(cube, ndim, nparams): efac = np.zeros(npsr) equad = np.zeros(npsr) Ared = np.zeros(npsr) gred = np.zeros(npsr) A = cube[0] gam = cube[1] for ii in range(npsr): efac[ii] = cube[ii + 2] equad[ii] = 10 ** cube[ii + 2 + npsr] Ared[ii] = 10 ** cube[ii + 2 + 2 * npsr] gred[ii] = cube[ii + 2 + 3 * npsr] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, Ared, gred, efac, equad, ORF) # print efac, rho, loglike return loglike
def myloglike(cube, ndim, nparams): efac = np.zeros(npsr) equad = np.zeros(npsr) Ared = np.zeros(npsr) gred = np.zeros(npsr) A = cube[0] gam = cube[1] for ii in range(npsr): efac[ii] = cube[ii + 2] equad[ii] = 10**cube[ii + 2 + npsr] Ared[ii] = 10**cube[ii + 2 + 2 * npsr] gred[ii] = cube[ii + 2 + 3 * npsr] loglike = PALLikelihoods.modelIndependentFullPTAPL(psr, F, Diag, f, A, gam, \ Ared, gred, efac, equad, ORF) #print efac, rho, loglike return loglike