highspy_basic_linear_programming_solver_quickstart.py

1import highspy
2import numpy as np
3
4# Create a HiGHS instance
5h = highspy.Highs()
6
7# Define the model:
8# Maximize x0 + x1
9# Subject to:
10# x0 + 2x1 <= 2
11# x0, x1 >= 0
12
13# Set objective sense
14h.changeObjectiveSense(highspy.ObjSense.kMaximize)
15
16# Add variables (columns)
17# addVar(lower_bound, upper_bound, objective_coefficient)
18h.addVar(0.0, np.inf, 1.0)
19h.addVar(0.0, np.inf, 1.0)
20
21# Add constraints (rows)
22# addRow(lower_bound, upper_bound, num_non_zero, indices, values)
23h.addRow(-np.inf, 2.0, 2, np.array([0, 1]), np.array([1.0, 2.0]))
24
25# Solve the model
26h.run()
27
28# Get the solution
29solution = h.getSolution()
30print(f"Status: {h.getModelStatus()}")
31print(f"Objective value: {h.getObjectiveValue()}")
32print(f"x0 = {solution.col_value[0]}")
33print(f"x1 = {solution.col_value[1]}")