Commit a4555703 authored by Roger.Grau.Haro's avatar Roger.Grau.Haro
Browse files

test

parent 80644400
Pipeline #6888 failed with stages
in 0 seconds
......@@ -131,20 +131,23 @@ for name in sorted(glob.glob(sys.argv[1])):
x_min //= 1000
x_max //= 1000
arr_y = np.zeros_like(x)
arr_y = np.zeros_like(x[ids])
x+=2200
x/=1000
model = models.Models()
for idx, i in enumerate(x):
xs = np.zeros_like(x[ids])
for idx, i in enumerate(x[ids]):
arr_y[idx] = ap.eval_full(list(ap.entry_data.keys())[0], i * u.km)
xs[idx] = i
print(arr_y[ids].shape, sr.shape)
print(arr_y.shape, sr.shape)
ans = curve_fit(model.quadratic_sin_model, x[ids], sr)
print(ans)
y_fit = model.quadratic_sin_model(x[ids], *ans[0])
print(*ans[0])
plt.plot(xs, arr_y)
plt.show()
plt.plot(x[ids], sr)
plt.plot(x[ids], y_fit)
plt.show()
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment