static void next_char(glp_data *data)
{ /* read next character */
int c;
if (data->c == XEOF)
glp_sdf_error(data, "unexpected end of file\n");
else if (data->c == '\n')
data->count++;
c = xfgetc(data->fp);
if (c < 0)
{ if (xferror(data->fp))
glp_sdf_error(data, "read error - %s\n", xerrmsg());
else if (data->c == '\n')
c = XEOF;
else
{ glp_sdf_warning(data, "missing final end of line\n");
c = '\n';
}
}
else if (c == '\n')
;
else if (isspace(c))
c = ' ';
else if (iscntrl(c))
glp_sdf_error(data, "invalid control character 0x%02X\n", c);
data->c = c;
return;
}
int glp_read_ipt(glp_prob *lp, const char *fname)
{ glp_data *data;
jmp_buf jump;
int i, j, k, ret = 0;
xprintf("Reading interior-point solution from '%s'...\n", fname);
data = glp_sdf_open_file(fname);
if (data == NULL)
{ ret = 1;
goto done;
}
if (setjmp(jump))
{ ret = 1;
goto done;
}
glp_sdf_set_jump(data, jump);
/* number of rows, number of columns */
k = glp_sdf_read_int(data);
if (k != lp->m)
glp_sdf_error(data, "wrong number of rows\n");
k = glp_sdf_read_int(data);
if (k != lp->n)
glp_sdf_error(data, "wrong number of columns\n");
/* solution status, objective value */
k = glp_sdf_read_int(data);
if (!(k == GLP_UNDEF || k == GLP_OPT))
glp_sdf_error(data, "invalid solution status\n");
lp->ipt_stat = k;
lp->ipt_obj = glp_sdf_read_num(data);
/* rows (auxiliary variables) */
for (i = 1; i <= lp->m; i++)
{ GLPROW *row = lp->row[i];
/* primal value, dual value */
row->pval = glp_sdf_read_num(data);
row->dval = glp_sdf_read_num(data);
}
/* columns (structural variables) */
for (j = 1; j <= lp->n; j++)
{ GLPCOL *col = lp->col[j];
/* primal value, dual value */
col->pval = glp_sdf_read_num(data);
col->dval = glp_sdf_read_num(data);
}
xprintf("%d lines were read\n", glp_sdf_line(data));
done:
if (ret) lp->ipt_stat = GLP_UNDEF;
if (data != NULL) glp_sdf_close_file(data);
return ret;
}
const char *glp_sdf_read_text(glp_data *data)
{ /* read text until end of line */
int c, len = 0;
for (;;)
{ c = data->c;
next_char(data);
if (c == ' ')
{ /* ignore initial spaces */
if (len == 0) continue;
/* and multiple ones */
if (data->item[len-1] == ' ') continue;
}
else if (c == '\n')
{ /* remove trailing space */
if (len > 0 && data->item[len-1] == ' ') len--;
/* and stop reading */
break;
}
/* add current character to the buffer */
data->item[len++] = (char)c;
if (len == sizeof(data->item))
glp_sdf_error(data, "line too long\n", data->item);
}
data->item[len] = '\0';
return data->item;
}
int glp_read_mip(glp_prob *mip, const char *fname)
{ glp_data *data;
jmp_buf jump;
int i, j, k, ret = 0;
xprintf("Reading MIP solution from '%s'...\n", fname);
data = glp_sdf_open_file(fname);
if (data == NULL)
{ ret = 1;
goto done;
}
if (setjmp(jump))
{ ret = 1;
goto done;
}
glp_sdf_set_jump(data, jump);
/* number of rows, number of columns */
k = glp_sdf_read_int(data);
if (k != mip->m)
glp_sdf_error(data, "wrong number of rows\n");
k = glp_sdf_read_int(data);
if (k != mip->n)
glp_sdf_error(data, "wrong number of columns\n");
/* solution status, objective value */
k = glp_sdf_read_int(data);
if (!(k == GLP_UNDEF || k == GLP_OPT || k == GLP_FEAS ||
k == GLP_NOFEAS))
glp_sdf_error(data, "invalid solution status\n");
mip->mip_stat = k;
mip->mip_obj = glp_sdf_read_num(data);
/* rows (auxiliary variables) */
for (i = 1; i <= mip->m; i++)
{ GLPROW *row = mip->row[i];
row->mipx = glp_sdf_read_num(data);
}
/* columns (structural variables) */
for (j = 1; j <= mip->n; j++)
{ GLPCOL *col = mip->col[j];
col->mipx = glp_sdf_read_num(data);
if (col->kind == GLP_IV && col->mipx != floor(col->mipx))
glp_sdf_error(data, "non-integer column value");
}
xprintf("%d lines were read\n", glp_sdf_line(data));
done:
if (ret) mip->mip_stat = GLP_UNDEF;
if (data != NULL) glp_sdf_close_file(data);
return ret;
}
double glp_sdf_read_num(glp_data *data)
{ /* read floating-point number */
double x;
next_item(data);
switch (str2num(data->item, &x))
{ case 0:
break;
case 1:
glp_sdf_error(data, "number `%s' out of range\n",
data->item);
case 2:
glp_sdf_error(data, "cannot convert `%s' to number\n",
data->item);
default:
xassert(data != data);
}
return x;
}
int glp_sdf_read_int(glp_data *data)
{ /* read integer number */
int x;
next_item(data);
switch (str2int(data->item, &x))
{ case 0:
break;
case 1:
glp_sdf_error(data, "integer `%s' out of range\n",
data->item);
case 2:
glp_sdf_error(data, "cannot convert `%s' to integer\n",
data->item);
default:
xassert(data != data);
}
return x;
}
int glp_read_graph(glp_graph *G, const char *fname)
{ glp_data *data;
jmp_buf jump;
int nv, na, i, j, k, ret;
glp_erase_graph(G, G->v_size, G->a_size);
xprintf("Reading graph from `%s'...\n", fname);
data = glp_sdf_open_file(fname);
if (data == NULL)
{ ret = 1;
goto done;
}
if (setjmp(jump))
{ ret = 1;
goto done;
}
glp_sdf_set_jump(data, jump);
nv = glp_sdf_read_int(data);
if (nv < 0)
glp_sdf_error(data, "invalid number of vertices\n");
na = glp_sdf_read_int(data);
if (na < 0)
glp_sdf_error(data, "invalid number of arcs\n");
xprintf("Graph has %d vert%s and %d arc%s\n",
nv, nv == 1 ? "ex" : "ices", na, na == 1 ? "" : "s");
if (nv > 0) glp_add_vertices(G, nv);
for (k = 1; k <= na; k++)
{ i = glp_sdf_read_int(data);
if (!(1 <= i && i <= nv))
glp_sdf_error(data, "tail vertex number out of range\n");
j = glp_sdf_read_int(data);
if (!(1 <= j && j <= nv))
glp_sdf_error(data, "head vertex number out of range\n");
glp_add_arc(G, i, j);
}
xprintf("%d lines were read\n", glp_sdf_line(data));
ret = 0;
done:
if (data != NULL) glp_sdf_close_file(data);
return ret;
}
static void next_item(glp_data *data)
{ /* read next item */
int len;
skip_pad(data);
len = 0;
while (!(data->c == ' ' || data->c == '\n'))
{ data->item[len++] = (char)data->c;
if (len == sizeof(data->item))
glp_sdf_error(data, "data item `%.31s...' too long\n",
data->item);
next_char(data);
}
data->item[len] = '\0';
return;
}
static void skip_pad(glp_data *data)
{ /* skip uninteresting characters and comments */
loop: while (data->c == ' ' || data->c == '\n')
next_char(data);
if (data->c == '/')
{ next_char(data);
if (data->c != '*')
glp_sdf_error(data, "invalid use of slash\n");
next_char(data);
for (;;)
{ if (data->c == '*')
{ next_char(data);
if (data->c == '/')
{ next_char(data);
break;
}
}
next_char(data);
}
goto loop;
}
return;
}
int glp_read_sol(glp_prob *lp, const char *fname)
{ glp_data *data;
jmp_buf jump;
int i, j, k, ret = 0;
xprintf("Reading basic solution from '%s'...\n", fname);
data = glp_sdf_open_file(fname);
if (data == NULL)
{ ret = 1;
goto done;
}
if (setjmp(jump))
{ ret = 1;
goto done;
}
glp_sdf_set_jump(data, jump);
/* number of rows, number of columns */
k = glp_sdf_read_int(data);
if (k != lp->m)
glp_sdf_error(data, "wrong number of rows\n");
k = glp_sdf_read_int(data);
if (k != lp->n)
glp_sdf_error(data, "wrong number of columns\n");
/* primal status, dual status, objective value */
k = glp_sdf_read_int(data);
if (!(k == GLP_UNDEF || k == GLP_FEAS || k == GLP_INFEAS ||
k == GLP_NOFEAS))
glp_sdf_error(data, "invalid primal status\n");
lp->pbs_stat = k;
k = glp_sdf_read_int(data);
if (!(k == GLP_UNDEF || k == GLP_FEAS || k == GLP_INFEAS ||
k == GLP_NOFEAS))
glp_sdf_error(data, "invalid dual status\n");
lp->dbs_stat = k;
lp->obj_val = glp_sdf_read_num(data);
/* rows (auxiliary variables) */
for (i = 1; i <= lp->m; i++)
{ GLPROW *row = lp->row[i];
/* status, primal value, dual value */
k = glp_sdf_read_int(data);
if (!(k == GLP_BS || k == GLP_NL || k == GLP_NU ||
k == GLP_NF || k == GLP_NS))
glp_sdf_error(data, "invalid row status\n");
glp_set_row_stat(lp, i, k);
row->prim = glp_sdf_read_num(data);
row->dual = glp_sdf_read_num(data);
}
/* columns (structural variables) */
for (j = 1; j <= lp->n; j++)
{ GLPCOL *col = lp->col[j];
/* status, primal value, dual value */
k = glp_sdf_read_int(data);
if (!(k == GLP_BS || k == GLP_NL || k == GLP_NU ||
k == GLP_NF || k == GLP_NS))
glp_sdf_error(data, "invalid column status\n");
glp_set_col_stat(lp, j, k);
col->prim = glp_sdf_read_num(data);
col->dual = glp_sdf_read_num(data);
}
xprintf("%d lines were read\n", glp_sdf_line(data));
done:
if (ret) lp->pbs_stat = lp->dbs_stat = GLP_UNDEF;
if (data != NULL) glp_sdf_close_file(data);
return ret;
}
