From 389f6ba0602c9ba86928da1ed6848eac9a090834 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?K=C3=BChle=2C=20Laura=20Christine=20=28lakue103=29?=
<laura.kuehle@uni-duesseldorf.de>
Date: Tue, 25 May 2021 23:08:07 +0200
Subject: [PATCH] Fixed bug in calculation of 'max_avg'.
---
Troubled_Cell_Detector.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/Troubled_Cell_Detector.py b/Troubled_Cell_Detector.py
index 521629a..915317d 100644
--- a/Troubled_Cell_Detector.py
+++ b/Troubled_Cell_Detector.py
@@ -3,6 +3,7 @@
@author: Laura C. Kühle, Soraya Terrab (sorayaterrab)
TODO: Fix cell averages and reconstructions to create data with an x-point stencil
+TODO: Fix bug for max_avg calculation -> Done
"""
import os
@@ -437,7 +438,7 @@ class Theoretical(WaveletDetector):
def _get_cells(self, multiwavelet_coeffs, projection):
troubled_cells = []
- max_avg = max(1, max(abs(projection[0][degree]) for degree in range(self._polynomial_degree+1)))
+ max_avg = np.sqrt(0.5) * max(1, max(abs(projection[0][cell+1]) for cell in range(self._num_coarse_grid_cells)))
for cell in range(self._num_coarse_grid_cells):
if self._is_troubled_cell(multiwavelet_coeffs, cell, max_avg):
--
GitLab