Add option for hipo_scaling

Author: filikat

docs/src/options/definitions.md

@@ -445,6 +445,11 @@
 - Type: boolean
 - Default: false
 
+## hipo\_scaling
+- HiPO scaling strategy: "curtis-reid" or "norm"
+- Type: string
+- Default: "curtis-reid"
+
 ## pdlp\_scaling
 - Scaling option for PDLP solver: Default = true
 - Type: boolean

highs/ipm/IpxWrapper.cpp

@@ -552,6 +552,17 @@ HighsStatus solveLpHipo(const HighsOptions& options, HighsTimer& timer,
     return HighsStatus::kError;
   }
 
+  if (options.hipo_scaling == kHipoCRscaling) {
+    hipo_options.scaling = hipo::kOptionCRscaling;
+  } else if (options.hipo_scaling == kHipoNormScaling) {
+    hipo_options.scaling = hipo::kOptionNormScaling;
+  } else {
+    highsLogUser(options.log_options, HighsLogType::kError,
+                 "Unknown value of option %s\n", kHipoScalingString.c_str());
+    model_status = HighsModelStatus::kSolveError;
+    return HighsStatus::kError;
+  }
+
   // block size option
   hipo_options.block_size = options.hipo_block_size;
 

highs/ipm/hipo/ipm/Options.h

@@ -35,11 +35,20 @@ enum OptionParallel {
   kOptionParallelDefault = kOptionParallelChoose
 };
 
+enum OptionScaling {
+  kOptionScalingMin = 0,
+  kOptionCRscaling = kOptionScalingMin,
+  kOptionNormScaling,
+  kOptionScalingMax = kOptionNormScaling,
+  kOptionScalingDefault = kOptionCRscaling
+};
+
 struct Options {
   // Solver options
   OptionNla nla = kOptionNlaDefault;
   OptionCrossover crossover = kOptionCrossoverDefault;
   OptionParallel parallel = kOptionParallelDefault;
+  OptionScaling scaling = kOptionScalingDefault;
 
   // Metis options
   bool metis_no2hop = false;

highs/lp_data/HighsOptions.h

@@ -307,6 +307,10 @@ const string kHipoTreeString = "tree";
 const string kHipoNodeString = "node";
 const string kHipoBothString = "both";
 
+const string kHipoScalingString = "hipo_scaling";
+const string kHipoCRscaling = "curtis-reid";
+const string kHipoNormScaling = "norm";
+
 struct HighsOptionsStruct {
   // Run-time options read from the command line
   std::string presolve;
@@ -373,6 +377,7 @@ struct HighsOptionsStruct {
   std::string hipo_parallel_type;
   HighsInt hipo_block_size;
   bool hipo_metis_no2hop;
+  std::string hipo_scaling;
 
   // Options for PDLP solver
   bool pdlp_scaling;
@@ -545,6 +550,7 @@ struct HighsOptionsStruct {
         hipo_parallel_type(""),
         hipo_block_size(0),
         hipo_metis_no2hop(false),
+        hipo_scaling(""),
         pdlp_scaling(false),
         pdlp_iteration_limit(0),
         pdlp_e_restart_method(0),
@@ -1272,6 +1278,13 @@ class HighsOptions : public HighsOptionsStruct {
                              advanced, &hipo_metis_no2hop, false);
     records.push_back(record_bool);
 
+    record_string =
+        new OptionRecordString(kHipoScalingString,
+                               "HiPO scaling option: \"curtis-reid\" "
+                               "or \"norm\".",
+                               advanced, &hipo_scaling, kHipoCRscaling);
+    records.push_back(record_string);
+
     record_bool = new OptionRecordBool(
         "pdlp_scaling", "Scaling option for PDLP solver: Default = true",
         advanced, &pdlp_scaling, true);