need to apply abs to test difference of mask with code

EdiGlacUQ · May 27, 2024 · a8accc9 · a8accc9
1 parent f6e3c3b
commit a8accc9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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