An LP matrix generator for GMPL (MathProg) and (by extension) AMPL.

Can be used in two ways:

• Compile an MPS file for solving elsewhere
• Solve directly with the HiGHS integration, no other binary needed

Currently a work-in-progress.

It works for a subset of GMPL (specifically the subset required to run Osemosys). See Known limitations section below.

There are a number of examples of varying complexity in the examples/ directory.

Developed by and for Climate Compatible Growth.

cargo install mosox

brew tap carderne/mosox-tap
brew install mosox

Binaries are built for a small set of systems and architectures. Available to download (compressed) from the Releases page. Please choose the appropriate archive.

Usage overview:

mosox help

MathProg Translation Kit

Usage: mosox <COMMAND>

Commands:
  compile    Load and output to MPS
  solve      Solve with HiGHS
  normalize  Normalize an MPS file for diffing
  compare    Compare two normalized MPS files with epsilon tolerance
  help       Print this message or the help of the given subcommand(s)

Options:
  -h, --help     Print help
  -V, --version  Print version

Compile MPS for one of the examples in this repo. Output will be written to to the provided file:

mosox compile examples/basic/model.mod -o output.mps

Or an example with a separate data file, piping MPS to a file:

mosox compile examples/sets/model.mod examples/sets/data.dat -o output.mps

Results will be saved to the output file:

mosox solve examples/basic/model.mod -o output.txt

# Loading model from model.mod
# Generating matrix
# Matrix compiled in 666.625µs
# Matrix: 3 rows, 2 cols, 6 nonzero
# Solving matrix with HiGHS
# Solved in 1.939083ms
# Objective value: 2200
# Results output to f.txt

Please install cargo-make:

cargo install cargo-make

The most useful dev commands are listed in Makefile.toml.

You can view available commands by running cargo make.

cargo make fmt
cargo make lint

In addition to other tests, this will run mosox against all the examples under examples/ and confirm that the output is identical to the existing MPS files.

cargo make test

This will additionally run a regression test against examples/osemosys_large if present.

cargo make testlarge

Run performance benchmarks across all examples (except osemosys_large):

cargo make bench

Include osemosys_large (requires the example to be present):

cargo make benchlarge

Run with glpsol comparison (requires glpsol on PATH):

cargo make benchglpsol

Note that some of the performance advantage over glpsol may be caused by the limitations enumerated below.

The matrix generator uses very roughly 2,000 bytes per non-zero. This is a significant overhead over the 12 bytes bytes per non-zero that would be needed in a format like CSC. However, as the table below shows, the matrix generator will rarely (if ever) use more memory than the solver.

This list of limitations is made with reference to the GNU MathProg Language Reference which can be viewed here or downloaded from the original here.

It is intended that all of these will ultimately be supported, and most of them are "trivial" to add.

• The following functions: abs, atan, ceil, cos, exp, floor, gmtime, length, log, prod, round, sin, sqrt, str2time, trunc, Irand224, Uniform, Normal.

• These arithmetic operators: less, div, mod
• These symbolic operators: & (string concatenation)
• These set expressions: conditional set expressions
• These set operators: diff, symdiff.
• These logical iterated expressions: forall, exists
• These logical operators: not

• within (parsed, not enforced)
• dimen (parsed, not enforced)

• Relational condition (parsed, not enforced)
• Superset expression (parsed, not enforced)
• Type specifier (integer, binary, symbolic) (parsed, not enforced)

• Bounds specified as expressions (currently only constant accepted)

• Multiple expressions (comma-separated)

Support for these statements is not planned (results should instead be parsed from the solver output): display, printf, for.