  Renormalization group flows of the $SU(N_f)\times SU(N_f)$ symmetric
Ginzburg-Landau potential are calculated for a general number of flavors,
$N_f$. Our approach does not rely on the $\epsilon$ expansion, but uses the
functional renormalization group, formulated directly in $d=3$ spatial
dimensions, with the inclusion of all possible (perturbatively) relevant and
marginal operators, whose number is considerably larger than those in $d=4$. We
find new, potentially infrared stable fixed points spanned throughout the
entire $N_f$ range. By conjecturing that the thermal chiral transition is
governed by these ``flavor continuous" fixed points, stability analyses show
that for $N_f\geq 5$ the chiral transition is of second-order, while for
$N_f=2,3,4$, it is of first-order. We argue that the $U_{\rm A}(1)$ anomaly
controls the strength of the first-order chiral transition for $N_f=2,3,4$, and
makes it almost indistinguishable from a second-order one, if it is
sufficiently weak at the critical point. This could open up a new strategy to
investigate the strength of the $U_{\rm A}(1)$ symmetry breaking around the
critical temperature.
