def METHOD_NAME( Ntot, old_array, new_array ):
    """
    Copy the `Ntot` elements of `old_array` into `new_array`, on GPU
    """
    # Loop over single particles
    ip = cuda.grid(1)
    if ip < Ntot:
        new_array[ip] = old_array[ip]