scs 0.5.1 → 0.5.2

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

@@ -0,0 +1,83 @@
+#include "glbopts.h"
+#include "linalg.h"
+#include "minunit.h"
+#include "problem_utils.h"
+#include "rw.h"
+#include "scs.h"
+#include "scs_matrix.h"
+#include "util.h"
+
+static const char *complex_PSD(void)
+{
+    ScsCone *k = (ScsCone *)scs_calloc(1, sizeof(ScsCone));
+    ScsData *d = (ScsData *)scs_calloc(1, sizeof(ScsData));
+    ScsSettings *stgs = (ScsSettings *)scs_calloc(1, sizeof(ScsSettings));
+    ScsSolution *sol = (ScsSolution *)scs_calloc(1, sizeof(ScsSolution));
+    ScsInfo info = {0};
+    scs_int exitflag;
+    scs_float perr, derr;
+    scs_int success;
+    const char *fail;
+
+    /* data */
+    scs_float Ax[] = {1.0, -1.0, -1.0, -1.0, -1.0, -1.0, 1.0, -1.0, -1.0, -1.0, 1.0, -1.0};
+    scs_int Ai[] = {0, 1, 2, 3, 4, 5, 0, 6, 7, 8, 0, 9};
+    scs_int Ap[] = {0, 2, 3, 4, 5, 6, 8, 9, 10, 12};
+    scs_float b[] = {1.0, 0., 0., 0., 0., 0., 0., 0., 0., 0.};
+    scs_float c[] = {1, 2 * sqrt(2), 3 * sqrt(2), 4 * sqrt(2), 5 * sqrt(2), 6,
+                     7 * sqrt(2), -8 * sqrt(2), 9};
+
+    scs_int m = 10;
+    scs_int n = 9;
+    scs_int cs[] = {3};
+    scs_int cssize = 1;
+    k->z = 1;
+    k->cs = cs;
+    k->cssize = cssize;
+
+    // computed offline
+    scs_float opt = -5.228930;
+    /* end data */
+
+    d->m = m;
+    d->n = n;
+    d->b = b;
+    d->c = c;
+
+    d->A = (ScsMatrix *)scs_calloc(1, sizeof(ScsMatrix));
+    d->A->m = m;
+    d->A->n = n;
+    d->A->x = Ax;
+    d->A->i = Ai;
+    d->A->p = Ap;
+
+    scs_set_default_settings(stgs);
+    stgs->eps_abs = 1e-7;
+    stgs->eps_rel = 1e-7;
+    stgs->eps_infeas = 1e-9;
+
+    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;
+
+    scs_printf("primal obj error %4e\n", perr);
+    scs_printf("dual obj error %4e\n", derr);
+
+    success = ABS(perr) < 1e-4 && ABS(derr) < 1e-4 && exitflag == SCS_SOLVED;
+
+    mu_assert("complex_PSD: SCS failed to produce outputflag SCS_SOLVED", success);
+
+    fail = verify_solution_correct(d, k, stgs, &info, sol, exitflag);
+    if (fail)
+        return fail;
+
+    /* kill data */
+    scs_free(d->A);
+    scs_free(k);
+    scs_free(stgs);
+    scs_free(d);
+    SCS(free_sol)(sol);
+
+    return fail;
+}

data/vendor/scs/test/rng.h

@@ -52,14 +52,14 @@ long ran_arr_buf[QUALITY];
 long ran_arr_dummy = -1, ran_arr_started = -1;
 long *ran_arr_ptr = &ran_arr_dummy; /* the next random number, or -1 */
 
-#define TT 70             /* guaranteed separation between streams */
-#define is_odd(x) ((x)&1) /* units bit of x */
+#define TT 70               /* guaranteed separation between streams */
+#define is_odd(x) ((x) & 1) /* units bit of x */
 
 #ifdef __STDC__
 void ran_start(long seed)
 #else
