def task_1(self):
        A = [[2, 1, -1, 0, 0, 1], [1, 0, 1, 1, 0, 0], [0, 1, 0, 0, 1, 0]]
        b = [2, 5, 0]
        c = [3, 2, 0, 3, -2, -4]
        d_lo = [0, -1, 2, 1, -1, 0]
        d_hi = [2, 4, 4, 3, 3, 5]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve([3, 4, 5])

        x_expected = [3 / 2, 1, 2, 3 / 2, -1, 0]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 3)
        self.assertAlmostEqual(f_val, 13, 5)
    def task_3(self):
        A = [[1, 0, 2, 2, -3, 3], [0, 1, 0, -1, 0, 1], [1, 0, 1, 3, 2, 1]]
        b = [15, 0, 13]
        c = [3, 0.5, 4, 4, 1, 5]
        d_lo = [0, 0, 0, 0, 0, 0]
        d_hi = [3, 5, 4, 3, 3, 4]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [3, 0, 4, 1.1818, 0.6364, 1.1818]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 4)
        self.assertAlmostEqual(f_val, 36.2727, 4)
    def task_2(self):
        A = [[1, -5, 3, 1, 0, 0], [4, -1, 1, 0, 1, 0], [2, 4, 2, 0, 0, 1]]
        b = [-7, 22, 30]
        c = [7, -2, 6, 0, 5, 2]
        d_lo = [2, 1, 0, 0, 1, 1]
        d_hi = [6, 6, 5, 2, 4, 6]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [5, 3, 1, 0, 4, 6]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 5)
        self.assertAlmostEqual(f_val, 67, 5)
    def task_4(self):
        A = [[1, 0, 0, 12, 1, -3, 4, -1], [0, 1, 0, 11, 12, 3, 5, 3],
             [0, 0, 1, 1, 0, 22, -2, 1]]
        b = [40, 107, 61]
        c = [2, 1, -2, -1, 4, -5, 5, 5]
        d_lo = [0, 0, 0, 0, 0, 0, 0, 0]
        d_hi = [3, 5, 5, 3, 4, 5, 6, 3]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [3, 5, 0, 1.8779, 2.7545, 3.0965, 6, 3]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 4)
        self.assertAlmostEqual(f_val, 49.6577, 4)
    def task_9(self):
        A = [[1, 3, 1, -1, 0, -3, 2, 1], [2, 1, 3, -1, 1, 4, 1, 1],
             [-1, 0, 2, -2, 2, 1, 1, 1]]
        b = [4, 12, 4]
        c = [2, -1, 2, 3, -2, 3, 4, 1]
        d_lo = [-1] * 8
        d_hi = [2, 3, 1, 4, 3, 2, 4, 4]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [-1, 0.4074, 1, 4, -0.3704, 1.7407, 4, 4]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 4)
        self.assertAlmostEqual(f_val, 37.5556, 4)
    def task_7(self):
        A = [[2, -1, 1, 0, 0, -1, 3], [0, 4, -1, 2, 3, -2, 2],
             [3, 1, 0, 1, 0, 1, 4]]
        b = [1.5, 9, 2]
        c = [0, 1, 2, 1, -3, 4, 7]
        d_lo = [0, 0, -3, 0, -1, 1, 0]
        d_hi = [3, 3, 4, 7, 5, 3, 2]

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [0, 1, 3.5, 0, 3.5, 1, 0]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 4)
        self.assertAlmostEqual(f_val, 1.5, 4)
    def task_5(self):
        A = [[1, -3, 2, 0, 1, -1, 4, -1, 0], [1, -1, 6, 1, 0, -2, 2, 2, 0],
             [2, 2, -1, 1, 0, -3, 8, -1, 1], [4, 1, 0, 0, 1, -1, 0, -1, 1],
             [1, 1, 1, 1, 1, 1, 1, 1, 1]]
        b = [3, 9, 9, 5, 9]
        c = [-1, 5, -2, 4, 3, 1, 2, 8, 3]
        d_lo = [0] * 9
        d_hi = [5] * 9

        task = DualTask(A, b, c, d_lo, d_hi)
        dual = DualSimplex(task)
        x, J, f_val = dual.solve()

        x_expected = [1.1579, 0.6942, 0, 0, 2.8797, 0, 1.0627, 3.2055, 0]
        for i, xi in enumerate(x):
            self.assertAlmostEqual(xi, x_expected[i], 4)
        self.assertAlmostEqual(f_val, 38.7218, 4)
Пример #8
0
    def solve(self):
        # print('\n==================================')

        if self.iteration == 10000:
            raise Exception("Достигнут лимит итераций")

        if len(self.queue) == 0:
            return

        task = self.queue.pop(0)
        dual = DualSimplex(task)

        # print('Solving task: ')
        # print(task.d_lo)
        # print(task.d_hi)

        try:
            x, _, f_val = dual.solve()
            x = np.around(x, decimals=4)
        except Exception as e:
            # print(e)
            return

        if f_val < self.mu:
            return

        if all_int(x):
            if self.mu < f_val:
                self.mu = f_val
                self.x = x
            return

        task1, task2 = self.__split_task(task, x)

        # print(x)
        # print('-----')

        self.queue.append(task1)
        self.queue.append(task2)
        self.solve()
        self.solve()