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
|
#include <algorithm>
#include <cassert>
#include "HCheckConfig.h"
#include "Highs.h"
#include "SpecialLps.h"
#include "catch.hpp"
const bool dev_run = false;
// No commas in test case name.
TEST_CASE("highs-model", "[highs_model]") {
std::string filename;
filename = std::string(HIGHS_DIR) + "/check/instances/adlittle.mps";
HighsStatus status;
Highs highs;
if (!dev_run) highs.setOptionValue("output_flag", false);
highs.readModel(filename);
const HighsLp& lp = highs.getLp();
HighsModel model;
model.lp_ = lp;
// Add an identity Hessian
HighsHessian& hessian = model.hessian_;
HighsInt dim = lp.num_col_;
hessian.dim_ = dim;
hessian.start_.resize(dim + 1);
hessian.index_.resize(dim);
hessian.value_.resize(dim);
for (int iCol = 0; iCol < dim; iCol++) {
hessian.start_[iCol] = iCol;
hessian.index_[iCol] = iCol;
hessian.value_[iCol] = 1.0;
}
hessian.start_[dim] = dim;
status = highs.passModel(model);
REQUIRE(status == HighsStatus::kOk);
status = highs.run();
REQUIRE(status == HighsStatus::kOk);
highs.clear();
model.clear();
SpecialLps special_lps;
HighsModelStatus require_model_status;
double optimal_objective;
special_lps.distillationLp(model.lp_, require_model_status,
optimal_objective);
dim = model.lp_.num_col_;
// A Hessian with dimesion but no nonzeros (ie identically zero)
// should be ignored OK
hessian.dim_ = dim;
hessian.start_.resize(dim + 1);
hessian.start_[0] = 0;
hessian.start_[dim] = 0;
if (!dev_run) highs.setOptionValue("output_flag", false);
status = highs.passModel(model);
REQUIRE(status == HighsStatus::kOk);
// A Hessian with small diagonal entries should cause an error. The
// first is eliminated when normalising the Hessian - due to being
// too small - and the second when checking the diagonal entries -
// for being less than small_matrix_value.
const double illegal_small_hessian_diagonal_entry = 1e-12;
const double illegal_negative_hessian_diagonal_entry = -1;
assert(dim == 2);
hessian.start_[1] = 1;
hessian.start_[2] = 2;
hessian.index_.resize(dim);
hessian.index_[0] = 0;
hessian.index_[1] = 1;
hessian.value_.resize(dim);
hessian.value_[0] = illegal_small_hessian_diagonal_entry;
hessian.value_[1] = illegal_negative_hessian_diagonal_entry;
status = highs.passModel(model);
REQUIRE(status == HighsStatus::kOk);
status = highs.run();
REQUIRE(status == HighsStatus::kError);
}
TEST_CASE("highs-integrality", "[highs_model]") {
HighsLp lp;
HighsModelStatus require_model_status;
double optimal_objective;
lp.model_name_ = "distillation";
lp.num_col_ = 2;
lp.num_row_ = 3;
lp.col_cost_ = {8, 10};
lp.col_lower_ = {0, 0};
lp.col_upper_ = {inf, inf};
lp.row_lower_ = {7, 12, 6};
lp.row_upper_ = {inf, inf, inf};
lp.a_matrix_.start_ = {0, 3, 6};
lp.a_matrix_.index_ = {0, 1, 2, 0, 1, 2};
lp.a_matrix_.value_ = {2, 3, 2, 2, 4, 1};
lp.sense_ = ObjSense::kMinimize;
lp.offset_ = 0;
lp.a_matrix_.format_ = MatrixFormat::kColwise;
require_model_status = HighsModelStatus::kOptimal;
optimal_objective = 31.2;
Highs highs;
highs.setOptionValue("output_flag", dev_run);
highs.passModel(lp);
HighsVarType integrality;
REQUIRE(highs.getColIntegrality(-1, integrality) == HighsStatus::kError);
REQUIRE(highs.getColIntegrality(0, integrality) == HighsStatus::kError);
REQUIRE(highs.getColIntegrality(lp.num_col_, integrality) ==
HighsStatus::kError);
lp.integrality_ = {HighsVarType::kContinuous, HighsVarType::kContinuous};
highs.passModel(lp);
REQUIRE(highs.getColIntegrality(-1, integrality) == HighsStatus::kError);
REQUIRE(highs.getColIntegrality(0, integrality) == HighsStatus::kOk);
REQUIRE(integrality == HighsVarType::kContinuous);
REQUIRE(highs.getColIntegrality(lp.num_col_, integrality) ==
HighsStatus::kError);
highs.run();
REQUIRE(highs.getModelStatus() == require_model_status);
REQUIRE(std::abs(highs.getInfo().objective_function_value -
optimal_objective) < 1e-5);
REQUIRE(highs.changeColIntegrality(0, HighsVarType::kInteger) ==
HighsStatus::kOk);
REQUIRE(highs.getColIntegrality(0, integrality) == HighsStatus::kOk);
REQUIRE(integrality == HighsVarType::kInteger);
}
|
