Couenne

Couenne (Convex Over and Under ENvelopes for Nonlinear Estimation) is a branch&bound algorithm to solve Mixed-Integer Nonlinear Programming (MINLP) problems of the form:

      min f₀(x,y)
             f_i(x,y) ≤ 0     i=1,2..., m
             x in Rⁿ, y in Z^p

where all f_i(x,y) are, in general, nonlinear functions.

Couenne aims at finding global optima of nonconvex MINLPs. It implements linearization, bound reduction, and branching methods within a branch-and-bound framework. Its main components are:

• an expression library;
• separation of linearization cuts;
• branching rules;
• bound tightening methods.

It is distributed on COIN-OR under the Eclipse Public License (EPL). The EPL is a license approved by the Open Source Initiative (OSI), thus Couenne is OSI Certified Open Source Software.

Download, installation and usage

Couenne is found on the COIN-OR project page. It can be downloaded with Subversion -- see also some instructions on using svn. Run the command:

svn co https://projects.coin-or.org/svn/Couenne/stable/0.5 Couenne

to get the source code of the stable version. Before building and installing Couenne, some third party packages are needed. These cannot be downloaded from COIN-OR, and have to be obtained independently. These packages are: ASL, Blas, Lapack, and HSL or MUMPS. The user is referred to the INSTALL file in each of the subdirectories of Couenne/ThirdParty for instructions on how to obtain them.

In general, the stable version of Couenne is subject to slight changes such as bug fixes. In order to be up-to-date with such changes, you may run the command svn update within the Couenne directory. All releases are also available as archive at https://www.coin-or.org/download/source/Couenne/. A release is not subject to change.

To install Couenne, we refer to general installation instructions for COIN-OR projects. We also suggest the excellent Ipopt compilation hints page for compiling on non-Linux systems, such as Mac and Windows. The impatient may want to issue the following commands:

cd Couenne
cd ThirdParty      # Read INSTALL.* file in each subdirectory and get third party software
cd ..
mkdir build
cd build
../configure -C
make
make install

The above commands place Couenne in the Couenne/build/bin/ directory, libraries in Couenne/build/lib/, and include files in Couenne/build/include/. An alternative directory can be specified with the --prefix option of configure. For instance, when replacing "../configure -C" above with

../configure -C --prefix=/usr/local

the Couenne executable will be installed in /usr/local/bin/, the libraries in /usr/local/lib/, and the include files in /usr/local/include/. Couenne is run as follows:

couenne instance.nl

where instance.nl is an AMPL stub (.nl) file. Such files can be generated from AMPL with the command "write gfilename;" (notice the "g" before the file name), for example.

You may also specify a set of options to tweak the performance of Couenne. These are found in the couenne.opt option file. A sample option file is given in the Couenne/src/ directory.

Documentation

A user manual is available, with explanations on most options available in Couenne. Doxygen documentation is also available, and it can be generated by running

make doxydoc

from the same build/ directory where you ran configure, make, and make install. Documentation in both html and LaTeX format can be found in the Doc/ subdirectory. Fire up your browser and take a look at Doc/html/index.html for documentation of Couenne.

Resources and links

Couenne is maintained by Pietro Belotti.

Web page: https://www.github.com/coin-or/Couenne

Dependencies: CoinUtils, Cbc, Cgl, Clp, Ipopt, and Osi (from COIN-OR); ASL (the Ampl Solver Library), Lapack, Blas, HSL, MUMPS, SCIP, and SoPlex.

External resources: COIN-OR, Eclipse Public License.

Report a bug, contribute to Couenne

As an open-source code, contributions to Couenne are welcome. To submit a contribution to Couenne, please follow the COIN-OR guidelines.

In order to report a bug, use the issue system.

In order to ensure that your issue is addressed in a timely fashion, please try to be as exhaustive as you can in the bug report, for instance by reporting what version of Couenne you have downloaded and what operating system you are using, and again by attaching the model/data files with which the crash occurred.

Contributors

