Fast nonasymptotic testing and support recovery for large sparse Toeplitz covariance matrices