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;");