Add support for a Williamson-type RK45 method.

This requires three registers of storage (the same as RK4) but
has a ~60% greater CFL limit, albeit at the cost of an extra
function evaluation.  However, for CFL-limited simulations a
~30% overall reduction in simulation time can be expected.

With effort this scheme can be reduced down to requiring only
two registers of storage.