numpy nditerの引数 op_flags["readwrite"]とは？

Question

def numerical_gradient(f, x):
    h = 1e-4 # 0.0001
    grad = np.zeros_like(x)

it = np.nditer(x, flags=['multi_index'], op_flags=['readwrite'])

    while not it.finished:
        idx = it.multi_index
        tmp_val = x[idx]
        x[idx] = float(tmp_val) + h
        fxh1 = f(x) # f(x+h)
        
        x[idx] = tmp_val - h 
        fxh2 = f(x) # f(x-h)
        grad[idx] = (fxh1 - fxh2) / (2*h)
        
        x[idx] = tmp_val # 値を元に戻す
        it.iternext()   
        
    return grad

ゼロから作るディープラーニングの
commonの中にある gradient.pyの関数numerical_gradienのコードです
https://github.com/oreilly-japan/deep-learning-from-scratch

このコードの大体の動作は、
http://www.aipacommander.com/entry/2017/05/14/172220
此方のサイトのおかげで理解することが出来たのですが、
以下の部分のop_flags=['readwrite']が何をしているのががわかりません

it = np.nditer(x, flags=['multi_index'], op_flags=['readwrite'])

これは何をしているのか、詳しいかたいればお願いします><

Answer 1

以下が参考になります。
Modifying Array Values

forループなどでイテレータで要素に参照する場合の要素への読み書きモードを指定するようです。

ただし提示されたコードではx[idx]のように元配列xをイテレータで取得した添え字idxで参照していますが、この場合はイテレータがreadonlyでも要素の書き換えはできるようです。

import numpy as np

# イテレータで要素を書き換え readwriteが必要
a = np.arange(6).reshape(2,3)
it = np.nditer( a, flags=['multi_index'], op_flags=['readwrite'])
for x in it:
    # To actually modify the element of the array,
    # x should be indexed with the ellipsis.
    x[...] = x * 2
print(a)

# 要素を添え字で参照して要素を書き換え イテレータはreadonlyでもよいようだ
a = np.arange(6).reshape(2,3)
it = np.nditer( a, flags=['multi_index'], op_flags=['readonly'])
while not it.finished:
    a[it.multi_index] = a[it.multi_index] * 2 # # 要素を添え字で参照
    it.iternext()
print(a)

# イテレータで要素を書き換え readonlyでは不可
a = np.arange(6).reshape(2,3)
it = np.nditer( a, flags=['multi_index'], op_flags=['readonly'])
for x in it:
    x[...] = x * 2 # assignment destination is read-only
print(a)

[[ 0  2  4]
 [ 6  8 10]]
[[ 0  2  4]
 [ 6  8 10]]
（assignment destination is read-only例外が発生）

Comments

[退会済みユーザー]

いつもありがとうございます！
x[...]という書き方が見慣れないのでなんだろうと思ったら、イテレーションする変数の要素までも書き換えるという事なんですね
x = x * 2 はイテレーションする変数の要素までは書き換えないのでそこが違うと

そしてそれをする際に、readwriteが必要になってくると

もう一つ聞きたいのですが、readonlyという事は、引数を省略しても問題はないという認識で合っているでしょうか？
試しにwhile not のイテレーションのコードで op_flagsを省略したら問題なく動作したので

[can110]

詳細な動作仕様は把握できていないのですが、結果からするとwhile not の例ではイテレータのモードに関係なく要素の値を書き換えられますね。