-void ran_start(seed)  /* do this before using ran_array */
-    long seed;        /* selector for different streams */
+void ran_start(seed) /* do this before using ran_array */
+    long seed;       /* selector for different streams */
 #endif
 {
   register int t, j;

data/vendor/scs/test/run_tests.c

@@ -33,11 +33,29 @@ _SKIP(test_validation)
 
 /* solve SDPs, requires blas / lapack */
 #if defined(USE_LAPACK) && NO_READ_WRITE == 0
+#include "problems/complex_PSD.h"
 #include "problems/random_prob.h"
 #include "problems/rob_gauss_cov_est.h"
 #else
 _SKIP(rob_gauss_cov_est)
 _SKIP(random_prob)
+_SKIP(complex_PSD)
+#endif
+
+#if defined(USE_SPECTRAL_CONES) && NO_READ_WRITE == 0
+#include "spectral_cones_problems/exp_design.h"
+#include "spectral_cones_problems/graph_partitioning.h"
+#include "spectral_cones_problems/robust_pca.h"
+#include "spectral_cones_problems/several_logdet_cones.h"
+#include "spectral_cones_problems/several_nuc_cone.h"
+#include "spectral_cones_problems/several_sum_largest.h"
+#else
+_SKIP(exp_design)
+_SKIP(robust_pca)
+_SKIP(graph_partitioning)
+_SKIP(several_sum_largest)
+_SKIP(several_nuc_cone)
+_SKIP(several_logdet_cones)
 #endif
 
 #if NO_READ_WRITE == 0 /* reads / writes */
@@ -56,6 +74,7 @@ static const char *all_tests(void) {
   mu_run_test(small_lp);
   mu_run_test(small_qp);
   mu_run_test(rob_gauss_cov_est);
+  mu_run_test(complex_PSD);
   mu_run_test(hs21_tiny_qp);
   mu_run_test(hs21_tiny_qp_rw);
   mu_run_test(qafiro_tiny_qp);
@@ -65,6 +84,12 @@ static const char *all_tests(void) {
   mu_run_test(max_ent);
   mu_run_test(mpc_bug);
   mu_run_test(test_exp_cone);
+  mu_run_test(exp_design);
+  mu_run_test(robust_pca);
+  mu_run_test(graph_partitioning);
+  mu_run_test(several_sum_largest);
+  mu_run_test(several_nuc_cone);
+  mu_run_test(several_logdet_cones);
   return 0;
 }
 int main(void) {

data/vendor/scs/test/spectral_cones_problems/exp_design.h

@@ -0,0 +1,141 @@
+#include "glbopts.h"
+#include "linalg.h"
+#include "minunit.h"
+#include "problem_utils.h"
+#include "rw.h"
+#include "scs.h"
+#include "scs_matrix.h"
+#include "util.h"
+
+// for SpectralSCS
+static const char *exp_design(void) {
+  ScsCone *k = (ScsCone *)scs_calloc(1, sizeof(ScsCone));
+  ScsData *d = (ScsData *)scs_calloc(1, sizeof(ScsData));
+  ScsSettings *stgs = (ScsSettings *)scs_calloc(1, sizeof(ScsSettings));
+  ScsSolution *sol = (ScsSolution *)scs_calloc(1, sizeof(ScsSolution));
+  ScsInfo info = {0};
+  scs_int exitflag;
+  scs_float perr, derr;
+  scs_int success;
+  const char *fail;
+
+  /* data */
+  scs_float Ax[] = {
+      -1.,         1.,          -1.,         3.24,        0.16,
+      1.,          4.84,        3.61,        1.,          -1.,
+      2.54558441,  -0.11313708, -0.14142136, 1.24450793,  0.26870058,
+      -2.12132034, -1.,         2.03646753,  0.05656854,  0.56568542,
+      0.93338095,  4.03050865,  0.28284271,  -1.,         1.,
+      0.04,        0.01,        0.16,        0.01,        2.25,
+      -1.,         1.13137085,  -0.02828427, -0.05656854, 0.16970563,
+      0.21213203,  -0.42426407, -1.,         0.64,        0.01,
+      0.16,        0.09,        2.25,        0.04,        -1.};
+  scs_int Ai[] = {7,  0,  8, 1, 2, 3, 4, 5,  6,  9, 1, 2, 3, 4, 5,
+                  6,  10, 1, 2, 3, 4, 5, 6,  11, 1, 2, 3, 4, 5, 6,
+                  12, 1,  2, 3, 4, 5, 6, 13, 1,  2, 3, 4, 5, 6, 14};
+  scs_int Ap[] = {0, 1, 3, 10, 17, 24, 31, 38, 45};
+
+  scs_float b[] = {1., 1., 1., 1., 1., 1., 1., 0., 0., 0., 0., 0., 0., 0., 0.};
+  scs_float c[] = {1., 0., 0., 0., 0., 0., 0., 0.};
+
+  scs_int m = 15;
+  scs_int n = 8;
+
+  scs_int z = 1;
+  scs_int l = 6;
+  scs_int *q = SCS_NULL;
+  scs_int qsize = 0;
+  scs_int *s = SCS_NULL;
+  scs_int ssize = 0;
+  scs_int ep = 0;
+  scs_int ed = 0;
+  scs_float *p = SCS_NULL;
+  scs_int psize = 0;
+  scs_int d_array[] = {3};
+  scs_int dsize = 1;
+  scs_int *nuc_m = SCS_NULL;
+  scs_int *nuc_n = SCS_NULL;
+  scs_int nucsize = 0;
+  scs_int *ell1 = SCS_NULL;
+  scs_int ell1_size = 0;
+  scs_int *sl_n = SCS_NULL;
+  scs_int *sl_k = SCS_NULL;
+  scs_int sl_size = 0;
+
+  // computed using mosek (the input of Ax is truncated, and mosek solved
+  // the problem with the non-truncated data)
+  scs_float opt = 3.0333290743428574;
+  /* end data */
+
+  d->m = m;
+  d->n = n;
+  d->b = b;
+  d->c = c;
+
+  d->A = (ScsMatrix *)scs_calloc(1, sizeof(ScsMatrix));
+
+  d->A->m = m;
+  d->A->n = n;
+  d->A->x = Ax;
+  d->A->i = Ai;
+  d->A->p = Ap;
+
+  k->z = z;
+  k->l = l;
+  k->q = q;
+  k->qsize = qsize;
+  k->s = s;
+  k->ssize = ssize;
+  k->ep = ep;
+  k->ed = ed;
+  k->p = p;
+  k->psize = psize;
+  k->d = d_array;
+  k->dsize = dsize;
+  k->nuc_m = nuc_m;
+  k->nuc_n = nuc_n;
+  k->nucsize = nucsize;
+  k->ell1 = ell1;
+  k->ell1_size = ell1_size;
+  k->sl_n = sl_n;
+  k->sl_k = sl_k;
+  k->sl_size = sl_size;
+
+  scs_set_default_settings(stgs);
+  stgs->eps_abs = 1e-7;
+  stgs->eps_rel = 1e-7;
+  stgs->eps_infeas = 1e-9;
+
+  stgs->log_csv_filename = "test_exp_design.csv";
+  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;
+
+  scs_printf("primal obj error %4e\n", perr);
+  scs_printf("dual obj error %4e\n", derr);
+
+  success = ABS(perr) < 1e-4 && ABS(derr) < 1e-4 && exitflag == SCS_SOLVED;
+
+  mu_assert("exp_design: SCS failed to produce outputflag SCS_SOLVED", success);
+
+  fail = 0;
+  // TODO: This test fails because of the complementary slackness check.
+  //       The complementary slackness tolerance is a bit too tight.
+  // fail = verify_solution_correct(d, k, stgs, &info, sol, exitflag);
+  // if (fail)
+  //   return fail;
+
+  mu_assert("exp_design: primal feas error: ", ABS(info.res_pri) < 1e-6);
+  mu_assert("exp_design: dual feas error: ", ABS(info.res_dual) < 1e-6);
+  mu_assert("exp_design: duality gap error: ", ABS(info.gap) < 1e-6);
+
+  /* kill data */
+  scs_free(d->A);
+  scs_free(k);
+  scs_free(stgs);
+  scs_free(d);
+  SCS(free_sol)(sol);
+
+  return fail;
+}

data/vendor/scs/test/spectral_cones_problems/graph_partitioning.h

@@ -0,0 +1,275 @@
+#include "glbopts.h"
+#include "linalg.h"
+#include "minunit.h"
+#include "problem_utils.h"
+#include "rw.h"
+#include "scs.h"
+#include "scs_matrix.h"
+#include "util.h"
+
+// for SpectralSCS
+static const char *graph_partitioning(void) {
+  ScsCone *k = (ScsCone *)scs_calloc(1, sizeof(ScsCone));
+  ScsData *d = (ScsData *)scs_calloc(1, sizeof(ScsData));
+  ScsSettings *stgs = (ScsSettings *)scs_calloc(1, sizeof(ScsSettings));
+  ScsSolution *sol = (ScsSolution *)scs_calloc(1, sizeof(ScsSolution));
+  ScsInfo info = {0};
+  scs_int exitflag;
+  scs_float perr, derr;
+  scs_int success;
+  const char *fail;
+
+  /* data */
+  scs_float Ax[] = {-1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1., 1., -1., 1.,
+                    -1., 1., -1., 1., -1., 1., -1., 1., -1.};
+  scs_int Ai[] = {1,   0, 2,   0, 42,  0, 81,  0, 119, 0, 156, 0, 192, 0,
+                  227, 0, 261, 0, 294, 0, 326, 0, 357, 0, 387, 0, 416, 0,
+                  444, 0, 471, 0, 497, 0, 522, 0, 546, 0, 569, 0, 591, 0,
+                  612, 0, 632, 0, 651, 0, 669, 0, 686, 0, 702, 0, 717, 0,
+                  731, 0, 744, 0, 756, 0, 767, 0, 777, 0, 786, 0, 794, 0,
+                  801, 0, 807, 0, 812, 0, 816, 0, 819, 0, 821};
+  scs_int Ap[] = {0,  1,  3,  5,  7,  9,  11, 13, 15, 17, 19, 21, 23, 25,
+                  27, 29, 31, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53,
+                  55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77, 79, 81};
+
+  scs_float b[] = {
+      0.,         0.,         -1.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        1.41421356, -0.,        -0.,        -0.,
+      -1.,        -0.,        1.41421356, -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -1.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -6.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        1.41421356, 1.41421356,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      1.41421356, -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -2.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        1.41421356, -0.,        -0.,        1.41421356,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -2.,        -0.,        1.41421356, -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -3.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      1.41421356, -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        1.41421356,
+      1.41421356, -0.,        -0.,        -6.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      1.41421356, -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        1.41421356, -0.,
+      -1.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -1.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -3.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      1.41421356, -0.,        1.41421356, -2.,        1.41421356, -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -3.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -4.,        1.41421356, -0.,        -0.,        1.41421356, -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -3.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -5.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        1.41421356, -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        1.41421356, -0.,        1.41421356, -0.,
+      -2.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        1.41421356,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -3.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        1.41421356, -0.,
+      1.41421356, -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -2.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        1.41421356,
+      -0.,        -0.,        1.41421356, -0.,        -0.,        -0.,
+      -2.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -1.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -2.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -2.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        1.41421356, -0.,        -0.,        1.41421356, -0.,
+      -0.,        -0.,        -3.,        -0.,        -0.,        -0.,
+      1.41421356, 1.41421356, -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -1.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -5.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        1.41421356, -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -2.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -4.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        1.41421356, -0.,        1.41421356, -0.,
+      -4.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -2.,
+      -0.,        -0.,        -0.,        -0.,        1.41421356, 1.41421356,
+      -0.,        -0.,        -0.,        -1.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -0.,        -0.,        -0.,
+      -3.,        -0.,        1.41421356, -0.,        -0.,        -0.,
+      -0.,        1.41421356, -4.,        -0.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -1.,        -0.,        -0.,
+      -0.,        -0.,        -0.,        -4.,        -0.,        -0.,
+      -0.,        1.41421356, -7.,        -0.,        -0.,        -0.,
+      -2.,        -0.,        -0.,        -4.,        -0.,        -3.};
+  scs_float c[] = {1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+                   0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+                   0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.};
+
+  scs_int m = 822;
+  scs_int n = 41;
+
+  scs_int z = 1;
+  scs_int l = 0;
+  scs_int *q = SCS_NULL;
+  scs_int qsize = 0;
+  scs_int *s = SCS_NULL;
+  scs_int ssize = 0;
+  scs_int ep = 0;
+  scs_int ed = 0;
+  scs_float *p = SCS_NULL;
+  scs_int psize = 0;
+  scs_int *d_array = SCS_NULL;
+  scs_int dsize = 0;
+  scs_int *nuc_m = SCS_NULL;
+  scs_int *nuc_n = SCS_NULL;
+  scs_int nucsize = 0;
+  scs_int *ell1 = SCS_NULL;
+  scs_int ell1_size = 0;
+  scs_int sl_n[] = {40};
+  scs_int sl_k[] = {3};
+  scs_int sl_size = 1;
+
+  // computed using mosek (the input of Ax is truncated, and mosek solved
+  // the problem with the non-truncated data)
+  scs_float opt = -0.7736762265822145;
+  /* end data */
+
+  d->m = m;
+  d->n = n;
+  d->b = b;
+  d->c = c;
+
+  d->A = (ScsMatrix *)scs_calloc(1, sizeof(ScsMatrix));
+
+  d->A->m = m;
+  d->A->n = n;
+  d->A->x = Ax;
+  d->A->i = Ai;
+  d->A->p = Ap;
+
+  k->z = z;
+  k->l = l;
+  k->q = q;
+  k->qsize = qsize;
+  k->s = s;
+  k->ssize = ssize;
+  k->ep = ep;
+  k->ed = ed;
+  k->p = p;
+  k->psize = psize;
+  k->d = d_array;
+  k->dsize = dsize;
+  k->nuc_m = nuc_m;
+  k->nuc_n = nuc_n;
+  k->nucsize = nucsize;
+  k->ell1 = ell1;
+  k->ell1_size = ell1_size;
+  k->sl_n = sl_n;
+  k->sl_k = sl_k;
+  k->sl_size = sl_size;
+
+  scs_set_default_settings(stgs);
+  stgs->eps_abs = 1e-7;
+  stgs->eps_rel = 1e-7;
+  stgs->eps_infeas = 1e-9;
+
+  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;
+
+  scs_printf("primal obj error %4e\n", perr);
+  scs_printf("dual obj error %4e\n", derr);
+
+  success = ABS(perr) < 1e-4 && ABS(derr) < 1e-4 && exitflag == SCS_SOLVED;
+
+  mu_assert("graph partitioning: SCS failed to produce outputflag SCS_SOLVED",
+            success);
+
+  fail = verify_solution_correct(d, k, stgs, &info, sol, exitflag);
+  if (fail)
+    return fail;
+
+  /* kill data */
+  scs_free(d->A);
+  scs_free(k);
+  scs_free(stgs);
+  scs_free(d);
+  SCS(free_sol)(sol);
+
+  return fail;
+}