diff --git a/app/cocal_methods.py b/app/cocal_methods.py index d49c55c0a3c6d327c52e481a3810e4d820b5002b..3a5e7371a0ee66d11e25de43476ae12dcb5b88c1 100644 --- a/app/cocal_methods.py +++ b/app/cocal_methods.py @@ -833,7 +833,7 @@ class CocalMethods: theta0=theta0, Nb=Nb, w_empirical_exp=w_empirical_exp, - weights=weights[mask], + weights=None, #weights[mask], ) umc_kwargs = { @@ -897,6 +897,16 @@ class CocalMethods: method="Powell", ) + # plotting code for debuggung purposes + # H_res = np.squeeze(self.freqz_core(res.x, Nb, w_empirical_exp)) + # H_0 = np.squeeze(self.freqz_core(theta0, Nb, w_empirical_exp)) + # plt.loglog(f, np.abs(H_empirical), label="empirical") + # plt.loglog(f, np.abs(H_res), label="fitted") + # plt.loglog(f, np.abs(H_0), label="init") + # plt.legend() + # plt.show() + + # alternative solvers # res = opt.minimize(self.evaluate_filter, x0=theta0, args=(Nb, H_empirical, w_empirical_exp)) # res = opt.minimize(self.evaluate_filter, x0=theta0, args=(Nb, H_empirical, w_empirical_exp), method="TNC")