• Pietro Belotti (Xpress Development Team, FICO, Birmingham UK)
• Timo Berthold (Xpress Development Team, FICO, Berlin)
• Pierre Bonami (IBM Cplex)
• Sonia Cafieri (École Nationale de l'Aviation Civile)
• François Margot (Carnegie Mellon University)
• Cameron Megaw (Clemson University)
• Stefan Vigerske (GAMS)
• Andreas Wächter (Northwestern University)

Acknowledgments

This project was initiated in 2006 within a collaboration between IBM and Carnegie Mellon University, aimed at developing algorithms for MINLP.

Credit should be given to our colleagues in this collaboration: Andreas, François, Pierre, Stefan, and Timo, who developed part of Couenne, and Larry T. Biegler, Gérard Cornuéjols, Ignacio E. Grossmann, and Jon Lee. Each has contributed an essential part of the development of Couenne.

Project Links

• Couenne executables provided by AMPL
• Couenne is available through the JuMP modeling language
• COIN-OR Initiative
• IBM CMU Open Source MINLP project
• Report a bug
• NEOS server for Optimization: You can submit a MINLP problem for solution with Couenne in AMPL format or GAMS format

Options for BonCouenne

Linearization options

convexification_cuts <num>

Specify the frequency (in terms of nodes) at which linearization cuts are generated. Default: 1. If 0, linearization cuts are never separated.

convexification_points <num>

Specify the number of points at which to convexify. Default: 1.

violated_cuts_only <yes|no>

If set to yes (default), only violated convexification cuts will be added.

art_lower <num>

Set artificial lower bound (for minimization problems), useful when a lower bound is known or for testing purposes. Default value is -10⁵⁰.

opt_window <num>

Multiplier for restricting variable bounds around known optimum (to be read from file with method CouenneProblem::readOptimum()). If the optimal value x,,i,, of the i-th variable is known, before starting Couenne its bounds will be intersected with interval [x_i-K(1+|x_i|),x_i+K(1+|x_i|)], where K is the value of the option. Default value is infinity.

use_quadratic <yes|no>

Use quadratic expressions and related exprQuad class. Still in testing, so default is "no".

Bound tightening options

feasibility_bt <yes|no>

Use feasibility-based bound tightening (strongly recommended). Default value is "yes".

optimality_bt

Optimality-based (expensive) bound tightening. Only recommended for problems with few variables and/or at the initial nodes of the B&B tree. Default is "no". If set to "yes", we recommend to couple it with a value of log_num_obbt_per_level of 0 (see below).

log_num_obbt_per_level <num>

Specify the frequency (in terms of nodes) for optimality-based bound tightening. Default is 0.

• If -1, apply at every node (expensive!).
• If 0, apply at root node only.
• If k>0, apply with probability 2^{(k - level)}, level being the current depth of the B&B tree.

aggressive_fbbt <yes|no"

Aggressive feasibility-based bound tightening (to use with NLP points). Default value is "yes". This is also computationally expensive.

log_num_abt_per_level <num>

Specify the frequency (in terms of nodes) for aggressive bound tightening (similar to log_num_obbt_per_level).

• If -1, apply at every node (expensive!);
• If 0, apply at root node only;
• If k>0, apply with probability 2^{(k - level)}, level being the current depth of the B&B tree.

Branching options

branch_fbbt <yes|no>

Apply bound tightening before branching. default: yes

branch_conv_cuts <yes|no>

Apply convexification cuts before branching (not active yet). Default: no.

branch_pt_select <string>

Chooses branching point selection strategy. Possible values are

• "lp-clamped": LP point clamped in [k,1-k] of the bound intervals (k defined by lp_clamp);
• "lp-central": LP point if within [k,1-k] of the bound intervals, middle point otherwise (k defined by branch_lp_clamp);
• "balanced": minimizes max distance from curve to convexification;
• "min-area": minimizes total area of the two convexifications;
• "mid-point": convex combination of current point and mid point;
• "no-branch": do not branch, return null infeasibility; for testing purposes only.

Default is mid-point.

branch_lp_clamp <num>

Defines a threshold for selecting an LP point as the branching point;