What follows couldn't have existed without the solution by S. Dolan and somewhat paraphrase it, because it comes from the "canonical factorization" of his quadratic. As it is too long for a comment, I write it as an answer.
Let
$$f(x,y):=3(x^2-x+1)(y^2-y+1)-2(x^2y^2 -xy+1)\tag{1}$$
We have to prove that, for all $x,y$,
$$f(x,y) \geq 0. \tag{2}$$
It is consequence of the (non-evident!!) identity
$$4(y^2 - 3y + 3)f(x,y)=\underbrace{(2xy^2- 6xy - 3y^2+ 6x + 5y- 3)^2+3(y^2-3y+1)^2}_{\text{sum of squares}}\tag{3}$$
Indeed, quadratic $y^2-3y+3>0$, having a negative discriminant, is always $>0$, proving assertion (2).
I am fully aware that finding a factorization such as (3) directly is not evident. But maybe some future reader will be able to say that it can be found by a direct reasoning...
Remark : there is a "sister identity" of (3) obtained by exchanging the roles of $x$ and $y$... Could it be a source of a better understanding ?