From 771f91954becc937ce3c5dcdb567587ae80f13c3 Mon Sep 17 00:00:00 2001
From: Eric Kooistra <kooistra@astron.nl>
Date: Thu, 20 Jun 2024 14:50:10 +0200
Subject: [PATCH] Support verbosity.
---
applications/lofar2/model/rtdsp/fourier.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/applications/lofar2/model/rtdsp/fourier.py b/applications/lofar2/model/rtdsp/fourier.py
index c2aad788ce..bfdc16b160 100644
--- a/applications/lofar2/model/rtdsp/fourier.py
+++ b/applications/lofar2/model/rtdsp/fourier.py
@@ -34,7 +34,7 @@ from .utilities import c_rtol, c_atol, ceil_pow2, is_even
# Time domain interpolation using DFT
###############################################################################
-def fourier_interpolate(HFfilter, Ncoefs):
+def fourier_interpolate(HFfilter, Ncoefs, verbosity=1):
"""Use Fourier interpolation to create final FIR filter coefs.
HF contains filter transfer function for N points, in order 0 to fs. The
@@ -82,7 +82,8 @@ def fourier_interpolate(HFfilter, Ncoefs):
# K + 1 + K = N values, because N is odd and K = N // 2
hInterpolated = np.fft.ifft(HFextended)
if np.allclose(hInterpolated.imag, np.zeros(Ncoefs), rtol=c_rtol, atol=c_atol):
- print('hInterpolated.imag ~= 0')
+ if verbosity:
+ print('hInterpolated.imag ~= 0')
else:
print('WARNING: hInterpolated.imag != 0 (max(abs) = %e)' % np.max(np.abs(hInterpolated.imag)))
return hInterpolated.real
--
GitLab