int
main(int argc, char **argv)
{
zvm_driver_t zvm;
int fence,
rc = 0;
openlog ("fence_zvm", LOG_CONS|LOG_PID, LOG_DAEMON);
memset(&zvm, 0, sizeof(zvm));
zvm.timeOut = DEFAULT_TIMEOUT;
zvm.delay = DEFAULT_DELAY;
if (argc > 1)
fence = get_options(argc, argv, &zvm);
else
fence = get_options_stdin(&zvm);
switch(fence) {
case ACT_OFFON : // OFFON
if ((rc = check_parm(&zvm)) == 0)
rc = zvm_smapi_imageRecycle(&zvm);
break;
case ACT_OFF : // OFF
if ((rc = check_parm(&zvm)) == 0)
rc = zvm_smapi_imageDeactivate(&zvm);
break;
case ACT_ON : // ON
if ((rc = check_parm(&zvm)) == 0)
rc = zvm_smapi_imageActivate(&zvm);
break;
case ACT_METADATA : // METADATA
rc = zvm_metadata();
break;
case ACT_STATUS : // STATUS
if ((rc = check_parm(&zvm)) == 0)
rc = zvm_smapi_imageQuery(&zvm);
break;
case ACT_MONITOR : // MONITOR
rc = 0;
break;
case ACT_LIST : // LIST
printf("N/A");
rc = 0;
break;
case ACT_HELP :
rc = usage();
}
closelog();
return (rc);
}
int glp_read_lp(glp_prob *P, const glp_cpxcp *parm, const char *fname)
{ /* read problem data in CPLEX LP format */
glp_cpxcp _parm;
struct csa _csa, *csa = &_csa;
int ret;
xprintf("Reading problem data from '%s'...\n", fname);
if (parm == NULL)
glp_init_cpxcp(&_parm), parm = &_parm;
/* check control parameters */
check_parm("glp_read_lp", parm);
/* initialize common storage area */
csa->P = P;
csa->parm = parm;
csa->fname = fname;
csa->fp = NULL;
if (setjmp(csa->jump))
{ ret = 1;
goto done;
}
csa->count = 0;
csa->c = '\n';
csa->token = T_EOF;
csa->image[0] = '\0';
csa->imlen = 0;
csa->value = 0.0;
csa->n_max = 100;
csa->ind = xcalloc(1+csa->n_max, sizeof(int));
csa->val = xcalloc(1+csa->n_max, sizeof(double));
csa->flag = xcalloc(1+csa->n_max, sizeof(char));
memset(&csa->flag[1], 0, csa->n_max * sizeof(char));
csa->lb = xcalloc(1+csa->n_max, sizeof(double));
csa->ub = xcalloc(1+csa->n_max, sizeof(double));
#if 1 /* 27/VII-2013 */
csa->lb_warn = csa->ub_warn = 0;
#endif
/* erase problem object */
glp_erase_prob(P);
glp_create_index(P);
/* open input CPLEX LP file */
csa->fp = glp_open(fname, "r");
if (csa->fp == NULL)
{ xprintf("Unable to open '%s' - %s\n", fname, get_err_msg());
ret = 1;
goto done;
}
/* scan very first token */
scan_token(csa);
/* parse definition of the objective function */
if (!(csa->token == T_MINIMIZE || csa->token == T_MAXIMIZE))
error(csa, "'minimize' or 'maximize' keyword missing\n");
parse_objective(csa);
/* parse constraints section */
if (csa->token != T_SUBJECT_TO)
error(csa, "constraints section missing\n");
parse_constraints(csa);
/* parse optional bounds section */
if (csa->token == T_BOUNDS) parse_bounds(csa);
/* parse optional general, integer, and binary sections */
while (csa->token == T_GENERAL ||
csa->token == T_INTEGER ||
csa->token == T_BINARY) parse_integer(csa);
/* check for the keyword 'end' */
if (csa->token == T_END)
scan_token(csa);
else if (csa->token == T_EOF)
warning(csa, "keyword 'end' missing\n");
else
error(csa, "symbol '%s' in wrong position\n", csa->image);
/* nothing must follow the keyword 'end' (except comments) */
if (csa->token != T_EOF)
error(csa, "extra symbol(s) detected beyond 'end'\n");
/* set bounds of variables */
{ int j, type;
double lb, ub;
for (j = 1; j <= P->n; j++)
{ lb = csa->lb[j];
ub = csa->ub[j];
if (lb == +DBL_MAX) lb = 0.0; /* default lb */
if (ub == -DBL_MAX) ub = +DBL_MAX; /* default ub */
if (lb == -DBL_MAX && ub == +DBL_MAX)
type = GLP_FR;
else if (ub == +DBL_MAX)
type = GLP_LO;
else if (lb == -DBL_MAX)
type = GLP_UP;
else if (lb != ub)
type = GLP_DB;
else
type = GLP_FX;
glp_set_col_bnds(csa->P, j, type, lb, ub);
}
}
/* print some statistics */
xprintf("%d row%s, %d column%s, %d non-zero%s\n",
P->m, P->m == 1 ? "" : "s", P->n, P->n == 1 ? "" : "s",
P->nnz, P->nnz == 1 ? "" : "s");
if (glp_get_num_int(P) > 0)
{ int ni = glp_get_num_int(P);
int nb = glp_get_num_bin(P);
if (ni == 1)
{ if (nb == 0)
xprintf("One variable is integer\n");
else
xprintf("One variable is binary\n");
}
else
{ xprintf("%d integer variables, ", ni);
if (nb == 0)
xprintf("none");
else if (nb == 1)
xprintf("one");
else if (nb == ni)
xprintf("all");
else
xprintf("%d", nb);
xprintf(" of which %s binary\n", nb == 1 ? "is" : "are");
}
}
xprintf("%d lines were read\n", csa->count);
/* problem data has been successfully read */
glp_delete_index(P);
glp_sort_matrix(P);
ret = 0;
done: if (csa->fp != NULL) glp_close(csa->fp);
xfree(csa->ind);
xfree(csa->val);
xfree(csa->flag);
xfree(csa->lb);
xfree(csa->ub);
if (ret != 0) glp_erase_prob(P);
return ret;
}
int glp_write_lp(glp_prob *P, const glp_cpxcp *parm, const char *fname)
{ /* write problem data in CPLEX LP format */
glp_cpxcp _parm;
struct csa _csa, *csa = &_csa;
glp_file *fp;
GLPROW *row;
GLPCOL *col;
GLPAIJ *aij;
int i, j, len, flag, count, ret;
char line[1000+1], term[500+1], name[255+1];
xprintf("Writing problem data to '%s'...\n", fname);
if (parm == NULL)
glp_init_cpxcp(&_parm), parm = &_parm;
/* check control parameters */
check_parm("glp_write_lp", parm);
/* initialize common storage area */
csa->P = P;
csa->parm = parm;
/* create output CPLEX LP file */
fp = glp_open(fname, "w"), count = 0;
if (fp == NULL)
{ xprintf("Unable to create '%s' - %s\n", fname, get_err_msg());
ret = 1;
goto done;
}
/* write problem name */
xfprintf(fp, "\\* Problem: %s *\\\n",
P->name == NULL ? "Unknown" : P->name), count++;
xfprintf(fp, "\n"), count++;
/* the problem should contain at least one row and one column */
if (!(P->m > 0 && P->n > 0))
{ xprintf("Warning: problem has no rows/columns\n");
xfprintf(fp, "\\* WARNING: PROBLEM HAS NO ROWS/COLUMNS *\\\n"),
count++;
xfprintf(fp, "\n"), count++;
goto skip;
}
/* write the objective function definition */
if (P->dir == GLP_MIN)
xfprintf(fp, "Minimize\n"), count++;
else if (P->dir == GLP_MAX)
xfprintf(fp, "Maximize\n"), count++;
else
xassert(P != P);
row_name(csa, 0, name);
sprintf(line, " %s:", name);
len = 0;
for (j = 1; j <= P->n; j++)
{ col = P->col[j];
if (col->coef != 0.0 || col->ptr == NULL)
{ len++;
col_name(csa, j, name);
if (col->coef == 0.0)
sprintf(term, " + 0 %s", name); /* empty column */
else if (col->coef == +1.0)
sprintf(term, " + %s", name);
else if (col->coef == -1.0)
sprintf(term, " - %s", name);
else if (col->coef > 0.0)
sprintf(term, " + %.*g %s", DBL_DIG, +col->coef, name);
else
sprintf(term, " - %.*g %s", DBL_DIG, -col->coef, name);
if (strlen(line) + strlen(term) > 72)
xfprintf(fp, "%s\n", line), line[0] = '\0', count++;
strcat(line, term);
}
}
if (len == 0)
{ /* empty objective */
sprintf(term, " 0 %s", col_name(csa, 1, name));
strcat(line, term);
}
xfprintf(fp, "%s\n", line), count++;
if (P->c0 != 0.0)
xfprintf(fp, "\\* constant term = %.*g *\\\n", DBL_DIG, P->c0),
count++;
xfprintf(fp, "\n"), count++;
/* write the constraints section */
xfprintf(fp, "Subject To\n"), count++;
for (i = 1; i <= P->m; i++)
{ row = P->row[i];
if (row->type == GLP_FR) continue; /* skip free row */
row_name(csa, i, name);
sprintf(line, " %s:", name);
/* linear form */
for (aij = row->ptr; aij != NULL; aij = aij->r_next)
{ col_name(csa, aij->col->j, name);
if (aij->val == +1.0)
sprintf(term, " + %s", name);
else if (aij->val == -1.0)
sprintf(term, " - %s", name);
else if (aij->val > 0.0)
sprintf(term, " + %.*g %s", DBL_DIG, +aij->val, name);
else
sprintf(term, " - %.*g %s", DBL_DIG, -aij->val, name);
if (strlen(line) + strlen(term) > 72)
xfprintf(fp, "%s\n", line), line[0] = '\0', count++;
strcat(line, term);
}
if (row->type == GLP_DB)
{ /* double-bounded (ranged) constraint */
sprintf(term, " - ~r_%d", i);
if (strlen(line) + strlen(term) > 72)
xfprintf(fp, "%s\n", line), line[0] = '\0', count++;
strcat(line, term);
}
else if (row->ptr == NULL)
{ /* empty constraint */
sprintf(term, " 0 %s", col_name(csa, 1, name));
strcat(line, term);
}
/* right hand-side */
if (row->type == GLP_LO)
sprintf(term, " >= %.*g", DBL_DIG, row->lb);
else if (row->type == GLP_UP)
sprintf(term, " <= %.*g", DBL_DIG, row->ub);
else if (row->type == GLP_DB || row->type == GLP_FX)
sprintf(term, " = %.*g", DBL_DIG, row->lb);
else
xassert(row != row);
if (strlen(line) + strlen(term) > 72)
xfprintf(fp, "%s\n", line), line[0] = '\0', count++;
strcat(line, term);
xfprintf(fp, "%s\n", line), count++;
}
xfprintf(fp, "\n"), count++;
/* write the bounds section */
flag = 0;
for (i = 1; i <= P->m; i++)
{ row = P->row[i];
if (row->type != GLP_DB) continue;
if (!flag)
xfprintf(fp, "Bounds\n"), flag = 1, count++;
xfprintf(fp, " 0 <= ~r_%d <= %.*g\n",
i, DBL_DIG, row->ub - row->lb), count++;
}
for (j = 1; j <= P->n; j++)
{ col = P->col[j];
if (col->type == GLP_LO && col->lb == 0.0) continue;
if (!flag)
xfprintf(fp, "Bounds\n"), flag = 1, count++;
col_name(csa, j, name);
if (col->type == GLP_FR)
xfprintf(fp, " %s free\n", name), count++;
else if (col->type == GLP_LO)
xfprintf(fp, " %s >= %.*g\n",
name, DBL_DIG, col->lb), count++;
else if (col->type == GLP_UP)
xfprintf(fp, " -Inf <= %s <= %.*g\n",
name, DBL_DIG, col->ub), count++;
else if (col->type == GLP_DB)
xfprintf(fp, " %.*g <= %s <= %.*g\n",
DBL_DIG, col->lb, name, DBL_DIG, col->ub), count++;
else if (col->type == GLP_FX)
xfprintf(fp, " %s = %.*g\n",
name, DBL_DIG, col->lb), count++;
else
xassert(col != col);
}
if (flag) xfprintf(fp, "\n"), count++;
/* write the integer section */
flag = 0;
for (j = 1; j <= P->n; j++)
{ col = P->col[j];
if (col->kind == GLP_CV) continue;
xassert(col->kind == GLP_IV);
if (!flag)
xfprintf(fp, "Generals\n"), flag = 1, count++;
xfprintf(fp, " %s\n", col_name(csa, j, name)), count++;
}
if (flag) xfprintf(fp, "\n"), count++;
skip: /* write the end keyword */
xfprintf(fp, "End\n"), count++;
#if 0 /* FIXME */
xfflush(fp);
#endif
if (glp_ioerr(fp))
{ xprintf("Write error on '%s' - %s\n", fname, get_err_msg());
ret = 1;
goto done;
}
/* problem data has been successfully written */
xprintf("%d lines were written\n", count);
ret = 0;
done: if (fp != NULL) glp_close(fp);
return ret;
}