$f$ is differentiable on the open set $\{(x,y)\in\mathbb R^2\mid x\neq y\}$ as a product of such functions. Let $x_0 \in \mathbb R$, we look at partial derivatives at $(x_0,x_0)$. We have $\lim_{h\to 0^+}\frac{f(x_0+h,x_0)-f(x_0,x_0)}h=\lim_{h\to 0^+}e^{x_0}\frac{e^h-1}h(2x_0-2)=e^{x_0}(2x_0-2)$ and $\lim_{h\to 0^-}\frac{f(x_0+h,x_0)-f(x_0,x_0)}h=\lim_{h\to 0^-}e^{x_0}\frac{1-e^h}h(2x_0-2)=-e^{x_0}(2x_0-2),$ so $f$ cannot have a partial derivative with respect to the first variable if $x_0\neq 1$. For $x_0=1$, we can see that the derivative at $(1,1)$ with respect to $x$ is $0$, and the derivative with respect to $y$ is also $0$. $f$ is differentiable at $(1,1)$ if and only if $\lim_{(x,y)\to (1,1)}\frac{f(x,y)-f(1,1)-\left(\frac{\partial f}{\partial x}(1,1)x+\frac{\partial f}{\partial y}(1,1)y\right)}{\sqrt{(x-1)^2+(y-1)^2}}=0.$ Since $\frac{\partial f}{\partial x}(1,1)=\frac{\partial f}{\partial y}(1,1)=0$, we have to see if $\lim_{(x,y)\to (1,1)}\frac{|e^x-e^y|(x+y-2)}{\sqrt{(x-1)^2+(y-1)^2}}=0$. But $\left|\frac{|e^x-e^y|(x+y-2)}{\sqrt{(x-1)^2+(y-1)^2}}\right|\leq |e^x-e^y|\frac{|x-1|+|y-1|}{\sqrt{(x-1)^2+(y-1)^2}}\leq 2|e^x-e^y|,$ so the above limit is $0$.
Conclusion: $f$ is differentiable at the points $(x,y), x\neq y$ and $(1,1)$.