  We construct efficient banks of templates suitable for searches of continuous
gravitational waves from isolated spinning neutron stars. We assume that the
search algorithm is based on the time-domain maximum-likelihood
$\mathcal{F}$-statistic and we consider narrow-band searches with several
spindown parameters included. Our template banks are suitable for both all-sky
and directed searches and they enable the usage of the fast Fourier transform
in the computation of the $\mathcal{F}$-statistic as well as, in the case of
all-sky searches, efficient resampling of data to barycentric time.
