I would like to simultaneously diagonalise the quadratic forms
$A=2x^2+3y^2+3z^2-2yz$,
and
$B=x^2+3y^2+3z^2+6xy+2yz-6zx$.
Of course there's a theorem saying this is possible and I followed the proof (diagonalise, via $P$, the positive definite $A$ to $I$ first and then find a unitary matrix to diagonalise $P^TBP$ etc.) but I ended up with really complicated matrices that made me doubt if the question intended me to do it that way. Could you please let me know if there are clever ways to solve this please. Preferably the clever way is also systematic.
Thanks!