'int' object is not callableの対処

Question

前提・実現したいこと

PythonでPython最急降下法を実装中、以下のエラーメッセージが発生しました。
変数名が問題なのでしょうか？

発生している問題・エラーメッセージ

TypeError                                 Traceback (most recent call last)
<ipython-input-88-1c0aba996c48> in <module>()
     31             pos_history.append(pos)
     32     return pos_history
---> 33 (H_1,H_2) = grad_descent((1.2, 1.2), 0.1, 0.1)
     34 plt.plot(H_1,H_2)
     35 plt.show()

<ipython-input-88-1c0aba996c48> in grad_descent(init_pos, a_0, p)
     25             break
     26         else:
---> 27             while (f(pos[0] + a * d[0],pos[1] + a * d[1]) > f(pos[1],pos[2]) - c * a * abs(np.linalg.norm(d))):#armijoの条件
     28                     a = a * p #↑をみたすように小さくしていく
     29             pos_new = [pos[0] + a * d[0],pos[1] + a * d[1]] #posの更新

<ipython-input-88-1c0aba996c48> in f(x_1, x_2)
      3 #最小化したい関数
      4 def f(x_1,x_2):
----> 5     return 100(x_2 - (x_1)**2)**2 + (1 - x_1)**2
      6 #再急降下方向
      7 def dsc_f(pos):

TypeError: 'int' object is not callable

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
#最小化したい関数
def f(x_1,x_2):
    return 100(x_2 - (x_1)**2)**2 + (1 - x_1)**2
#再急降下方向
def dsc_f(pos):
    x_1 = pos[0]
    x_2 = pos[1]
    return   (-(400 * x_1**3 - 400 * x_1 * x_2 + 2 * x_1 - 2),-((-200) * x_1 ** 2 + 200 * x_2))

#a_0:learning rateの初期値
#p:aの動かし幅
def grad_descent(init_pos ,a_0,p):
    eps = 1e-10
    iteration_max = 500 #爆発しないように
    a = a_0
    init_pos = np.array(init_pos)
    pos = init_pos
    pos_history = [init_pos]
    
    for i in range(iteration_max):
        d = np.array(dsc_f(pos))
        if abs(np.linalg.norm(d)) < eps:
            break
        else:
            while (f(pos[0] + a * d[0],pos[1] + a * d[1]) > f(pos[1],pos[2]) - c * a * abs(np.linalg.norm(d))):#armijoの条件
                    a = a * p #↑をみたすように小さくしていく
            pos_new = [pos[0] + a * d[0],pos[1] + a * d[1]] #posの更新
            pos = pos_new
            pos_history.append(pos)
    return pos_history
(H_1,H_2) = grad_descent((1.2, 1.2), 0.1, 0.1)
plt.plot(H_1,H_2)
plt.show()

試したこと

ここに問題に対して試したことを記載してください。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

おそらく
return 100(x_2 - (x_1)**2)**2 + (1 - x_1)**2
の100の後に演算子が抜けているのでは？

Comments

[aaa-1]

素早い回答ありがとうございます！
助かりました。

Answer 2

return 100(x_2 - (x_1)**2)**2 + (1 - x_1)**2

この部分が、「100という関数を呼び出す」と認識されてエラーとなっているようです。100 * (x-2...のようにしましょう。

Comments

[aaa-1]

素早い解答ありがとうございます！
助かりました。