maxntimes removed; it's not necessary as that's not the limiting factor
(covmat is symmetric with dims = (nfreq, nfreq))