Efficient high-order spectral element discretizations for building block operators of CFD

•Efficient implementation of tensor-product operators for high-order spectralelement methods.•Exploitation of parametrization, tailored unroll-jam and blocking of loops.•50% of the peak performance attained over a wide range of polynomial degrees.•Benefit of new operator implementation demonstrated on combustion problem with 1: 72 · 109 mesh points.