When a quadratic form is not chi-square

Suppose $\uline{X}_{3\times1}\sim\text{Normal}_3(\uline{\mu},\Sigma)$ and

$$Y=X_1^{2}-X_2^{2}+X_3^{2}+0.5X_1X_3 .$$

Then the probability distribution of $Y$ is: