Let $\hat{\varrho}, \tilde{\varrho}$ be probability density functions on $\mathbb R^d$ where $\tilde{\varrho} \in L^{\infty} (\mathbb R^d)$. For $\varepsilon \in [0, 1]$, we define $\varrho_{\varepsilon}:=(1-\varepsilon) \hat{\varrho}+\varepsilon \tilde{\varrho}$. Set $M=\|\tilde{\varrho}\|_{L^{\infty}}$ and look at $\varepsilon \mapsto \int \varrho_{\varepsilon} \log \varrho_{\varepsilon}$.
It is mentioned in the proof of Proposition 8.7 in Santambrogio's Optimal transport for applied mathematicians that
The integrand can be differentiated in $\varepsilon$ pointwisely, thus getting $\left(1+\log \varrho_{\varepsilon}\right)(\tilde{\varrho}-\hat{\varrho})$. For $\varepsilon<1 / 2$, we can check that these functions are dominated by $(\hat{\varrho}+M)(|\log \hat{\varrho}|+\log M+1)$.
Could you explain how we get the inequality $$ |\left(1+\log \varrho_{\varepsilon}\right)(\tilde{\varrho}-\hat{\varrho})| \le (\hat{\varrho}+M)(|\log \hat{\varrho}|+\log M+1), \quad \forall \varepsilon \in (0, \frac{1}{2}) $$ ?
Thank you so much for your elaboration!