def check(f, x0, minx):
tracex = minimize(f, x0, LEARNING_RATE, h, PRECISION)
start = tracex[0]
stop = tracex[-1]
assert round(stop*1000) == round(minx*1000)
def check(f, x0, minx):
tracex = minimize(f, x0, LEARNING_RATE, h, PRECISION)
start = tracex[0]
stop = tracex[-1]
assert round(stop * 10000) == round(minx * 10000)
(cosf, 1.2, 0.98865134),
(cosf, 0.05, 0.29691298),
(x2, 2.0, 2.0),
(x2, -1.0, 2.0),
(x2, -9.3, 2.0),
(x3, 2.0, 1 / 3.0),
(x3, 0.5, 1 / 3.0),
(x3, 1 / 3.0, 1 / 3.0),
(x3, 1.0, 1 / 3.0),
]
print """def cosf(x): return np.cos(3 * np.pi * x) / x
def x2(x): return (x-2)**2 + 1
def x3(x): return 5*x**3 + 2*x**2 - 3*x
"""
for t in tests:
f = t[0]
x0 = t[1]
minx = t[2]
tracex = minimize(f, x0, LEARNING_RATE, h, PRECISION)
start = tracex[0]
stop = tracex[-1]
# print "Start f(x=%2.8f)=%2.8f" % (start, f(start))
print "x0 = %1.2f; min = %s(x=%2.3f)=%2.3f" % (x0, f.__name__, stop, f(stop))
