We propose a fast and accurate numerical method for pricing European swaptions in multi-factor Gaussian term structure models. Our method can be used to accelerate the calibration of such models to the volatility surface. The pricing of an interest rate option in such a model involves evaluating...