10 March 2023
In this article we continue the Python Production mix series. Specifically, we build Model 11 using the SciPy library.
SciPy is an open source Python package for a wide range of scientific computing applications, including optimization, integration, interpolation, eigenvalue problems, algebraic equations, differential equations, statistics and many other classes of problems.
Our objective is to compare a linear programming model built using SciPy with the same model built using Pyomo.
Articles in this series
Articles in the Python Production mix series:
- Python optimization Rosetta Stone
- Production mix - Model 1, Pyomo concrete
- Production mix - Model 2, Pyomo separate data
- Production mix - Model 3, Pyomo external data
- Production mix - Model 4, Pyomo json file
- Production mix - Model 5, Pyomo using def
- Production mix - Model 6, Pyomo abstract
- Production mix - Model 7, PuLP
- Production mix - Model 8, OR-Tools
- Production mix - Model 9, Gekko
- Production mix - Model 10, CVXPY
- Production mix - Model 11, SciPy
- Production mix - Conclusions
Download the model
The model is available on GitHub.
Formulation for Model 11
For this model, we're using the same general formulation that we used for previous models in this series, as shown in Figure 1.
Model 11 Python code
The first task is to import the libraries that are needed for our program. As shown in Figure 2, we import the SciPy library, which we've previously installed, along with some other libraries that we'll use.
The data for Model 11 is shown in Figure 3. We're using the json format that we used in some previous models.
As our engine we specify the HiGHS solver, which is a fast, open source linear solver that is packaged with SciPy.
We import the data from the json file using the code shown in Figure 4. This code is the same as the code we used for previous json files, apart from the filename.
As shown in Figure 5, we assign our data to suitable data structures. Unlike most other models in this series, we don't assign the data to a
Model object. Instead, we'll define the objective function and constraints using lists, as expected by SciPy's
linprog function. Even so, the general structure of this code is similar to previous models.
Define the model
Like our CVXPY model, SciPy defines a linear program using of matrix notation. Note that the variables are implicit, rather than being explicitly declared. We will be able to access the variable values later, using the built-in name
We start by defining our parameters as empty numpy arrays (i.e., arrays of zeros, where the array length is the number of products). Then we populate the arrays using the data we loaded from the json file.
Note that SciPy always minimizes, while we want to maximize, so we need to change the sign of our objective function coefficients.
Finally, we create the objective function, constraints' left-hand side, and constraints' right-hand side lists. Note that all of our constraints are in the form ≤, which we'll specify using the
b_ub notation in the next section.
As shown in Figure 7, we specify and solve the model in a single step. We also include options like the solver engine and time limit.
linprog function expects that a linear program is expressed in the form shown in Figure 8.
- Objective function coefficients (transposed) multiplied by the
- Inequality constraints expressed as upper bounds.
- Equality constraints.
- Lower and upper bounds on the variables.
Therefore, we specify the
bounds parameters in our call to the
linprog function. We have no equality constraints in this model, so we omit the
The code for processing the solver result, as shown in Figure 9, is similar to the code for previous models.
linprog returns a Boolean that indicates whether the model solved successfully.
The code for writing the output, as shown in Figure 10, is similar to our previous models. The main difference is the syntax for extracting the variable values, along with the slack and dual values, from the solver.
When we find an optimal solution, the output is as shown in Figure 11. This output is similar to previous models, except for the model's name.
Note that, since SciPy always minimizes, we need to change the sign of the objective function to represent our maximization objective. Similarly, the signs of the slack and dual values are inverted (though we have printed them as reported).
Evaluation of this model
This SciPy model has more in common with our CVXPY model than with our Pyomo models. That is, both CVXPY and SciPy use a matrix notation that reflects the standard mathematical definition of a linear program. Although some mathematical purists prefer matrix notation, from a modelling perspective it can be less intuitive. Even so, despite the syntax of SciPy being markedly different to Pyomo, the general structure of the Python program is familiar.
SciPy is a capable modelling library, with a scope that is much broader than optimization. Consequently, because SciPy is not focussed solely on optimization modelling, it lacks the specific modelling support that makes Pyomo easier to use. Therefore, for a linear programming model, there is no reason to use SciPy compared with Pyomo, though SciPy may be suitable for other types of modelling.
In the next article we'll conclude the Production mix series with a summary of the libraries we've used and our experience of working with the libraries.
In this article we built the Production mix model using the SciPy library. Compared with the Pyomo models, the code is broadly similar though the model specification is markedly different to Pyomo.
For some types of modelling, SciPy would be a capable and appropriate choice. Here we're solving a straightforward linear programming model, for which SciPy offers no significant advantage compared with Pyomo. Therefore, we tend to prefer Pyomo over SciPy for this type of model.
In the next article, we'll summarize and compare the optimization libraries that we've used in this series of articles.
If you would like to know more about this model, or you want help with your own models, then please contact us.