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
125
126
127
128
129
130
131
132
133
134
135
|
/***************************************************************************
* Copyright (C) 2009 by BUI Quang Minh *
* minh.bui@univie.ac.at *
* *
* This program 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 2 of the License, or *
* (at your option) any later version. *
* *
* This program 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 this program; if not, write to the *
* Free Software Foundation, Inc., *
* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. *
***************************************************************************/
#include <string.h>
#include <sstream>
#include "utils/tools.h"
#include "gurobiwrapper.h"
#define tolerance 0.000001
/**
interface to call GUROBI LP solver
@param filename name of input lp file
@param ntaxa number of taxa
@param score (OUT) returned optimal score
@param variables (OUT) array of returned solution
@param verbose_mode verbose mode
@return
-1 if gurobi was not installed properly or does not exist at all
0 if everything works file,
5 if solution is not optimal,
6 if some variable has wrong name,
7 if returned solution is not binary. In this case, one should run the solver
again with strict binary variable constraint.
*/
int gurobi_solve(char *filename, int ntaxa, double *score, double *variables, int verbose_mode, int num_threads) {
int ret = 0;
*score = -1;
string command;
ostringstream ss;
ss << "gurobi_cl Threads=" << num_threads << " ResultFile=" << filename
<< ".sol MIPGap=0 "<< filename << " >" << filename << ".log ";
command = ss.str();
if (verbose_mode >= VB_MED)
cout << command << endl;
int sys_ret = system(command.c_str());
if (sys_ret != 0) {
cout << "gurobi_cl could not be executed. Make sure it was installed with proper license." << endl;
cout << command << endl;
return -1;
}
command = filename;
command += ".sol";
try {
ifstream in;
in.exceptions(ios::failbit | ios::badbit);
in.open(command.c_str());
string str;
while (!in.eof()) {
// remove the failbit
in.exceptions(ios::badbit);
if(!(in >> str)) break;
// set the failbit again
in.exceptions(ios::failbit | ios::badbit);
if (str[0] != 'x') continue;
int index = convert_int(str.substr(1).c_str());
if (index < 0 || index >= ntaxa) {
cout << "Index x_" << index << " is not in the range!" << endl;
ret = 6;
break;
}
double value;
in >> value;
if (value > tolerance && (1.0 - value) > tolerance) {
if (verbose_mode >= VB_MED) cout << endl << str << " = " << value;
ret = 7;
if (!verbose_mode) break;
}
variables[index] = value;
}
in.clear();
// set the failbit again
in.exceptions(ios::failbit | ios::badbit);
in.close();
} catch(ios::failure) {
outError(ERR_READ_INPUT);
} catch (string str) {
outError(str);
}
command = filename;
command += ".log";
try {
ifstream in;
in.exceptions(ios::failbit | ios::badbit);
in.open(command.c_str());
string str;
while (!in.eof()) {
in.exceptions(ios::badbit);
if(!(in >> str)) break;
// set the failbit again
in.exceptions(ios::failbit | ios::badbit);
if (str != "Best" && str != "Optimal") continue;
in >> str;
if (str != "objective") continue;
in >> str;
// remove the ending comma ,
if (*str.rbegin() == ',') str.erase(str.length()-1);
*score = convert_double(str.c_str());
break;
}
in.clear();
// set the failbit again
in.exceptions(ios::failbit | ios::badbit);
in.close();
} catch(ios::failure) {
outError(ERR_READ_INPUT);
} catch (string str) {
outError(str);
}
return ret;
}
|
