Filter B_flow in the highest resolution expressions (#333)

TulipaEnergy · Dec 5, 2023 · feaf33c · feaf33c
1 parent 3e91b75
commit feaf33c
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/model.jl b/src/model.jl
@@ -209,8 +209,9 @@ function create_model(
         model,
         incoming_flow_highest_resolution[a ∈ A, rp ∈ RP, B ∈ Ph[(a, rp)]],
         sum(
-            duration(B, B_flow, rp) * flow[(u, a), rp, B_flow] for u in inneighbor_labels(graph, a),
-            B_flow ∈ graph[u, a].partitions[rp]
+            duration(B, B_flow, rp) * flow[(u, a), rp, B_flow] for
+            u in inneighbor_labels(graph, a), B_flow ∈ graph[u, a].partitions[rp] if
+            B_flow[end] ≥ B[1] && B[end] ≥ B_flow[1]
         )
     )
 
@@ -219,7 +220,8 @@ function create_model(
         outgoing_flow_highest_resolution[a ∈ A, rp ∈ RP, B ∈ Ph[(a, rp)]],
         sum(
             duration(B, B_flow, rp) * flow[(a, v), rp, B_flow] for
-            v in outneighbor_labels(graph, a), B_flow ∈ graph[a, v].partitions[rp]
+            v in outneighbor_labels(graph, a), B_flow ∈ graph[a, v].partitions[rp] if
+            B_flow[end] ≥ B[1] && B[end] ≥ B_flow[1]
         )
     )