We introduce a new numerical method called the complex Fourier series (CFS) method proposed by Chan (2017) to price options with an early-exercise feature-American, Bermudan and discretely monitored barrier options-under exponential Lévy asset dynamics. This new method allows us to quickly and accurately compute the values of early-exercise options and their Greeks. We also provide an error analysis to demonstrate that, in many cases, we can achieve an exponential convergence rate in the pricing method as long as we choose the correct truncated computational interval. Our numerical analysis indicates that the CFS method is computationally more comparable or favourable than the methods currently available. Finally, the superiority of the CFS method is illustrated with real financial data by considering Standard & Poor's depositary receipts (SPDR) exchange-traded fund (ETF) on the S&P 500 R index options, which are American options traded