  We propose a time-adaptive, high-order compact finite difference scheme for
option pricing in a family of stochastic volatility models. We employ a
semi-discrete high-order compact finite difference method for the spatial
discretisation, and combine this with an adaptive time discretisation,
extending ideas from [LSRHF02] to fourth-order multistep methods in time.
