From 02f25498287111dce0057d2f257abcd515c4f2b4 Mon Sep 17 00:00:00 2001
From: Abdoulbari Zakir <abdoulbari.zakir@gmail.com>
Date: Mon, 16 Dec 2024 12:28:20 +0100
Subject: [PATCH] compute rhs

---
 ...opt_gestion_second_membre_cas_lineaire.cpp | 27 ++++++++++++++++++-
 1 file changed, 26 insertions(+), 1 deletion(-)

diff --git a/src/solver/optimisation/opt_gestion_second_membre_cas_lineaire.cpp b/src/solver/optimisation/opt_gestion_second_membre_cas_lineaire.cpp
index 36bb303964..407074bd5e 100644
--- a/src/solver/optimisation/opt_gestion_second_membre_cas_lineaire.cpp
+++ b/src/solver/optimisation/opt_gestion_second_membre_cas_lineaire.cpp
@@ -43,6 +43,28 @@ static void shortTermStorageLevelsRHS(
     }
 }
 
+static void shortTermStorageCumulationRHS(
+        const std::vector<::ShortTermStorage::AREA_INPUT>& shortTermStorageInput,
+        int numberOfAreas,
+        std::vector<double>& SecondMembre,
+        const CORRESPONDANCES_DES_CONTRAINTES_HEBDOMADAIRES&
+        CorrespondancesDesContraintesHebdomadaires)
+{
+    for (int areaIndex = 0; areaIndex < numberOfAreas; areaIndex++)
+    {
+        for (auto& storage: shortTermStorageInput[areaIndex])
+        {
+            for (const auto& constraint: storage.additional_constraints)
+            {
+                const int clusterGlobalIndex = storage.clusterGlobalIndex;
+                int cnt = CorrespondancesDesContraintesHebdomadaires.ShortTermStorageCumulation
+                        [constraint.globalIndex];
+                SecondMembre[cnt] = constraint.rhs;
+            }
+        }
+    }
+}
+
 void OPT_InitialiserLeSecondMembreDuProblemeLineaire(PROBLEME_HEBDO* problemeHebdo,
                                                      int PremierPdtDeLIntervalle,
                                                      int DernierPdtDeLIntervalle,
@@ -145,7 +167,6 @@ void OPT_InitialiserLeSecondMembreDuProblemeLineaire(PROBLEME_HEBDO* problemeHeb
                                   ProblemeAResoudre->SecondMembre,
                                   CorrespondanceCntNativesCntOptim,
                                   hourInTheYear);
-
         for (uint32_t interco = 0; interco < problemeHebdo->NombreDInterconnexions; interco++)
         {
             if (const COUTS_DE_TRANSPORT& CoutDeTransport = problemeHebdo->CoutDeTransport[interco];
@@ -377,6 +398,10 @@ void OPT_InitialiserLeSecondMembreDuProblemeLineaire(PROBLEME_HEBDO* problemeHeb
         }
     }
 
+    shortTermStorageCumulationRHS(problemeHebdo->ShortTermStorage,
+                                  problemeHebdo->NombreDePays,
+                                  ProblemeAResoudre->SecondMembre,
+                                  problemeHebdo->CorrespondanceCntNativesCntOptimHebdomadaires);
     if (problemeHebdo->OptimisationAvecCoutsDeDemarrage)
     {
         OPT_InitialiserLeSecondMembreDuProblemeLineaireCoutsDeDemarrage(problemeHebdo,