numpy.delete()
を使ってみました。
Z = [np.delete(x, i) for i, x in enumerate([np.arange(N)]*N)]
実行速度は neco さんが書かれたコードよりほんの少しは速そうですが、timeit モジュールのオプション(-n
と -r
)に指定する値で実行時間が結構変動します。まぁ、同程度といったところでしょうか。
$ sudo sh -c "echo 3 > /proc/sys/vm/drop_caches"
$ python -m timeit -v -n 10 -r 10 -s 'import numpy as np;N=10000' \
'Z=[np.delete(x, i) for i, x in enumerate([np.arange(N)]*N)]'
raw times: 1.42 1.07 1.03 1.03 1.03 1.04 1.03 1.03 1.02 1.02
10 loops, best of 10: 102 msec per loop
$ sudo sh -c "echo 3 > /proc/sys/vm/drop_caches"
$ python -m timeit -v -n 10 -r 10 -s 'import numpy as np;N=10000' \
'Z=[np.concatenate((np.arange(0, i), np.arange(i+1, N))) for i in range(N)]'
raw times: 3.54 1.5 1.31 1.12 1.33 1.17 1.18 1.13 1.28 1.31
10 loops, best of 10: 112 msec per loop
この様な方法もありますよ、という事で参考までにどうぞ。
追記
他の回答者の方々のスクリプトも計測してみました。
$ python -m timeit -v -n 10 -r 10 -s 'import numpy as np;N=10000' \
'Z=np.tile(np.arange(N), (N, 1))[np.eye(N, N) == 0].reshape(N, N - 1)'
raw times: 7.59 7.58 7.59 7.59 7.59 7.59 7.58 7.58 7.58 7.59
10 loops, best of 10: 758 msec per loop
$ python -m timeit -v -n 10 -r 10 -s 'import numpy as np;N=10000' \
'Z=np.tile(np.arange(N), N)[np.tile([False] + [True] * N, N)[:N*N]].reshape(N, N-1)'
raw times: 6.35 6.35 6.35 6.35 6.34 6.35 6.34 6.34 6.34 6.34
10 loops, best of 10: 634 msec per loop
$ python -m timeit -v -n 10 -r 10 \
-s "$(cat <<'EOS'
import numpy as np
N=10000
def by_np_slice():
nums = np.arange(0, N, dtype=np.int)
result = np.empty((N, N-1), dtype=np.int)
for i, a in enumerate(result):
a[:i] = nums[:i]
a[i:] = nums[i+1:]
return result
EOS
)" 'Z=by_np_slice()'
raw times: 2.66 2.66 2.66 2.65 2.65 2.66 2.66 2.65 2.66 2.65
10 loops, best of 10: 265 msec per loop
$ python -m timeit -v -n 10 -r 10 \
-s "$(cat <<'EOS'
import numpy as np
N=10000
def by_np_arange(N):
result = np.empty((N, N-1),dtype=np.int)
for i in range(N):
a = np.take(result,i)
a = np.concatenate((np.arange(0,i), np.arange(i+1,N)))
return result
EOS
)" 'Z=by_np_arange(N)'
raw times: 0.978 0.976 0.977 0.975 0.979 0.976 0.974 0.976 0.974 0.974
10 loops, best of 10: 97.4 msec per loop
Z = [range(0,i) + range(i+1,N) for i in range(N)]
とした場合、N=10000
で 3.79秒 -> 1.14秒になりました(timeit
で計測)。 – metropolis 18年8月18日 0:18