Update comments. Remove unneeded import

Author: Opt-Mucca

highs/mip/HighsCutGeneration.h

@@ -20,7 +20,6 @@
 #include "util/HighsCDouble.h"
 #include "util/HighsInt.h"
 #include "util/HighsRandom.h"
-#include "util/HighsSparseVectorSum.h"
 
 class HighsLpRelaxation;
 class HighsTransformedLp;
@@ -119,20 +118,20 @@ class HighsCutGeneration {
 
   /// Single Node Flow Relaxation for flow cover cuts
   struct SNFRelaxation {
-    HighsInt numNnzs;                    // number of nonzeros
-    std::vector<HighsInt> coef;          // coefficients of cols in SNFR
-    std::vector<double> vubCoef;         // coefficients in vub of cols in SNFR
-    std::vector<double> binSolval;       // vub bin col sol in SNFR
-    std::vector<double> contSolval;      // real sol in SNFR
-    std::vector<HighsInt> origBinCols;   // orig bin col used in SNFR
-    std::vector<HighsInt> origContCols;  // orig cont cols used in SNFR
-    std::vector<double> aggrBinCoef;     // aggr coef of orignal bin-col in SNFR
-    std::vector<double> aggrContCoef;  // aggr coef of original cont-col in SNFR
-    std::vector<double> aggrConstant;  // aggr original constant in SNFR
-
-    std::vector<HighsInt> flowCoverStatus;  // (+1) in fcover (-1) not in fcover
-    double rhs;
-    double lambda;
+    HighsInt numNnzs; // |N-| + |N+|
+    std::vector<HighsInt> coef; // (+-1) coefficient of col in SNFR
+    std::vector<double> vubCoef; // u_j in y'_j <= u_j x_j in SNFR
+    std::vector<double> binSolval; // lp[x_j], y'_j <= u_j x_j in SNFR
+    std::vector<double> contSolval; // lp[y'_j] in y'_j <= u_j x_j in SNFR
+    std::vector<HighsInt> origBinCols; // orig x_i, y'_j <= u_j x_j in SNFR
+    std::vector<HighsInt> origContCols; // orig y_i used to make y'_j in SNFR
+    std::vector<double> aggrBinCoef; // c_i row coef of x_i in orig aggrrow
+    std::vector<double> aggrContCoef; // a_i row coef of y_i in orig aggrrow
+    std::vector<double> aggrConstant; // constant shift used in SNFR transform
+
+    std::vector<HighsInt> flowCoverStatus; // (+1) in f-cover (-1) notin f-cover
+    double rhs; // in \sum_{j \in N+} y'_j - \sum_{j \in N-} y'_j <= b
+    double lambda; // in sum_{j in C+} u_j - sum_{j in C-} u_j = b + lambda
   };
 
  private:

highs/mip/HighsTransformedLp.h

@@ -23,7 +23,6 @@
 #include "util/HighsInt.h"
 #include "util/HighsSparseVectorSum.h"
 
-class HighsCutGeneration;
 class HighsLpRelaxation;
 
 /// Helper class to compute single-row relaxations from the current LP