From d851b6f743c75f5689e8ea6b3a8e60279aece5b7 Mon Sep 17 00:00:00 2001 From: landmanbester Date: Mon, 25 Mar 2024 12:17:28 +0200 Subject: [PATCH] use abs(model) for hard-thresholding in model2comps --- pfb/workers/model2comps.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/pfb/workers/model2comps.py b/pfb/workers/model2comps.py index 09eb40204..11c3cc011 100644 --- a/pfb/workers/model2comps.py +++ b/pfb/workers/model2comps.py @@ -121,7 +121,7 @@ def _model2comps(**kw): wsums[...] = 1.0 if opts.min_val is not None: - model = np.where(model >= opts.min_val, model, 0.0) + model = np.where(np.abs(model) >= opts.min_val, model, 0.0) if not np.any(model): raise ValueError(f'Model is empty or has no components above {opts.min_val}')