Let $$ f(x) = \log(1+e^{2x+1}) - 2\log(1 + e^{2x}) + \log(1 + e^{2x-1}). $$ According to Wolfram Alpha, $$ \int_{-\infty}^\infty f(x)\,dx = \frac 12.\tag{$*$} $$ $f(x)$ is a "bump function" built out of the softplus function $y=\log(1 + e^x)$, that I want to use as a kernel for nonparametric regression. I want to normalize $f(x)$ to have integral $1$, and the definite integral (*) gives the normalizing factor.
But I can't figure out how to do the integral by hand. Substitution gives (and Wolfram Alpha confirms) $$ \int f(x)\,dx = \operatorname{Li}_2(-e^{2x}) - \frac12\operatorname{Li}_2(-e^{2x-1}) - \frac12\operatorname{Li}_2(-e^{2x+1}) + C, $$ where $\operatorname{Li}_2(x)$ is the dilogarithm function. I haven't had any success in applying dilogarithm identities to the above, though.
Can anyone give a derivation of $(*)$?