scs 0.4.1 → 0.4.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 9005db60b3ebc21bfe74e18fbd90db752801cce6232c5dcbbf9746ed6a3772b9
-  data.tar.gz: 8d1eade7f1501cbf46e820a2e7b20133f008e14d2bcc098629643a9b750b55b3
+  metadata.gz: 82ca3e4eb97e0a3d55e0755babc5ba5cca8d0eb0c171a0733c38916726697bd5
+  data.tar.gz: f91bdb81c22052f788c615dd9bc598941741f1fdba59d741bcaa4b6633c1eae5
 SHA512:
-  metadata.gz: 98ecf4f8beb7215e0ca44a57e3345b257054948c311e8e749f31d698e27c7f6c66ccbaaa71bbb091a92d4a23809b9b98a5d65767ee8d58e78cc7a8f08a32e042
-  data.tar.gz: 4701d88bad2ae68509fba757a922019400a78f8d8c2bdf0e33f35f7e409c49b0e8306d0cc3281a542508e3a73176042bf16a6039aee5101430ff378689791686
+  metadata.gz: 835cbd3b7a6fcd085663ef49c54d10c2bfecee3c2435b2b873cfeab01efa90dba7f7a0f59f4487dad15135570399c75c002864dadc4031dda83af1de0abb509a
+  data.tar.gz: a9401bac5dd33e24e140061d66cc52fe02c2982bc63b96ebf44f5c13affc8e6f8477bf4907043902278b90e673aca848e77bcad8fc90f2f94550729b25a5de81

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## 0.4.2 (2024-07-10)
+
+- Updated SCS to 3.2.6
+
 ## 0.4.1 (2023-04-09)
 
 - Updated SCS to 3.2.3

data/LICENSE.txt

@@ -1,7 +1,7 @@
 The MIT License (MIT)
 
 Copyright (c) 2012 Brendan O'Donoghue (bodonoghue85@gmail.com)
-Copyright (c) 2019-2023 Andrew Kane
+Copyright (c) 2019-2024 Andrew Kane
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal

data/lib/scs/ffi.rb

@@ -36,7 +36,7 @@ module SCS
           "ScsMatrix *a",
           "ScsMatrix *p",
           "scs_float *b",
-          "scs_float *c",
+          "scs_float *c"
         ]
 
         m::Cone = struct [
@@ -52,7 +52,7 @@ module SCS
           "scs_int ep",
           "scs_int ed",
           "scs_float *p",
-          "scs_int psize",
+          "scs_int psize"
         ]
 
         m::Solution = struct [

data/lib/scs/version.rb

@@ -1,3 +1,3 @@
 module SCS
-  VERSION = "0.4.1"
+  VERSION = "0.4.2"
 end

data/lib/scs.rb

@@ -2,9 +2,9 @@
 require "fiddle/import"
 
 # modules
-require "scs/matrix"
-require "scs/solver"
-require "scs/version"
+require_relative "scs/matrix"
+require_relative "scs/solver"
+require_relative "scs/version"
 
 module SCS
   class Error < StandardError; end

data/vendor/scs/CITATION.cff

@@ -10,7 +10,7 @@ authors:
 - family-names: "Boyd"
   given-names: "Stephen"
 title: "SCS: Spltting Conic Solver"
-version: 3.2.3
+version: 3.2.6
 date-released: 2023
 url: "https://github.com/cvxgrp/scs"
 

data/vendor/scs/CMakeLists.txt

@@ -6,7 +6,7 @@ cmake_minimum_required(VERSION 3.5)
 project(
   scs
   LANGUAGES C
-  VERSION 3.2.3)
+  VERSION 3.2.6)
 
 # Defines the CMAKE_INSTALL_LIBDIR, CMAKE_INSTALL_BINDIR and many other useful
 # macros. See https://cmake.org/cmake/help/latest/module/GNUInstallDirs.html
@@ -121,7 +121,15 @@ message(STATUS "Printing is NOT ${NO_PRINTING}")
 option(NO_READ_WRITE "Disables all read/write functionality" OFF)
 message(STATUS "Read/write functionality is NOT ${NO_READ_WRITE}")
 
