From a8accc97e8f9821ac6ae5bcbc28e79dcb1c457cc Mon Sep 17 00:00:00 2001 From: dngoldberg <31406058+dngoldberg@users.noreply.github.com> Date: Mon, 27 May 2024 13:54:18 +0100 Subject: [PATCH] need to apply abs to test difference of mask with code --- fenics_ice/solver.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/fenics_ice/solver.py b/fenics_ice/solver.py index b8d112d..9c9b77e 100644 --- a/fenics_ice/solver.py +++ b/fenics_ice/solver.py @@ -1460,7 +1460,7 @@ def comp_Q_vaf(self, verbose=False): if self.params.error_prop.qoi_apply_vaf_mask: if self.params.error_prop.qoi_vaf_mask_usecode: code = float(self.params.error_prop.qoi_vaf_mask_code) - msk_ex = conditional((self.vaf_mask-code)<1.e-10, 1.0, 0.0) + msk_ex = conditional(abs(self.vaf_mask-code)<1.e-10, 1.0, 0.0) else: msk_ex = conditional(self.vaf_mask > 0.0, 1.0, 0.0) Q_vaf = msk_ex * HAF * dx