I know by brutal calculation this identity holds always:
$$(u × v) \cdot (x × y) = \begin{vmatrix} u \cdot x & v \cdot x \\ u \cdot y & v \cdot y \\ \end{vmatrix}$$
for arbitrary vectors $u$, $v$, $x$, $y$.
I'd like to know where it come from naturally.
Thank you.
I suppose you are talking about dot products and cross products in $\mathbb R^3$. The identity can be seen as a special case of Cauchy-Binet formula. Alternatively, note that $$ \pmatrix{x^T\\ y^T\\ (u\times v)^T}\pmatrix{u&v&x\times y} =\pmatrix{u\cdot x&v\cdot x&0\\ u\cdot y&v\cdot y&0\\ 0&0&(u\times v)\cdot(x\times y)}.\tag{1} $$ Therefore the RHS has determinant $$ \det\pmatrix{u\cdot x&v\cdot x\\ u\cdot y&v\cdot y}[(u\times v)\cdot(x\times y)].\tag{2} $$ On the other hand, since $\det\pmatrix{p&q&r}\equiv(p\times q)\cdot r$ (this actually is the definition of cross product), the LHS of $(1)$ has determinant $$ [(u\times v)\cdot(x\times y)]^2.\tag{3} $$ The result now follows by equating $(2)$ and $(3)$ with the use of a continuity argument.