Reverse indexation.

Author: SebKrantz

src/models/solve_allocation_by_duality.jl

@@ -181,7 +181,7 @@ function hessian_duality(
         ind = 0
 
         # https://stackoverflow.com/questions/38901275/inbounds-propagation-rules-in-julia
-        @inbounds for (jdnd, jn) in zip(hess_str[1], hess_str[2])
+        @inbounds for (jdnd, jn) in zip(hess_str[2], hess_str[1])
             ind += 1
             # First get the indices of the element and the respective derivative 
             j = (jn-1) % J + 1

src/models/solve_allocation_by_duality_cgc.jl

@@ -168,7 +168,7 @@ function hessian_duality_cgc(
         ind = 0
 
         # https://stackoverflow.com/questions/38901275/inbounds-propagation-rules-in-julia
-        @inbounds for (jdnd, jn) in zip(hess_str[1], hess_str[2])
+        @inbounds for (jdnd, jn) in zip(hess_str[2], hess_str[1])
             ind += 1
             # First get the indices of the element and the respective derivative 
             j = (jn-1) % J + 1