def eq_pnn_corr_ff_pnd_corr_ff(k, w, pf, part=p.P_DFLT):
eq0_tmp = p.part(eq0, part, k, w)
pi_pnn_corr_ff_tmp = p.part(pi_pnn_corr_ff, part, k, w, pf)
pi_pnd_corr_ff_tmp = p.part(pi_pnd_corr_ff, part, k, w, pf)
return eq0_tmp - pi_pnn_corr_ff_tmp - pi_pnd_corr_ff_tmp
def eq_pnn_corr(k, w, pf, part=p.P_DFLT):
pi_pnn_corr_tmp = p.part(pi_pnn_corr, part, k, w, pf)
eq0_tmp = p.part(eq0, part, k, w)
# return eq0_tmp - pi_pnn_corr_tmp
return eq0(k, w) - pi_pnn_corr(k, w, pf)
def eq_pnn_pnd(k, w, pf, part=p.P_DFLT):
'''Equation with all interactions.'''
pi_pnn_tmp = p.part(pi_pnn, part, k, w, pf)
pi_pnd_tmp = p.part(pi_pnd, part, k, w, pf)
return eq0(k, w) - pi_pnn_tmp - pi_pnd_tmp
def eq_pnd(k, w, pf, pi_pnd_part=p.P_DFLT):
'''Equation for pND interaction.'''
pi_pnd_tmp = p.part(pi_pnd, pi_pnd_part, k, w, pf)
return eq0(k, w) - pi_pnd(k, w, pf)
def eq_pnn(k, w, pf, pi_pnn_part=p.P_DFLT):
'''Equation for pNN interaction.'''
pi_pnn_tmp = p.part(pi_pnn, pi_pnn_part, k, w, pf)
eq0_tmp = p.part(eq0, pi_pnn_part, k, w)
return eq0_tmp - pi_pnn_tmp
def s(k, w, pf, pi_pnd_dw_part=p.P_DFLT):
'''One term of S(k).'''
pi_pnd_dw_tmp = p.part(pi_pnd_dw, pi_pnd_dw_part, k, w, pf)
return 2 * w / (2 * w - pi_pnd_dw_tmp)