Trigonometric Quadrature Fourier Features for Scalable Gaussian Process Regression