Non-perturbative renormalization of Phi-derivable approximations in theories with multiple fields

We provide a renormalization procedure for Phi-derivable approximations in theories coupling different types of fields. We illustrate our approach on a scalar phi^4 theory coupled to fermions via a Yukawa-like interaction. The non-perturbative renormalization amounts to fixing the scalar coupling via a set of nested Bethe-Salpeter equations coupling fermions to scalars.