Add option for hipo_scaling

Author: filikat

docs/src/options/definitions.md

@@ -445,6 +445,11 @@
 - Range: {0, 2147483647}
 - Default: 128
 
+## 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

@@ -564,6 +564,17 @@ HighsStatus solveLpHipo(const HighsOptions& options, HighsTimer& timer,
   }
   hipo_options.ordering = options.hipo_ordering;
 
+  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,12 +35,21 @@ 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;
   std::string ordering = kHighsChooseString;
+  OptionScaling scaling = kOptionScalingDefault;
 
   // Ipm parameters
   Int max_iter = kMaxIterDefault;

highs/lp_data/HighsOptions.h

@@ -313,6 +313,10 @@ const string kHipoMetisString = "metis";
 const string kHipoAmdString = "amd";
 const string kHipoRcmString = "rcm";
 
+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;
@@ -379,6 +383,7 @@ struct HighsOptionsStruct {
   std::string hipo_parallel_type;
   std::string hipo_ordering;
   HighsInt hipo_block_size;
+  std::string hipo_scaling;
 
   // Options for PDLP solver
   bool pdlp_scaling;
@@ -551,6 +556,7 @@ struct HighsOptionsStruct {
         hipo_parallel_type(""),
         hipo_ordering(""),
         hipo_block_size(0),
+        hipo_scaling(""),
         pdlp_scaling(false),
         pdlp_iteration_limit(0),
         pdlp_e_restart_method(0),
@@ -1289,6 +1295,13 @@ class HighsOptions : public HighsOptionsStruct {
         advanced, &hipo_block_size, 0, 128, kHighsIInf);
     records.push_back(record_int);
 
+    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);