Use the nodecomposition option which is available in th_autocovariances.m: for large models with the HP filter option, this may spare a lot of time, if only the unconditional moments are needed on output;
Merge pull request #702 from rattoma/master

Use the nodecomposition option which is available in th_autocovarianc...