-set(COMPILER_OPTS "-DUSE_LAPACK -DCTRLC")
+option(USE_LAPACK "Whether to use BLAS/LAPACK" ON)
+message(STATUS "BLAS/LAPACK usage is ${USE_LAPACK}")
+
+# Enable OpenMP support
+option(USE_OPENMP "Compile with OpenMP support" OFF)
+message(STATUS "OpenMP parallelization is ${USE_OPENMP}")
+
+set(COMPILER_OPTS "-DCTRLC")
+set(LAPACK_LINK_LIBRARIES "")
 
 # Primitive types
 if(SFLOAT)
@@ -146,7 +154,17 @@ if(NO_READ_WRITE)
   set(COMPILER_OPTS "-DNO_READ_WRITE=1 ${COMPILER_OPTS}")
 endif()
 
+if (USE_LAPACK)
+  set(COMPILER_OPTS "-DUSE_LAPACK ${COMPILER_OPTS}")
+  list(APPEND LAPACK_LINK_LIBRARIES "blas" "lapack")
+endif()
+
+if(USE_OPENMP)
+  set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS} -fopenmp")
+endif()
+
 message(STATUS "COMPILER_OPTS = ${COMPILER_OPTS}")
+message(STATUS "CMAKE_C_FLAGS = ${CMAKE_C_FLAGS}")
 
 # TODO this is a hack that overwrites the scs_types.h file, we should find a way
 # to do this that doesn't pollute the master directory.
