From f4b66ee2c0a446e78d505c7d8801af447f323f44 Mon Sep 17 00:00:00 2001 From: Harry Ho Date: Tue, 14 Dec 2021 12:34:18 +0800 Subject: [PATCH] analyze_sayma_data: Add max abs dev --- analyze_sayma_data.py | 3 +++ 1 file changed, 3 insertions(+) diff --git a/analyze_sayma_data.py b/analyze_sayma_data.py index f00d753..95cb8e5 100644 --- a/analyze_sayma_data.py +++ b/analyze_sayma_data.py @@ -53,6 +53,7 @@ def report(*args, **kwargs): rpt.append(" Minimum: {:>10.4f} ps".format(prm["phase_ps_min"])) rpt.append(" Maximum: {:>10.4f} ps".format(prm["phase_ps_max"])) rpt.append(" Standard Deviation: {:>10.4f} ps".format(prm["phase_ps_std"])) + rpt.append(" Max Absolute Deviation: {:>10.4f} ps".format(prm["phase_ps_maxabsdev"])) rpt.append(" Mean Absolute Deviation: {:>10.4f} ps".format(prm["phase_ps_meanabsdev"])) rpt.append("") # TODO: Use jinja2 to produce a report @@ -110,6 +111,8 @@ def main(): rpt_params["phase_ps_min"] = data_phase_ps.min() rpt_params["phase_ps_max"] = data_phase_ps.max() rpt_params["phase_ps_std"] = data_phase_ps.std() + rpt_params["phase_ps_maxabsdev"] = \ + np.absolute(data_phase_ps - data_phase_ps.max()).mean() rpt_params["phase_ps_meanabsdev"] = \ np.absolute(data_phase_ps - data_phase_ps.mean()).mean() # Generate the report