for文を使ったインデックスについて

Question

print(b）では（x[0]、インデックスの０番目ということで）3.0001、
print(c)では（x[1]、インデックスの１番目ということで）4.0001の出力を予想していたのですが、
以下のとおりになります。どう解釈したらよいのでしょう。よろしくお願いいたします。

import numpy as np
x=np.array([3.0,4.0])
h=1e-4

for i in range(x.size):
    a=x[i]
    x[i]=a+h
    a=x[i]
    print(a)
    #3.0001   4.0001
    b=x[0]
    print(b)
    #3.0001   3.0001
    c=x[1]
    print(c)
    # 4.0   4.0001

Comments

[meg_]

> どう解釈したらよいのでしょう。
出力が2回ずつあるため、何回目のどこのことを質問されているのか分かりません。

[okahijiki]

分かりにくい文章にて失礼しました。どうぞ、ご容赦くださいませ。

Answer 1

aは配列ではなくスカラーの浮動小数だからです。

以下を見て考えてみてください。

>>> for i in range(x.size):
...     a=x[i]
...     print(type(a))
...     x[i]=a+h
...     a=x[i]
...     print(a)
...     #3.0001   4.0001
...     b=x[0]
...     print(b)
...     #3.0001   3.0001
...     c=x[1]
...     print(c)
...     # 4.0   4.0001
...
<class 'numpy.float64'>
3.0001
3.0001
4.0
<class 'numpy.float64'>
4.0001
3.0001
4.0001

期待している結果を得るには以下のような方法があります。

>>> import numpy as np
>>> x=np.array([3.0,4.0])
>>> h=1e-4
>>>
>>> a = x + h
>>> print(a)
[3.0001 4.0001]
>>> #3.0001   4.0001
>>> b =np.full_like(x,x[0]+h)
>>> print(b)
[3.0001 3.0001]
>>> #3.0001   3.0001
>>> c = np.stack([x, x + h])
>>> print(c[:,1])
[4.     4.0001]
>>> # 4.0   4.0001

Comments

[okahijiki]

ご回答ありがとうございます。
たしかに、type(a)としても、
listにはなりませんでしたー。

[okahijiki]

こちらから質問しておきながら、申し訳ありません。
ほんじつ、ちょっと用がありまして。明日、じっくり考えてみたいと存じます。
ありがとうございます！

[okahijiki]

np.full_likeにて、xと同じ配列の要素数とし、
np.stackにより、xとx+hを結合し、そのうえで[:,1]により２列目を抽出するのですね。
勉強になりました。ありがとうございました〜