highs-sys 1.14.2

Rust binding for the HiGHS linear programming solver. See http://highs.dev.
Documentation 
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

#include "HCheckConfig.h"
#include "catch.hpp"
#include "util/HighsCDouble.h"
#include "util/HighsRandom.h"

void testCeil(HighsCDouble x) {
  double ceil_x;
  double double_x;
  ceil_x = double(ceil(x));
  double_x = double(x);
  REQUIRE(ceil_x >= double_x);
  REQUIRE(ceil(x) >= x);
}

void testFloor(HighsCDouble x) {
  double floor_x;
  double double_x;
  floor_x = double(floor(x));
  double_x = double(x);
  REQUIRE(floor_x <= double_x);
  REQUIRE(floor(x) <= x);
}

TEST_CASE("HighsCDouble-ceil", "[util]") {
  HighsCDouble x;
  x = -1e-34;
  testCeil(x);
  x = -1e-32;
  testCeil(x);
  x = -1e-30;
  testCeil(x);
  x = -1e-23;
  testCeil(x);
  x = -1e-12;
  testCeil(x);
  x = -1e-1;
  testCeil(x);
  x = -0.99;
  testCeil(x);

  x = 0.99;
  testCeil(x);
  x = 1e-1;
  testCeil(x);
  x = 1e-12;
  testCeil(x);
  // This and rest failed in #2041
  x = 1e-23;
  testCeil(x);
  x = 1e-30;
  testCeil(x);
  x = 1e-32;
  testCeil(x);
  x = 1e-34;
  testCeil(x);

  HighsRandom rand;
  for (HighsInt k = 0; k < 1000; k++) {
    double man = rand.fraction();
    HighsInt power = 2 - rand.integer(5);
    double exp = std::pow(10, power);
    x = man * exp;
    testCeil(x);
  }
}

TEST_CASE("HighsCDouble-floor", "[util]") {
  HighsCDouble x;

  x = 1e-34;
  testFloor(x);
  x = 1e-32;
  testFloor(x);
  x = 1e-30;
  testFloor(x);
  x = 1e-23;
  testFloor(x);
  x = 1e-12;
  testFloor(x);
  x = 1e-1;
  testFloor(x);
  x = 0.99;
  testFloor(x);

  x = -0.99;
  testFloor(x);
  x = -1e-1;
  testFloor(x);
  x = -1e-12;
  testFloor(x);
  // This and rest failed in #2041
  x = -1e-23;
  testFloor(x);
  x = -1e-30;
  testFloor(x);
  x = -1e-32;
  testFloor(x);
  x = -1e-34;
  testFloor(x);

  HighsRandom rand;
  for (HighsInt k = 0; k < 1000; k++) {
    double man = rand.fraction();
    HighsInt power = 2 - rand.integer(5);
    double exp = std::pow(10, power);
    x = man * exp;
    testFloor(x);
  }
}