Variance reduction in computations of neoclassical transport in stellarators using a δf method