Quasi maximum likelihood estimation and inference in multivariate volatility models remains a challenging computational task if, for example, the dimension is high. One of the reasons is that typically numerical procedures are used to compute the score and the Hessian, and often they are...