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

...es.m: I have experienced that for large models with the HP filter option, this may spare a lot of time, if only the unconditional moments are needed on output;