multiply gn return by 2 for consistency

EdiGlacUQ · Jul 3, 2021 · eac0786 · eac0786
1 parent f16d344
commit eac0786
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/runs/run_eigendec.py b/runs/run_eigendec.py
@@ -117,7 +117,7 @@ def ghep_action(x):
             """Hessian action w/o preconditioning"""
             ddJ_val = slvr.H_GN.action(cntrl, x)
             # reg_op.inv_action(ddJ_val.vector(), xg.vector()) <- gnhep_prior
-            return function_get_values(ddJ_val)
+            return 2.0 * function_get_values(ddJ_val)
         else:
             """Hessian action w/o preconditioning"""
             _, _, ddJ_val = slvr.ddJ.action(cntrl, x)