Let us consider two numbers of the form $a^2 - ab + b^2$ and $c^2 - cd + d^2$ which are not both divisible by $3$ and such that $(a, b) = 1$ and $(c,d) = 1$. Running some computations it seems that the product $$(a^2 -ab + b^2)(c^2 - cd + d^2) $$ is still of the form $h^2 - hk + k^2$ for some suitable coprime integers $h,k$. Is this true?
I tried to prove it by writing down explicitly the product and looking for patterns, but I had no luck. Any help would be appreciated!