Remove constant factors: seems to work without.

Author: SebKrantz

src/models/solve_allocation_by_duality.jl

@@ -236,7 +236,7 @@ function hessian_duality(
             end
 
             # Assign result
-            values[ind] = -obj_factor * term + 1e-5 # Somehow need this increment to make it work
+            values[ind] = -obj_factor * term # + 1e-5 # Somehow need this increment to make it work
         end
     end
     return

src/models/solve_allocation_by_duality_cgc.jl

@@ -230,7 +230,7 @@ function hessian_duality_cgc(
             end
 
             # Assign result
-            values[ind] = -obj_factor * term + 1e-6
+            values[ind] = -obj_factor * term # + 1e-6
         end
     end
     return