Just to add that I already checked the (excellent) notebooks on implementing AR/MA models given at this address and which were very helpful in building the code I shared above for the univariate AR(1) stochastic volatility model. I also checked some of the test codes available here but none used scan to implement multivariate time series models.