1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
|
/* wript.c (write interior-point solution in GLPK format) */
/***********************************************************************
* This code is part of GLPK (GNU Linear Programming Kit).
*
* Copyright (C) 2010-2016 Andrew Makhorin, Department for Applied
* Informatics, Moscow Aviation Institute, Moscow, Russia. All rights
* reserved. E-mail: <mao@gnu.org>.
*
* GLPK is free software: you can redistribute it and/or modify it
* under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* GLPK is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
* or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
* License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GLPK. If not, see <http://www.gnu.org/licenses/>.
***********************************************************************/
#include "env.h"
#include "prob.h"
/***********************************************************************
* NAME
*
* glp_write_ipt - write interior-point solution in GLPK format
*
* SYNOPSIS
*
* int glp_write_ipt(glp_prob *P, const char *fname);
*
* DESCRIPTION
*
* The routine glp_write_ipt writes interior-point solution to a text
* file in GLPK format.
*
* RETURNS
*
* If the operation was successful, the routine returns zero. Otherwise
* it prints an error message and returns non-zero. */
int glp_write_ipt(glp_prob *P, const char *fname)
{ glp_file *fp;
GLPROW *row;
GLPCOL *col;
int i, j, count, ret = 1;
char *s;
#if 0 /* 04/IV-2016 */
if (P == NULL || P->magic != GLP_PROB_MAGIC)
xerror("glp_write_ipt: P = %p; invalid problem object\n", P);
#endif
if (fname == NULL)
xerror("glp_write_ipt: fname = %d; invalid parameter\n", fname)
;
xprintf("Writing interior-point solution to '%s'...\n", fname);
fp = glp_open(fname, "w"), count = 0;
if (fp == NULL)
{ xprintf("Unable to create '%s' - %s\n", fname, get_err_msg());
goto done;
}
/* write comment lines */
glp_format(fp, "c %-12s%s\n", "Problem:",
P->name == NULL ? "" : P->name), count++;
glp_format(fp, "c %-12s%d\n", "Rows:", P->m), count++;
glp_format(fp, "c %-12s%d\n", "Columns:", P->n), count++;
glp_format(fp, "c %-12s%d\n", "Non-zeros:", P->nnz), count++;
switch (P->ipt_stat)
{ case GLP_OPT: s = "OPTIMAL"; break;
case GLP_INFEAS: s = "INFEASIBLE (INTERMEDIATE)"; break;
case GLP_NOFEAS: s = "INFEASIBLE (FINAL)"; break;
case GLP_UNDEF: s = "UNDEFINED"; break;
default: s = "???"; break;
}
glp_format(fp, "c %-12s%s\n", "Status:", s), count++;
switch (P->dir)
{ case GLP_MIN: s = "MINimum"; break;
case GLP_MAX: s = "MAXimum"; break;
default: s = "???"; break;
}
glp_format(fp, "c %-12s%s%s%.10g (%s)\n", "Objective:",
P->obj == NULL ? "" : P->obj,
P->obj == NULL ? "" : " = ", P->ipt_obj, s), count++;
glp_format(fp, "c\n"), count++;
/* write solution line */
glp_format(fp, "s ipt %d %d ", P->m, P->n), count++;
switch (P->ipt_stat)
{ case GLP_OPT: glp_format(fp, "o"); break;
case GLP_INFEAS: glp_format(fp, "i"); break;
case GLP_NOFEAS: glp_format(fp, "n"); break;
case GLP_UNDEF: glp_format(fp, "u"); break;
default: glp_format(fp, "?"); break;
}
glp_format(fp, " %.*g\n", DBL_DIG, P->ipt_obj);
/* write row solution descriptor lines */
for (i = 1; i <= P->m; i++)
{ row = P->row[i];
glp_format(fp, "i %d %.*g %.*g\n", i, DBL_DIG, row->pval,
DBL_DIG, row->dval), count++;
}
/* write column solution descriptor lines */
for (j = 1; j <= P->n; j++)
{ col = P->col[j];
glp_format(fp, "j %d %.*g %.*g\n", j, DBL_DIG, col->pval,
DBL_DIG, col->dval), count++;
}
/* write end line */
glp_format(fp, "e o f\n"), count++;
if (glp_ioerr(fp))
{ xprintf("Write error on '%s' - %s\n", fname, get_err_msg());
goto done;
}
/* interior-point solution has been successfully written */
xprintf("%d lines were written\n", count);
ret = 0;
done: if (fp != NULL)
glp_close(fp);
return ret;
}
/* eof */
|
