Consider the partial Loewner order $\le_L$ for symmetric matrices: let $A,B$ be symmetric matrices of the same dimension, we say $A\le_L B$ if $B-A$ is positive definite. Now let $f:\mathbb{R}^n\to \mathbb{R}^n$ be a Lipschitz function, i.e., there exists $L\ge 0$ such that $ |f(x)-f(y)|\le L|x-y|$. I was wondering whether the following statement holds: there exists $C\ge 0$ such that for all $\epsilon\in (0,1]$ and $x\in \mathbb{R}^n$, $$ f(x)x^\top+xf(x)^\top \le_L \epsilon I+\frac{C}{\epsilon} xx^\top. $$
The claim holds for $n=1$. In this case, $\le_L$ is equivalent to the Euclidean order. By the Lipschitz continuity of $f$, $ f(x)\le |f(0)|+L|x|$ for all $x\in \mathbb{R}$. Then by Young's inequality, for all $\epsilon\in (0,1]$, \begin{align*} xf(x) \le |f(0)||x|+L|x|^2\le \epsilon+ \left(L+\frac{|f(0)|^2}{4\epsilon}\right)x^2 \le \epsilon+ \frac{4L+|f(0)|^2}{4\epsilon}x^2. \end{align*} It is not clear to me how to extend the argument to a multidimensional setting.