X - UWSpace - University of Waterloo

* Use advanced basis at each iteration */
status = CPXsetintparam(env, CPX-PARAM-ADVIND, 1);
if ( status 1 {
fprintf (stderr, "Failed to scale.\nm);
goto TERMINATE;
/* For barrier method, use option 1 to solve subproblems */
status = CPXsetintparam (env, CPX-PARAM-BARALG, Il ;
if ( scatus {
fprintf (scderr, "Failed to scale. \nm
goto TERMINATE;
1
/* Create the problern- */
Ip = CPXcreateprob (env, &scatus, "3-regionmi ;
1
if ( lp == NULL i I
1 ;
fprintf (stderr, 'Failed co create LP. \na
goto TERMINATE;
/* Now cûpy the probiem data into CPLEX */
status=CPXcopylp (env,Ip, ( inc) LpSub [cur-block] ->n, ( int 1 LpSub [CU-blockl ->m,
-l,objcoeffs, LpSub[cur-blockl->rhs,sençe,(int*) (LpSub[=ur-blockl-clpnts),
maccnt, (inc*) LpSublcur~blockl-~rInnmbs,LpSub~cur~block]-~coeffs,
LpSub[cur,blockj->la-bndç. LpSub[cur,block]->up-bnds, NaL) ;
if ( status ) E
fprintf (stderr, "Failed co copy problem data.\nm);
goto TERMINATE;
1
/* Delete the first row because the original vector has always the objective
function coefficients in row zero '/
status = CPXdelrows (env, Ip, O, 0);
if ( status i I
I
fprintf (stderr, "Failed CO delete row O\nW) ;
goto TERMINATE;
Add an artificial variable to each linking constraints */
for (k=L; k 0)
status = CPX.newcols(env, lp, (înt) LSPbmlkl, LpSubfkl->arti,obj, NWLL.
NUU, NULL,
NtJtLl:
i
for ( F=O; im-LOWSPb-m[cur,blockW; im-l;i++)
stacus=CPXchgccef (ew, lp, (int) i, (=nt) (LpSub[cu,Slock] ->n) +
LSPb~mI~~5lockl+i-LpS~~~~r~'olock~-~~+LOWSPb~m[~~block]l
, arti-coef);
1 ;