@@ -226,8 +244,8 @@ target_include_directories(
 # Compiled with blas and lapack, can solve LPs, SOCPs, SDPs, ECPs, and PCPs
 target_compile_definitions(${${PROJECT_NAME}_DIRECT} PRIVATE ${COMPILER_OPTS})
 
-# The library depends on math (m) blas and lapack
-target_link_libraries(${${PROJECT_NAME}_DIRECT} PRIVATE m blas lapack)
+# The library depends on math (m) and (optionally) blas and lapack
+target_link_libraries(${${PROJECT_NAME}_DIRECT} PRIVATE m ${LAPACK_LINK_LIBRARIES})
 
 # Set some properties
 set_target_properties(
@@ -274,8 +292,8 @@ target_include_directories(
 target_compile_definitions(${${PROJECT_NAME}_INDIRECT} PRIVATE ${COMPILER_OPTS}
                                                                -DINDIRECT)
 
-# The library depends on math (m) blas and lapack
-target_link_libraries(${${PROJECT_NAME}_INDIRECT} PUBLIC m blas lapack)
+# The library depends on math (m) and (optionally) blas and lapack
+target_link_libraries(${${PROJECT_NAME}_INDIRECT} PUBLIC m ${LAPACK_LINK_LIBRARIES})
 
 # Set some properties
 set_target_properties(

data/vendor/scs/Makefile

@@ -151,7 +151,7 @@ $(GPUINDIR)/private.o: $(GPUINDIR)/private.c
 # 	$(ARCHIVE) $@ $^
 # 	- $(RANLIB) $@
 
-$(OUT)/libscsgpuindir.$(SHARED): $(SCS_O) $(SCS_OBJECTS) $(GPUINDIR)/private.o $(LINSYS)/scs_matrix.o $(LINSYS)/csparse.o $(LINSYS)/gpu/gpu.o
+$(OUT)/libscsgpuindir.$(SHARED): $(SCS_INDIR_O) $(SCS_OBJECTS) $(GPUINDIR)/private.o $(LINSYS)/scs_matrix.o $(LINSYS)/csparse.o $(LINSYS)/gpu/gpu.o
 	mkdir -p $(OUT)
 	$(CC) $(CFLAGS) -shared -Wl,$(SONAME),$(@:$(OUT)/%=%) -o $@ $^ $(LDFLAGS) $(BLASLDFLAGS) $(CULDFLAGS)
 

data/vendor/scs/README.md

@@ -8,7 +8,7 @@
 
 
 SCS (`splitting conic solver`) is a numerical optimization package for solving
-large-scale convex cone problems. The current version is `3.2.3`.
+large-scale convex cone problems. The current version is `3.2.6`.
 
 The full documentation is available [here](https://www.cvxgrp.org/scs/).
 

data/vendor/scs/include/glbopts.h

@@ -14,7 +14,7 @@ extern "C" {
 
 /* SCS VERSION NUMBER ----------------------------------------------     */
 /* string literals automatically null-terminated */
-#define SCS_VERSION ("3.2.3")
+#define SCS_VERSION ("3.2.6")
 
 /* verbosity level */
 #ifndef VERBOSITY
@@ -39,7 +39,8 @@ extern "C" {
 #define WRITE_DATA_FILENAME (0)
 #define LOG_CSV_FILENAME (0)
 #define TIME_LIMIT_SECS (0.)
-
+/* Tolerance to check negativity condition for infeasibility */
+#define INFEAS_NEGATIVITY_TOL (1e-9)
 /* redefine printfs as needed */
 #if NO_PRINTING > 0     /* Disable all printing */
 #define scs_printf(...) /* No-op */

data/vendor/scs/include/scs_blas.h

@@ -17,8 +17,10 @@ extern "C" {
 /* single or double precision */
 #ifndef SFLOAT
 #define BLAS(x) d##x
+#define BLASI(x) id##x
 #else
 #define BLAS(x) s##x
+#define BLASI(x) is##x
 #endif
 #else
 /* this extra indirection is needed for BLASSUFFIX to work correctly as a
@@ -28,8 +30,10 @@ extern "C" {
 /* single or double precision */
 #ifndef SFLOAT
 #define BLAS(x) stitch__(d, x, BLASSUFFIX)
+#define BLASI(x) stitch__(id, x, BLASSUFFIX)
 #else
 #define BLAS(x) stitch__(s, x, BLASSUFFIX)
+#define BLASI(x) stitch__(is, x, BLASSUFFIX)
 #endif
 #endif
 

data/vendor/scs/linsys/cpu/direct/private.c

@@ -1,6 +1,6 @@
 #include "private.h"
 
-const char *scs_get_lin_sys_method() {
+const char *scs_get_lin_sys_method(void) {
   return "sparse-direct-amd-qdldl";
 }
 

data/vendor/scs/linsys/cpu/indirect/private.c

@@ -1,6 +1,6 @@
 #include "private.h"
 
-const char *scs_get_lin_sys_method() {
+const char *scs_get_lin_sys_method(void) {
   return "sparse-indirect-scs";
 }
 

data/vendor/scs/scs.mk

@@ -143,7 +143,7 @@ MKLFLAGS = -L$(MKLROOT) -L$(MKLROOT)/lib -Wl,--no-as-needed -lmkl_rt -lmkl_gnu_t
 USE_OPENMP = 0
 ifneq ($(USE_OPENMP), 0)
   override CFLAGS += -fopenmp
-  LDFLAGS += -lgomp
+  LDFLAGS += -fopenmp
 endif
 
 ############ SDPS: BLAS + LAPACK ############

data/vendor/scs/src/cones.c

@@ -331,7 +331,7 @@ char *SCS(get_cone_header)(const ScsCone *k) {
 static scs_int set_up_sd_cone_work_space(ScsConeWork *c, const ScsCone *k) {
   scs_int i;
 #ifdef USE_LAPACK
-  blas_int n_max = 0;
+  blas_int n_max = 1;
   blas_int neg_one = -1;
   blas_int info = 0;
   scs_float wkopt = 0.0;
@@ -356,7 +356,7 @@ static scs_int set_up_sd_cone_work_space(ScsConeWork *c, const ScsCone *k) {
    &info);
 
   if (info != 0) {
-    scs_printf("FATAL: syev failure, info = %li\n", (long)info);
+    scs_printf("FATAL: syev workspace query failure, info = %li\n", (long)info);
     return -1;
   }
   c->lwork = (blas_int)(wkopt + 1); /* +1 for int casting safety */

data/vendor/scs/src/linalg.c

@@ -102,12 +102,16 @@ extern "C" {
 scs_float BLAS(nrm2)(blas_int *n, const scs_float *x, blas_int *incx);
 scs_float BLAS(dot)(const blas_int *n, const scs_float *x, const blas_int *incx,
                     const scs_float *y, const blas_int *incy);
-scs_float BLAS(lange)(const char *norm, const blas_int *m, const blas_int *n,
-                      const scs_float *a, blas_int *lda, scs_float *work);
 void BLAS(axpy)(blas_int *n, const scs_float *a, const scs_float *x,
                 blas_int *incx, scs_float *y, blas_int *incy);
 void BLAS(scal)(const blas_int *n, const scs_float *sa, scs_float *sx,
                 const blas_int *incx);
+blas_int BLASI(amax)(blas_int *n, const scs_float *x, blas_int *incx);
+
+/* Possibly not working correctly on all platforms.
+scs_float BLAS(lange)(const char *norm, const blas_int *m, const blas_int *n,
+                      const scs_float *a, blas_int *lda, scs_float *work);
+*/
 
 #ifdef __cplusplus
 }
@@ -140,10 +144,20 @@ scs_float SCS(norm_2)(const scs_float *v, scs_int len) {
   return BLAS(nrm2)(&blen, v, &bone);
 }
 
+/* Possibly not working correctly on all platforms.
+scs_float SCS(norm_inf)(const scs_float *a, scs_int len) {
+  blas_int bone = 1;
+  blas_int blen = (blas_int)len;
+  return BLAS(lange)("Max", &blen, &bone, a, &blen, SCS_NULL);
+}
+*/
+
 scs_float SCS(norm_inf)(const scs_float *a, scs_int len) {
   blas_int bone = 1;
   blas_int blen = (blas_int)len;
-  return BLAS(lange)("Max", &blen, &bone, a, &bone, SCS_NULL);
+  scs_int idx = (scs_int)BLASI(amax)(&blen, a, &bone);
+  /* Returned idx is 1-based. */
+  return ABS(a[idx - 1]);
 }
 
 /* axpy a += sc*b */

data/vendor/scs/src/rw.c

@@ -60,7 +60,7 @@ static size_t read_int(scs_int *dest, size_t file_int_sz, size_t nitems,
   }
   void *ptr = scs_calloc(nitems, file_int_sz);
   size_t val = fread(ptr, file_int_sz, nitems, fin);
-  scs_int i;
+  size_t i;
   switch (file_int_sz) {
   case 4:
     for (i = 0; i < nitems; ++i) {
@@ -328,8 +328,15 @@ void SCS(log_data_to_csv)(const ScsCone *k, const ScsSettings *stgs,
                   "tau_normalized,"
                   "kap_normalized,"
                   "ax_nrm_inf,"
+                  "ax_s_nrm_inf"
                   "px_nrm_inf,"
                   "aty_nrm_inf,"
+                  "xt_p_x,"
+                  "xt_p_x_tau,"
+                  "ctx,"
+                  "ctx_tau,"
+                  "bty,"
+                  "bty_tau,"
                   "b_nrm_inf,"
                   "c_nrm_inf,"
                   "scale,"
@@ -385,8 +392,15 @@ void SCS(log_data_to_csv)(const ScsCone *k, const ScsSettings *stgs,
   fprintf(fout, "%.16e,", r_n->tau);
   fprintf(fout, "%.16e,", r_n->kap);
   fprintf(fout, "%.16e,", SCS(norm_inf)(r->ax, w->d->m));
+  fprintf(fout, "%.16e,", SCS(norm_inf)(r->ax_s, w->d->m));
   fprintf(fout, "%.16e,", SCS(norm_inf)(r->px, w->d->n));
   fprintf(fout, "%.16e,", SCS(norm_inf)(r->aty, w->d->n));
+  fprintf(fout, "%.16e,", r->xt_p_x);
+  fprintf(fout, "%.16e,", r->xt_p_x_tau);
+  fprintf(fout, "%.16e,", r->ctx);
+  fprintf(fout, "%.16e,", r->ctx_tau);
+  fprintf(fout, "%.16e,", r->bty);
+  fprintf(fout, "%.16e,", r->bty_tau);
   fprintf(fout, "%.16e,", SCS(norm_inf)(w->b_orig, w->d->m));
   fprintf(fout, "%.16e,", SCS(norm_inf)(w->c_orig, w->d->n));
   fprintf(fout, "%.16e,", w->stgs->scale);

data/vendor/scs/src/scs.c

@@ -112,6 +112,9 @@ static void print_init_header(const ScsData *d, const ScsCone *k,
   if (stgs->time_limit_secs) {
     scs_printf("\t  time_limit_secs: %.2e\n", stgs->time_limit_secs);
   }
+#ifdef _OPENMP
+  scs_printf("\t  compiled with openmp parallelization enabled\n");
+#endif
   if (lin_sys_method) {
     scs_printf("lin-sys:  %s\n\t  nnz(A): %li, nnz(P): %li\n", lin_sys_method,
                (long)d->A->p[d->A->n], d->P ? (long)d->P->p[d->P->n] : 0l);
@@ -204,13 +207,13 @@ static void compute_residuals(ScsResiduals *r, scs_int m, scs_int n) {
   r->res_unbdd_a = NAN;
   r->res_unbdd_p = NAN;
   r->res_infeas = NAN;
-  if (r->ctx_tau < 0) {
+  if (r->ctx_tau < -INFEAS_NEGATIVITY_TOL) {
     nm_ax_s = NORM(r->ax_s, m);
     nm_px = NORM(r->px, n);
     r->res_unbdd_a = SAFEDIV_POS(nm_ax_s, -r->ctx_tau);
     r->res_unbdd_p = SAFEDIV_POS(nm_px, -r->ctx_tau);
   }
-  if (r->bty_tau < 0) {
+  if (r->bty_tau < -INFEAS_NEGATIVITY_TOL) {
     nm_aty = NORM(r->aty, n);
     r->res_infeas = SAFEDIV_POS(nm_aty, -r->bty_tau);
   }
@@ -584,22 +587,22 @@ static void print_summary(ScsWork *w, scs_int i, SCS(timer) * solve_timer) {
   scs_printf("\n");
 
 #if VERBOSITY > 0
-  scs_printf("Norm u = %4f, ", SCS(norm_2)(w->u, w->d->n + w->d->m + 1));
-  scs_printf("Norm u_t = %4f, ", SCS(norm_2)(w->u_t, w->d->n + w->d->m + 1));
-  scs_printf("Norm v = %4f, ", SCS(norm_2)(w->v, w->d->n + w->d->m + 1));
-  scs_printf("Norm rsk = %4f, ", SCS(norm_2)(w->rsk, w->d->n + w->d->m + 1));
-  scs_printf("Norm x = %4f, ", SCS(norm_2)(w->xys_orig->x, w->d->n));
-  scs_printf("Norm y = %4f, ", SCS(norm_2)(w->xys_orig->y, w->d->m));
-  scs_printf("Norm s = %4f, ", SCS(norm_2)(w->xys_orig->s, w->d->m));
-  scs_printf("Norm |Ax + s| = %1.2e, ", SCS(norm_2)(r->ax_s, w->d->m));
-  scs_printf("tau = %4f, ", w->u[w->d->n + w->d->m]);
-  scs_printf("kappa = %4f, ", w->rsk[w->d->n + w->d->m]);
-  scs_printf("|u - u_t| = %1.2e, ",
+  scs_printf("Norm u = %1.6e, ", SCS(norm_2)(w->u, w->d->n + w->d->m + 1));
+  scs_printf("Norm u_t = %1.6e, ", SCS(norm_2)(w->u_t, w->d->n + w->d->m + 1));
+  scs_printf("Norm v = %1.6e, ", SCS(norm_2)(w->v, w->d->n + w->d->m + 1));
+  scs_printf("Norm rsk = %1.6e, ", SCS(norm_2)(w->rsk, w->d->n + w->d->m + 1));
+  scs_printf("Norm x = %1.6e, ", SCS(norm_2)(w->xys_orig->x, w->d->n));
+  scs_printf("Norm y = %1.6e, ", SCS(norm_2)(w->xys_orig->y, w->d->m));
+  scs_printf("Norm s = %1.6e, ", SCS(norm_2)(w->xys_orig->s, w->d->m));
+  scs_printf("Norm |Ax + s| = %1.6e, ", SCS(norm_2)(r->ax_s, w->d->m));
+  scs_printf("tau = %1.6e, ", w->u[w->d->n + w->d->m]);
+  scs_printf("kappa = %1.6e, ", w->rsk[w->d->n + w->d->m]);
+  scs_printf("|u - u_t| = %1.6e, ",
              SCS(norm_diff)(w->u, w->u_t, w->d->n + w->d->m + 1));
-  scs_printf("res_infeas = %1.2e, ", r->res_infeas);
-  scs_printf("res_unbdd_a = %1.2e, ", r->res_unbdd_a);
-  scs_printf("res_unbdd_p = %1.2e, ", r->res_unbdd_p);
-  scs_printf("ctx_tau = %1.2e, ", r->ctx_tau);
+  scs_printf("res_infeas = %1.6e, ", r->res_infeas);
+  scs_printf("res_unbdd_a = %1.6e, ", r->res_unbdd_a);
+  scs_printf("res_unbdd_p = %1.6e, ", r->res_unbdd_p);
+  scs_printf("ctx_tau = %1.6e, ", r->ctx_tau);
   scs_printf("bty_tau = %1.2e\n", r->bty_tau);
 #endif
 

data/vendor/scs/test/problems/mpc_bug.h

@@ -0,0 +1,19 @@
+#include "glbopts.h"
+#include "problems/test_prob_from_data_file.h"
+#include "scs.h"
+
+static const char *mpc_bug(void) {
+  const char *fail;
+  scs_float OPT1 = -0.473957794500; /* from scs */
+  scs_float OPT2 = -0.029336830816; /* from scs */
+  scs_float OPT3 = -0.002215217478; /* from scs */
+  fail = _test_prob_from_data("test/problems/mpc_bug1", OPT1);
+  if (fail) {
+    return fail;
+  }
+  fail = _test_prob_from_data("test/problems/mpc_bug2", OPT2);
+  if (fail) {
+    return fail;
+  }
+  return _test_prob_from_data("test/problems/mpc_bug3", OPT3);
+}

data/vendor/scs/test/problems/test_prob_from_data_file.h

@@ -19,6 +19,8 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
   scs_float perr, derr;
   scs_int success;
   const char *fail;
+  scs_float xt_p_x;
+  scs_float *px = SCS_NULL;
 
   read_status = SCS(read_data)(file, &d, &k, &stgs);
 
@@ -28,12 +30,24 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
 
   stgs->eps_abs = 1e-6;
   stgs->eps_rel = 1e-6;
+  /* Force verbosity for the test */
+  stgs->verbose = 1;
 
   sol = (ScsSolution *)scs_calloc(1, sizeof(ScsSolution));
   exitflag = scs(d, k, stgs, sol, &info);
 
-  perr = SCS(dot)(d->c, sol->x, d->n) - OPT;
-  derr = -SCS(dot)(d->b, sol->y, d->m) - OPT;
+  if (d->P) {
+    /* px = Px */
+    px = (scs_float *)scs_calloc(d->n, sizeof(scs_float));
+    memset(px, 0, d->n * sizeof(scs_float));
+    SCS(accum_by_p)(d->P, sol->x, px);
+    xt_p_x = SCS(dot)(px, sol->x, d->n);
+  } else {
+    xt_p_x = 0.;
+  }
+
+  perr = 0.5 * xt_p_x + SCS(dot)(d->c, sol->x, d->n) - OPT;
+  derr = -0.5 * xt_p_x - SCS(dot)(d->b, sol->y, d->m) - OPT;
   scs_printf("primal obj error %4e\n", perr);
   scs_printf("dual obj error %4e\n", derr);
 
@@ -47,7 +61,9 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
   SCS(free_cone)(k);
   SCS(free_sol)(sol);
   scs_free(stgs);
-
+  if (px) {
+    scs_free(px);
+  }
   if (fail) {
     scs_printf("%s: FAILED\n", file);
   }

data/vendor/scs/test/run_tests.c

@@ -43,9 +43,11 @@ _SKIP(random_prob)
 #if NO_READ_WRITE == 0 /* reads / writes */
 #include "problems/hs21_tiny_qp_rw.h"
 #include "problems/max_ent.h"
+#include "problems/mpc_bug.h"
 #else
 _SKIP(hs21_tiny_qp_rw)
 _SKIP(max_ent)
+_SKIP(mpc_bug)
 #endif
 
 static const char *all_tests(void) {
@@ -61,6 +63,7 @@ static const char *all_tests(void) {
   mu_run_test(unbounded_tiny_qp);
   mu_run_test(random_prob);
   mu_run_test(max_ent);
+  mu_run_test(mpc_bug);
   mu_run_test(test_exp_cone);
   return 0;
 }

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: scs
 version: !ruby/object:Gem::Version
-  version: 0.4.1
+  version: 0.4.2
 platform: ruby
 authors:
 - Andrew Kane
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2023-04-09 00:00:00.000000000 Z
+date: 2024-07-11 00:00:00.000000000 Z
 dependencies: []
 description:
 email: andrew@ankane.org
@@ -102,6 +102,10 @@ files:
 - vendor/scs/test/problems/infeasible_tiny_qp.h
 - vendor/scs/test/problems/max_ent
 - vendor/scs/test/problems/max_ent.h
+- vendor/scs/test/problems/mpc_bug.h
+- vendor/scs/test/problems/mpc_bug1
+- vendor/scs/test/problems/mpc_bug2
+- vendor/scs/test/problems/mpc_bug3
 - vendor/scs/test/problems/qafiro_tiny_qp.h
 - vendor/scs/test/problems/random_prob
 - vendor/scs/test/problems/random_prob.h
@@ -135,7 +139,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.4.10
+rubygems_version: 3.5.11
 signing_key:
 specification_version: 4
 summary: SCS - the splitting conic solver - for Ruby