/* Remove FILEs, honoring options specified via X.
Return RM_OK if successful. */
enum RM_status
rm (char *const *file, struct rm_options const *x)
{
enum RM_status rm_status = RM_OK;
if (*file)
{
int bit_flags = (FTS_CWDFD
| FTS_NOSTAT
| FTS_PHYSICAL);
if (x->one_file_system)
bit_flags |= FTS_XDEV;
FTS *fts = xfts_open (file, bit_flags, NULL);
while (1)
{
FTSENT *ent;
ent = fts_read (fts);
if (ent == NULL)
{
if (errno != 0)
{
error (0, errno, _("fts_read failed"));
rm_status = RM_ERROR;
}
break;
}
enum RM_status s = rm_fts (fts, ent, x);
assert (VALID_STATUS (s));
UPDATE_STATUS (rm_status, s);
}
if (fts_close (fts) != 0)
{
error (0, errno, _("fts_close failed"));
rm_status = RM_ERROR;
}
}
return rm_status;
}
static bool process_inactive (lpprob_struct *orig_lp, int oxkndx)
/*
This routine handles the data structure updates for an inactive variable
x<k>. We need to have a look at the bounds l<k> and u<k>, and perhaps
update the status kept in dy_origvars. We need to add the contribution
c<k>l<k> or c<k>u<k> to the objective function. Finally, if we've reloaded
b & blow due to a bound or rhs change, we need to walk the column a<k>
and adjust b<i> (and perhaps blow<i>) for each nonzero a<ik> in the active
system.
Parameters:
orig_lp: the original lp problem
oxkndx: index of x<k> in orig_sys
Returns: TRUE if the update is made without incident, FALSE otherwise.
*/
{ int oaindx,aindx,ndx ;
double xk,lk,uk,ck ;
pkvec_struct *ak ;
pkcoeff_struct *aik ;
consys_struct *orig_sys ;
flags xkstatus ;
const char *rtnnme = "process_inactive" ;
orig_sys = orig_lp->consys ;
xkstatus = getflg(orig_lp->status[oxkndx],vstatSTATUS) ;
# ifdef DYLP_PARANOIA
/*
Any inactive variable should be nonbasic, and the paranoid check is looking
to make sure of this.
*/
if (!VALID_STATUS(xkstatus))
{ errmsg(300,rtnnme,(int) xkstatus,
consys_nme(orig_sys,'v',oxkndx,FALSE,NULL),oxkndx) ;
return (FALSE) ; }
if (flgoff(xkstatus,vstatNONBASIC|vstatNBFR))
{ errmsg(433,rtnnme,
dy_sys->nme,dy_prtlpphase(dy_lp->phase,TRUE),dy_lp->tot.iters,
"inactive",consys_nme(orig_sys,'v',oxkndx,TRUE,NULL),oxkndx,
dy_prtvstat(xkstatus)) ;
return (FALSE) ; }
# endif
/*
The bounds can change arbitrarily, and the client may not be maintaining
the status vector, but we're limited in what we can do --- bounds and status
are our only clues to the value of an inactive variable. (Contrast with the
equivalent section in process_active.)
*/
lk = orig_sys->vlb[oxkndx] ;
uk = orig_sys->vub[oxkndx] ;
ck = orig_sys->obj[oxkndx] ;
/*
Start with the case that both bounds are finite. Use a previous status of
NBLB or NBUB. Otherwise, guess from the sign of the objective coefficient.
`Dirty' fixed variables are marked as unloadable.
*/
if (lk > -dy_tols->inf && uk < dy_tols->inf)
{ if (atbnd(lk,uk) && lk != uk)
{ if (flgon(xkstatus,vstatNBLB|vstatNBUB))
{ setflg(xkstatus,vstatNOLOAD) ; }
else
{ if (ck < 0)
{ xkstatus = vstatNBUB|vstatNOLOAD ; }
else
{ xkstatus = vstatNBLB|vstatNOLOAD ; } }
# ifndef DYLP_NDEBUG
if (dy_opts->print.setup >= 3)
{ dyio_outfmt(dy_logchn,dy_gtxecho,"\n\tDirty fixed variable %s (%d)",
consys_nme(orig_sys,'v',oxkndx,0,0),oxkndx) ;
dyio_outfmt(dy_logchn,dy_gtxecho,
" assigned status %s.",dy_prtvstat(xkstatus)) ;
dyio_outfmt(dy_logchn,dy_gtxecho,
"\n\t original lb = %g, ub = %g, diff = %g, tol = %g",
lk,uk,uk-lk,dy_tols->pfeas) ; }
# endif
}
else
if (lk == uk)
{ xkstatus = vstatNBFX|vstatNOLOAD ; }
else
if (flgon(xkstatus,vstatNBLB|vstatNBUB))
{ xkstatus = orig_lp->status[oxkndx] ; }
else
{ if (ck < 0)
{ xkstatus = vstatNBUB ; }
else
{ xkstatus = vstatNBLB ; } } }
/*
Variables with one bound, or no bounds. No choices here.
*/
else
if (lk > -dy_tols->inf)
{ xkstatus = vstatNBLB ; }
else
if (uk < dy_tols->inf)
{ xkstatus = vstatNBUB ; }
else
{ xkstatus = vstatNBFR ; }
/*
Determine the variable's value and set up the status entries.
The default case in the switch below should never execute, but it serves
for paranoia and lets gcc conclude xk will always have a value.
Consider whether it's really a good idea to change orig_lp->status.
*/
switch (getflg(xkstatus,vstatSTATUS))
{ case vstatNBLB:
case vstatNBFX:
{ xk = lk ;
break ; }
case vstatNBUB:
{ xk = uk ;
break ; }
case vstatNBFR:
{ xk = 0 ;
break ; }
default:
{ xk = 0 ;
errmsg(1,rtnnme,__LINE__) ;
return (FALSE) ; } }
orig_lp->status[oxkndx] = xkstatus ;
dy_origvars[oxkndx] = -((int) xkstatus) ;
/*
Note any contribution to the objective and constraint rhs & rhslow values.
*/
dy_lp->inactzcorr += xk*orig_sys->obj[oxkndx] ;
if (flgon(orig_lp->ctlopts,lpctlRHSCHG|lpctlLBNDCHG|lpctlUBNDCHG))
{ ak = NULL ;
if (consys_getcol_pk(orig_sys,oxkndx,&ak) == FALSE)
{ errmsg(122,rtnnme,orig_sys->nme,"variable",
consys_nme(orig_sys,'v',oxkndx,TRUE,NULL),oxkndx) ;
if (ak != NULL) pkvec_free(ak) ;
return (FALSE) ; }
for (ndx = 0, aik = &ak->coeffs[0] ; ndx < ak->cnt ; ndx++, aik++)
{ oaindx = aik->ndx ;
if (ACTIVE_CON(oaindx))
{ aindx = dy_origcons[oaindx] ;
dy_sys->rhs[aindx] -= aik->val*xk ;
if (dy_sys->ctyp[aindx] == contypRNG)
dy_sys->rhslow[aindx] -= aik->val*xk ; } }
pkvec_free(ak) ; }
/*
And we're done. Print some information and return.
*/
# ifndef DYLP_NDEBUG
if (dy_opts->print.crash >= 4)
{ dyio_outfmt(dy_logchn,dy_gtxecho,"\n\t %s (%d) %s inactive with value ",
consys_nme(orig_sys,'v',oxkndx,FALSE,NULL),oxkndx,
dy_prtvstat(xkstatus)) ;
switch (getflg(xkstatus,vstatSTATUS))
{ case vstatNBFX:
case vstatNBLB:
case vstatNBUB:
case vstatNBFR:
{ dyio_outfmt(dy_logchn,dy_gtxecho,"%g.",xk) ;
break ; }
default:
{ dyio_outfmt(dy_logchn,dy_gtxecho,"??.") ;
break ; } } }
# endif
return (TRUE) ; }
