From eae6047c5a676bf7366f379e9fcb62bf9f75a2c8 Mon Sep 17 00:00:00 2001
From: Maximilian Gruber <maximilian.gruber@ptb.de>
Date: Fri, 31 May 2024 15:17:56 +0200
Subject: [PATCH] undo use of weights in fitting process

---
 app/cocal_methods.py | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/app/cocal_methods.py b/app/cocal_methods.py
index d49c55c..3a5e737 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")
-- 
GitLab