Merge pull request ERGO-Code#1773 from ERGO-Code/fix-1694

Fix 1694

Author: jajhall

check/TestMipSolver.cpp

@@ -665,7 +665,8 @@ TEST_CASE("IP-with-fract-bounds-no-presolve", "[highs_test_mip_solver]") {
   lp.col_cost_ = {1, -2, 3};
   lp.col_lower_ = {2.5, 2.5, 2.5};
   lp.col_upper_ = {6.5, 5.5, 7.5};
-  lp.integrality_ = {HighsVarType::kInteger, HighsVarType::kInteger, HighsVarType::kInteger};
+  lp.integrality_ = {HighsVarType::kInteger, HighsVarType::kInteger,
+                     HighsVarType::kInteger};
 
   // Solve
   highs.passModel(lp);

src/pdlp/CupdlpWrapper.cpp

@@ -73,8 +73,18 @@ HighsStatus solveLpCupdlp(const HighsOptions& options, HighsTimer& timer,
   void* presolvedmodel = NULL;
   void* model2solve = NULL;
 
+  // WIP on zbook?
+  //
+  //  HighsInt size_of_CUPDLPscaling = sizeof(CUPDLPscaling);
+  //
   CUPDLPscaling* scaling = (CUPDLPscaling*)cupdlp_malloc(sizeof(CUPDLPscaling));
 
+  // WIP on zbook?
+  //
+  //  printf("size_of_CUPDLPscaling = %d\n", size_of_CUPDLPscaling);
+  //  scaling->ifRuizScaling = 1;
+  //  printf("scaling->ifRuizScaling = %d\n", scaling->ifRuizScaling);
+
   // claim solvers variables
   // prepare pointers
   CUPDLP_MATRIX_FORMAT src_matrix_format = CSC;
@@ -103,7 +113,7 @@ HighsStatus solveLpCupdlp(const HighsOptions& options, HighsTimer& timer,
   const cupdlp_int local_log_level = getCupdlpLogLevel(options);
   if (local_log_level) cupdlp_printf("Solving with cuPDLP-C\n");
 
-  Init_Scaling(local_log_level, scaling, nCols, nRows, cost, rhs);
+  H_Init_Scaling(local_log_level, scaling, nCols, nRows, cost, rhs);
   cupdlp_int ifScaling = 1;
 
   CUPDLPwork* w = cupdlp_NULL;
@@ -125,8 +135,8 @@ HighsStatus solveLpCupdlp(const HighsOptions& options, HighsTimer& timer,
   memcpy(csc_cpu->colMatElem, csc_val, nnz * sizeof(double));
 
   cupdlp_float scaling_time = getTimeStamp();
-  PDHG_Scale_Data_cuda(local_log_level, csc_cpu, ifScaling, scaling, cost,
-                       lower, upper, rhs);
+  H_PDHG_Scale_Data_cuda(local_log_level, csc_cpu, ifScaling, scaling, cost,
+                         lower, upper, rhs);
   scaling_time = getTimeStamp() - scaling_time;
 
   cupdlp_float alloc_matrix_time = 0.0;

src/pdlp/cupdlp/README.md

@@ -101,7 +101,21 @@ The `cuPDLP-c` iteration count is held in `pdhg->timers->nIter`, but `pdhg` is d
 
 - Make CupldlpWrapper.cpp look more like C++ than C
 
-
-
-
-
++## Using a GPU
++
++### Install CUDA
++
++* sudo apt update && sudo apt upgrade
++* sudo apt autoremove nvidia* --purge
++* sudo apt update && sudo apt upgrade
++* nvcc --version
++* sudo apt install nvidia-cuda-toolkit
++
++### Building PDLP
++
++export HIGHS_HOME=/home/jajhall/install
++export CUDA_HOME=/usr/lib/cuda
++
++
++
++

src/pdlp/cupdlp/cupdlp_scaling_cuda.c

@@ -228,11 +228,11 @@ cupdlp_retcode cupdlp_pc_scaling_cuda(CUPDLPcsc *csc, cupdlp_float *cost,
   return retcode;
 }
 
-cupdlp_retcode PDHG_Scale_Data_cuda(cupdlp_int log_level,
-				    CUPDLPcsc *csc, cupdlp_int ifScaling,
-                                    CUPDLPscaling *scaling, cupdlp_float *cost,
-                                    cupdlp_float *lower, cupdlp_float *upper,
-                                    cupdlp_float *rhs) {
+cupdlp_retcode H_PDHG_Scale_Data_cuda(cupdlp_int log_level,
+				      CUPDLPcsc *csc, cupdlp_int ifScaling,
+				      CUPDLPscaling *scaling, cupdlp_float *cost,
+				      cupdlp_float *lower, cupdlp_float *upper,
+				      cupdlp_float *rhs) {
   cupdlp_retcode retcode = RETCODE_OK;
   // scaling->dObjScale = 1.0;
 
@@ -395,10 +395,10 @@ cupdlp_retcode PDHG_Scale_Data_cuda(cupdlp_int log_level,
   return retcode;
 }
 
-cupdlp_retcode Init_Scaling(cupdlp_int log_level,
-			    CUPDLPscaling *scaling, cupdlp_int ncols,
-                            cupdlp_int nrows, cupdlp_float *cost,
-                            cupdlp_float *rhs) {
+cupdlp_retcode H_Init_Scaling(cupdlp_int log_level,
+			      CUPDLPscaling *scaling, cupdlp_int ncols,
+			      cupdlp_int nrows, cupdlp_float *cost,
+			      cupdlp_float *rhs) {
   cupdlp_retcode retcode = RETCODE_OK;
 
   scaling->ifRuizScaling = 1;

src/pdlp/cupdlp/cupdlp_scaling_cuda.h

@@ -11,17 +11,17 @@
 extern "C" {
 #endif
 
-cupdlp_retcode PDHG_Scale_Data_cuda(cupdlp_int log_level,
-				    CUPDLPcsc *csc, cupdlp_int ifScaling,
-                                    CUPDLPscaling *scaling, cupdlp_float *cost,
-                                    cupdlp_float *lower, cupdlp_float *upper,
-                                    cupdlp_float *rhs);
-
-cupdlp_retcode Init_Scaling(cupdlp_int log_level,
-			    CUPDLPscaling *scaling, cupdlp_int ncols,
-                            cupdlp_int nrows, cupdlp_float *cost,
-                            cupdlp_float *rhs);
+cupdlp_retcode H_PDHG_Scale_Data_cuda(cupdlp_int log_level,
+				      CUPDLPcsc *csc, cupdlp_int ifScaling,
+				      CUPDLPscaling *scaling, cupdlp_float *cost,
+				      cupdlp_float *lower, cupdlp_float *upper,
+				      cupdlp_float *rhs);
 
+cupdlp_retcode H_Init_Scaling(cupdlp_int log_level,
+			      CUPDLPscaling *scaling, cupdlp_int ncols,
+			      cupdlp_int nrows, cupdlp_float *cost,
+			      cupdlp_float *rhs);
+  
 #ifdef __cplusplus
 }
 #endif