code: def apply_fseries_time_shift(htilde, dt, kmin=0, copy=True):
      
      if htilde.precision != 'single':
          raise NotImplementedError("CUDA version of apply_fseries_time_shift only supports single precision")
  
      if copy:
          out = htilde.copy()
      else:
          out = htilde
  
      kmin = numpy.int32(kmin)
      kmax = numpy.int32(len(htilde))
      nb = int(numpy.ceil(kmax / nt_float))
      if nb > 1024:
          raise ValueError("More than 1024 blocks not supported yet")
  
      phi = numpy.float32(-2 * numpy.pi * dt * htilde.delta_f)
      fseries_ts_fn.prepared_call((nb, 1), (nt, 1, 1), out.data.gpudata, phi, kmin, kmax)
      if copy:
          htilde = FrequencySeries(out, delta_f=htilde.delta_f, epoch=htilde.epoch,
                                   copy=False)
      return htilde
docs: Shifts a frequency domain waveform in time. The waveform is assumed to
      be sampled at equal frequency intervals.