  We present a method for calculating the kinetic energy of localised functions
represented on a regular real space grid. This method uses fast Fourier
transforms applied to restricted regions commensurate with the simulation cell
and is applicable to grids of any symmetry. In the limit of large systems it
scales linearly with system size. Comparison with the finite difference
approach shows that our method offers significant improvements in accuracy
without loss of efficiency.
