Now calling ekk_instance_.computeBasisCondition(); but LP needs to be…

… passed
ERGO-Code · Sep 18, 2024 · 091ab68 · 091ab68
1 parent ccef75c
commit 091ab68
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/src/lp_data/Highs.cpp b/src/lp_data/Highs.cpp
@@ -1948,7 +1948,13 @@ HighsStatus Highs::getReducedColumn(const HighsInt col, double* col_vector,
 }
 
 HighsStatus Highs::getKappa(double& kappa) {
-  return HighsStatus::kError;
+  printf("Highs::getKappa basis_.valid = %d, ekk_instance_.status_.has_invert = %d\n",
+	 int(basis_.valid),
+	 int(ekk_instance_.status_.has_invert));
+  if (!ekk_instance_.status_.has_invert)
+    return invertRequirementError("getBasisInverseRow");
+  kappa = ekk_instance_.computeBasisCondition();
+  return HighsStatus::kOk;
 }
 
 HighsStatus Highs::setSolution(const HighsSolution& solution) {

diff --git a/src/simplex/HEkk.h b/src/simplex/HEkk.h
@@ -134,6 +134,7 @@ class HEkk {
   HighsBasis getHighsBasis(HighsLp& use_lp) const;
 
   const SimplexBasis& getSimplexBasis() { return basis_; }
+  double computeBasisCondition();
 
   HighsStatus initialiseSimplexLpBasisAndFactor(
       const bool only_from_known_basis = false);
@@ -333,7 +334,6 @@ class HEkk {
   HighsStatus returnFromEkkSolve(const HighsStatus return_status);
   HighsStatus returnFromSolve(const HighsStatus return_status);
 
-  double computeBasisCondition();
   void initialiseAnalysis();
   std::string rebuildReason(const HighsInt rebuild_reason);