A new method is proposed for efficient computation of auto- and cross-spectral densities in the stochastic modelling of ocean waves and wave loads. The commonly used cos-2s type directional distribution is given special attention. The new method is compared with traditional methods that use numerical integration. The new method avoids spurious hat errors and improves computational efficiency.