def test_fit_engine_list():
a = 1
b = 2
c = 3
m = quadratic_model('',
a, 'free', [-1, 1],
b, 'free', [-1, 1],
c, 'free', [-1, 1])
x = np.arange(-1, 1, 0.01)
y = x**2 + 1
datav = [(x, y), (x, y+2)]
out = fit_engine_list(m, datav)
assert_equal(len(out), 2)
def test_quadratic_fit():
x = np.arange(-1, 1, .01)
a = 1
b = 2
c = 3
true_val = [a, b, c]
y = a * x*x + b * x + c
m = quadratic_model('',
a, 'free', [-1, 1],
b, 'free', [-1, 1],
c, 'free', [-1, 1])
result, yfit = fit_engine(m, x, y)
out = result.values
fitted_val = (out['a'], out['b'], out['c'])
assert_array_almost_equal(true_val, fitted_val)
