vendor/scs/src/scs.c in scs-0.4.1 vs vendor/scs/src/scs.c in scs-0.4.2
- old
+ new
@@ -110,10 +110,13 @@
(int)acceleration_lookback, (int)acceleration_interval);
}
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);
}
@@ -202,17 +205,17 @@
r->res_pri = SAFEDIV_POS(nm_ax_s_btau, r->tau);
r->res_dual = SAFEDIV_POS(nm_px_aty_ctau, r->tau);
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);
}
}
@@ -582,25 +585,25 @@
scs_printf("%*.2e ", (int)HSPACE,
(SCS(tocq)(solve_timer) + w->setup_time) / 1e3);
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
#ifdef MATLAB_MEX_FILE
mexEvalString("drawnow;");