Gams Users Guide PDF
Gams Users Guide PDF
c 2007 GAMS Development Corporation, Washington, DC, USA
Table of Contents
1 Introduction 1
1.1 Motivation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Basic Features of GAMS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2.1 General Principles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2.2 Documentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.3 Portability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.4 User Interface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.5 Model Library . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Organization of the Book . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
3 GAMS Programs 27
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2 The Structure of GAMS Programs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2.1 Format of GAMS Input . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2.2 Classification of GAMS Statements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
4 TABLE OF CONTENTS
4 Set Definitions 35
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2 Simple Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2.1 The Syntax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2.2 Set Names . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.2.3 Set Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2.4 Associated Text . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
4.2.5 Sequences as Set Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.2.6 Declarations for Multiple Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.3 The Alias Statement: Multiple Names for a Set . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4.4 Subsets and Domain Checking . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4.5 Multi-dimensional Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.5.1 One-to-one Mapping . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.5.2 Many-to-many Mapping . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.6 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
7 Variables 59
7.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.2 Variable Declarations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.2.1 The Syntax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.2.2 Variable Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
7.2.3 Styles for Variable Declaration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
7.3 Variable Attributes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
7.3.1 Bounds on Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
7.3.2 Fixing Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
7.3.3 Activity Levels of Variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.4 Variables in Display and Assignment Statements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.4.1 Assigning Values to Variable Attributes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.4.2 Variable Attributes in Assignments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
7.4.3 Displaying Variable Attributes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
7.5 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
8 Equations 65
8.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
8.2 Equation Declarations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
8.2.1 The Syntax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
8.2.2 An Illustrative Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
8.3 Equation Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
8.3.1 The Syntax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
8.3.2 An Illustrative Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
8.3.3 Scalar Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
8.3.4 Indexed Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
8.3.5 Using Labels Explicitly in Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
8.4 Expressions in Equation Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
8.4.1 Arithmetic Operators in Equation Definitions . . . . . . . . . . . . . . . . . . . . . . . . . 68
8.4.2 Functions in Equation Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
8.4.3 Preventing Undefined Operations in Equations . . . . . . . . . . . . . . . . . . . . . . . . 68
8.5 Data Handling Aspects of Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
10 GAMS Output 79
10.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
10.2 The Illustrative Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
10.3 Compilation Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
10.3.1 Echo Print of the Input File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
10.3.2 The Symbol Reference Map . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
10.3.3 The Symbol Listing Map . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
10.3.4 The Unique Element Listing - Map . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
10.3.5 Useful Dollar Control Directives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
10.4 Execution Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
10.5 Output Produced by a Solve Statement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
10.5.1 The Equation Listing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
10.5.2 The Column Listing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
10.5.3 The Model Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
10.5.4 The Solve Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
10.5.5 Solver Report . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
10.5.6 The Solution Listing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
10.5.7 Report Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
10.5.8 File Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
10.6 Error Reporting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
10.6.1 Compilation Errors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
10.6.2 Compilation Time Errors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
10.6.3 Execution Errors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10.6.4 Solve Errors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10.7 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
Appendix 156
A Glossary 157
Index 271
List of Tables
2.1 Data for the transportation problem (adapted from Dantzig, 1963) . . . . . . . . . . . . . . . . . . 6
2.2 The basic components of a GAMS model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.3 Permissible variable types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Introduction
1.1 Motivation
Substantial progress was made in the 1950s and 1960s with the development of algorithms and computer codes
to solve large mathematical programming problems. The number of applications of these tools in the 1970s was
less then expected, however, because the solution procedures formed only a small part of the overall modeling
effort. A large part of the time required to develop a model involved data preparation and transformation and
report preparation. Each model required many hours of analyst and programming time to organize the data
and write the programs that would transform the data into the form required by the mathematical programming
optimizers. Furthermore, it was difficult to detect and eliminate errors because the programs that performed the
data operations were only accessible to the specialist who wrote them and not to the analysts in charge of the
project.
GAMS was developed to improve on this situation by:
Providing a high-level language for the compact representation of large and complex models
Allowing changes to be made in model specifications simply and safely
Allowing unambiguous statements of algebraic relationships
Permitting model descriptions that are independent of solution algorithms
The design of GAMS has incorporated ideas drawn from relational database theory and mathematical program-
ming and has attempted to merge these ideas to suit the needs of strategic modelers. Relational database theory
provides a structured framework for developing general data organization and transformation capabilities. Math-
ematical programming provides a way of describing a problem and a variety of methods for solving it. The
following principles were used in designing the system:
1. All existing algorithmic methods should be available without changing the users model representation.
Introduction of new methods, or of new implementations of existing methods, should be possible
without requiring changes in existing models. Linear, nonlinear, mixed integer, mixed integer nonlinear
optimizations and mixed complementarity problems can currently be accommodated.
2. The optimization problem should be expressible independently of the data it uses. This separation of
logic and data allows a problem to be increased in size without causing an increase in the complexity of
the representation.
2 Introduction
3. The use of the relational data model requires that the allocation of computer resources be automated.
This means that large and complex models can be constructed without the user having to worry about
details such as array sizes and scratch storage.
1.2.2 Documentation
The GAMS model representation is in a form that can be easily read by people and by computers. This means
that the GAMS program itself is the documentation of the model, and that the separate description required in
the past (which was a burden to maintain, and which was seldom up-to-date) is no longer needed. Moreover, the
design of GAMS incorporates the following features that specifically address the users documentation needs:
A GAMS model representation is concise, and makes full use of the elegance of the mathematical
representation.
All data transformations are specified concisely and algebraically. This means that all data can be
entered in their most elemental form and that all transformations made in constructing the model and
in reporting are available for inspection.
Explanatory text can be made part of the definition of all symbols and is reproduced whenever
associated values are displayed.
All information needed to understand the model is in one document.
Of course some discipline is needed to take full advantage of these design features, but the aim is to make models
more accessible, more understandable, more verifiable, and hence more credible.
1.2.3 Portability
The GAMS system is designed so that models can be solved on different types of computers with no change. A
model developed on a small personal computer can later be solved on a large mainframe. One person can develop
a model that is later used by others, who may be physically distant from the original developer. In contrast to
previous approaches, only one document need be moved the GAMS statement of the model. It contains all
the data and logical specifications needed to solve the model.
2.1 Introduction
The introductory part of this book ends with a detailed example of the use of GAMS for formulating, solving,
and analyzing a small and simple optimization problem. Richard E. Rosenthal of the Naval Postgraduate School
in Monterey, California wrote it. The example is a quick but complete overview of GAMS and its features. Many
references are made to other parts of the book, but they are only to tell you where to look for more details; the
material here can be read profitably without reference to the rest of the book.
The example is an instance of the transportation problem of linear programming, which has historically served
as a laboratory animal in the development of optimization technology. [See, for example, Dantzig (1963)1 .] It is
a good choice for illustrating the power of algebraic modeling languages like GAMS because the transportation
problem, no matter how large the instance at hand, possesses a simple, exploitable algebraic structure. You will
see that almost all of the statements in the GAMS input file we are about to present would remain unchanged if
a much larger transportation problem were considered.
In the familiar transportation problem, we are given the supplies at several plants and the demands at several
markets for a single commodity, and we are given the unit costs of shipping the commodity from plants to
markets. The economic question is: how much shipment should there be between each plant and each market so
as to minimize total transport cost?
The algebraic representation of this problem is usually presented in a format similar to the following.
Indices:
i = plants
j = markets
Given Data:
ai = supply of commodity of plant i (in cases)
bi = demand for commodity at market j
cij = cost per unit shipment between plant i and market j ($/case)
Decision Variables:
xij = amount of commodity to ship from plant i to market j (cases),
where xij 0, for all i, j
Constraints: P
Observe supply limit at plant i: x aj for all i (cases)
Pj ij
Satisfy demand at market j: i ij bP
x j P for all j (cases)
Objective Function: Minimize i j cij xij ($K)
Note that this simple example reveals some modeling practices that we regard as good habits in general and
that are consistent with the design of GAMS. First, all the entities of the model are identified (and grouped) by
type. Second, the ordering of entities is chosen so that no symbol is referred to before it is defined. Third, the
1 Dantzig, George B. (1963). Linear Programming and Extensions, Princeton University Press, Princeton N.J.
6 A GAMS Tutorial by Richard E. Rosenthal
units of all entities are specified, and, fourth, the units are chosen to a scale such that the numerical values to be
encountered by the optimizer have relatively small absolute orders of magnitude. (The symbol $K here means
thousands of dollars.)
The names of the types of entities may differ among modelers. For example, economists use the terms exogenous
variable and endogenous variable for given data and decision variable, respectively. In GAMS, the terminology
adopted is as follows: indices are called sets, given data are called parameters, decision variables are called
variables, and constraints and the objective function are called equations.
The GAMS representation of the transportation problem closely resembles the algebraic representation above.
The most important difference, however, is that the GAMS version can be read and processed by the computer.
Table 2.1: Data for the transportation problem (adapted from Dantzig, 1963)
As an instance of the transportation problem, suppose there are two canning plants and three markets, with the
data given in table 2.1. Shipping distances are in thousands of miles, and shipping costs are assumed to be $90.00
per case per thousand miles. The GAMS representation of this problem is as follows:
Sets
i canning plants / seattle, san-diego /
j markets / new-york, chicago, topeka / ;
Parameters
Variables
x(i,j) shipment quantities in cases
z total transportation costs in thousands of dollars ;
Positive Variable x ;
Equations
cost define objective function
supply(i) observe supply limit at plant i
demand(j) satisfy demand at market j ;
If you submit a file containing the statements above as input to the GAMS program, the transportation model
will be formulated and solved. Details vary on how to invoke GAMS on different of computers, but the simplest
(no frills) way to call GAMS is to enter the word GAMS followed by the input files name. You will see a number
of terse lines describing the progress GAMS is making, including the name of the file onto which the output is
being written. When GAMS has finished, examine this file, and if all has gone well the optimal shipments will
be displayed at the bottom as follows.
You will also receive the marginal costs (simplex multipliers) below.
chicago topeka
seattle 0.036
san-diego 0.009
These results indicate, for example, that it is optimal to send nothing from Seattle to Topeka, but if you insist
on sending one case it will add .036 $K (or $36.00) to the optimal cost. (Can you prove that this figure is correct
from the optimal shipments and the given data?)
2. GAMS statements may be laid out typographically in almost any style that is appealing to the user.
Multiple lines per statement, embedded blank lines, and multiple statements per line are allowed. You
will get a good idea of what is allowed from the examples in this tutorial, but precise rules of the road
are given in the next Chapter.
3. When you are a beginning GAMS user, you should terminate every statement with a semicolon, as in
our examples. The GAMS compiler does not distinguish between upper-and lowercase letters, so you
are free to use either.
4. Documentation is crucial to the usefulness of mathematical models. It is more useful (and most likely
to be accurate) if it is embedded within the model itself rather than written up separately. There are
at least two ways to insert documentation within a GAMS model. First, any line that starts with an
asterisk in column 1 is disregarded as a comment line by the GAMS compiler. Second, perhaps more
important, documentary text can be inserted within specific GAMS statements. All the lowercase
words in the transportation model are examples of the second form of documentation.
8 A GAMS Tutorial by Richard E. Rosenthal
Inputs: Outputs:
Sets Echo Print
Declaration
Reference Maps
Assignment of members
Equation Listings
Data (Parameters, Tables, Scalars) Status Reports
Declaration
Assignment of values Results
Variables
Declaration
Assignment of type
5. As you can see from the list of input components above, the creation of GAMS entities involves two
steps: a declaration and an assignment or definition. Declaration means declaring the existence of
something and giving it a name. Assignment or definition means giving something a specific value or
form. In the case of equations, you must make the declaration and definition in separate GAMS
statements. For all other GAMS entities, however, you have the option of making declarations and
assignments in the same statement or separately.
6. The names given to the entities of the model must start with a letter and can be followed by up to
thirty more letters or digits.
2.3 Sets
Sets are the basic building blocks of a GAMS model, corresponding exactly to the indices in the algebraic
representations of models. The Transportation example above contains just one Set statement:
Sets
i canning plants / seattle, san-diego /
j markets / new-york, chicago, topeka / ;
The effect of this statement is probably self-evident. We declared two sets and gave them the names i and j. We
also assigned members to the sets as follows:
i = {Seattle, San Diego}
j = {New York, Chicago, Topeka}.
You should note the typographical differences between the GAMS format and the usual mathematical format
for listing the elements of a set. GAMS uses slashes / rather than curly braces {} to delineate the set simply
2.4 Data 9
because not all computer keyboards have keys for curly braces. Note also that multiword names like New York
are not allowed, so hyphens are inserted.
The lowercase words in the sets statement above are called text. Text is optional. It is there only for internal
documentation, serving no formal purpose in the model. The GAMS compiler makes no attempt to interpret the
text, but it saves the text and parrots it back to you at various times for your convenience.
It was not necessary to combine the creation of sets i and j in one statement. We could have put them into
separate statements as follows:
The placement of blank spaces and lines (as well as the choice of upper- or lowercase) is up to you. Each GAMS
user tends to develop individual stylistic conventions. (The use of the singular set is also up to you. Using set
in a statement that makes a single declaration and sets in one that makes several is good English, but GAMS
treats the singular and plural synonymously.)
A convenient feature to use when you are assigning members to a set is the asterisk. It applies to cases when the
elements follow a sequence. For example, the following are valid set statements in GAMS.
Alias (t,tp);
the name tp is like a t0 in mathematical notation. It is useful in models that are concerned with the interactions
of elements within the same set.
The sets i, j, t, and m in the statements above are examples of static sets, i.e., they are assigned their members
directly by the user and do not change. GAMS has several capabilities for creating dynamic sets, which acquire
their members through the execution of set-theoretic and logical operations. Dynamic sets are discussed in Chap-
ter 12, page107. Another valuable advanced feature is multidimensional sets, which are discussed in Section 4.5,
page 39.
2.4 Data
The GAMS model of the transportation problem demonstrates all of the three fundamentally different formats
that are allowable for entering data. The three formats are:
Lists
Tables
Direct assignments
The next three sub-sections will discuss each of these formats in turn.
10 A GAMS Tutorial by Richard E. Rosenthal
The first format is illustrated by the first Parameters statement of the example, which is repeated below.
Parameters
This statement has several effects. Again, they may be self-evident, but it is worthwhile to analyze them in
detail. The statement declares the existence of two parameters, gives them the names a and b, and declares their
domains to be i and j, respectively. (A domain is the set, or tuple of sets, over which a parameter, variable, or
equation is defined.) The statement also gives documentary text for each parameter and assigns values of a(i)
and b(j) for each element of i and j. It is perfectly acceptable to break this one statement into two, if you
prefer, as follows.
Here are some points to remember when using the list format.
1. The list of domain elements and their respective parameter values can be laid out in almost any way
you like. The only rules are that the entire list must be enclosed in slashes and that the element-value
pairs must be separated by commas or entered on separate lines.
2. There is no semicolon separating the element-value list from the name, domain, and text that precede
it. This is because the same statement is being used for declaration and assignment when you use the
list format. (An element-value list by itself is not interpretable by GAMS and will result in an error
message.)
3. The GAMS compiler has an unusual feature called domain checking, which verifies that each domain
element in the list is in fact a member of the appropriate set. For example, if you were to spell Seattle
correctly in the statement declaring Set i but misspell it as Seatle in a subsequent element-value list,
the GAMS compiler would give you an error message that the element Seatle does not belong to the
set i.
4. Zero is the default value for all parameters. Therefore, you only need to include the nonzero entries in
the element-value list, and these can be entered in any order .
5. A scalar is regarded as a parameter that has no domain. It can be declared and assigned with a Scalar
statement containing a degenerate list of only one value, as in the following statement from the
transportation model.
Scalar f freight in dollars per case per thousand miles /90/;
If a parameters domain has two or more dimensions, it can still have its values entered by the list format. This
is very useful for entering arrays that are sparse (having few non-zeros) and super-sparse (having few distinct
non-zeros).
2.4 Data 11
The effect of this statement is to declare the parameter d and to specify its domain as the set of ordered pairs in
the Cartesian product of i and j. The values of d are also given in this statement under the appropriate heading.
If there are blank entries in the table, they are interpreted as zeroes.
As in the list format, GAMS will perform domain checking to make sure that the row and column names of the
table are members of the appropriate sets. Formats for entering tables with more columns than you can fit on
one line and for entering tables with more than two dimensions are given in Chapter 5, page 43.
It is important to emphasize the presence of the semicolon at the end of the first line. Without it, the GAMS
compiler would attempt to interpret both lines as parts of the same statement. (GAMS would fail to discern a
valid interpretation, so it would send you a terse but helpful error message.)
The effects of the first statement above are to declare the parameter c, to specify the domain (i,j), and to provide
some documentary text. The second statement assigns to c(i,j) the product of the values of the parameters
f and d(i,j). Naturally, this is legal in GAMS only if you have already assigned values to f and d(i,j) in
previous statements.
The direct assignment above applies to all (i,j) pairs in the domain of c. If you wish to make assignments for
specific elements in the domain, you enclose the element names in quotes. For example,
c(Seattle,New-York) = 0.40;
csquared = sqr(c);
12 A GAMS Tutorial by Richard E. Rosenthal
e = m*csquared;
w = l/lamda;
eoq(i) = sqrt( 2*demand(i)*ordcost(i)/holdcost(i));
t(i) = min(p(i), q(i)/r(i), log(s(i)));
euclidean(i,j) = qrt(sqr(xi(i) - xi(j) + sqr(x2(i) - x2(j)));
present(j) = future(j)*exp(-interest*time(j));
The summation and product operators to be introduced later can also be used in direct assignments.
2.5 Variables
The decision variables (or endogenous variables ) of a GAMS-expressed model must be declared with a Variables
statement. Each variable is given a name, a domain if appropriate, and (optionally) text. The transportation
model contains the following example of a Variables statement.
Variables
x(i,j) shipment quantities in cases
z total transportation costs in thousands of dollars ;
This statement results in the declaration of a shipment variable for each (i,j) pair. (You will see in Chapter 8,
page 65, how GAMS can handle the typical real-world situation in which only a subset of the (i,j) pairs is
allowable for shipment.)
The z variable is declared without a domain because it is a scalar quantity. Every GAMS optimization model
must contain one such variable to serve as the quantity to be minimized or maximized.
Once declared, every variable must be assigned a type. The permissible types are given in table 2.3.
The variable that serves as the quantity to be optimized must be a scalar and must be of the free type. In our
transportation example, z is kept free by default, but x(i,j) is constrained to non-negativity by the following
statement.
Positive variable x ;
Note that the domain of x should not be repeated in the type assignment. All entries in the domain automatically
have the same variable type.
Section 2.10 describes how to assign lower bounds, upper bounds, and initial values to variables.
2.6 Equations
The power of algebraic modeling languages like GAMS is most apparent in the creation of the equations and
inequalities that comprise the model under construction. This is because whenever a group of equations or
inequalities has the same algebraic structure, all the members of the group are created simultaneously, not
individually.
2.6 Equations 13
Equations must be declared and defined in separate statements. The format of the declaration is the same as
for other GAMS entities. First comes the keyword, Equations in this case, followed by the name, domain and
text of one or more groups of equations or inequalities being declared. Our transportation model contains the
following equation declaration:
Equations
cost define objective function
supply(i) observe supply limit at plant i
demand(j) satisfy demand at market j ;
Keep in mind that the word Equation has a broad meaning in GAMS. It encompasses both equality and inequality
relationships, and a GAMS equation with a single name can refer to one or several of these relationships. For
example, cost has no domain so it is a single equation, but supply refers to a set of inequalities defined over the
domain i.
Before going into equation definition we describe the summation notation in GAMS. Remember that GAMS is
designed for standard keyboards and line-by-line input readers, so it is not possible (nor would it be convenient
for the user) to employ the standard mathematical notation for summations.
The summation notation in GAMS can be used for simple and complex expressions. The format is based on
the idea of always thinking of a summation as an operator with two arguments: Sum(index of summation,
summand) A comma separates the two arguments, and if the first argument requires a comma then it should be
in parentheses. The second argument can be any mathematical expression including another summation.
As a simple example, the transportation problem contains the expression
Sum(j, x(i,j))
P
that is equivalent to j xij .
A slightly more complex summation is used in the following example:
Sum((i,j), c(i,j)*x(i,j))
P P
that is equivalent to i j cij xij .
The last expression could also have been written as a nested summation as follows:
In Section 11.3, page 100, we describe how to use the dollar operator to impose restrictions on the summation
operator so that only the elements of i and j that satisfy specified conditions are included in the summation.
Products are defined in GAMS using exactly the same format as summations, replacing Sum by Prod. For example,
2. The domain
4. The symbol ..
5. Left-hand-side expression
7. Right-hand-side expression
The key idea here is that the definition of the demand constraints is exactly the same whether we are
solving the toy-sized example above or a 20,000-node real-world problem. In either case, the user enters
just one generic equation algebraically, and GAMS creates the specific equations that are appropriate
for the model instance at hand. (Using some other optimization packages, something like the extract
above would be part of the input, not the output.)
In many real-world problems, some of the members of an equation domain need to be omitted or
differentiated from the pattern of the others because of an exception of some kind. GAMS can readily
accommodate this loss of structure using a powerful feature known as the dollar or such-that
operator, which is not illustrated here. The domain restriction feature can be absolutely essential for
keeping the size of a real-world model within the range of solvability.
The relational operators have the following meanings:
=l= less than or equal to
=g= greater than or equal to
=e= equal to
It is important to understand the difference between the symbols = and =e=. The = symbol is used
only in direct assignments, and the =e= symbol is used only in equation definitions. These two
contexts are very different. A direct assignment gives a desired value to a parameter before the solver
is called. An equation definition also describes a desired relationship, but it cannot be satisfied until
after the solver is called. It follows that equation definitions must contain variables and direct
assignments must not.
2.7 Objective Function 15
Variables can appear on the left or right-hand side of an equation or both. The same variable can
appear in an equation more than once. The GAMS processor will automatically convert the equation to
its equivalent standard form (variables on the left, no duplicate appearances) before calling the solver.
An equation definition can appear anywhere in the GAMS input, provided the equation and all
variables and parameters to which it refers are previously declared. (Note that it is permissible for a
parameter appearing in the equation to be assigned or reassigned a value after the definition. This is
useful when doing multiple model runs with one GAMS input.) The equations need not be defined in
the same order in which they are declared.
This statement may seem superfluous, but it is useful to advanced users who may create several models in one
GAMS run. If we were to use the explicit list rather than the shortcut /all/, the statement would be written as
The domains are omitted from the list since they are not part of the equation name. The list option is used when
only a subset of the existing equations comprises a specific model (or sub-model) being generated.
Once a model has been declared and assigned equations, we are ready to call the solver. This is done with a solve
statement, which in our example is written as
that calls for a printout of the final levels, x.l, and marginal (or reduced costs), x.m, of the shipment variables,
x(i,j). GAMS will automatically format this printout in to dimensional tables with appropriate headings.
The format for referencing these quantities is the variable or equations name followed by the fields name, followed
(if necessary) by the domain (or an element of the domain).
GAMS allows the user complete read-and write-access to the database. This may not seem remarkable to you
now, but it can become a greatly appreciated feature in advanced use. Some examples of use of the database
follow.
x.up(i,j) = capacity(i,j) ;
x.lo(i,j) = 10.0 ;
x.up(seattle,new-york) = 1.2*capacity(seattle,new-york) ;
2.10 The .lo, .l, .up, .m Database 17
It is assumed in the first and third examples that capacity(i,j) is a parameter that was previously declared and
assigned values. These statements must appear after the variable declaration and before the Solve statement.
All the mathematical expressions available for direct assignments are usable on the right-hand side.
In nonlinear programming it is very important for the modeler to help the solver by specifying as narrow a range
as possible between lower and upper bound. It is also very helpful to specify an initial solution from which
the solver can start searching for the optimum. For example, in a constrained inventory model, the variables
are quantity(i), and it is known that the optimal solution to the unconstrained version of the problem is a
parameter called eoq(i). As a guess for the optimum of the constrained problem we enter
quantity.l(i) = 0.5*eoq(i) ;
(The default initial level is zero unless zero is not within the bounded range, in which case it is the bound closest
to zero.)
It is important to understand that the .lo and .up fields are entirely under the control of the GAMS user. The
.l and .m fields, in contrast, can be initialized by the user but are then controlled by the solver.
Appending these commands to the original transportation problem input results in the following output:
For an example involving marginal, we briefly consider the ratio constraints that commonly appear in blending
and refining problems. These linear programming models are concerned with determining the optimal amount of
each of several available raw materials to put into each of several desired finished products. Let y(i,j) be the
variable for the number of tons of raw material i put into finished product j. Suppose the ratio constraint is that
no product can consist of more than 25 percent of one ingredient, that is,
Unfortunately, this relaxed constraint has no realistic significance. The constraint we are interested in relaxing
(or tightening) is the nonlinear form of the ration constraint. For example, we would like to know the marginal
benefit arising from changing the ratio constraint to
We can in fact obtain the desired marginals by entering the following transformation on the undesired marginals:
Notice that the assignment statement for amr accesses both .m and .l records from the database. The idea behind
the transformation is to notice that
is equivalent to
A great deal of unnecessary anxiety has been caused by textbooks and users manuals that give the reader the false
impression that flawless use of advanced software should be easy for anyone with a positive pulse rate. GAMS
is designed with the understanding that even the most experienced users will make errors. GAMS attempts to
catch the errors as soon as possible and to minimize their consequences.
3 Sets
4 i canning plants / seattle, san-diego /
5 j markets / new-york, chicago, topeka / ;
6
7 Parameters
8
9 a(i) capacity of plant i in cases
10 / seattle 350
11 san-diego 600 /
12
13 b(j) demand at market j in cases
14 / new-york 325
15 chicago 300
16 topeka 275 / ;
17
18 Table d(i,j) distance in thousands of miles
19 new-york chicago topeka
20 seattle 2.5 1.7 1.8
2.11 GAMS Output 19
The reason this echo print starts with line number 3 rather than line number 1 is because the input file contains
two dollar-print-control statements. This type of instruction controls the output printing, but since it has nothing
to do with defining the optimization model, it is omitted from the echo. The dollar print controls must start in
column 1.
The $title statement causes the subsequent text to be printed at the top of each page of output. The $offupper
statement is needed for the echo to contain mixed upper- and lowercase. Other available instructions are given
in Appendix D, page 193.
In this case, the GAMS compiler indicates that something is wrong with the set element sum. At the bottom of
the echo print, we see the interpretation of error code 160:
20 A GAMS Tutorial by Richard E. Rosenthal
Error Message
160 UNIQUE ELEMENT EXPECTED
The problem is that sum is a reserved word denoting summation, so our set element must have a unique name like
summer. This is a common beginners error. The complete list of reserved words is shown in the next chapter.
Example 2: Another common error is the omission of a semicolon preceding a direct assignment or equation
definition. In our transportation example, suppose we omit the semicolon prior to the assignment of c(i,j), as
follows.
It is not uncommon for one little offense like our missing semicolon to generate five intimidating error messages.
The lesson here is: concentrate on fixing the first error and ignore the other! The first error detected (in line 17),
code 97, indicate that GAMS thinks the symbols in line 17 are a continuation of the documentary text at the
end of line 16 rather than a direct assignment as we intended. The error message also appropriately advises us
to check the preceding line for a missing semicolon.
Unfortunately, you cannot always expect error messages to be so accurate in their advice. The compiler cannot
read your mind. It will at times fail to comprehend your intentions, so learn to detect the causes of errors by
picking up the clues that abound in the GAMS output. For example, the missing semicolon could have been
detected by looking up the c entry in the cross-reference list (to be explained in the next section) and noticing
that it was never assigned.
Example 3: Many errors are caused merely by spelling mistakes and are caught before they can be damaging.
For example, with Seattle spelled in the table differently from the way it was introduced in the set declaration,
we get the following error message.
4 sets
5 i canning plants /seattle, san-diego /
6 j markets /new-york, chicago, topeka / ;
7
8 table d(i,j) distance in thousand of miles
9 new-york chicago topeka
10 seatle 2.5 1.7 1.8
**** $170
11 san-diego 2.5 1.8 1.4 ;
Error Message
170 DOMAIN VIOLATION FOR ELEMENT
Example 4: Similarly, if we mistakenly enter dem(j) instead of b(j) as the right-hand side of the demand
constraint, the result is
2.11 GAMS Output 21
Example 5: The next example is a mathematical error, which is sometimes committed by novice modelers and
which GAMS is adept at catching. The following is mathematically inconsistent and, hence, is not an interpretable
statement. X
For all i, xij = 100
i
There are two errors in this equation, both having to do with the control of indices. Index i is over-controlled
and index j is under-controlled.
You should see that index i is getting conflicting orders. By appearing in the quantifier for all i, it is supposed
to remain fixed for each instance of the equation. Yet, by appearing as an index of summation, it is supposed to
vary. It cant do both. On the other hand, index j is not controlled in any way, so we have no way of knowing
which of its possible values to use.
If we enter this meaningless equation into GAMS, both errors are correctly diagnosed.
A great deal more information about error reporting is given in Section 10.6, page 92. Comprehensive error
detection and well-designed error messages are a big help in getting models implemented quickly and correctly.
For example, the cross-reference list tells us that the symbol A is a parameter that was declared in line 10, defined
(assigned value) in line 11, and referenced in line 43. The symbol I has a more complicated entry in the cross-
reference list. It is shown to be a set that was declared and defined in line 5. It is referenced once in lines 10, 19,
26, 28, 31, 38, 45 and referenced twice in lines 41 and 43. Set I is also used as a controlling index in a summation,
equation definition or direct parameter assignment in lines 28, 41, 43 and 45.
For the GAMS novice, the detailed analysis of the cross-reference list may not be important. Perhaps the most
likely benefit he or she will get from the reference maps will be the discovery of an unwanted entity that mistakenly
entered the model owing to a punctuation or syntax error.
The second part of the reference map is a list of model entities grouped by type and listed with their associated
documentary text. For example, this list is as follows.
sets
i canning plants
j markets
parameters
a capacity of plant i in cases
b demand at market j in cases
c transport cost in 1000s of dollars per case
d distance in thousands of miles
f freight in dollars per case per thousand miles
variables
x shipment quantities in cases
z total transportation costs in 1000s of dollars
equations
cost define objective function
demand satisfy demand at market j
supply observe supply limit at plant i
models
transport
The default output is a maximum of three specific equations for each generic equation. To change the default,
insert an input statement prior to the solve statement:
option limrow = r ;
2.11 GAMS Output 23
where c is the desired number of columns. (Setting limrow and limcol to 0 is a good way to save paper after
your model has been debugged.)
In nonlinear models, the GAMS equation listing shows first-order Taylor approximations of the nonlinear equa-
tions. The approximations are taken at the starting values of the variables.
MODEL STATISTICS
The BLOCK counts refer to the number of generic equations and variables. The SINGLE counts refer to individual
rows and columns in the specific model instance being generated. For nonlinear models, some other statistics are
given to describe the degree of non-linearity in the problem.
S O L V E S U M M A R Y
The status reports are preceded by the same **** string as an error message, so you should probably develop
the habit of searching for all occurrences of this string whenever you look at an output file for the first time.
The desired solver status is 1 NORMAL COMPLETION, but there are other possibilities, documented in Section 10.5,
page 85, which relate to various types of errors and mishaps.
There are eleven possible model statuss, including the usual linear programming termination states (1 OPTIMAL,
3 UNBOUNDED, 4 INFEASIBLE), and others relating to nonlinear and integer programming. In nonlinear pro-
gramming, the status to look for is 2 LOCALLY OPTIMAL. The most the software can guarantee for nonlinear
programming is a local optimum. The user is responsible for analyzing the convexity of the problem to determine
whether local optimality is sufficient for global optimality.
24 A GAMS Tutorial by Richard E. Rosenthal
In integer programming, the status to look for is 8 INTEGER SOLUTION. This means that a feasible integer solution
has been found. More detail follows as to whether the solution meets the relative and absolute optimality tolerances
that the user specifies.
The single dots . in the output represent zeroes. The entry EPS, which stands for epsilon, mean very small but
nonzero. In this case, EPS indicates degeneracy. (The slack variable for the Seattle supply constraint is in the
basis at zero level. The marginal is marked with EPS rather than zero to facilitate restarting the optimizer from
the old basis.)
If the solvers results contain either infeasibilities or marginal costs of the wrong sign, then the offending entries
are marked with INFES or NOPT, respectively. If the problem terminates unbounded, then the rows and columns
corresponding to extreme rays are marked UNBND.
At the end of the solvers solution report is a very important report summary, which gives a tally of the total
number of non-optimal, infeasible, and unbounded rows and columns. For our example, the report summary
shows all zero tallies as desired.
After the solvers report is written, control is returned from the solver back to GAMS. All the levels and marginals
obtained by the solver are entered into the GAMS database in the .l and .m fields. These values can then be
transformed and displayed in any desired report. As noted earlier, the user merely lists the quantities to be
displayed, and GAMS automatically formats and labels an appropriate array. For example, the input statement.
2.12 Summary 25
chicago topeka
seattle 0.036
san-diego 0.009
As seen in reference maps, equation listings, solution reports, and optional displays, GAMS saves the documentary
text and parrots it back throughout the output to help keep the model well documented.
2.12 Summary
This tutorial has demonstrated several of the design features of GAMS that enable you to build practical opti-
mization models quickly and effectively. The following discussion summarizes the advantages of using an algebraic
modeling language such as GAMS versus a matrix generator or conversational solver.
By using an algebra-based notation, you can describe an optimization model to a computer nearly as
easily as you can describe it to another mathematically trained person.
Because an algebraic description of a problem has generality, most of the statements in a GAMS model
are reusable when new instances of the same or related problems arise. This is especially important in
environments where models are constantly changing.
You save time and reduce generation errors by creating whole sets of closely related constraints in one
statement.
You can save time and reduce input errors by providing formulae for calculating the data rather than
entering them explicitly.
The model is self-documenting. Since the tasks of model development and model documentation can
be done simultaneously, the modeler is much more likely to be conscientious about keeping the
documentation accurate and up to date.
The output of GAMS is easy to read and use. The solution report from the solver is automatically
reformatted so that related equations and variables are grouped together and appropriately labeled.
Also, the display command allows you to modify and tabulate results very easily.
If you are teaching or learning modeling, you can benefit from the insistence of the GAMS compiler
that every equation be mathematically consistent. Even if you are an experienced modeler, the
hundreds of ways in which errors are detected should greatly reduce development time.
By using the dollar operator and other advanced features not covered in this tutorial, one can efficiently
implement large-scale models. Specific applications of the dollar operator include the following:
1. It can enforce logical restrictions on the allowable combinations of indices for the variables and
equations to be included in the model. You can thereby screen out unnecessary rows and
columns and keep the size of the problem within the range of solvability.
2. It can be used to build complex summations and products, which can then be used in
equations or customized reports.
3. It can be used for issuing warning messages or for terminating prematurely conditioned upon
context-specific data edits.
26 A GAMS Tutorial by Richard E. Rosenthal
3
GAMS Programs
3.1 Introduction
This chapter provides a look at the structure of the GAMS language and its components. It should be emphasized
again that GAMS is a programming language, and that programs must be written in the language to use it. A
GAMS program is contained in a disk file, which is normally constructed with a text editor of choice. When
GAMS is run, the file containing the program (the input file) is submitted to be processed. After this processing
has finished the results, which are in the output file(s), can be inspected with a text editor. On many machines
a few terse lines appear on the screen while GAMS runs, keeping the user informed about progress and error
detection. But it is the responsibility of the user to inspect the output file carefully to see the results and to
diagnose any errors.
The first time or casual reader can skip this chapter: the discussion of specific parts of the language in the next
Chapters does not assume an understanding of this chapter.
statement;
statement;
statement; statement; statement;
the words that you are now reading is an example of a very
long statement which is stretched over two lines;
Blanks and end-of-lines can generally be used freely between individual symbols or words. GAMS is not case
sensitive, meaning that lower and upper case letters can be mixed freely but are treated identically. Up to
255 characters can be placed on a line and completely blank lines can be inserted for easier reading.
Not all lines are a part of the GAMS language. Two special symbols, the asterisk * and the dollar symbol $
can be used in the first position on a line to indicate a non-language input line. An asterisk in column one means
28 GAMS Programs
that the line will not be processed, but treated as a comment. A dollar symbol in the same position indicates
that compiler options are contained in the rest of the line.
Multiple files can be used as input through the use of the $include facility described in Appendix D. In short,
the statement,
$include file1
inserts the contents of the specified file (file1 in this case) at the location of the call. A more complex versions of
this is the $batinclude which is described in Appendix D.
A declaration statement describes the class of symbol. Often initial values are provided in a declaration, and
then it may be called a definition. The specification of symbolic relationships for an equation is a definition. The
declaration and definition statements are:
acronym parameter equation declaration
set scalar equation definition
alias table model
variable
Execution statements are instructions to carry out actions such as data transformation, model solution, and
report generation. The execution statements are:
option display solve
assignment abort loop
for while repeat
execute
Although there is great freedom about the order in which statements can be placed in a GAMS program, certain
orders are commonly used. The two most common arrangements are discussed in the next sub-section.
set c "crops" ;
parameter yield "crop yield" ;
Style 1 Style 2
Data: Model:
Set declarations and definitions Set declarations
Parameter declarations and definitions Parameter declarations
Assignments Variable declarations
Displays Equation declaration
Model: Equation definition
Variable declarations Model definition
Equation declaration Data:
Equation definition Set definitions
Model definition Parameter definitions
Solution: Assignments
Solve Displays
Displays Solution:
Solve
Displays
The first statement declares that the identifier c is a set and the second defines the elements in the set
Sets and parameters used in the equations must be declared before the equations are specified; they
can be defined, however, after the equation specifications but before a specific equation is used in a
solve statement. This gives GAMS programs substantial organizational flexibility.
Scalars and tables are not separate data types but are a shorthand way to declare a symbol to be a parameter,
and to use a particular format for initializing the numeric data.
Definitions have common characteristics, for example:
parameter a (i,j) input-output matrix
data-type-keyword identifier domain list text
The domain list and the text are always optional characteristics. Other examples are:
In the last example a number of identifiers (separated by commas) are declared in one statement.
As noted previously, we can use any mix of lower and upper case. GAMS makes no distinction
between upper and lower case.
3.4.1 Characters
A few characters are not allowed in a GAMS program because they are illegal or ambiguous on some machines.
Generally all unprintable and control characters are illegal. The only place where any character is legal is in an
$ontext-$offtext block as illustrated in the section on comments below. For completeness the full set of legal
characters are listed in table 3.2. Most of the uncommon punctuation characters are not part of the language,
but can be used freely in text or comments.
3.4.3 Identifiers
Identifiers are the names given to sets, parameters, variables, models, etc. GAMS requires an identifier to start
with a letter followed by more letters or digits. The length of an identifier is currently limited to 31 characters.
3.4 Language Items 31
Identifiers can only contain alphanumeric characters (letters or numbers). Examples of legal identifiers are:
A name used for one data type cannot be reused for another.
3.4.4 Labels
Labels are set elements. They may be up to 31 characters long and can be used in quoted or unquoted form.
The unquoted form is simpler to use but places restrictions on characters used, in that any unquoted label must
start with a letter or digit and can only be followed by letters, digits, or the sign characters + and -. Examples
of unquoted labels are:
In quoted labels, quotes are used to delimit the label, which may begin with and/or include any legal character.
Either single or double quotes can be used but the closing quote has to match the opening one. A label quoted
with double quotes can contain a single quote (and vice versa). Most experienced users avoid quoted labels
because they can be tedious to enter and confusing to read. There are a couple of special circumstances. If one
wants to make a label stand out, then one can, for instance, put asterisks in it and indent it. A more subtle
example is that GAMS keywords can be used as labels if they are quoted. If one needs to use labels like no, ne
or sum then they will have to be quoted. Examples of quoted labels are:
Labels do not have a value. The label 1986 does not have the numerical value 1986 and the label
01 is different from the label 1.
The rules for constructing identifiers and labels are shown in table 3.4.
3.4.5 Text
Identifiers and simple labels can also be associated with a line of descriptive text. This text is more than a
comment: it is retained by GAMS and is displayed whenever results are written for the identifier.
Text can be quoted or unquoted. Quoted text can contain any character except the quote character used. Single
or double quotes can be used but must match. Text has to fit on one line and cannot exceed 80 characters in
length. Text used in unquoted form must follow a number of mild restrictions. Unquoted text cannot start with
a reserved word, .. or = and must not include semicolon ;, commas ,, or slashes /. End of lines terminate
a text. These restrictions are a direct consequence of the GAMS syntax and are usually followed naturally by the
user. Some examples are:
32 GAMS Programs
this is text
final product shipment (tpy)
"quoted text containing otherwise illegal characters ; /,"
use single quotes to put a "double" quote in text
3.4.6 Numbers
Numeric values are entered in a style similar to that used in other computer languages
The common distinction between real and integer data types does not exist in GAMS. If a number is
used without a decimal point it is still stored as a real number.
In addition, GAMS uses an extended range arithmetic that contains special symbols for infinity (INF), negative
infinity (-INF), undefined (UNDF), epsilon (EPS), and not available (NA). One cannot enter UNDF; it is only produced
by an operation that does not have a proper result, such as division by zero. All the other special symbols can
be entered and used as if they were ordinary numbers.
The following example shows various legal ways of entering numbers:
0 156.70 -135 .095 1.
2e10 2e+10 15.e+10 .314e5 +1.7
0.0 .0 0. INF -INF
EPS NA
The letter e denotes the well-known scientific notation allowing convenient representation of very large or small
numbers.
For example:
1e-5 = 1 105 = 0.00001 3.56e6 = 3.56 106 = 3, 560, 000
GAMS uses a smaller range of numbers than many computers are able to handle. This has been
done to ensure that GAMS programs will behave in the same way on a wide variety of machines,
including personal computers. A good general rule is to avoid using or creating numbers with
absolute values greater than 1.0e+20.
A number can be entered with up to ten significant digits on all machines, and more on some.
3.4.7 Delimiters
As mentioned before, statements are separated by a semicolon ;. However, if the next statement begins with a
reserved word (often called keyword in succeeding chapters), then GAMS does not require that the semicolon be
used.
The characters comma , and slash / are used as delimiters in data lists, to be introduced later. The comma
terminates a data element (as does an end-of-line) and the slash terminates a data list.
3.4.8 Comments
A comment is an explanatory text that is not processed or retained by the computer. There are three ways to
include comments in a GAMS program.
The first, already mentioned above, is to start a line with an asterisk * in the first character position. The
remaining characters on the line are ignored but printed on the output file.
The second is to use special block delimiters that cause GAMS to ignore an entire section of the program. The
$ symbol must be in the first character position. The choice between the two ways is a matter of individual taste
or utility. The example below illustrates the use of the block comment.
3.5 Summary 33
$ontext
Following a $ontext directive in column 1 all lines are ignored by GAMS but
printed on the output file until the matching $offtext is encountered, also
in column 1. This facility is often used to logically remove parts of programs
that are not used every time, such as statements producing voluminous reports.
Every $ontext must have a matching $offtext in the same file
$offtext
The third style of comment allows embedding a comment within a line. It must be enabled with the compiler
option $inlinecom or $eolcom as in the following example.
$eolcom #
$inlinecom {}
x = 1 ; # this is a comment
y = 2 ; { this is also a comment } z = 3 ;
3.5 Summary
This completes the discussion of the components of the GAMS language. Many unfamiliar terms used in this
chapter have been further explained in the Glossary.
34 GAMS Programs
4
Set Definitions
4.1 Introduction
Sets are fundamental building blocks in any GAMS model. They allow the model to be succinctly stated and
easily read. In this chapter we will discuss how sets are declared and initialized. There are some more advanced
set concepts, such as assignments to sets as well as lag and lead operations, but these are not introduced until
much later in the book. However the topics covered in this chapter will be enough to provide a good start on
most models.
In GAMS notation, because of character set limitations, the same set must be written as
set S /a, b, c/
The set statement begins with the keyword set (or sets). S is the name of the set, and its members are a, b,
and c. They are labels, but are often referred to as elements or members.
set name is the internal name of the set (also called an identifier) in GAMS. The accompanying text is used to
describe the set or element immediately preceding it.
In quoted labels, quotes are used to delimit the label, which may begin with and/or include any legal character.
Either single or double quotes can be used but the closing quote has to match the opening one. A label quoted
with double quotes can contain a single quote (and vice versa). Most experienced users avoid quoted labels
because they can be tedious to enter and confusing to read. There are a couple of special circumstances. If one
wants to make a label stand out, then to put asterisks in it and indent it, as below, is common. A more subtle
example is that it is possible to use GAMS keywords as labels if they are quoted. If one need to use labels like
no, ne or sum then they will have to be quoted.
Examples of quoted labels are:
Labels do not have a value. The label 1986 does not have the numerical value 1986 and the label
01 is different from the label 1.
Each element in a set must be separated from other elements by a comma or by an end-of-line. In contrast, each
element is separated from any associated text by a blank.
Consider the following example from the Egyptian fertilizer model [FERTS], where the set of fertilizer nutrients
could be written as
or as
set cq "nutrients" / N
P2O5 / ;
The order in which the set members are listed is normally not important. However, if the members represent, for
example, time periods, then it may be useful to refer to next or previous member. There are special operations
to do this, and they will be discussed in Chapter 13. For now, it is enough to remember that the order in which
set elements are specified is not relevant, unless and until some operation implying order is used. At that time,
the rules change, and the set becomes what we will later call an ordered set.
Notice that text may have embedded blanks, and may include special characters such as parentheses. There
are, however, restrictions on special characters in text. Include slashes, commas or semicolons only if the text is
enclosed in quotes. A set definition like
will cause errors since the slash between dollars and ounce will signal the beginning of the set declaration, and
the GAMS compiler will treat ounce as the name of the first element. Further, the slash before gold-price will be
treated as the end of the set definition, and gold-price will be treated as a new set. However, by enclosing the
explanatory text in quotes, this problem is avoided. The following text is valid:
which means that the set includes the ten elements 1991, 1992,...,2000. GAMS builds up these label lists by
looking at the differences between the two labels. If the only characters that differ are digits, and if the number
(say L) formed by these digits in the left one is less than that from the right one (R), then a label is constructed
for every integer in the sequence L to R. Any non-numeric differences or other inconsistencies cause errors.
The following example illustrates the most general form of the asterisked definition:
although the sets, which have 20 members each, have 11 members in common. As a last example, the following
are all illegal because they are not consistent with the rule given above for making lists:
Note one last time that set elements (often referred to as labels) can contain the sign characters - and + as
well as letters and numbers.
at the beginning of the program. When this is done the set keyword need only be used once. If you prefer to
intermingle set declarations with other statements, you will have to use a new set statement for each additional
group of sets.
The following example below shows how two sets can be declared together. Note that the semicolon is used only
after the last set is declared.
sets
s "Sector" / manuf
agri
services
government /
r "regions" / north
eastcoast
midwest
sunbelt / ;
and a second name for the set c is established with either of the following statements
where cp is the new set that can be used instead of the original set c.
The newly introduced set can be used as an alternative name for the original set, and will always
contain only the same elements as the original set.
The alias statement can be used to introduce more than one new name for the original set.
where the new sets cp, cpp, cppp are all new names for the original set c.
The order of the sets in the alias statement does not matter. The only restriction set by GAMS is
that exactly one of the sets in the statement be defined earlier. All the other sets are introduced by
the alias statement.
We will not demonstrate the use of set aliases until later. Just remember they are used for cases when a set has
to be referred to by more than one name.
set
i "all sectors" / light-ind, food+agr, heavy-ind, services /
t(i) "traded sectors" / light-ind, food+agr, heavy-ind /
nt "non-traded sectors" / services / ;
Some types of economic activity, for example exporting and importing, may be logically restricted to a subset
of all sectors. In order to model the trade balance, for example, we need to know which sectors are traded, and
one obvious way is to list them explicitly, as in the definition of the set t above. The specification t(i) means
that each member of the set t must also be a member of the set i. GAMS will enforce this relationship, which
is called domain checking. Obviously the order of declaration is important: the membership of i must be known
before t is declared for checking to be done. There will be much more on this topic in succeeding chapters. For
now it is important to note that domain checking will find any spelling errors that might be made in establishing
the members of the set t. These would cause errors in the model if they went undetected.
It is legal but unwise to define a subset without reference to the larger set, as is done above for the set nt. If
services were misspelled no error would be marked, but the model would give incorrect results. So we urge you
to use domain checking whenever possible. It catches errors and allows you to write models that are conceptually
cleaner because logical relationships are made explicit.
This completes the discussion of sets in which the elements are simple. This is sufficient for most GAMS appli-
cations; however, there are a variety of problems for which it is useful to have sets that are defined in terms of
two or more other sets.
The next two sub-sections explain how to express one-to-one and many-to-many mappings between sets.
set c "countries"
/ jamaica
haiti
guyana
brazil / ;
set p "ports"
/ kingston
s-domingo
georgetown
belem / ;
Then a set can be created to associate each port with its country, viz.,
/ kingston .jamaica
s-domingo .haiti
georgetown .guyana
belem .brazil /;
The dot between kingston and jamaica is used to create one such pair. Blanks may be used freely around
the dot for readability. The set ptoc has four elements, and each element consists of a port-country pair. The
notation (p,c) after the set name ptoc indicates that the first member of each pair must be a member of the set
p of ports, and that the second must be in the set c of countries. This is a second example of domain checking.
GAMS will check the set elements to ensure that all members belong to the appropriate sets.
set i / a, b /
j / c, d, e /
ij1(i,j) /a.c, a.d/
ij2(i,j) /a.c, b.c/
ij3(i,j) /a.c, b.c, a.d, b.d/ ;
ij1 represents a one-to-many mapping where one element of i maps onto many elements of j.
ij2 represents a many-to-one mapping where many elements of i map onto one element of j.
ij3 is the most general case where many elements of i map on to many elements of j.
set i / a, b /
j / c, d, e /
ij1(i,j) /a.(c,d)/
ij2(i,j) /(a,b).c/
ij3(I,j) /(a,b).(c,d)/ ;
The parenthesis provides a list of elements that can be expanded when creating pairs.
When complex sets like this are created, it is important to check that the desired set has been
obtained. The checking can be done by using a display statement.
The concepts may be generalized to set with more than two labels per set element. Mathematically these are
called 3-tuples, 4-tuples, or more generally, n-tuples.
This section ends with some examples to illustrate definitions of multi-label set elements. Some examples of the
compact representation of sets of n-tuples using combinations of dots, parentheses and commas are shown in
table 4.1.
Construct Result
(a,b).c.d a.c.d, b.c.d
(a,b).(c,d) .e a.c.e, b.c.e, a.d.e, b.d.e
(a.1*3).c (a.1, a.2, a.3).c or a.1.c, a.2.c, a.3.c
1*3. 1*3. 1*3 1.1.1, 1.1.2, 1.1.3, ..., 3.3.3
Note that the asterisk can also be used in conjunction with the dot. Recall that the elements of the list 1*4 are
{1, 2, 3, 4} when examining the examples in table 4.1.
4.6 Summary 41
4.6 Summary
In GAMS, a simple set consists of a set name and the elements of the set. Both the name and the elements may
have associated text that explains the name or the elements in more detail. More complex sets have elements that
are pairs or even n-tuples. These sets with pairs and n-tuples are ideal for establishing relationships between the
elements in different sets. GAMS also uses a domain checking capability to help catch labeling inconsistencies
and typographical errors made during the definition of related sets.
The discussion here has been limited to sets whose members are all specified as the set is being declared. For
many models this is all you need to know about sets. Later we will discuss more complicated concepts, such as
sets whose membership changes in different parts of the model (assignment to sets) and other set operations such
as unions, complements and intersections.
42 Set Definitions
5
5.1 Introduction
One of the basic design paradigms of the GAMS language has been to use data in its most basic form, which
may be scalar, list oriented, or tables of two or more dimensions. Based on this criterion, three data types are
introduced in this chapter.
Scalar Single (scalar) data entry.
Parameter List oriented data.
Table Table oriented data. Must involve two or more dimensions.
Each of these data types will be explained in detail in the following sections.
Initialization of data can only be done once for parameters; thereafter data must be modified with
assignment statements.
5.2 Scalars
The scalar statement is used to declare and (optionally) initialize a GAMS parameter of dimensionality zero.
That means there are no associated sets, and that there is therefore exactly one number associated with the
parameter.
Scalar name is the internal name of the scalar (also called an identifier) in GAMS. The accompanying text is
used to describe the element immediately preceding it. Signed num is a signed number and is assigned to be the
value of scalar name.
As with all identifiers, scalar name has to start with a letter followed by more letters or digits. It can only contain
alphanumeric characters, and can be up to 31 characters long. Explanatory text must not exceed 254 characters
and must all be contained on the same line as the identifier or label it describes.
44 Data Entry: Parameters, Scalars & Tables
The statement above initializes rho and life, but not irr. Later on another scalar statement can be used to
initialize irr, or, (looking ahead to a notion that will be developed later), an assignment statement could be used
to provide the value:
irr = 0.07;
5.3 Parameters
While parameter is a data type that encompasses scalars and tables, the discussion in this chapter will focus on
the use of parameters in data entry. List oriented data can be read into GAMS using the parameter statement.
Param name is the internal name of the parameter (also called an identifier) in GAMS. The accompanying text is
used to describe the parameter immediately preceding it. Signed num is a signed number and is declared to be
the value of the entry associated with the corresponding element.
As with all identifiers, param name has to start with a letter followed by more letters or digits. It can only contain
alphanumeric characters, and can be up to 31 long. Explanatory text must not exceed 254 characters and must
all be contained on the same line as the identifier or label it describes.
A parameter may be indexed over one or more sets (the maximum number being 10). The elements in the data
should belong to the set that the parameter is indexed over.
Parameter initialization requires a list of data elements, each consisting of a label and a value. Slashes must be
used at the beginning and end of the list, and commas must separate data elements entered more than one to a
line. An equals sign or a blank may be used to separate the label-tuple from its associated value. A parameter
can be defined in a similar syntax to that used for a set.
The domain checking specification for dd means that there will be a vector of data associated with it, one number
corresponding to every member of the set j listed. The numbers are specified along with the declaration in a
format very reminiscent of the way we specified sets: in this simple case a label followed by a blank separator and
then a value. Any of the legal number entry formats are allowable for the value. The default data value is zero.
Since monterrey has been left out of the data list, then the value associated with dd(monterrey), the market
share in monterrey, would be zero.
We can also put several data elements on a line, separated by commas:
parameter salaries(employee,manager,department)
/anderson .murphy .toy = 6000
hendry .smith .toy = 9000
hoffman .morgan .cosmetics = 8000 / ;
All the mechanisms using asterisks and parenthesized lists that we introduced in our discussion of sets are available
here as well. Below is an artificial example, in which a very small fraction of the total data points are initialized.
GAMS will mark an error if the same label combination (or label-tuple) appears more than once in a data list.
In this example, the twelve elements row1.col2 to row1.col7 and row4.col2 to row4.col7 are all initialized
at 12, the single element row10.col10 at 17, and the seven elements rows1.col10 to row7.col10 at 33. The
other 80 elements (out of a total of 100) remain at their default value, which is 0. This example shows the ability
of GAMS to provide a concise initialization, or definition, for a sparse data structure.
5.4 Tables
Tabular data can be declared and initialized in GAMS using a table statement. For 2- and higher-dimensional
parameters this provides a more concise and easier method of data entry than the list based approach, since each
label appears only once (at least in small tables).
Table name is the internal name of the table (also called an identifier) in GAMS. The accompanying text is used
to describe the parameter immediately preceding it. Signed num is a signed number and is declared to be the
value of the entry associated with the corresponding element.
The table statement is the only statement in the GAMS language that is not free format.
The rules for forming simple tables are straightforward. The components of the header line are the by now
familiar keyword-identifier-domain list-text sequence, the domain-list and text being optional. Labels
are used on the top and the left to map out a rectangular grid that contains the data values. The order of labels
is unimportant, but if domain checking has been specified each label must match one in the associated set. Labels
must not be repeated, but can be left out if the corresponding numbers are all zero or not needed. At least
one blank must separate all labels and data entries. Blank entries imply that the default value (zero) will be
associated with that label combination.
Notice also that, in contrast to the set, scalar, and parameter statements, only one identifier can
be declared and initialized in a table statement.
sets i "plants"
/ inchon,ulsan,yosu /
m "productive units"
atmos-dist "atmospheric distillation unit"
steam-cr "steam cracker"
aromatics "aromatics unit"
hydrodeal "hydrodealkylator" / ;
table ka(m,i) "initial cap. of productive units (100 tons per yr)"
inchon ulsan yosu
atmos-dist 3702 12910 9875
steam-cr 517 1207
aromatics 181 148
hydrodeal 180 ;
In the example above, the row labels are drawn from the set m, and those on the column from the set i. Note
that the data for each row is aligned under the corresponding column headings.
If there is any uncertainty about which data column a number goes with, GAMS will protest with an
error message and mark the ambiguous entry.
table ka(m,i) initial cap. of productive units (100 tons per yr)
inchon ulsan
atmos-dist 3702 12910
steam-cr 517
aromatics 181
hydrodeal 180
+ yosu
atmos-dist 9875
steam-cr 1207
aromatics 148 ;
The crucial item is the plus + sign above the row labels and to the left of the column labels in the continued
part of the table. The row labels have been duplicated, except that hydroreal has been left out, not having
associated data. Tables can be continued as many times as necessary.
table upgrade(strat,size,tech)
small.tech1 small.tech2 medium.tech1 medium.tech2
strategy-1 .05 .05 .05 .05
strategy-2 .2 .2 .2 .2
strategy-3 .2 .2 .2 .2
strategy-4 .2 .2
Here we encounter the display statement again. It causes the data associated with upgrade and upgradex to
be listed on the output file.
5.5 Acronyms
An acronym is a special data type that allows the use of strings as values.
Acronym name is an identifier and follows the same naming convention as other identifiers like names of sets,
parameters, or tables.
m-2 wednesday
m-3 friday
m-4 monday
m-5 thursday
/ ;
In the example above, data entries are in the form of strings like monday and tuesday. By declaring each of
those character strings as acronyms, this kind of data entry can be used by GAMS. Sections 6.2.7 and 11.2.5 will
explain the further use of acronyms once entered in this form.
5.6 Summary
In this chapter, the declaration and initialization of parameters with the parameter, scalar and the table
statement have been discussed. The next chapter will describe how this data can be changed with assignment
statements.
50 Data Entry: Parameters, Scalars & Tables
6
6.1 Introduction
Data once initialized may require manipulation in order to bring it to the form required in the model. The first
part of this chapter will deal explicitly with parameter manipulation. The rest of the chapter will be devoted to
explaining the ramifications: indexed assignment functions, index operations.
scalar x / 1.5/ ;
x = 1.2;
x = x + 2;
The scalar x is initialized to be 1.5. The second statement changes the value to 1.2, and the third changes it
to 3.2. The second and third statement assignments have the effect of replacing the previous value of x, if any,
with a new one.
Note that the same symbol can be used on the left and right of the = sign. The new value is not available until
the calculation is complete, and the operation gives the expected result.
An assignment cannot start with a reserved word. A semicolon is therefore required as a delimiter
before all assignments.
This means that for every member of the set d, a value is assigned to DJ. This can be written in GAMS as
follows,
dj(d) = 2.75*da(d) ;
This assignment is known technically as an indexed assignment and set d will be referred to as the controlling
index or controlling set.
The index sets on the left hand side of the assignment are together called the controlling domain of
the assignment
The extension to two or more controlling indices should be obvious. There will be an assignment made for each
label combination that can be constructed using the indices inside the parenthesis. Consider the following example
of an assignment to all 100 data elements of a.
The calculation in the last statement is carried out for each of the 100 unique two-label combinations that can
be formed from the elements of row and col. The first of these is, explicitly,
a(r-7,c-4) = -2.36 ;
This statement assigns a constant value to one element of a. All other elements of a remain unchanged. Either
single or double quotes can be used around the labels.
The number of controlling indices on the left of the = sign should be at least as many as the number
of indices on the right. There should be no index on the right hand side of the assignment that is not
present on the left unless it is operated on by an indexed operator
a(row,col-2) = 22 - c(col) 1;
6.3 Expressions 53
GAMS will flag this statement as an error since col is an index on the right hand side of the equation but not
on the left.
Each set is counted only once to determine the number of controlling indices. If the same set appears
more than once within the controlling domain, the second and higher occurrences of the set should
be aliases of the original set in order for the number of controlling indices to be equal to the
number of indices.
Consider the following statement as an illustration,
This statement has only one controlling index (row). If one steps through the elements of row one at a time
assignments will be made only to the diagonal entries in b. This will assign exactly 10 values! None of the
off-diagonal elements of b will be filled.
If an additional name is provided for row and used in the second index position, then there will be two controlling
indices and GAMS will make assignments over the full Cartesian product, all 100 values. Consider the following
example,
alias(row,rowp) ;
b(row,rowp) = 7.7 - r(row) + r(rowp) ;
Extended range arithmetic will be discussed later in this Section. The values most often used are NA in incomplete
tables, and INF for variable bounds.
Acronyms contain no numeric value, and are treated as character strings only.
6.3 Expressions
An expression is an arbitrarily complicated specification for a calculation, with parentheses nested as needed for
clarity and intent. In this section, the discussion of parameter assignments will continue by showing in more detail
the expressions that can be used on the right of the = sign. All numerical facilities available in both standard and
extended arithmetic will be covered.
They are listed above in precedence order, which determines the order of evaluation in an expression without
parentheses.
Consider, for example:
x = 5 + 4*3**2 :
x = 5 + (4*(3**2)) ;
It is better to use parentheses than to rely on the precedence of operators, since it prevents errors
and clarifies intentions.
Expressions may be freely continued over many lines: an end-of-line is permissible at any point
where a blank may be used. Blanks may be used for readability around identifiers, parentheses and
operator symbols. Blanks are not allowed within identifiers or numbers, and are significant inside the
quote marks used to delimit labels.
x**n is calculated inside GAMS as exp[n*log(x)]. This operation is not defined if x has a negative
value, and an error will result. If the possibility of negative values for x is to be admitted and the
exponent is known to be an integer, then a function call, power(x,n), is available.
Three additional capabilities are available to add power and flexibility of expression calculations. They are indexed
operations, functions and extended range arithmetic.
These four operations are performed over one or more controlling indices. The syntax in GAMS for these
operations is,
If there is only one controlling index, the parentheses around it can be removed. The most common of these is
sum, which is used to calculate totals over the domain of a set. Consider the following simple example adapted
from [ANDEAN] for illustration.
P
This would be written, using normal mathematical representation, as totCm = i Cim .
The index over which the summation is done, i, is separated from the reserved word sum by a left parenthesis
and from the data term capacity(i,m) by a comma. i is again called the controlling index for this operation.
The scope of the control is the pair of parentheses () that starts immediately after the sum. It is not likely to be
useful to have two independent index operations controlled by the same index.
It is also possible to sum simultaneously over the domain of two or more sets, in which case more parentheses are
needed. Also, of course, an arithmetic expression may be used instead of an identifier;
The smin and smax operations are used to find the largest and smallest values over the domain of the index set
or sets. The index for the smin and smax operators is specified in the same manner as in the index for the sum
operator. Consider the following example to find the largest capacity,
lrgunit = smax((i,m),capacity(i,m));
6.3.3 Functions
GAMS provides commonly used standard functions such as exponentiation, and logarithmic, and trigonometric
functions. The complete list of available functions is given in table 6.1. There are cautions to be taken when
functions appear in equations; these are dealt with in Section 8.4, page 68.
Function Description
errorf(x) Integral of the standard normal distribution from to x
exp(x) Exponential, ex
log(x) Natural logarithm, loge x
log10(x) Common logarithm, log10 x
normal(x,y) Random number normally distributed with mean x and standard deviation y
uniform(x,y) Random number with uniform distribution between x and y
abs(x) Absolute Value of x, i.e. |x|
ceil(x) Ceiling of x. Smallest integer x
floor(x) Floor of x. Largest integer x
mapval(x) Mapping function. Assigns unique numbers to special values.
max(x,y,..) Largest value among all arguments.
min(x,y,..) Smallest value among all arguments
mod(x,y) Remainder. x y trunc(x/y)
power(x,y) Integer power. xy , where y must be an integer
round(x) round x to the nearest integer
round(x,y) Rounds x to y decimal places right (+) or left (-) to the decimal point
sign(x) Returns 1 if x > 0, 1 if x < 0, and 0 if x = 0
sqr(x) Square of x. x2
sqrt(x) Square root of x. x
trunc(x) sign(x)*floor(abs(x))
arctan(x) tan1 x. Result in radians
cos(x) cos x; x in radians
sin(x) sin x; x in radians
x(j) = log(y(j)) ;
which replaces the current value of x with the natural logarithm of y over the domain of the index set j.
GAMS has defined the results of all arithmetic operations and all function values using these special values.
The results can be inspected by running the model library problem [CRAZY] . As one would expect, 1+INF
evaluates to INF, and 1-EPS to 1.
The mapval function should be used in comparisons involving extended range arithmetic. Only the
extended range arithmetic shown in the table above give non-zero values for mapval. For example,
mapval(a) takes a value of 6 if a is inf. All regular numbers result in a mapval of 0.
The following table shows a selection of results for exponentiation and division for a variety of input parameters.
Value Operations
a b a**b power(a,b) a/b
2 2 4 4 1
-2 2 undf 4 -1
2 2.1 4.28 undf .952
na 2.5 na na na
3 0 1 1 undf
inf 2 inf inf inf
2 inf undf undf 0
One should avoid creating or using numbers with absolute values larger than 1.0E20. If a number is
too large, it may be treated by GAMS as undefined (UNDF), and all values derived from it in a model
may be unusable. Always use INF (or -INF) explicitly for arbitrarily large numbers
When an attempted arithmetic operation is illegal or has undefined results because of the value of arguments
(division by zero is the normal example), an error is reported and the result is set to undefined (UNDF).
From there on, UNDF is treated as a proper data value and does not trigger additional error messages.
GAMS will not solve a model if an error has been detected, but will terminate with an error
condition.
It is thus always necessary to anticipate conditions that will cause errors, such as divide by zero. This is most
easily done with the dollar control, and will be discussed in the next section.
6.4 Summary 57
6.4 Summary
GAMS provides powerful facilities for data manipulation with parallel assignment statements, built-in functions
and extended range arithmetic.
58 Data Manipulations with Parameters
7
Variables
7.1 Introduction
This chapter covers the declaration and manipulation of GAMS variables. Many of the concepts covered in the
previous Chapters are directly applicable here.
A variable is the GAMS name for what are called endogenous variables by economists, columns or activities
by linear programming experts, and decision variables by industrial Operations Research practitioners. They are
the entities whose values are generally unknown until after a model has been solved. A crucial difference between
GAMS variables and columns in traditional mathematical programming terminology is that one GAMS variable
is likely to be associated with many columns in the traditional formulation.
Var type is the optional variable type that is explained in detail later. Var name is the internal name of the
variable (also called an identifier) in GAMS. An identifier has to start with a letter followed by more letters or
digits. It can only contain alphanumeric characters, and can be up to 31 characters long. The accompanying text
is used to describe the set or element immediately preceding it. This must not exceed 254 characters and must
all be contained on the same line as the identifier it describes.
One important difference between variable and parameter declarations is that values cannot be initialized in a
variable declaration.
A typical variable statement, adapted from [RAMSEY], is shown below for illustration:
The declaration of k above implies, as usual, that references to k are restricted to the domain of the set t. A model
that includes k will probably have several corresponding variables in the associated mathematical programming
60 Variables
problem: most likely one for each member of t. In this way, very large models can be constructed using a small
number of variables. (It is quite unusual for a model to have as many as 50 distinct variables.) It is still unclear
from the declaration whether utility is not domain checked or whether it is a scalar variable, i.e., one without
associated sets. Later references will be used to settle the issue.
It is important that variable declarations include explanatory text and that this be as descriptive as possible,
since the text is used to annotate the solution output. Note the use of per instead of / in the text above:
slashes are illegal in all unquoted text.
The default type is free, which means that if the type of the variable is not specified, it will not be bounded
at all. The most frequently used types are free and positive, for descriptions of variables for which negative
values are meaningless, such as capacities, quantities or prices.
Four additional, although more exotic, variable types - sos1, sos2, semicont and semiint are available in GAMS.
These are explained in Section 17.2.1.
variables
u(c,i) "purchase of domestic materials (mill units per yr)"
v(c.j) "imports (mill tpy)"
e(c,i) "exports (mill tpy)"
phi "total cost (mill us$)"
phipsi "raw material cost (mill us$)" ;
positive variables u, v, e ;
It is possible to declare an identifier more than once, but that the second and any subsequent
declarations should only add new information that does not contradict what has already been
entered.
The second popular way of declaring variables is to list them in groups by type. We rewrite the example above
using this second method:
7.3 Variable Attributes 61
free variables
phi "total cost (mill us$)"
phipsi "raw material cost (mill us$)"
positive variables
u(c,i) "purchase of domestic materials (mill units per yr)"
v(c,j) "imports (mill typ)"
e(c,i) "exports (mill typ)" ;
The user distinguishes between these suffix numbers when necessary by appending a suffix to the variable name.
For binary and integer variable types, the consequences of the type declaration cannot be completely
undone.
Bounds on variables are the responsibility of the user. After variables have been declared, default bounds have
already been assigned: for many purposes, especially in linear models, the default bounds are sufficient. In
nonlinear models, on the other hand, bounds play a far more important role. It may be necessary to provide
bounds to prevent undefined operations, such as division by zero.
It is also often necessary to define a reasonable solution space that will help to make the nonlinear programming
problem be solved more efficiently.
The lower bound cannot be greater than the upper: if you happen to impose such a condition,
GAMS will exit with an error condition.
Note that, in the first statement, the index set covering the domain of x appears after the suffix. The first
assignment puts an upper bound on all variables associated with the identifier x. The statement on the second
line bounds one particular entry. The statement on the last line sets the level values of the variables in c to four
times the values in the parameter cinit.
Remember that the order is important in assignments, and notice that the two pairs of statements below produce
very different results. In the first case, the lower bound for c(1985) will be 0.01, but in the second, the lower
bound is 1.
Everything works as described in the previous chapter, including the various mechanisms described there of
indexed operations, dollar operations, subset assignments and so on.
As with parameters, a variable must have some non-default data values associated with it before one can use
it in a display statement or on the right hand side of an assignment statement. After a solve statement (to be
discussed later) has been processed or if non-default values have been set with an assignment statement, this
condition is satisfied.
The .fx suffix is really just a shorthand for .lo and .up and can therefore only be used only on the
left-hand side of an assignment statement.
The output looks similar, except that (of course) the listing shows which of the values is being displayed. Because
zeroes, and especially all zero rows or columns, are suppressed, the patterns seen in the level and marginal displays
will be quite different, since non-zero marginal values are often associated with activity levels of zero.
We should mention here a clarification of our previous discussion of displays. It is actually the default values
that are suppressed on display output. For parameters and variable levels, the default is zero, and so zero entries
are not shown. For bounds, however, the defaults can be non-zero. The default value for the upper bound of a
positive variable is +INF, and if above you also would display v.up, for example, you will see:
If any of the bounds have been changed from the default value, then only the entries for the changed elements will
be shown. This sounds confusing, but since few users display bounds it has not proved troublesome in practice.
7.5 Summary
Remember that wherever a parameter can appear in a display or an assignment statement, a variable can also
appear - provided that it is qualified with one of the four suffixes. The only places where a variable name can
appear without a suffix is in a variable declaration, as shown here, or in an equation definition, which is discussed
in the next chapter.
64 Variables
8
Equations
8.1 Introduction
Equations are the GAMS names for the symbolic algebraic relationships that will be used to generate the
constraints in the model. As with variables, one GAMS equation will map into arbitrarily many individual
constraints, depending on the membership of the defining sets.
Eqn name is the internal name of the equation (an identifier) in GAMS. An identifier has to start with a letter
followed by more letters or digits. It can only contain alphanumeric characters, and can be up to 31 characters
long. The accompanying text is used to describe the set or element immediately preceding it. This must not
exceed 254 characters and must all be contained on the same line as the identifier it describes.
There are no modifying keywords as there are with variables, and no initializing data list as there may be with
parameters or sets.
equations
cost total cost definition
invb(q) inventory balance
sbal(q,s) shift employment balance ;
66 Equations
The declaration of the first equation follows the keyword equations. This declaration begins with the name of
the equation, in this case cost, and is followed by the text, namely Total cost definition. The equation
cost above is a scalar equation, which will produce at most one equation in the associated optimization problem.
By contrast, the equation sbal is declared over the sets q (4 members) and s (2 members), and is thus likely to
produce eight individual equations, one for each unique combination of labels. The circumstances under which
less than eight equations might be produced will be discussed in later chapters. It is certainly true, however, that
no more than eight equations will be produced.
Eqn name is the name of the equation as in the equation declaration. The two dots .. are always required
between the equation name and start of the algebra. The expressions in the equation definition can be of the
forms discussed in the Chapters before, but can involve variables as well. Eqn type refers to the symbol between
the two expressions that form the equation, and can be of the following types,
=e= Equality: rhs must equal lhs
=g= Greater than: lhs must be greater than or equal to rhs
=l= Less than: lhs must be less than or equal to rhs
=n= No relationships enforced between lhs and rhs. This equation type is rarely used.
As with the assignment statement, equation definitions can be carried over as many lines of input as
needed. Blanks can be inserted to improve readability, and expressions can be arbitrarily
complicated.
An equation, once defined, can not be altered or re-defined. If one needs to change the logic, a new
equation with a new name will have to be defined. It is possible, however, to change the meaning of
an equation by changing the data it uses, or by using exception handling mechanisms (dollar
operations) built into the definition
Obj is the name of the equation being defined. The =e= symbol means that this is an equality. Any of the
following forms of the equation are mathematically equivalent,
8.3 Equation Definitions 67
The arrangement of the terms in the equation is a matter of choice, but often a particular one is
chosen because it makes the model easier to understand.
Domain checking ensures that the domain over which an equation is defined must be the set or a
subset of the set over which the equation is declared.
Consider the following example of a singly indexed equation, meaning one that produces a separate constraint
for each member of the driving (or controlling) set.
As t has three members, three constraints will be generated , each one specifying separately for each member of
t, the dependence of g on m. Mew and xsi are parameters: the data associated with them are used in building up
the individual constraints. These data do not have to be known when the equation is defined, but do have to be
when a model containing the equation is solved.
The extension to two or more index positions on the left of the .. should be obvious. There will be one constraint
generated for each label combination that can constructed using the indices inside the parenthesis. Here are two
examples from [AIRCRAFT], a scheduling model.
The domain of definition of both equations is the Cartesian product of j and h: constraints will be generated for
every label pair that can be constructed from the membership of the two sets.
All the mechanisms that may be used to evaluate expressions in assignments are also available in
equations.
Consider the following example adapted from [CHENERY] showing parentheses and exponentiation,
The occurrence of any function with endogenous arguments implies that the model is not linear.
Functions with endogenous arguments can be further classified into types listed in table 8.1.
Smooth functions can be used routinely in nonlinear models, but non-smooth ones may cause numerical problems
and should be used only if unavoidable, and only in a special model type called dnlp. However, the use of the
dnlp model type is strongly discouraged and the use of binary variables is recommended to model non-smooth
functions. Discontinuous functions are not allowed at all with variable arguments.
A fuller discussion is given in Chapter 9, page 71. For convenience, all the available functions are classified in
table 8.2.
Classification is for the first argument only. The second must be an integer and is normally a constant.
c.lo(t) = 0.01 ;
util .. utility =e= sum(t, beta(t)*log(c(t))) ;
The bounding on c(t) away from 0 prevents the log function from being undefined.
9.1 Introduction
This chapter brings together all the concepts discussed in previous chapters by explaining how to specify a model
and solve it.
Model name is the internal name of the model (also called an identifier) in GAMS. The accompanying text is
used to describe the set or element immediately preceding it. Eqn name is the name of an equation that has been
declared prior to the model statement.
As with all identifiers, model name has to start with a letter followed by more letters or digits. It can only contain
alphanumeric characters, and can be up to 10 characters long. Explanatory text must not exceed 80 characters
and must all be contained on the same line as the identifier or label it describes.
An example of a model definition in GAMS is shown below.
The model is called transport and the keyword all is a shorthand for all known (declared) equations.
Several models can be declared (and defined) in one model statement. This is useful when experimenting with
different ways of writing a model, or if one has different models that draw on the same data. Consider the
following example, adapted from [PROLOG], in which different groups of the equations are used in alternative
versions of the problem. Three versions are solved the linear, nonlinear, and expenditure versions. The model
statement to define all three is
72 Model and Solve Statements
where cb, rc, etc. are the names of the equations. We will describe below how to obtain the solution to each of
the three models.
This sets the solver an upper limit of 60 seconds to attempt to solve the problem to optimality.
It is important to remember that GAMS itself does not solve your problem, but passes the problem
definition to one of a number of separate solver programs.
Model name is the name of the model as defined by a model statement. Var name is the name of the objective
variable that is being optimized. Model type is one of the model types described before. An example of a solve
statement in GAMS is shown below.
Solve and using are reserved words. Transport is the name of the model, lp is the model type, minimizing
is the direction of optimization, and cost is the objective variable. The opposite of minimizing is maximizing,
both reserved words. Note that an objective variable is used instead of an objective row or function
The objective variable must be scalar and type free, and must appear in the least one of the
equations in the model.
The next two sub-sections will describe briefly below what happens when a solve statement is processed, and more
details on how the resulting output is to be interpreted will be given in the next chapter. After that sequences of
solve statements will be discussed. The final section will describe options that are important in controlling solve
statements.
3. Each equation fits into the specified problem class (linearity for lp, continuous derivatives for nlp, as
we outlined above)
2. Debugging and comprehension aids are produced and written to the output file (EQUATION LISTING,
etc).
3. GAMS verifies that there are no inconsistent bounds or unacceptable values (for example NA or UNDF)
in the problem.
4. Any errors detected at this stage cause termination with as much explanation as possible, using the
GAMS names for the identifiers causing the trouble.
5. GAMS passes control to the solution subsystem and waits while the problem is solved.
9.4 Programs with Several Solve Statements 75
6. GAMS reports on the status of the solution process and loads solution values back into the GAMS
database. This causes new values to be assigned to the .l and .m fields for all individual equations and
variables in the model. A row by row and column by column listing of the solution is provided by
default. Any apparent difficulty with the solution process will cause explanatory messages to be
displayed. Errors caused by forbidden nonlinear operations are reported at this stage.
The outputs from these steps, including any possible error messages, are discussed in detail in the next chapter.
When there is more than one solve statement in your program, GAMS uses as much information as possible
form the previous solution to provide a starting point in the search for the next solution.
qs(upper,fuel-oil,sulfur) = 3.5 ;
solve oil using lp maximizing phi;
report(cr,p,base) = z.l(cr,p) ;
report(sulfur,limit,base) = qs(upper,fuel-oil,sulfur);
qs (upper,fuel-oil,sulfur) = 3.4 ;
solve oil using lp maximizing phi ;
report(cr,p,one) = z.l(cr,p) ;
76 Model and Solve Statements
report(sulfur,limit,one) = qs (upper,fuel-oil,sulfur);
qs(upper,fuel-oil,sulfur) = 5.0 ;
solve oil using lp maximizing phi ;
report(cr,p,two) = z.l(cr,p) ;
report(sulfur,limit,two) = qs(upper,fuel-oil,sulfur);
display report ;
This example shows not only how simply sensitivity analysis can be done, but also how the associated multi-case
reporting can be handled. The parameter qs is used to set the upper bound on the sulfur content in the fuel
oil, and the value is retrieved for the report.
The output from the display is shown below. Notice that there is no production at all if the permissible sulfur
content is lowered. The case attributes have been listed in the row SULFUR.LIMIT. The wild card domain is useful
when generating reports: otherwise it would be necessary to provide special sets containing the labels used in the
report. Any mistakes made in spelling labels used only in the report should be immediately apparent, and their
effects should be limited to the report. Section 14.5.1, page 124, contains more detail on how to arrange reports
in a variety of ways.
Another use of multiple solve statements is to permit iterative solution of different blocks of equations, solution
values from the first are used as data in the next. These decomposition methods are useful for certain classes
of problems because the sub-problems being solved are smaller, and therefore more tractable. One of the most
common examples of such a method is the Generalized Benders Decomposition method.
An example of a problem that is solved in this way is an input-output system with endogenous prices, described
in Henaff (1980)1 . The model consists of two groups of equations. The first group uses a given final demand
vector to determine the output level in each sector. The second group uses some exogenous process and input-
output data to compute sectoral price levels. Then the resulting prices are used to compute a new vector of final
demands, and the two block of equations are solved again. This iterative procedure is repeated until satisfactory
convergence is obtained. Henaff has used GAMS statements to perform this kind of calculation. The statements
that solve the system for the first time and the next iteration are shown below:
pbar(ta) = (sum(ipd.l(i,ta))/4.);
d(i,t) = (db(i)*g(t))/(pd.l(i,t)/pbar(t)) ;
1 Henaff, Patrick (1980). An Input-Output Model of the French Economy, Masters Thesis, Department of Economics, University
of Maryland.
9.5 Making New Solvers Available with GAMS 77
Mb is a set of material balance (input-output) equations, and output is a total output equation. Dual is a group
of price equations, and totp is an equation that sums all the sectoral prices. The domestic prices pd used in the
calculation of the average price pbar are divided by four because there are four sectors in this example. Also
the .l is appended to pd to indicate that this is the level of the variable in the solution of the model namely in
dualmodel. Thus the iterative procedure uses solution values from one iteration to obtain parameter values for
the next one. In particular, both pbar and pd are used to compute the demand d for the i-th product in time
period t, d(i,t). Also, the base year demand db and the growth factor g are used in that calculation. Then
when the new final demand vector d is calculated, the two blocks of equations are solved again.
This completes the discussion of the model and solve statements. In the next chapter the various components of
GAMS output are described in some detail.
78 Model and Solve Statements
10
GAMS Output
10.1 Introduction
The output from GAMS contains many aids for checking and comprehending a model. In this chapter the contents
of the output file are discussed. Ways by which the amount of diagnostic output produced can be controlled will
also be discussed, although complete lists of all these controls are not given until later. A small nonlinear model,
[ALAN] by Alan S. Manne, is used to illustrate the output file, and list it piece by piece as we discuss the
various components. The possibilities for extension to large models with voluminous output (which is when the
diagnostics are really useful) should be apparent.
The output from a GAMS run is produced on one file, which can be read using any text editor. The default name
of this output file depends on the operating system, but Appendix C describes how this default can be changed.
The display statement, described in detail in Chapter 14, can be used to export information from the GAMS
program to the listing file.
$Offtext
* This model has been modified for use in the documentation
Positive Variable x;
Note that the first line number shown is 9. If the lines on the input are counted, it can be seen that this comment
line shown above appears after 8 lines of dollar directives and comments.
The line starting $title has caused text of the users choice to be put on the page header, replacing the default
tile, which just announces GAMS. The following $- directives are used to display more information in the output
file and we be discussed. The text within the $ontext-$offtext pair is listed without line numbers, whereas
comments starting with asterisks have line numbers shown. Line numbers always refer to the physical line number
in your input file.
Dollar control directives are only listed if a directive to list them is enabled, or if they contain errors.
10
11 Set i securities /hardware,software,show-biz,t-bills/; alias (i,j);
10.3 Compilation Output 81
12
13 Scalar target target mean annual return on portfolio % /10/,
14 lowyield yield of lowest yielding security,
15 highrisk variance of highest security risk ;
16
17 Parameters mean(i) mean annual returns on individual securities (%)
18
19 / hardware 8
20 software 9
21 show-biz 12
22 t-bills 7 /
23
24 Table v(i,j) variance-covariance array (%-squared annual return)
25
26 hardware software show-biz t-bills
27
28 hardware 4 3 -1 0
29 software 3 6 1 0
30 show-biz -1 1 10 0
31 t-bills 0 0 0 0 ;
32
33 lowyield = smin(i, mean(i)) ;
34 highrisk = smax(i, v(i,i)) ;
35 display lowyield, highrisk ;
36
37 Variables x(i) fraction of portfolio invested in asset i
38 variance variance of portfolio
39
40 Positive Variable x;
41
42 Equations fsum fractions must add to 1.0
43 dmean definition of mean return on portfolio
44 dvar definition of variance;
45
46 fsum.. sum(i, x(i)) =e= 1.0 ;
47 dmean.. sum(i, mean(i)*x(i)) =e= target;
48 dvar.. sum(i, x(i)*sum(j,v(i,j)*x(j))) =e= variance;
49
50 Model portfolio / fsum, dmean, dvar / ;
51
52 Solve portfolio using nlp minimizing variance;
That is the end of the echo of the input file. If errors had been detected, the explanatory messages would be
found in this section of the listing file. All discussion of error messages have been grouped in the section 10.6,
page 92.
Symbol Listing
REF 50
HIGHRISK PARAM DECLARED 15 ASSIGNED 34 REF 35
I SET DECLARED 11 DEFINED 11 REF 11
17 24 33 2*34 37 46
2*47 2*48 CONTROL 33 34 46
47 48
J SET DECLARED 11 REF 24 2*48
CONTROL 48
LOWYIELD PARAM DECLARED 14 ASSIGNED 33 REF 35
MEAN PARAM DECLARED 17 DEFINED 19 REF 33
47
PORTFOLIO MODEL DECLARED 50 DEFINED 50 IMPL-ASN 52
REF 52
TARGET PARAM DECLARED 13 DEFINED 13 REF 47
V PARAM DECLARED 24 DEFINED 24 REF 34
48
VARIANCE VAR DECLARED 38 IMPL-ASN 52 REF 48
52
X VAR DECLARED 37 IMPL-ASN 52 REF 40
46 47 2*48
For each symbol, the name and type of the symbol are first provided. For example, the last symbol listed is X
which is defined to be of type VAR. The complete list of data types are given in table 10.1.
Then comes a list of references to the symbol, grouped by reference type and identified by the line number in the
output file. The actual reference can then be found by referring to the echo print of the program, which has line
numbers on it. In the case of the symbol X in the example above, the list of references as shown in the symbol
reference map are as follows,
DECLARED 37
IMPL-ASN 52
REF 40 46 47 2*48
This means that X is declared on line 37, implicitly assigned through a solve statement on line 52, and referenced
on lines 40, 46, and 47. The entry 2*48 means that there are two references to X on line 48 of the input file .
The complete list of reference types is given below.
DECLARED This is where the identifier is declared as to type. This must be the first appearance of the
identifier.
DEFINED This is the line number where an initialization (a table or a data list between slashes) or symbolic
definition (equation) starts for the symbol.
ASSIGNED This is when values are replaced because the identifier appears on the left of an assignment
statement.
IMPL-ASN This is an implicit assignment: an equation or variable will be updated as a result of being referred
to implicitly in a solve statement.
CONTROL This refers to the use of a set as the driving index in an assignment, equation, loop or other
indexed operation (sum, prod, smin or smax).
REF This is a reference: the symbol has been referenced on the right of an assignment, in a display, in
an equation, or in a model or solve statement.
10.3 Compilation Output 83
Symbol Listing
SETS
I securities
J Aliased with I
PARAMETERS
VARIABLES
EQUATIONS
MODELS
PORTFOLIO
ELEMENT REFERENCES
The dollar control directives are compiler directives that can be put in the input file to control the appearance
and amount of detail in the output produced by the GAMS compiler.
$offlisting, $onlisting
This directive stops the echo print of the input file. $onlisting restores the default.
$offupper, $onupper
This directive causes the echo print of the portion of the GAMS program following the directive to appear
on the output file in the case that it has been entered in. This is the default on newer GAMS systems. It
is necessary if case conventions have been used in the program, for example to distinguish between
variables and equations. $onupper, will cause all echo print to be in upper case.
$ontext, $offtext
$ontext-$offtext pairs are used to create block comments that are ignored by GAMS. Every $ontext
must have a matching $offtext in the same file. The $offtext must be on a line by itself.
$title text
The text can have up to 80 characters. This causes every page of the output to have the title specified.
In all dollar control directives, the $ symbol must be in the first character position on the line.
Dollar control directives, are dynamic: they affect only what happens after they are encountered,
and they can be set and reset wherever appropriate.
They are remembered in continued compilations started from work files The directives that do not have following
text can be entered many to a line, as shown below for the map controls.
If errors are detected because of illegal data operations, a brief message indicating the cause and the line number
of the offending statement will appear.
10.5 Output Produced by a Solve Statement 85
(LHS = 0)
The equation listing is an extremely useful debugging aid. It shows the variables that appear in each
constraint, and what the individual coefficients and right-hand-side value evaluate to after the data
manipulations have been done.
Most of the listing is self-explanatory. The name, text, and type of constraints are shown. The four dashes are
useful for mechanical searching.
All the terms that depend on variables are collected on the left, and all the constant terms are
combined into one number on the right, any necessary sign changes being made.
Four places of decimals are shown if necessary, but trailing zeroes following the decimal point are suppressed.
E-format is used to prevent small numbers being displayed as zero.
The nonlinear equations are treated differently. If the coefficient of a variable in the equation listing
is enclosed in parentheses, then the corresponding constraint is nonlinear, and the value of the
coefficient depends on the activity levels of one or more of the variables. The listing is not algebraic,
but shows the partial derivative of each variable evaluated at their current level values.
Note that, in the equation listing from our example, the equation dvar is nonlinear. A simpler example will help
to clarify the point. Consider the following equation and associated level values.
The coefficient of x is determined by first differentiating the equation above with respect to x. This results in
2*(2*x.l)*power(y.l,3)+ 5, which evaluates to 221. Similarly the coefficient of y is obtained by differenti-
ating the equation above with respect to y which results in 2*(sqr(x.l)*3*sqr(y.l) + 1.5/sqr(y.l), giving
216.1667. Notice that the coefficient of y could not have been determined if its level had been left at zero. The
attempted division by zero would have produced an error and premature termination.
The result of evaluating the left-hand-side of the equation at the initial point is shown at the end of each individual
equation listing. In the example above it is 225.5, and the three asterisks (***) are a warning that the constraint
is infeasible at the starting point.
The order in which the equations are listed depends on how the model was defined. If it was defined
with a list of equation names, then the listing will be in the order in that list. If it was defined as
/all/, then the list will be in the order of declaration of the equations. The order of the entries for
the individual constraints is determined by the label entry order.
X(hardware)
(.LO, .L, .UP = 0, 0, +INF)
1 FSUM
8 DMEAN
(0) DVAR
X(software)
(.LO, .L, .UP = 0, 0, +INF)
1 FSUM
9 DMEAN
(0) DVAR
X(show-biz)
(.LO, .L, .UP = 0, 0, +INF)
1 FSUM
12 DMEAN
(0) DVAR
VARIANCE
(.LO, .L, .UP = -INF, 0, +INF)
-1 DVAR
The order in which the variables appear is the order in which they were declared.
MODEL STATISTICS
The BLOCK counts refer to GAMS equations and variables, the SINGLE counts to individual rows and columns
in the problem generated. The NON ZERO ELEMENTS entry refers to the number of non-zero coefficients in the
problem matrix.
There are four entries that provide additional information about nonlinear models. The NON LINEAR N-Z entry
refers to the number of nonlinear matrix entries in the model.
All forms of nonlinearity do not have the same level of complexity. For example, x*y is a simpler form of
nonlinearity than exp(x*y). So, even though both these terms count as 1 nonlinear entry in the matrix, additional
information is required to provide the user with a feel for the complexity of the nonlinearity. GAMS provides
the CODE LENGTH entry as a good yardstick for this purpose. There are two other entries - DERIVATIVE POOL and
CONSTANT POOL that provide some more information about the nonlinearity. In general, the more nonlinear a
problem is, the more difficult it is to solve.
The times that follow statistics are also useful. The GENERATION TIME is the time used since the syntax check
finished. This includes the time spent in generating the model. The measurement units are given, and represent
ordinary clock time on personal computers, or central processor usage (CPU) time on other machines.
S O L V E S U M M A R Y
The common part of the solve summary is shown above. It can be found mechanically by searching for four
asterisks. The explanation for the information provided in this section follows.
MODEL PORTFOLIO
TYPE NLP
88 GAMS Output
SOLVER MINOS5
This provides the name of the solver used to solve the model.
OBJECTIVE VARIANCE
DIRECTION MINIMIZE
These provide the solver status and model status for the problem, and are discussed in greater detail at the end
of this subsection.
This provides the value of the objective function at the termination of the solve. If the Solver and Model have
the right status, this value is the optimum value for the problem.
These two entries provide the amount of CPU time (in seconds) taken by the solver, as well as the upper limit
allowed for the solver. The solver will stop as soon as the limit on time usage has been reached. The default
limit on time usage is 1000 seconds. This limit can be changed by entering a line containing the statement
option reslim = xx ; in the program before the solve statement, where xx is the required limit on CPU time
in seconds.
These two entries provide the number of iterations used by the solver, as well as the upper limit allowed for the
solver. The solver will stop as soon as this limit is reached. The default limit on iterations used is 1000. This
limit can be changed by entering a line containing the statement option iterlim = nn ; in the program before
the solve statement, where nn is the required limit on the iterations used.
EVALUATION ERRORS 0 0
These two entries provide the number of numerical errors encountered by the solver, as well as the upper limit
allowed for the solver. These errors result due to numerical problems like division by 0. This is suppressed for
LP, RMIP, and MIP models since evaluation errors are not applicable for these model types. The default limit on
evaluation errors used is 0. This limit can be changed by entering a line containing the statement option domlim =
nn ; in the program before the solve statement, where nn is the required limit on the evaluation errors allowed.
The SOLVER STATUS and MODEL STATUS require special explanation. The status for the solver (the state of the
program) and the model (what the solution looks like) are characterized, and a complete list of possible MODEL
STATUS and SOLVER STATUS messages is given below.
Here is a list of possible MODEL STATUS messages:
1 OPTIMAL
This means that the solution is optimal. It only applies to linear problems or relaxed mixed integer
problems (RMIP).
2 LOCALLY OPTIMAL
This message means that a local optimum has been found. This is the message to look for if the problem
is nonlinear, since all we can guarantee for general nonlinear problems is a local optimum.
10.5 Output Produced by a Solve Statement 89
3 UNBOUNDED
This means that the solution is unbounded. This message is reliable if the problem is linear, but
occasionally it appears for difficult nonlinear problems that are not truly unbounded, but that lack some
strategically placed bounds to limit the variables to sensible values.
4 INFEASIBLE
This means that the linear problem is infeasible. Something is probably misspecified in the logic or the
data.
5 LOCALLY INFEASIBLE
This message means that no feasible point could be found for the nonlinear problem from the given
starting point. It does not necessarily mean that no feasible point exists.
6 INTERMEDIATE INFEASIBLE
This means that the current solution is not feasible, but that the solver program stopped, either because
of a limit (iteration or resource), or because of some sort of difficulty. Check the solver status for more
information.
7 INTERMEDIATE NONOPTIMAL
This is again an incomplete solution, but it appears to be feasible.
8 INTEGER SOLUTION
An integer solution has been found to a MIP (mixed integer problem). There is more detail following about
whether this solution satisfies the termination criteria (set by options optcr or optca).
9 INTERMEDIATE NON-INTEGER
This is an incomplete solution to a MIP. An integer solution has not yet been found.
10 INTEGER INFEASIBLE
There is no integer solution to a MIP. This message should be reliable.
11 LIC PROBLEM - NO SOLUTION
The solver cannot find the appropriate license key needed to use a specific subsolver.
ERROR UNKNOWN
After a solver error, the model status is unknown.
ERROR NO SOLUTI0N
An error occurred and no solution has been returned. No solution will be returned to GAMS because of
errors in the solution process.
14 NO SOLUTION RETURNED
A solution is not expected for this solve. For example, the convert solver only reformats the model but
does not give a solution.
15 SOLVED UNIQUE
This indicates a unique solution to a CNS model.
16 SOLVED
A CNS model has been solved but multiple solutions may exist.
17 SOLVED SINGULAR
A CNS model has been solved but the point is singular.
18 UNBOUNDED - NO SOLUTION
The model is unbounded and no solution can be provided.
19 INFEASIBLE - NO SOLUTION
The model is infeasible and no solution can be provided.
2 ITERATION INTERRUPT
This means that the solver was interrupted because it used too many iterations. Use option iterlim to
increase the iteration limit if everything seems normal.
3 RESOURCE INTERRUPT
This means that the solver was interrupted because it used too much time. Use option reslim to increase
the time limit if everything seems normal.
4 TERMINATED BY SOLVER
This means that the solver encountered difficulty and was unable to continue. More detail will appear
following the message.
5 EVALUATION ERROR LIMIT
Too many evaluations of nonlinear terms at undefined values. You should use bounds to prevent forbidden
operations, such as division by zero. The rows in which the errors occur are listed just before the solution.
6 CAPABILITY PROBLEMS
The solver does not have the capability required by the model, for example, BARON has a more limited
set of functions than other solvers.
7 LICENSING PROBLEMS
The solver cannot find the appropriate license key needed to use a specific subsolver.
8 USER INTERRUPT
The user has sent a message to interrupt the solver via the interrupt button in the IDE or sending a
control C from a command line.
9 ERROR SETUP FAILURE
The solver encountered a fatal failure during problem set-up time.
ERROR SOLVER FAILURE
The solver encountered a fatal error.
ERROR INTERNAL SOLVER FAILURE
The solver encountered an internal fatal error.
12 SOLVE PROCESSING SKIPPED
The entire solve step has been skipped. This happens if execution errors were encountered and the GAMS
parameter ExeErr has been set to a nonzero value, or the property MaxExecError has a nonzero value.
ERROR SYSTEM FAILURE
This indicates a completely unknown or unexpected error condition.
GAMS/MINOS
B. A. Murtagh, University of New South Wales
and
P. E. Gill, W. Murray, M. A. Saunders and M. H. Wright
Systems Optimization Laboratory, Stanford University.
The line work space allocated -- 0.04 MB provides the amount of memory used by the solver for the prob-
lem. If the amount of memory the solver estimates that it needs is not available, GAMS will return a message
saying that not enough memory was allocated. GAMS will also return the maximum amount of memory available
on the machine. The user can direct the solver to use less memory by entering a line containing the statement
mymodel.workspace = xx; were mymodel is the name of the model being solved as specified by the model state-
ment, and xx is the amount of memory in Megabytes. Note that the solver will attempt to solve the problem
with xx MB of memory, however it is not guaranteed to succeed since the problem may require more memory.
More information can be obtained for a successful run by entering a line containing the statement option sysout
= on ; in the program above the solve statement.
The order of the equations and variables are the same as in the symbol listing described before and will be
described later
The four columns associated with each entry have the following meaning,
LOWER lower bound (.lo)
LEVEL level value (.l)
UPPER upper bound (.up)
MARGINAL marginal (.m)
92 GAMS Output
For variables the values in the LOWER and UPPER columns refer to the lower and upper bounds. For equations
they are obtained from the (constant) right-hand-side value and from the relational type of the equation. These
relationships were described in Chapter 8.
The LEVEL and MARGINAL values have been determined by the solver, and the values shown are used
to update the GAMS values. In the list they are shown with fixed precision, but the values are
returned to GAMS with full machine accuracy. The single dots . on the list represent zero.
EPS is the GAMS extended value that means very close to but different from zero. It is common to see a marginal
value given as EPS, since GAMS uses the convention that marginal are zero for basic variables, and not zero for
others.
EPS is used with non-basic variables whose marginal values are very close to, or actually, zero, or in
nonlinear problems with superbasic variables whose marginals are zero or very close to it. A
superbasic variable is one between its bounds at the final point but not in the basis.
There are brief explanations of technical terms used in this section in the Glossary. For models that are not solved
to optimality, some constraints may additionally be marked with certain flags. The list of these flags and their
description is given below.
INFES The row or column is infeasible. This mark is made for any entry whose level value is not between
the upper and lower bounds.
NOPT The row or column is non-optimal. This mark is made for any non-basic entries for which the
marginal sign is incorrect, or superbasic ones for which the marginal value is too large.
UNBND The row or column that appears to cause the problem to be unbounded.
If our example had display output for reporting, it would come here.
Effective error detection and recovery are important parts of any modeling system. GAMS is designed around
the assumption that the error State is the normal state of modeling. Experience shows that most compilations
during the early stages of development will produce errors. Not to Worry! The computer is much better at
checking details that the human mind and should be able to provide positive feedback and suggestions about how
to correct errors or avoid ambiguities. Developing a model is like writing a paper or an essay ; many drafts and
rewrites are required until the arguments are presented in the most effective way for the reader and meet all the
requirements of proper English. GAMS acts like a personal assistant with knowledge of mathematical modeling
and of the syntactic and semantic details of the language.
Errors are detected at various stages in the modeling process. Most of them are caught at the compilation stage,
which behaves like the proofreading stage of the modeling process. Once a problem has passed through the
rigorous test of this stage, the error rate drops almost to zero. Most of the execution runs, which are much
more expensive than compilation, proceed without difficulties because GAMS knows about modeling and has
anticipated problems. Many of the typical errors made with conventional programming languages are associated
with concepts that do not exist in GAMS. Those error sources such as address calculations, storage assignment,
subroutine linkages, input-output and flow control create problems at execution time, are difficult to locate,
often lead to long and frustrating searches, and leave the computer user intimidated. GAMS takes a radically
different approach. Errors are spotted as early as possible, are reported in a way understandable to the user,
including clear suggestions for how to correct the problem, and a presentation of the source of the error in terms
of the users problem.
All errors are marked with four asterisks **** at the beginning of a line in the output listing.
As soon as an error is detected, processing will be stopped at the next convenient opportunity. A model will
never be solved after an error has been detected. The only remedy is to fix the error and repeat the run.
Errors are grouped into the three phases of GAMS modeling! compilation, execution and model generation (which
includes the solution that follows). The following three sub-sections discuss these types of errors.
Compilation errors were discussed in some detail in Chapter 2. There is some overlap between the material in
those sections and this one. Several hundred different types of errors can be detected during compilation and
can often be traced back to just one specific symbol in the GAMS input. Most of the errors will be caused by
simple mistakes: forgetting to declare an identifier, putting indices in the wrong order, leaving out a necessary
semicolon, or misspelling a label. For errors that are not caused by mistakes, the explanatory error message text
will help you diagnose the problem and correct it.
When a compilation error is discovered, a $-symbol and error number are printed below the
offending symbol (usually to the right) on a separate line that begins with the four asterisks.
If more than one error is encountered on a line (possibly because the first error caused a series of other spurious
errors) the $-signs may be suppressed and error number squeezed. GAMS will not list more than 10 errors on
any one line.
At the end of the echo print of the program, a list of all error numbers encountered, together with a
description of the probable cause of each error, will be printed. The error messages are
self-explanatory and will not be listed here.
It is worth noting that it is easy to produce a model that does not do what you want it to do, but does not
contain errors in the sense that the term is being used in this section. The best precaution is to check your work
carefully and build in as many automatic consistency checks as possible.
One mistake that may cause confusion is if a GAMS reserved word is used for a label or an identifier. In this
case, it is impossible to provide helpful messages for technical reasons.
94 GAMS Output
In some cases, an error may not be detected until the statement following its occurrence, where it
may produce a number of error conditions whose explanations seem quite silly. Always check
carefully for the cause of the first error is such a group, and look at the previous statement (and
especially for missing semicolons) if nothing seems obvious.
The following example illustrates the general reporting format for compiler errors.
The reporting format for errors found while analyzing solve statements is more complicated than for normal
compilation errors, mainly because many things must be checked. All identifiers referenced must be defined or
assigned, the mathematics in the equations must match the model class, and so on. More elaborate reporting is
required to accurately describe any problems found. The solve statement is only checked if the model has been
found to be error free up to this point. This is not only because the check is comparatively expensive, but also
because many erroneous and confusing messages can be produced while checking a solve in a program containing
other errors.
Solve error messages are reported in two places and in two formats. First, they are shown
immediately below the solve statement with a short text including the name of any offending
identifier and the type of model involved. This will be sufficient in most cases. Second, a longer
message with some hints appears with the rest of the error messages at the end of the compilation.
The example below illustrates how the general reporting format for compiler errors associated with a solve
statement.
1 variables x,y, z ;
2 equations eq1 , eq2;
3
4 eq1.. x**2 - y =e= z ;
5 eq2.. min(x,y) =l= 20 ;
6
7 model silly / all / ;
8 solve silly using lp maximizing z ;
**** $54,51,256
**** THE FOLLOWING LP ERRORS WERE DETECTED IN MODEL SILLY:
**** 54 IN EQUATION EQ1 .. ENDOG OPERANDS FOR **
**** 51 IN EQUATION EQ2 .. ENDOG ARGUMENT(S) IN FUNCTION
9
Error Messages
1 variable x;
2 equation eq1;
3
4 eq1.. x =l= 10 ;
5 x.lo = 10 ;
6 x.up = 5 ;
7 model wrong /eq1/;
8 solve wrong using lp maximizing x ;
9
Some solve statement require the evaluation of nonlinear functions and the computation of derivatives. Since
these calculations are not carried out by GAMS but by other subsystems not under its direct control, errors
associated with these calculations are reported in the solution report. Unless reset with the domlim option the
subsystems will interrupt the solution process if arithmetic exceptions are encountered. They are then reported
on the listing as shown in the following example:
1 variable x, y;
2 equation one;
3
4 one.. y =e= sqrt(10/x);
5 x.l = 10;
6 x.lo = 0;
7
8 model divide / all / ;
9 solve divide maximizing y using nlp;
S O L V E S U M M A R Y
Note that the solver status returned with a value of 5, meaning that the solver has been interrupted because more
than domlim evaluation errors have been encountered. The type of evaluation error and the equation causing the
error are also reported.
If the solver returns an intermediate solution because of evaluation errors, the following solve will still be at-
tempted. The only fatal GAMS error that can be caused by a solver program is the failure to return any solution
at all. If this happens, as mentioned above, all possible information is listed on the GAMS output file and any
solves following will not be attempted.
10.7 Summary
This is the end of teh sequentila discussion of the basic features of the GAMS language. All further chapters are
geared towards more advanced use of GAMS.
11
11.1 Introduction
This chapter deals with the way in which conditional assignments, expressions and equations are made in GAMS.
The index operations already described are very powerful, but it is necessary to allow for exceptions of one sort
or another. For example, heavy trucks may not be able use a particular route because of a weak bridge, or some
sectors in an economy may not produce exportable product. The use of a subset in an indexed expression has
already been shown to provide some ability to handle exceptions.
Numerical expressions can also serve as logical conditions - a result of zero is treated as a logical
value of False, and a non-zero result is treated as a logical value of True.
2*a - 4
This expression results in a logical value of False when a is 2 because the expression numerically evaluates to 0.
For all other values of a, the expression results in a non-zero value, and therefore is equivalent to a logical value
of True.
The following example of a numerical relationship illustrates its use as a logical condition.b
(sqr(a) > a)
This condition evaluates to False if 1 a 1. For all other values of a, this condition evaluates to True. Note
that the same expression can also be written as (sqr(a) gt a).
The truth table generated by these logical operators is given in table 11.1.
Operands Results
a b a and b a or b a xor b not a
0 0 0 0 0 1
0 non-zero 0 1 1 1
non-zero 0 0 1 1 0
non-zero non-zero 1 1 0 0
set i /1*10/
subi(i) /1*3/ ;
The set subi(i) results in a logical value of True for all elements that belong to subi and False for all elements
of i that do not belong to subi.
The use of set membership as a logical condition is an extremely powerful feature of GAMS and while its use
will be illustrated later on in this chapter, its full power becomes clear when considered with the description of
dynamic sets later.
dayofweek = wednesday
dayofweek <> thursday
GAMS follows the convention that the result of a relational operation is zero if the assertion is False,
and one if True.
Consider the following example for illustration,
x = (1 < 2) + (2 < 3)
The expression to the right of the assignment evaluates to 2 since both logical conditions within parenthesis are
true and therefore assume a value of 1. Note that this is different from the assignment below,
x = (1 < 2) or (2 < 3)
Operation Operator
Exponentiation **
Numerical Operators
- Multiplication, Division *, /
- Unary operators - Plus, Minus +, -
- Binary operators - addition, subtraction +, -
Numerical Relationship operators <, <=, =, <>, >=, >
Logical Operators
- not not
- and and
- or, xor or, xor
Note that in the case of operators with the same precedence, the order in which the operator appears in the
expression is used as the precedence criterion, with the order reducing from left to right.
It is always advisable to use parentheses rather than relying on the precedence order of operators. It
prevents errors and makes the intention clear.
x - 5*y and z - 5
is treated equivalent to (x - (5*y)) and (z-5). However, note that the use of parenthesis does make the
expression clearer to understand.
100 Conditional Expressions, Assignments and Equations
The dollar logical conditions cannot contain variables. Variable attributes (like .l and .m) are
permitted however.
The dollar operator is used to model conditional assignments, expressions, and equations. The following subsection
provides an example that will clarify its use. The next section will deal individually with the topic of using dollar
conditions to model conditional assignments, expressions, and equations respectively.
11.3.1 An Example
Consider the following simple condition,
If the condition is not satisfied, no assignment is made. Note that one can read the $ as such that to clarify the
meaning: a, such that b is greater than 1.5, equals 2 .
For nested dollar conditions, all succeeding expressions after the dollar must be enclosed in
parentheses.
u(k)$(s(k)$t(k)) = a(k) ;
11.4 Conditional Assignments 101
K, s(k), t(k), and i are sets, while u(k) and a(i) are parameters. The assignment will be made only for those
members of k that are also members of both s and t. Note the position of the parenthesis in the dollar condition.
The statement above can be rewritten as
To assist with the readability of statements, it is strongly recommended to use the logical and
operator instead of nesting dollar operators.
The effect of the dollar condition is significantly different depending on which side of the assignment
it is in.
In many cases, it may be possible to use either of the two forms of the dollar condition to describe an
assignment. In such a case, clarity of logic should be used as the criterion for choice.
The next two subsections describe the use of the dollar condition on each side of the assignment.
For an assignment statement with a dollar condition on the left-hand side, no assignment is made
unless the logical condition is satisfied. This means that the previous contents of the parameter on
the left will remain unchanged for labels that do not satisfy the condition.
If the parameter on the left-hand side of the assignment has not previously been initialized or
assigned any values, zeroes will be used for any label for which the assignment was suppressed.
rho(i)$(sig(i) ne 0) = (1./sig(i)) - 1. ;
The parameter sig(i) has been previously defined in the model and the statement above is used to calculate
rho(i). The dollar condition on the statement protects against dividing by zero. If any of the values associated
with sig(i) turn out to be zero, no assignment is made and the previous values of rho(i) remain. As it happens,
rho(i) was previously not initialized, and therefore all the labels for which sig(i) is 0 will result in a value of 0.
Now recall the convention, explained in Section 11.2.1 that non zero implies True and zero implies False. The
assignment above could therefore be written as
rho(i)$sig(i) = (1./sig(i)) - 1. ;
For an assignment statement with a dollar condition on the right hand side, an assignment is always
made. If the logical condition is not satisfied, the corresponding term that the logical dollar
condition is operating on evaluates to 0.
Consider the following example, which is a slight modification to the one described in Section 11.3.1,
102 Conditional Expressions, Assignments and Equations
Therefore an if-then-else type of construct is implied, but the else operation is predefined and never made explicit.
Notice that the statement in the illustrative example above can be re-written with an explicit if-then-else and
equivalent meaning as
This use of this feature is more apparent for instances when an else condition needs to be made explicit. Consider
the next example adapted from [FERTD]. The set i is the set of plants, and are calculating mur(i), the cost of
transporting imported raw materials. In some cases a barge trip must be followed by a road trip because the
plant is not alongside the river and we must combine the separate costs. The assignment is:
This means that if the entry in the distance table is not zero, then the cost of shipping using that link, which has a
fixed and a variable components, is added to the total cost,. If there is no distance entry, there is no contribution
to the cost, presumably because that mode is not used.
where shipped is the amount of material shipped from i to j, and totcost is the total cost of all shipment. The
equation above can be written as,
Index j appears separately from i in congestfac(j). The equation then needs to be simplified as,
Note that the presence of j separately in the indexed expression necessitated the use of ij(i,j) rather than ij.
u(k)$s(k) = a(k) ;
where k and s(k) are sets, while u and a are parameters. This can be rewritten as,
u(s) = a(s) ;
Note that the assignment has been filtered through the conditionality without the use of the dollar operator. This
is a cleaner and more understable representation of the assignment. This feature gets more useful when dealing
with tuples (sets with multiple indices).
Consider the following example for calculating the travel cost for a fictional parcel delivery service between
collection sites (i) and regional transportation hubs (j),
set i /miami,boston,chicago,houston,sandiego,phoenix,baltimore/
j /newyork,detroit,losangeles,atlanta/ ;
set ij(i,j) /
boston.newyork
baltimore.newyork
miami.atlanta
houston.atlanta
chicago.detroit
sandiego.losangeles
phoenix.losangeles / ;
The set ij denotes the regional transportation hub for each collection site. Factor is the cost estimate per unit
mile. The cost of transporting parcels (shipcost) from a local collection site (i) to a regional hub(j) is then
provided by the following assignment,
shipcost(i,j)$ij(i,j) = factor*distance(i,j) ;
Note that i and j do not appear separately in the assignment for shipcost. The assignment can then be simply
written as,
shipcost(ij) = factor*distance(ij) ;
If i or j appear separately in any assignment, the above simplification cannot be made. For example, consider
that travelcost depended not only on factor and the distance between collection sites and regional hubs but also
the load on the regional hub.
Parameter congestfac(j) /
newyork 1.5
detroit 0.7
losangeles 1.2
atlanta 0.9/ ;
Congestfac is a parameter used to model the congestion at each regional hub. The unit cost of shipment is then
computed as follows:
shipcost(i,j)$ij(i,j) = factor*congestfac(j)*distance(i,j) ;
shipcost(ij) = factor*congestfac(j)*distance(ij) ;
The above representation has the index j on the right hand side, but not on the left hand side. As explained
before, GAMS will flag this assignment as an error. However, the following representation will work:
shipcost(ij(i,j)) = factor*congestfac(j)*distance(ij) ;
In the above assignment ij is specifically denoted as a tuple of i and j which then appear on the left hand side.
A common use of dollar controlled index operations is where the control is itself a set. This
importance of this concept will become apparent with the discussion of dynamic sets.
A set has been used to define the mapping between mines and ports in Chapter 4. Another typical example is
a set-to-set mapping defining the relationship between states and regions, used for aggregating data obtained by
state to the models requirements (by region).
sets r / west,east /
s / florida,texas,vermont,maine) /
corr(r,s) / north.(vermont,maine)
south.(florida,texas) /
parameter y(r)
income (s) "income of each state"
/ florida 4.5, vermont 4.2
texas 6.4, maine 4.1 / ;
The set corr provides a correspondence to show which states belong to which regions. The parameter income is
the income of each state. Y(r) can be calculated with this assignment statement:
For each region r, the summation over s is only over those pairs of (r,s) for which corr(r,s) exists. Concep-
tually, set existence is analogous to the Boolean value True or the arithmetic value not zero. The effect is that
only the contributions for vermont and maine are included in the total for north, and south includes only
texas and florida.
Note that the summation above can also be written as sum(s, income(s)$corr(r,s)) but this form is not as
easy to read as controlling the index of summation.
where shipped is the amount of material shipped from i to j, and totcost is the total cost of all shipment. The
equation above can be written as,
Index j appears separately from i in congestfac(j). The equation then needs to be simplified as,
Note that the presence of j separately in the indexed expression necessitated the use of ij(i,j) rather than ij.
The term is added to the right hand side of the equation only for those elements of i that belong to t(i).
Controlling indexing operations using the dollar condition can also be done as with any assignment. Consider the
following supply balance (sb) equation from [GTM],
The purpose of the dollar control over the domain of definition of equations is to restrict the number
of constraints generated to less than that implied by the domain of the defining sets.
cc(m,i)$mpos(m,i)..
sum(p$ppos(p,i), b(m,p)*z(p,i)) =l= util*k(m,i);
106 Conditional Expressions, Assignments and Equations
Cc is a capacity constraint defined for all units (m) and locations (i).
Not all types of units exist at all locations, however, and the mapping set mpos(m,i) is used to restrict the number
of constraints actually generated. The control of the summation over p with ppos(p,i) is an additional one, and
is required because not all processes (p) are possible at all locations (i).
parameter bigM(i,j) ;
variable shipped(i,j) ;
binary variable bin(i,j) ;
equation logical(i,j) ;
logical(i,j)$ij(i,j).. shipped(i,j) =l= bigM(i,j)*bin(i,j) ;
The equation logical relates the continuous variable shipped(i,j) to the binary variable bin(i,j). This can
be simplified as follows:
Note that if the right hand side of the equation contained any term that was indexed over i or j separately, then
the equation logical(i,j)$ij(i,j) would have to be simplified as logical(ij(i,j)).
12
Dynamic Sets
12.1 Introduction
All the sets that have been discussed so far had their membership declared as the set itself was declared, and
the membership was never changed. In this chapter we will discuss changing the membership of sets. A set
whose membership can change is called a dynamic set to contrast it with a static set whose membership will never
change. The distinction is important and will be discussed in detail in this chapter. This is a topic that has been
avoided until now because of a potential confusion for new users. Advanced Users will, however, find it useful.
Set name is the internal name of the set (also called an identifier) in GAMS. Yes and no are keywords used in
GAMS to denote membership or absence respectively from the assigned set.
The most important principle to follow is that a dynamic set should always be domain checked at
declaration time to be a subset of a static set (or sets).
It is of course possible to use dynamic sets that are not domain checked, and this provides additional
power, flexibility, lack of intelligibility, and danger. Any label is legal as long as the dimensionally,
once established, is preserved.
Note that the sets subitem1 and subitem2 are declared like any other set. The two sets become dynamic because
of assignments. They are also domain checked: the only members they will ever be able to have must also be
members of item. And item is a static set and henceforth its membership is frozen. The first two assignments
each add one new element to subitem1. The third is an example of the familiar indexed assignment: subitem2 is
assigned all the members of item. The output caused by the display statement, that will reveal the membership
of the sets, is shown below for verification.
The elements are displayed in the order specified in the declaration of item.
Dynamic sets, like static sets, can have up to 10 dimensions. The following example illustrates assignments for
multi-dimensional sets.
supply(pencil,bic) = yes ;
supply(pen,sup) = yes ;
All the mechanisms using asterisks and parenthesized lists that we introduced in the discussion on static sets in
chapter 4 are available for dynamic sets as well.
One can make an assignment over the domain of a dynamic set because dynamic sets are known to be proper
subsets of static sets. This is not the same as doing domain checking using a dynamic set.
The following example, adapted from the Section 12.2.2 illustrates the use of dynamic sets as domains:
subitem1(item) = no
subitem1(subitem2) = yes;
The first assignment ensures that subitem1 is empty. Note that this can also be done with parameters. For
example,
parameter inventory(item) ;
inventory(subitem1) = 25 ;
12.3 Using Dollar Controls with Dynamic Sets 109
The trick is to declare the equation over the entire domain but define it over the dynamic set.
To repeat the important point: the equation is declared over allr but referenced over r. Then arbitrary assign-
ments can be made to r within the membership of allr.
12.3.1 Assignments
Dynamic sets can be used inside dollar conditions within assignments defining other dynamic sets or parameters.
As an illustration of its use in defining other dynamic sets, the two statements in the example from Section 12.2.4
can be written with equivalent effect as
subitem1(item) = yes$subitem2(item) ;
The value used in the implied else that goes with dollar on the right is no in a set assignment,
rather than zero which is used with normal data.
The second example from Section 12.2.4 can be rewritten as follows to illustrate the use of dynamic sets in defining
parameters,
inventory(item)$subitem1(item) = 25 ;1
110 Dynamic Sets
This example has been shown only for illustration. Note that the second statement above can also be rewritten
tersely as
totinv = sum(subitem1,inventory(subitem1)) ;
This is not always possible. Consider the following artificially created example,
parameter totsales(dep) ;
totsales(dep) = sum(item$supply(item,bic), sales(dep,item)) ;
The assignment above is used to find the total sales of all departments that sell items supplied by bic. Note that
the dynamic set is used to limit the domain of summation to those for which supply(item,bic) is true.
12.3.3 Equations
Dynamic sets can be used inside dollar conditions in equations both as part of the equation algebra, or while
defining the domain of the equation. The first case is similar to the case of assignments discussed in Section 12.3.1.
The latter case is used to restrict the equation over the domain of a dynamic set. The equation defined in the
example from Section 12.2.5 can be rewritten with equivalent effect as follows,
The domain of definition of equation prodbal is restricted to those elements that belong to the dynamic set r.
inventory(item)$subitem1(item) = 25 ;
prodbal(allr)$r(allr).. activity(allr)*price =e= revenue(allr) ;
inventory(subitem1) = 25 ;
prodbal(r).. activity(r)*price =e= revenue(r) ;
The membership of subitem3 is set equal to all the elements of subitem1 and all the elements of subitem2. The
operation above is equivalent to the following longer way of representation,
The membership of subitem3 is set equal to only those present in both subitem1 and subitem2. The operation
above is equivalent to the following longer way of representation,
The membership of subitem3 is set equal to all those in item but not in subitem1. The operation above is
equivalent to the following longer way of representation,
subitem3(item)=yes; subitem3(subitem1)=no;
The membership of subitem3 is set equal to all elements that are members of subitem1 but subitem2. The
operation above is equivalent to the following longer way of representation,
subitem3(item)=yes$(subitem1(item)); subitem3(subitem2)=no;
12.5 Summary
The purpose of set assignments is to make calculations based on given data (the static sets) for use in exception
handling. It is one more example of the principle of entering a small amount of data and building a model up
from the most elemental information.
112 Dynamic Sets
13
13.1 Introduction
In our original discussion of sets in Chapter 4, we said that unless there is a special need to do things differently,
a one-dimensional set should be regarded as an unordered collection of labels. In this chapter we will discuss
special features that can be used when you need to be able to deal with a set as if it were a sequence.
For example, in economic models that explicitly represent conditions in different time periods, it is necessary to
refer to the next or previous time period, because there must be links between the periods. As another example,
stocks of capital are normally tracked through such models by equations of the form stocks at the end of period
n are equal to stocks at the end of period n 1 plus net gains during period n. Location problems, where the
formulation may require a representation of contiguous areas, as in a grid representation of a city, and scheduling
problems are other classes of problems in which sets must also have the properties of sequences.
Models involving sequences of time periods are often called dynamic models, because they describe
how conditions change over time. This use of the word dynamic unfortunately has a different
meaning from that used in connection with sets, but this is unavoidable.
Ordered sets must be static sets. In other words, no order is possible for dynamic sets.
GAMS maintains one list of unique elements - the labels that are used as elements in one or more
sets. The order of the elements in any one set is the same as the order of those elements in that
unique element list. This means that the order of a set may not be what it appears to be if some of
the labels were used in an earlier definition.
The map of your labels in the GAMS order can be seen by putting the compiler directive
$onuellist somewhere before the first set declaration.
A good rule of thumb is that if the labels in a set one wants to be ordered have not been used
already, then they will be ordered.
The map is shown with the other compiler maps after the listing of your program. In the example below we show
ordered and unordered sets and the map showing the order. The input is:
114 Sets as Sequences: Ordered Sets
$onuellist
set t1 / 1987, 1988, 1989, 1990, 1991 /
t2 / 1983, 1984, 1985, 1986, 1987 /
t3 / 1987, 1989, 1991, 1983, 1985 / ;
The map below shows the entry order (the important one) and the sorted order, obtained by sorting the labels
into dictionary order. The single digits on the left are the sequence numbers of the first label on that line.
G e n e r a l A l g e b r a i c M o d e l i n g S y s t e m
Unique Element Listing
A set can always be made ordered by moving its declaration closer to the beginning of the program. With these
restrictions in mind, we move on the operations that are used in dealing with sets as sequences.
As a result of the statements above, the value of val(1985) will be 1, val(1986) will be 2 and so on.
A common use of ord is in setting up vectors that represent quantities growing in some analytically specified way.
For example, suppose a country has 56 million people in the base period and population is growing at the rate of
1.5 percent per year. Then the population in succeeding years can be calculated by using:
It is often useful to simulate general matrix operations in GAMS. The first index on a two dimensional parameter
can conveniently represent the rows, and the second the columns, and order is necessary. The example below
shows how to set the upper triangle of a matrix equal to the row index plus the column index, and the diagonal
and lower triangle to zero.
which fixes the variable for the last member only: no assignment is made for other members of t. The advantage
of this way of fixing c is that the membership of t can be changed safely and this statement will always fix c for
the last one.
The difference between these two types of operators involves the handling of endpoints in the sequence. GAMS
provides some built in facilities to deal with this issue, but in any work involving sequences the user must think
carefully about the treatment of endpoints, and all models will need special exception handling logic to deal with
them.
In the linear case, the members of the set that are endpoints are left hanging. In other words, there are no
members preceding the first member or following the last one. This may cause the use of non-existent elements.
The next section will describe how this is handled in GAMS. This form of the lag and lead operators is useful for
modeling time periods that do not repeat. A set of years (say 1990 to 1997) is an example. The operators are +
and -.
GAMS is able to distinguish linear lag and lead operators (+,-) from arithmetic operators by
context.
In the circular case, the first and last members of the set are assumed to be adjacent, so as to form a circular
sequence of members. The notion is that first - 1 is a reference to last and last + 2 is the same as first
+ 1 and so on. All references and assignments are defined. This is useful for modeling time periods that repeat,
such as months of the year or hours in the day. It is quite natural to think of January as the month following
December. Agricultural farm budget models and workforce scheduling models are examples of applications where
circular leads occur naturally. The operators are ++ and --.
The next two sections will describe the use of these lag and lead operators in assignment statements and in
equations respectively.
involves the definition of a domain of the assignment. The concepts behind reference and assignment are equally
valid for the linear and circular forms of the lag and lead operator. However, the importance of the distinction
between reference and assignment is not pronounced for circular lag and lead operators because non-existent
elements are not used in this case.
A reference to a non-existent element causes the default value (zero in this case) to be used, whereas
an attempt to assign to a non-existent element results in no assignment being made.
The next two sub-sections provide examples illustrating the use of the linear form of the lag and lead operators for
reference and assignment. Section 13.5.3 will illustrate the use of the circular form of the lag and lead operator.
The option statement suppresses the decimal places from the display. The results are shown below.
---- 6 PARAMETER A
Y-1987 1987, Y-1988 1988, Y-1989 1989, Y-1990 1990, Y-1991 1991
---- 6 PARAMETER B
Y-1988 1987, Y-1989 1988, Y-1990 1989, Y-1991 1990
For a, as expected, the values 1987, 1988 up to 1991 are obtained corresponding to the labels y-1987, y-1988
and so on. b is initialized to 1.
For b, the assignment is done over all members of t, and for each, the value of a from the previous period is
assigned to the current member of b. If no previous period, as with y-1987, zero is used, and b(y-1987)
becomes zero, replacing the previous value of 1.
---- 6 PARAMETER A
Y-1987 1987, Y-1988 1988, Y-1989 1989, Y-1990 1990, Y-1991 1991
---- 6 PARAMETER C
Y-1987 -1, Y-1988 -1, Y-1989 1987, Y-1990 1988, Y-1991 1989
The assignment to a is explained in Section 13.5.1. The assignment to c is different. It is best to spell it out in
words. For each member of t in sequence, find the member of c associated with t+2. If it exists, replace its value
with that of a(t). If not (as with y-1990 and y-1991) make no assignment. The first member of t is y+1987,
13.6 Lags and Leads in Equations 117
and therefore the first assignment is made to c(y-1989) which takes the value of a(y-1987), viz., 1987. No
assignments at all are made to c(y-1987) or c(y-1988): these two retain their previous values of 1.
The lag (or lead) value does not have to be an explicit constant: it can be arbitrary expression, provided that it
evaluates to an integer. If it does not, error messages will be produced. A negative result causes a switch in sense
(from lag to lead, for example). The following is guaranteed to set d(t) to all zero:
The parameter lagval2 is used for reference while lagval1 if used for assignment. Notice that the case of
circular lag and lead operators does not lead to any non-existent elements. The difference between reference and
assignment is therefore not important. Note that the following two statements from the example above,
lagval2(s) = val(s--2) ;
leadval(s++1) = val(s) ;
are equivalent to
lagval2(s++2) = val(s) ;
leadval(s) = val(s--1) ;
The use of reference and assignment have been reversed with no difference in effect.
The next two sub-sections provide examples illustrating the use of the linear form of the lag and lead operators
in equations for reference and to modify the domain of its definition. Section 13.5.3 will illustrate the use of the
circular form of the lag and lead operator in equations.
118 Sets as Sequences: Ordered Sets
The declaration of t is included, as are a couple of dynamic sets that are used to handle the first and last periods
(terminal conditions) in a clean way.
The interesting equation is kk, the capital balance. The set t contains members 1990 to 2000, and so there will
be a capital stock constraint for 1991 to 2000. Spelling out the constraint for 1991,
The lead operator on the domain of definition has restricted the number of constraints generated so that there are
no references to non-existent variables: the generated problem will have 10 kk constraints defining the relationship
between the 11 k capital values.
The other interesting point in the [RAMSEY] excerpt is that the constraint tc is explicitly defined only for the
final period because of the assignment to the set tlast. Notice the use of dynamic sets to control the domain of
the two equations. The set tfirst is also used in other parts of the model to set initial conditions, particularly
the capital stock in the first period, k(1990).
The dollar condition will cause one of the individual equations to be suppressed.
However, note that using lags and leads in the equation domain will always cause one or more individual equations
to be suppressed, and this may not be desirable in every case. Consider the following modified set of constraints
to the one discussed in the previous example. It is expressed with the lag and lead operators being used to control
the domain of the equation definition.
Here, the important boundary is the one at the beginning of the set rather than at the end. This can be expressed
more compactly as
In general, the choice between using lag and lead operators as reference or in domain control is often a matter of
taste.
13.7 Summary 119
equation matbal(s) ;
In this example, four individual examples are generated. They are listed below.
13.7 Summary
This chapter introduced the concept of ordering in sets. All the features in GAMS that dealt with this issue
including the ord and card functions, as well as the linear and circular forms of the lag and lead operators were
described in detail.
120 Sets as Sequences: Ordered Sets
14
14.1 Introduction
In this chapter we will provide more detail about display statements, including the controls that a user has
over the layout and appearance of the output. These controls are a compromise to provide some flexibility. The
display statement will not provide a publication quality reporting function, but is instead aimed for functionality
that is easy to use, and provides graceful defaults. The execution of the display statement allows the data to be
written into the listing file only.
Ident-ref means the name without domain lists or driving indices of a set or parameter, or a sub-field of an
equation or variable. The identifier references and the text can be mixed and matched in any order, and the
whole statement can be continued over several lines.
The output produced by a display consists of labels and data. For sets, the character string yes (indicating
existence) is used instead of values.
Only the non-default values are displayed for all data types.
The default value is generally zero, except for the .lo and .up subtypes of variables and equations. The
default values for these are shown in table 14.1.
14.3 An Example
An example of a display statement is given below.
The resulting listing file will contain the following section that corresponds to the display statement.
122 The Display Statement
.lo .up
Variable
positive 0 +INF
free -INF +INF
negative -INF 0
integer 0 100
binary 0 1
Equation
=g= 0 +INF
=n= -INF +INF
=c= 0 +INF
=l= -INF 0
=e= 0 0
S1 0.330, S3 0.670
T5 T7
S1 0.109 0.221
S3 0.221 0.449
Note that the only the non-zero values are displayed. In the case of multi-dimensional identifiers, the data is
reported in a tabular form that is easy to read.
For 7 to 10 indices, the natural progression is followed. The labels vary slowest for the first index position, and
14.5 Display Controls 123
quickest for the highest. Within each index position the order is the GAMS entry order of the labels.
The order of the indices is always as in the declaration statement for the symbol. One can declare them in the
order that is found appealing, or make an assignment to a new identifier with a different order.
The only way to change the order in which the labels for each index position appear on display
output is to change the order of appearance of the labels in the GAMS program. This is most easily
done by declaring a set whose only purpose is to list all the labels in the order that is needed. Make
this set the very first declaration in the GAMS program.
14.4.1 Example
Consider the following example. X has four dimensions or index positions. It is initialized using parameter format
and then displayed as shown below:
INDEX 1 = first
i ii
one .b 5.636
two .a -2.939 0.029
two .b 10.346
three.b -6.316
INDEX 1 = second
i ii
Notice that there are two sub-tables, one for each label in the first index position. Note that the zero in the
list for x(first,one,a,ii) has vanished, since zero values are suppressed in each sub-table separately.
The order of the labels is not the same as in the input data list.
GAMS has rounded or converted numbers to E-format where necessary and the output is as follows:
INDEX 1 = first
i ii
one .b 5.6
two .a -2.9 2.873000E-2
two .b 10.3
three.b -6.3
INDEX 1 = second
i ii
option ident:d-value:
Ident represent the name of a parameter, variable or equation, and d-value must be (as before) in the range
0 and 8 . Exactly d-value places of decimals will be shown on all displays of ident that follow. This form can
be extended to control layout of the data. The general form is:
option:ident:d-value:r-value:c-value ;
Here r-value means the number of index positions that are combined to form the row label and c-value means
the number on the column headers.
The example discussed in the previous section is further extended in order to illustrate the local display control.
i ii
14.5 Display Controls 125
Five places of decimals are shown, and three labels are used to mark the rows and one on the column. Since this
is a four-dimensional structure, there are no remaining indices to be used as sub-table labels (on the plane), and
we now have the results in one piece. The option statement is checked for consistency against the dimensionality
of the identifier, and error messages issued if necessary. Here is an example that puts two indices on each of the
row and column labels, and retains five decimal places:
The output is :
This output nicely illustrates the label order used. The first index varies the slowest, the last the fastest, and
each one runs from beginning to end before the next one to the left advances. This ordering scheme is also used
on equation and column lists and on the solution report, all produced by the solve statement.
126 The Display Statement
15
15.1 Introduction
In this chapter, the put writing facility of the GAMS language is introduced. The purpose of this writing facility
is to output individual items under format control onto different files. Unlike the display statement, the entire
set of values for indexed identifiers cannot be output using a single put statement (identifiers are the names
given to data entities such as the names for parameters, sets, variables, equations, models, etc). While its
structure is more complex and requires more programming than is required for the display statement, there is
much greater flexibility and control over the output of individual items.
In this chapter, the working of the put writing facility is described as well as the syntax for accessing files and
globally formatting documents using file suffixes for various attributes of a file. The put writing facility enables
one to generate structured documents using information that is stored by the GAMS system. This information is
available using numerous suffixes connected with identifiers, models, and the system. Formatting of the document
can be facilitated by the use of file suffixes and control characters.
The put writing facility generates documents automatically when GAMS is executed. A document is written to
an external file sequentially, a single page at a time. The current page is stored in a buffer, which is automatically
written to an external file whenever the page length attribute is exceeded. Consequently, the put writing facility
only has control of the current page and does not have the ability to go back into the file to alter former pages
of the document. However, while a particular page is current, information placed on it can be overwritten or
removed at will.
The basic structure of the put writing facility in its simplest form is:
file fname(s);
put fname;
put item(s);
where fname represents the name used inside the GAMS model to refer to an external file. Items are any type
of output such as explanatory text, labels, parameters, variable or equation values. In the basic structure shown
above, the first line defines the one or more files which you intend to write to. The second line assigns one of
these defined files as the current file, that is the file to be written to. Lastly, the third line represents the actual
writing of output items to the current file.
128 The Put Writing Facility
15.3 An Example
It is instructive to use a small example to introduce the basics of the put writing facility. The example will be
based on the transportation model [TRNSPORT]. The following program segment could be placed at the end of
the transportation model to create a report:
put results;
put Transportation Model Results// ;
loop((i,j), put i.tl, @12, j.tl, @24, x.l(i,j):8:4 /);
In the first line, the internal file names factors and results are defined and connected to the external file names
factors.dat and results.dat. These internal file names are used inside the model to reference files, which are
external to the model. The second line of this example assigns the file factors.dat as the current file, that is
the file which is currently available to be written to.
In the third line of the example, writing to the document begins using a put statement with the textual item
Transportation Model Factors. Notice that the text is quoted. The slashes following the quoted text
represent carriage returns. The example continues with another textual item followed by the scalar f. Notice
that these output items are separated with commas. Blanks, commas, and slashes serve as delimiters for separating
different output items. As mentioned above, the slash is used as a carriage return. Commas and blank spaces
serve as item delimiters. These delimiters leave the cursor at the next column position in the document following
the last item written. In most cases, the blank and the comma can be used interchangeably; however, the comma
is the stronger form and will eliminate any ambiguities.
In the fifth line of the program above, the cursor is repositioned to the first column of the sixth row of the output
file where another textual item is written. The cursor control characters # and @ serve to reposition the cursor
to a specific row or column as designated by the row or column number following the cursor control character.
Lastly, the put statement is terminated with a semicolon.
Next, the parameters a and b are written along with their corresponding set labels. Only one element of the
index set can be written using a put. To write the entire contents of the parameters a and b, the put statement
is embedded inside a loop which iterates over the index set. In the example above, the set element labels are
identified using their set identifier and the suffix .tl. As can be seen, the set element labels are located starting
in the third column and the parameter a at column 15. The example continues with the display of another quoted
textual item followed by the parameter b. When executed, the factors.dat file will look like:
Plant capacity
seattle 350.00
san-diego 600.00
Market demand
new-york 325.00
chicago 300.00
topeka 275.00
This output has been formatted using the default file format values. The methods to change these defaults will
be described later in this chapter.
In the last two lines of the example, the file results.dat is made current and the values associated with the
variable x along with their corresponding set element index labels are written line by line. The output results of
15.4 Output Files 129
the variable x are formatted by specifying a field width of eight spaces with four of these spaces reserved for the
decimal. Notice that the local formatting options are delimited with colons. The results.dat file will look like:
With just this brief introduction to the put writing facility, it is easy to envision its many uses such as report
writing, providing output to a file for use by another computer program, or simply the display of intermediate
calculations. But, the surface of the put writing facility has just barely been scratched. In the sections that follow,
the many features and structure of the put writing facility are described in more detail, along with examples.
where file is the keyword used to define files. Fname is the internal file name and is used inside the GAMS model
to refer to an external file. External files are the actual files that output is written to. During file declaration,
the external file name and explanatory text are optional. When the external file name is omitted, GAMS will
provide a system specific default external file name, often fname.put. Note that multiple files can be defined
using a single file statement. Consider the following example:
file class1
class2 this defines a specific external file /report.txt/
con this defines access to the console (screen) for PC systems;
The first output file is recognized in the model by the name class1 and corresponds to the default file class1.put
for a PC system. The second output file is recognized in the model by the name class2 and it corresponds to
the defined external file report.txt. Lastly, the special internal file name con is defined to write output to the
console (screen) for a PC systems. Writing to the screen can be useful to advise the user of various aspects of the
model during the models execution.
As indicated by this syntax, multiple files can be sequentially written using a single put statement. Note that
only one file is current at a time. After the output items following an internal file name are written, the current
file is reassigned based on the next internal file name in the statement. The last internal file name used in a put
statement remains as the current file until a subsequent put statement uses an internal file name.
130 The Put Writing Facility
where myfile is the internal name of the file to be closed, and item(s) are the final entries into the file before it
is closed. If the internal file name is omitted from the putclose statement, the current put file is closed. Note
that after using the putclose command, the file does not have to be redefined in order to use it again. Simply
make the file current and use put statements as would be done normally. Of course, the existing file will either
be overwritten or appended to depending on the value of the append file suffix.
One application where this is useful is to write the solver option file from within the GAMS model.
Option file statements can be written using put and the file closed with a putclose prior to the solve
statement. This makes the option file available for use by the solver.
The following example shows the creation and closing of an option file for the MINOS solver:
This program segment would be placed inside the GAMS model prior to the solve statement.
class2.ap = 1;
Any items put into report.txt will from that point on be added to the end of the existing file contents. If the
file had not existed, the file would be created.
4 Formatted output; Non-numeric output is quoted, and each item is delimited with a blank space.
5 Formatted output; Non-numeric output is quoted, and each item is delimited with commas.
6 Formatted output; Non-numeric output is quoted, and each item is delimited with tabs.
page size (.ps) Used to specify the number of rows (lines) which can be placed on a page of the document.
Can be reset by the user at any place in the program. However, an error will result if set to a value less
than the number of rows which have already been written to the current page. Maximum value is 130.
The default value is 60
page width (.pw) Used to specify the number of columns (characters) which can be placed on a single row of
the page. Can be reset by the user at any place in the program. However, an error will result if set to a
value less than the number of rows or columns which have already been written to the current page. The
default value is 132.
top margin (.tm) Number of blank lines to be placed at the top margin of the page. These lines are in
addition to the number of lines specified in the .ps file suffix. Default value is 0.
bottom margin (.bm) Number of blank lines to be placed in the bottom margin of the page. These lines are
in addition to the number of lines specified in the .ps file suffix. This is functional with .pc option 0 only.
Default value is 0.
alphabetic case (.case) Used to specify the case in which alphabetic characters are displayed in the output
file.
0 (default) Causes mixed case to be displayed.
1 Causes the output to be displayed in upper case regardless of the case used for the input.
To illustrate the use of these file suffixes, the following example involves formatting report.txt so that the pages
are 72 spaces wide with 58 lines of output, an additional top margin of 6 lines, using ASCII page control characters
(inserted every 64 lines), and with the output displayed in upper case.
Using a value of 4, 5, or 6 for the print control suffix (.pc) will cause data to be squeezed and
therefore will ignore spacing information provided by the user through the @ character. However,
these values can be used to pass data on to be read by spreadsheets.
Title Block
Header Block
Window
132 The Put Writing Facility
In the illustrative example described in Section 15.3, all the data was written to the window. A title block might
have been included, if more elaboration were needed, to provide the model name along with the page number.
In addition, a header block might have been used to display a disclaimer or an instruction, which we wanted
consistently, repeated on every page. Once this information is placed in the title or header blocks, it is displayed
on each page thereafter unless modified. This could be especially useful for a long document covering many pages.
The size of any area within a given page is based entirely on the number of lines put into it. Note that the total
number of lines for all areas must fit within the specified page size. If the total number of lines written to the
title and header block equals or exceeds the page size, an overflow error will be displayed in the program listing.
When this occurs, this means there is no room remaining on the page to write to the window.
As mentioned above, the syntax for writing an output item to any of the three possible writing areas of the page
is basically the same, the only difference being the choice of put keyword. This is illustrated by writing to the
title block of our report.dat file:
In this case, the text GAMS Put Example has been placed in the first column of the first row of the title block.
Any subsequent pages in the report.dat file will now start with this information.
If the title block was modified or the header block was started after the window of the current page
has been written to, these modifications would appear in the next page and not the current page.
15.6.2 Paging
Paging occurs automatically whenever a page is full. However, note that the window must be used in order for
the page to be written to the output file. When a page has no output in its window, the page is not written to
file regardless of whether there are output items in the title or header blocks. To force a page that has an empty
window out to file, simply write something innocuous to the window such as:
put ;
Now the window of the page has been initiated and it will be written.
As an illustration, consider the example discussed in the previous section. One can add page numbers to the title
of the report file by modifying the puttl statement to read
This causes the word page followed by the page number to appear on the title of every page starting at column 65.
The methods for identifying and using each of these different types of output items are described in the following
sub-sections.
set element text (.te(index)) Displays the text associated with an element of a set. Notice that the .te
suffix requires a driving index. This driving index controls the set, which will be displayed and does not
necessarily have to be the same as the controlled set. Often a subset of indices of the controlled set is
used.
text fill (.tf) Used to control the display of missing text for set elements.
0 no fill
1 fill existing only
2 (default) fill always
subset of sites
i3 San Francisco
i4 Los Angeles
i5 i5
In this example, the symbol text for the identifier of the subset j is written first. This is followed with the labels
for the subset j and the associated element text found in its domain, that is, the set i. Notice the driving set j
is used for the element text specification of the set i. Since there was no set element text associated with the i5
element of set i, the set element label was displayed again. By placing the following before the last line:
out.tf = 0;
The missing element text is now no longer replaced with the label text. The resulting file out.put file would now
look like:
subset of sites
i3 San Francisco
i4 Los Angeles
i5
subset of sites
i1 NO Seattle
i2 NO Portland
i3 YES San Francisco
i4 YES Los Angeles
i5 YES
The second columns represents whether the element belongs to set j or not.
The field width is specified with the number of spaces to be allocated to the field. Variable length field widths
are possible by using a suffix value of 0. This forces the field width to match the exact size of the item being
displayed. If a textual output item does not fit within the specified field, truncation occurs to the right. For
numeric output items, the decimal portion of a number is rounded or scientific notation used to fit the number
within the given field. If a number is still too large, asterisks replace the value in the output file.
As an example, to set the global numeric field width to four spaces from its default of 12 in the file out.put, we
would use the following statement:
out.nw = 4;
136 The Put Writing Facility
item:{<>}width:decimals;
The item is followed by a justification symbol, the field width, and the number of decimals to be displayed.
The specification of the number of decimals is only valid for numeric output. The following local justification
symbols are applicable:
> right justified
< left justified
<> center justified
Omitting any of the components causes their corresponding global format settings to be used. As with global
formatting, when the field width is given a value of 0, the field width is variable in size. The item, width, and
decimals are delimited with colons as shown above. The use of the local format feature as well as the inclusion
any of the components for justification, field width, or the number of decimals is entirely optional.
The following example shows some examples of the local formatting feature:
numeric round format (.nr) Allows one to display a numeric value in scientific notation, which would
otherwise be displayed as zero because of being smaller than the number of decimals allowed by the .nd
suffix. This situation occurs when a number is smaller than the .nd specification, but is larger than the
zero tolerance level set by .nz. In many situations, it is important to know that these small values exist.
The default is 1.
0 displayed in F or E format
1 rounded to fit fields
2 displayed in scientific notation
numeric zero tolerance (.nz) Sets the tolerance level for which a number will be rounded to zero for display
purposes. When it is set equal to zero, rounding is determined by the field width. Default value is 1.0e-5.
The maximum size of a displayed number must fit within 20 spaces using at most 10 significant digits. The
remaining 10 spaces are used for the sign, exponential notation, or padding with zeros.
15.13 Cursor Control 137
The following illustrative example shows the results of different combinations of these numeric file suffixes. The
example uses five combinations of the numeric file suffixes .nd, .nz, .nr, and .nw. Four number values, each of
which is shifted by three decimal places from its predecessor, are used with these suffix combinations. The com-
binations are chosen to show various format results when these suffix values are used together in put statements:
For readability, the numeric values have purposely been made left justified using the .nj suffix since the numeric
field width is changed as the model goes through the suffix combinations. The following is the resulting file
out.put, which shows the value/suffix combinations:
Notice that in comb1, the display of values switch to exponential notation when a value becomes smaller than the
number of decimal places allowed. This is triggered by the suffix .nr being set to zero. Of particular interest is
value3 for comb2 and comb3. Value3 is greater than the zero tolerance level in .nz, but smaller than the number
of decimals allowed by .nd. In comb2, since .nr is set to zero, the value is displayed in exponential format. In
comb3, .nr is set to 1, so this small value is rounded to 0. In comb5, value1 is rounded to an integer because of
.nd being set to 0.
The convention for updating the values stored for the .cc suffix is that it are updated at the
conclusion of a put statement. Consequently, the .cc value remains constant throughout the writing
of items for the next put statement, even if multiple items are displayed.
The following example illustrates the updating of the cursor control suffixes and the use of cursor control charac-
ters. The example is trivial but instructive:
out.cc = 1 x
y
z out.cc = 23
Initially, the scalar lmarg is set to a specific value to use as an alignment tab. Symbols which hold common
alignment values such as margins or tabs are often useful for large structured documents. The first put statement
uses the current column cursor control character to relocate the cursor. In this example, the cursor is moved to
column 8where out.cc and its value is displayed.
The second put statement illustrates the updating of the cursor control suffixes by writing the letters x, y, and
z on three different lines. Each is preceded by the cursor being moved to the out.cc value. Initially, the value
for the cursor control suffice is 20. Since a single put statement is used for these three items, the out.cc value
remains constant and consequently the letters end up in the same column. Following this put statement, the
out.cc value is updated to 23, which is the location of the cursor at the end of the second put statement (note
the additional blank spaces displayed with the letter z).
The convention for updating the values stored for the .cr suffix is that it are updated at the conclusion of a
put statement. Consequently, the .cr value remains constant throughout the writing of items for the next put
statement, even if multiple items are displayed. Its behavior is similar to that of .cl.
15.14 Paging Control 139
Unlike the row and column control, the last line suffix is updated continuously. Last line suffixes are especially
useful for modifying the various writing areas of a page.
The .tlll and .hdll suffixes may not hold values applicable to the current page because when the
title or header blocks are modified, they correspond to the title or header blocks of the next page
whenever the window has been written to on the current page.
Not only can this suffix be used to determine the last line used in a writing area, but it can also be
used to delete lines within this area.
In the following example, the header section will be completely deleted by resetting the .hdll suffix to 0.
file out;
puthd out This header statement will be eliminated;
out.hdll = 0;
In this example, a header is initially written. By changing the .hdll suffix to 0, the cursor is reset to the top
of the header block. Consequently the header will not be written unless something new is added to the header
block.
putpage out This text is placed in window and the page ends;
Two additional file suffixes that can help the user in determining when to page a file are:
last page (.lp) Indicates the number of pages that are already in the document. Note that setting this to 0
does not erase the pages that have previously been written to the file.
window size (.ws) Shows the number of rows, which can be placed in the window, considering the number of
lines that are in the title and header blocks of the current page and the existing page size. The .ws file
suffix value is calculated by GAMS and is not changeable by the user. This suffix is useful for manual
pagination when used in conjunction with the .ll file suffix.
140 The Put Writing Facility
To illustrate its use, the following statement could be inserted at any point of a program to detect the number of
errors, which have occurred up to its location. The choice of output file could be the same file, a different file, or
the console as appropriate:
In this example it is assumed that the files out.put and error.put have previously been defined with a file
statement. With this statement, the number of put errors that occur in the file out.put are displayed in the file
error.put. Using putpage would allow the immediate display to the screen of a PC system at the location of
this statement if the console had been the output device.
be prepared with comma delimiters and text items in quotes. This is implemented by using .pc suffix value 5.
Delimited files are different than normal put files. All output items are written with variable field widths and
separated by delimiters. Consequently, all global and local format specifications for field widths and justification
are ignored by GAMS. Note that the number of decimals for numeric items can still be specified with the .nd
file suffix. Each item is written immediately following the previous delimiter on the same line unless the cursor is
reset.
Avoid horizontal cursor relocations in a program, which creates a delimited file. Horizontally
relocating the cursor in a delimited file is potentially damaging since a delimiter could be overwritten.
While the comma is the most common delimiting character for spreadsheets, other delimiters like blank space
and tab characters can also be used.
15.17.1 An Example
In the following example, the capacity sub-table of the [MEXSS] report program is prepared as a delimited file.
The following program segment demonstrates .pc suffix value 5. The program segment could be placed at the
end of the original [MEXSS] model:
The first line of this program segment creates the file out.put as the delimited file. Notice that in the remainder
of this program, field widths, justifications, and horizontal cursor relocations are completely avoided. All text
items are quoted. The following is the resulting output file:
"CAPACITY (tons)","AHMSA","FUNDIDORA","SICARTSA","HYLSA","HYLSAP"
"BLAST FURNACES",3.25,1.40,1.10,0.00,0.00
"oPEN HEARTH FURNACES",1.50,0.85,0.00,0.00,0.00
"BASIC OXYGEN CONVERTERS",2.07,1.50,1.30,0.00,0.00
"DIRECT REDUCTION UNITS",0.00,0.00,0.00,0.98,1.00
"ELECTRIC ARC FURNACES",0.00,0.00,0.00,1.13,0.56
Notice that each item is delimited with a comma and that textual output is quoted.
142 The Put Writing Facility
16
16.1 Introduction
The previous chapters have focused on the ability of GAMS to describe models. This chapter will describe the
various programming features available in GAMS to help the advanced user. The various programming flow
control features discussed in this chapter are
Loop Statement If-Else Statement
For Statement While Statement
loop(controlling_domain[$(condition)],
statement {; statement}
) ;
If the controlling domain consists of more than one set, then parentheses are required around it.
The loop statement causes GAMS to execute the statements within the scope of the loop for each member of
the driving set(s) in turn. The order of evaluation is the entry order of the labels. A loop is thus another, more
general, type of indexed operation. The loop set may be dollar-controlled and does not need to be static or
nested. Loops may be controlled by more than one set.
It is illegal to modify any controlling set inside the body of the loop.
144 Programming Flow Control Features
16.2.2 Examples
Consider a hypothetical case when a growth rate is empirical:
set t / 1985*1990 /
parameter pop(t) / 1985 3456 /
growth(t) / 1985 25.3, 1986 27.3, 1987 26.2
1988 27.1, 1989 26.6, 1990 26.6 /;
in an iterative rather than a parallel way. In this example there is one statement in the scope of the loop, and
one driving, or controlling, set.
A loop is often used to perform iterative calculations. Consider the following example, which uses finds square
roots by Newtons method. This example is purely for illustration - in practice, the function sqrt should be used.
Newtons method is the assertion that if x is an approximation to the square root of v, then (x + v/x)/2 is a
better one
scalars
target "number whose square root is needed" /23.456 /
sqrtval "final approximation to sqrt(target)"
curacc "accuracy of current approximation"
reltol "required relative accuracy" / 1.0e-06 / ;
if (condition,
statements;
{elseif condition, statements }
[else statements;]
);
16.3.2 Examples
Consider the following set of statements
p(i)$(f <= 0) = -1 ;
p(i)$((f > 0) and (f < 1)) = p(i)**2 ;
p(i)$(f > 1) = p(i)**3 ;
q(j)$(f <= 0) = -1 ;
q(j)$((f > 0) and (f < 1)) = q(j)**2 ;
q(j)$(f > 1) = q(j)**3 ;
if (f <= 0,
p(i) = -1 ;
q(j) = -1 ;
elseif ((f > 0) and (f < 1)),
p(i) = p(i)**2 ;
q(j) = q(j)**2 ;
else
p(i) = p(i)**3 ;
q(j) = q(j)**3 ;
) ;
The body of the if statement can contain solve statements. For instance, consider the following bit of GAMS
code:
if ((ml.modelstat eq 4),
* model ml was infeasible
* relax bounds on x and solve again
x.up(j) = 2*x.up(j) ;
solve ml using lp minimizing lp ;
else
if ((ml.modelstat ne 1),
abort "error solving model ml ;
);
);
The following GAMS code is illegal since one cannot define equations inside an if statement.
if (s gt 0,
eq.. sum(i,x(i)) =g= 2 ;
);
The following GAMS code is illegal since one cannot make declarations inside an if statement.
if (s gt 0,
scalar y ; y = 5 ;
);
146 Programming Flow Control Features
while(condition,
statements;
);
16.4.2 Examples
One can use while statements to control the solve statement. For instance, consider the following bit of GAMS
code that randomly searches for a global optimum of a non-convex model:
In this example, a non-convex model is solved from 1000 random starting points, and the global solution is tracked.
The model [PRIME] from the model library illustrates the use of the while statement through an example where
the set of prime numbers less than 200 are generated
The following GAMS code is illegal since one cannot define equations inside a while statement.
while (s gt 0,
eq.. sum(i,x(i)) =g= 2 ;
);
The following GAMS code is illegal since one cannot make declarations inside a while statement.
while(s gt 0,
scalar y ; y = 5 ;
);
Note that i is not a set but a parameter. Start and end are the start and end, and incr is the increment by
which i is changed after every pass of the loop.
The values of start, end and incr need not be integer. The start and end values can be positive or
negative real numbers. The value of incr has to be a positive real number.
16.5.2 Examples
One can use for statements to control the solve statement. For instance, consider the following bit of GAMS
code that randomly searches for a global optimum of a non-convex model:
scalar i ;
scalar globmin ; globmin = inf ;
option bratio = 1 ;
for (i = 1 to 1000,
x.l(j) = uniform(0,1) ;
solve ml using nlp minimizing obj ;
if (obj.l le globmin,
globmin = obj.l ;
globinit(j) = x.l(j) ;
);) ;
In this example, a non-convex model is solved from 1000 random starting points, and the global solution is tracked.
The use of real numbers as start, end and incr can be understood from the following example,
Notice that the value of s was incremented by 1.4 with each pass of the loop as long as it did not exceed 0.3.
The following GAMS code is illegal since one cannot define equations inside a for statement.
for (s = 1 to 5 by 1,
eq.. sum(i,x(i)) =g= 2 ;
);
The following GAMS code is illegal since one cannot make declarations inside a for statement.
for (s=1 to 5 by 1,
scalar y ; y = 5 ;
);
148 Programming Flow Control Features
17
17.1 Introduction
This chapter introduces special features in GAMS that do not translate across solvers, or are specific to certain
model types. These features can be extremely useful for relevant models, and are among the most widely used.
integer variables These can take on integer values between. The default lower and upper bounds are
0 and 100 respectively.
In addition to these two, two new types of discrete variables that are introduced in this section. Both these
variables exploit special structures in MIP models during the solution phase. These are the following
Special Ordered Sets (SOS) The precise definition of special ordered sets differ from one solver to another
and the development of these features has been driven more by internal algorithmic consideration than by
broader modeling concepts. GAMS offers sos1 and sos2 variables as two types of compromise features
that model special ordered sets. Sections 17.2.2 and 17.2.3 discuss these two types of variables in greater
detail.
Semi-continuous variables GAMS offers semicont and semiint variables to model this class of variables.
These are explained in Sections 17.2.3 and 17.2.4.
The presence of any of the above types of discrete variables requires a mixed integer model and all the discreteness
is handled by the branch and bound algorithm in the same way as binary and general integer variables are handled.
The members of the innermost (the right-most) index belongs to the same set. For example, in the sets defined
above, s1 represents one special ordered set of type 1 with i elements, t1 defines k sets of j elements each, and
w1 defines (i,j) sets with k elements each.
The default bounds for SOS1 variables are 0 to +. As with any other variable, the user may set
these bounds to whatever is required.
The user can, in addition, explicitly provide whatever convexity row that the problem may need
through an equation that requires the members of the SOS set to be less than a certain value. Any
such convexity row would implicitly define bounds on each of the variables.
The equation defsoss1 implicitly defines the non-zero value that one of the elements of the SOS1 variable s1 can
take.
A special case of SOS1 variables is when exactly one of the elements of the set have to be non-zero. In this case,
the defsoss1 equation will be
A common use of the use of this set is for the case where the non-zero value is 1. In such cases, the SOS1 variable
behaves like a binary variable. It is only treated differently by the solver at the level of the branch and bound
algorithm. For example, consider the following example to model the case where at most one out of n options
can be selected. This is expressed as
The variable x can be made binary without any change in meaning and the solution provided by the solver will
be indistinguishable from the SOS1 case.
The use of special ordered sets may not always improve the performance of the branch and bound algorithm.
If there is no natural order the use of binary variables may be a better choice. A good example of this is the
assignment problem.
Not all MIP solvers allow SOS1 variables. Furthermore, among the solvers that allow their use, the
precise definition can vary from solver to solver. Any model can contains these variables may not be
transferable among solvers. Please verify how the solver you are interested in handles SOS1 variables
by checking the relevant section of the Solver Manual.
The default bounds for SOS2 variables are 0 to +. As with any other variable, the user may set
these bounds to whatever is required.
The members of the innermost (the right-most) index belongs to the same set. For example, in the sets defined
above, s2 represents one special ordered set of type 2 with i elements, t2 defines k sets of j elements each, and
w2 defines (i,j) sets with k elements each.
[PRODSCHX] shows SOS type formulations with binary, SOS1 and SOS2 sets. The default bounds for SOS
variables are 0 to +. As with any other variable, the user may set these bounds to whatever is required.
Not all MIP solvers allow SOS2 variables. Furthermore, among the solvers that allow their use, the
precise definition can vary from solver to solver. Any model can contains these variables may not be
transferable among solvers. Please verify how the solver you are interested in handles SOS2 variables
by checking the relevant section of the Solver Manual.
semicont variable x ;
x.lo = 1.5 ; x.up = 23.1 ;
The above slice of code declares the variable x to be semi-continuous variable that can either be 0, or can behave
as a continuous variable between 1.5 and 23.1.
Not all MIP solvers allow semi-continuous variables. Please verify that the solver you are interested
in can handle semi-continuous variables by checking the relevant section of the Solver Manual.
The lower bound has to be less than the upper bound, and both bounds have to be greater than 0.
GAMS will flag an error if it finds that this is not the case.
semiint variable x ;
x.lo = 2 ; x.up = 25 ;
The above slice of code declares the variable x to be semi-continuous variable that can either be 0, or can take
any integer value between 2 and 25.
Not all MIP solvers allow semi-integer variables. Please verify that the solver you are interested in
can handle semi-integer variables by checking the relevant section of the Solver Manual.
The lower bound (L) has to be less than the upper bound (U ), and both bounds have to be greater
than 0. GAMS will flag an error during model generation if it finds that this is not the case.
The bounds for semiint variables have to take integer values. GAMS will flag an error during model
generation if it finds that this is not the case.
152 Special Language Features
mymodel.prioropt = 1 ;
where mymodel is the name of the model specified in the model statement. The default value is 0 in which case
priorities will not be used.
Using the .prior suffix sets the priorities of the individual variables. Note that there is one prior value for each
individual component of a multidimensional variable. Priorities can be set to any real value. The default value
is 1. As a general rule of thumb, the most important variables should be given the highest priority.
The following example illustrates its use,
z.prior(i,small) = 3 ;
z.prior(i,medium) = 2 ;
z.prior(i,large) = 1 ;
In the above example, z(i,large) variables are branched on before z(i, small) variables.
The lower the value given to the .prior suffix, the higher the priority for branching.
All members of any SOS1 or SOS2 set should be given the same priority value since it is the set itself
which is branched upon rather than the individual members of the set.
The statement should be inserted somewhere after the model statement and before the solve statement. In order
to turn scaling off again, set the model.scaleopt parameter to 0 before the next solve.
The scale factor of a variable or an equation is referenced with the suffix .scale, i.e. the scale factor of variable
x(i) is referenced as x.scale(i). Note that there is one scale value for each individual component of a multi-
dimensional variable or equation. Scale factors can be defined using assignment statements. The default scale
factor is always 1.
GAMS scaling is in most respects hidden from the user. The solution values reported back from a solution
algorithm are always reported in the users notation. The algorithms versions of the equations and variables are
17.3 Model Scaling - The Scale Option 153
only reflected in the derivatives in the equation and column listings in the GAMS output if the options limrow
and limcol are positive, and the debugging output from the solution algorithm generated with sysout option
set to on.
By setting x1.scale to 0.01 and x2.scale to 10, the model seen by the solver is,
Note that the solver does not see the variables x1 or x2, but rather the scaled (and better-behaved) variables
xprime1 and xprime2.
Upper and lower bounds on variables are automatically scaled in the same way as the variable itself.
The user may have to perform a combination of equation and variable scaling until a well-scaled
model is obtained.
Consider the following example,
x1.scale = 0.1 ;
eq1.scale = 5 ;
eq2.scale = 5 ;
will result in the solver seeing the following well scaled model,
Glossary
acronym A GAMS data type used to give logical classifications to data points.
algorithm This term may be used in two ways. It is either a prescription for how to solve a problem, or a
particular solver system.
basic A classification of a row or column that is in the basis maintained by solution methods that use linear
programming iterations.
binding An inequality constraint is binding when the value of the associated slack is zero.
bounds Upper and lower limits on the possible values that a column may assume in a feasible solution. May be
infinite, meaning that no limit is imposed.
column An individual decision variable in the model seen by a solver program. Many may be associated with
one GAMS variable
compilation The initial phase of GAMS processing, when the program is being checked for syntax and
consistency.
constant set A set is constant if it remains unchanged. It has to be initialized with a set definition statement
and cannot be changed using assignment statement. Sets used in domain definitions must be constant.
Sets used in lag operations must be ordered as well. Sometimes the word static is used instead of constant.
constraint A relationship between columns that must hold in a feasible solution. There may be many
constraints associated with one GAMS equation.
continuous There are two contexts. First a classification of a function. A plot of the function values will be a
line without breaks in it. Second, a classification of variables. A continuous variable may assume any value
within its bounds.
data types Each symbol or identifier has to be declared to be one of the seven data types, which are set,
parameter, variable, equation, model, file and acronym. The keywords scalar and table do not
introduce separate data types but rather comprise a shorthand way to declare a symbol to be a parameter
that will use a particular format for specifying initial values.
declaration The entry of a symbol and the specification of its data type. A declaration may include the
specification of initial values, and then it is more properly called a definition.
default The value used, or the action taken, if the user provides no information.
158 Glossary
direction Either maximization or minimiz depeding on whether the use is interested in the llatargest or
the smalle possible vue f the obe function.
discontiuouio]TJET1001116.919599.726cm0g0G1001-116.919-599.726cmBT/F89.963Tf121.9599.726Td[(A)-333(c)
dollar operator An ope used for excep handling in aison1(ign)atmet statements and in equlatation
definiti
domai checking The check that ensu that only legal label cmbinati are used on every assignment to,
or reference of an identi
domai defin The label cmbinationion334(who)1(s)-1(e)-334(d)1(ata)-333(w)-1(i)1(ll)-333(b)-28(e)-333(up)-27(dated)-333
domai restrict condit The alte to the domain of defin caused when a dollar operator is
used on the left (of the = in an assignment or of the .. in an equlatation deniti
dr set The seton334(t)1(hat)-333(determine)-333(the)-334(d)1(omain)-333(of)-333(d)1(e)-atfiniti or that control and index
dynamic set A seton333(is)-334(d)1(ynamic)-334(i)1(f)-333(it)-333(has)-333(b)-28(ee)-atn change with an aison1(ign)atmet sta
used with lag ope or in domain defin.
endogenouio]TJET1001106.99339.622cm0g0G1001-106.99-339.622cmBT/F89.963Tf111.971339.622Td[(Data)-333(
equation The GAMS data type use tospe re reonshiplation334(b)-27(e)-attween activity leels of varlatiables.
execution The secon phlatase of GAMS processing, whe GAMS ilation334(acon1(t)1(ually)-333(carr)1(yin)1(g)-334(ou)1(t)-333
execution sttio]TJET1001154.95256.906cm0g0G1001-154.95-256.906cmBT/F89.963Tf159.931256.906Td[(Instru)1
gene. Someon334(e)-1(x)atampare theon334(aison1(i)atgnme)-atnt andl-333(th)1(e)]TJ/F229.963Tf273.9730Td[(option)]TJ
exogenouio]TJET1001100.307213.394cm0g0G1001-100.307-213.394cmBT/F89.963Tf105.288213.394Td[(Data)-328
feasible solution A solution to a model in which all column activity levels are within the bounds and all the
constraints are satisfied.
GAMS coordinator The person who looks after the administration of a GAMS system, and who will know
what solvers are available and can tell you who to approach for help with GAMS problems. Unlikely to
apply to personal computer versions.
index position(s) Another way of describing the set(s) that must be used when referencing a symbol of
dimensionality one or more (i.e., a vector or a matrix).
inequality constraint A constraint in which the imposed relationship between the columns is not fixed, but
must be either greater than or equal to, or less than or equal to, a constant. The GAMS symbols =g= and
=l= are used in equation definitions to specify these relationships.
infeasible Used to describe either a model that has no feasible solution, or an intermediate solution that is not
feasible (although feasible solutions may exist). See feasible, above.
initialization Associating initial values with sets or parameters using lists as part of the declaration or
definition, or (for parameters only) using table statements.
list One of the ways of specifying initial values. Used with sets or parameters, most often for one-dimensional
but also for two and higher dimensional data structures.
list format One of the ways in which sets and parameters, can be initialized and all symbol classes having data
can be displayed. Each unique label combination is specified in full, with the associated non-default value
alongside.
marginal Often called reduced costs or dual values. The values, which are meaningful only for non-basic rows
or columns in optimal solutions, contain information about the rate at which the objective value will
change of if the associated bound or right hand side is changed.
model generation The initial phase of processing a solve statement: preparing a problem description for the
solver.
nonbasic A column that is not basic and (in nonlinear problems) not superbasic. Its value will be the same as
the one of the finite bounds (or zero if there are no finite bounds) if the solution is feasible.
nonlinear nonzero In a linear programming problem, the nonzero elements are constant. In a nonlinear
problem, one or more of them vary because their values depend on that of one or more columns. The ratio
of nonlinear (varying) to linear (constant) non linear zero elements is a good indicator of the pervasiveness
of non-linearities in the problem.
nonoptimal There are two contexts. First, describing a variable: a non-basic variable that would improve the
objective value if made basic. The sign of the marginal value is normally used to test for non-optimality.
Second, for a solution: other solutions exists with better objective values.
nonsmooth A classification of function that does not have continuous first derivatives, but has continuous
function values. A plot of the function values will be a line with kinks in it.
nonzero element The coefficient of a particular column in a particular row if it is not zero. Most
mathematical programming problems are sparse meaning that only a small proportion of the entries in the
full tableau of dimensions number of rows by number of columns is different from zero.
objective row (or function) Solver system require the specification of a row on (for nonlinear systems) a
function whose value will be maximized or minimized. GAMS users, in contrast, must specify a scalar
variable.
160 Glossary
objective value The current value of the objective row or of the objective variable.
optimal A feasible solution in which the objective value is the best possible.
option The statement that allows users to change the default actions or values in many different parts of the
system.
ordered set A set is ordered if its content has been initialized with a set definition statement and the entry
order of the individual elements of the set has the same partial order as the chronological order of the
labels. A set name alone on the left-hand side of an assignment statement destroys the ordered property.
Lag and Ord operations rely on the relative position of the individual elements and therefore require
ordered sets. Ordered sets are by definition constant.
output file A disk file containing output. A GAMS task produces one such file that can be inspected.
parameter A constant or group of constants that may be a scalar, a vector, or a matrix of two or more
dimensions. Of the six data types in GAMS.
problem type A model class that is dependent on functional form and specification. Examples are linear,
nonlinear, and mixed integer programs.
program A GAMS input file that provides a representation of the model or models.
relational operator This term may be used in two ways. First, in an equation definition it describes the type
of relationships the equation specifies, for example equality, as specified with the =e= symbol. Second, in
a logical expression, the symbols eq, ne, lt and so on are also called relational operators, and are used to
specify a required relationship between two values.
scalar One of the forms of parameter inputs. Used for single elements.
set A collection of elements (labels). The set statement is used to declare and define a set.
simplex method The standard algorithm used to solve linear programming problems.
slack variable An artificial column introduced by a solver into a linear programming problem. Makes the
implementation of simplex method much easier.
solver A computer code used to solve a given problem type. An example is GAMS/MINOS, which is used to
solve either linear or nonlinear programming problems.
statements Sometimes called units. The fundamental building block of GAMS programs. Statements or
sentences that define data structures, initial values, data modifications, and symbolic relationships.
Examples are table, parameter, variable, model, assignment and display statements.
superbasic In nonlinear programming, a variable that it is not in the basis but whose value is between the
bounds. Nonlinear algorithms often search in the space defined by the superbasic variables.
symbol An identifier.
table One of the ways of initializing parameters. Used for two and higher dimensional data structures.
Professor Paul Samuelson is fond of saying that he hopes each generation economists will be able to stand on
the shoulders of the previous generation. The library of models included with the GAMS system is a reflection
of this desire. We believe that the quality of modeling will be greatly improved and the productivity of modelers
enhanced if each generation can stand on the shoulders of the previous generation by beginning with the previous
models and enhancing and improving them. Thus the GAMS systems includes a large library, collectively called
GAMSLIB .
The models included have been selected not only because they collectively provide strong shoulders for new users
to stand on, but also because they represent interesting and sometimes classic problems. For example the trade-
off between consumption and investment is richly illustrated in the Ramsey problem, which can be solved using
nonlinear programming methods. Examples of other problems included in the library are production and shipment
by firms, investment planning in time and space, cropping patterns in agriculture, operation of oil refineries and
petrochemical plants, macroeconomics stabilization, applied general equilibrium, international trade in aluminum
and in copper, water distribution networks, and relational databases.
Another criterion for including models in the library is that they illustrate the modeling capabilities GAMS offers.
For example, the mathematical specification of cropping patterns can be represented handily in GAMS. Another
example of the systems capability is the style for specifying initial solutions as staring points in the search for
the optimal solution of dynamic nonlinear optimization problems.
Finally, some models have been selected for inclusion because they have been used in other modeling systems.
Examples are network problems and production planning models. These models permit the user to compare how
problems are set up and solved in different modeling systems.
Most of the models have been contributed by GAMS users. The submission of new models is encouraged. If you
would like to see your model in a future release of the library, please send the model and associated documents
and reports to GAMS Development Corporation.
The most convenient way (Windows only) to access the model library is from within the GAMS IDE by going
through: File Model Library Open GAMS Model Library. A window will pop up and give you access to all
models.
Another way to access the library is through the gamslib command. This command copies a model from the
library directory into the current directory. If you enter gamslib without any parameters, the command syntax
will be displayed as shown below:
>gamslib modelname [target]
or
>gamslib modelnum [target]
where modelname is the modelname, modelnum is the model sequence number, and target is the target file name.
If the target file name is not provided, the default is modelname.gms. For example, the [TRNSPORT] model
could be copied in any of the following ways
164 The GAMS Model Library
The full and annotated list of the models of the GAMS Model Library is available at:
http://www.gams.com/modlib/modlib.htm
Chemical Engineering
Name Description Type No. Author
alkyl Simplified Alkylation Process NLP 165 Berna, T
batchdes Optimal Design for Chemical Batch Processing MINLP 119 Kocis, G R
catmix Catalyst Mixing COPS 2.0 #14 NLP 242 Dolan, E D
chem Chemical Equilibrium Problem NLP 21 Bracken, J
csched Cyclic Scheduling of Continuous Parallel Units MINLP 222 Jain, V
feedtray Optimum Feed Plate Location MINLP 122 Viswanathan, J
gasoil Catalytic cracking of gas oil COPS 2.0 #12 NLP 240 Ascher, U M
haverly Haverlys pooling problem example NLP 214 Adhya, N
hda Synthesis: Hydrodealkylation of Toluene MINLP 123 Kocis, G R
kport Product Portfolio Optimization MINLP 270 Kallrath, J
methanol Methanol to hydrocarbons COPS 2.0 #13 NLP 241 Ascher, U M
minlphi Heat Integrated Distillation Sequences NLP 118 Floudas, C A
minlphix Heat Integrated Distillation Sequences MINLP 227 Floudas, C A
nonsharp Synthesis of General Distillation Sequences MIP 120 Aggarwal, A
nsharpx Synthesis of General Distillation Sequences MINLP 226 Aggarwal, A
pinene Isometrization of alpha-pinene COPS 2.0 #5 NLP 233 Ascher, U M
pool Pooling problem NLP 254 Adhya, N
process Alkylation Process Optimization NLP 20 Bracken, J
procsel Structural Optimization of Process Flowsheets MINLP 116 Kocis, G R
reaction Logical Inference for Reaction path synthesis MIP 121 Raman, R
synheat Simultaneous Optimization for Hen Synthesis MINLP 117 Yee, T F
wall Chemical Equilibrium Problem NLP 76 Wall, T W
wallmcp Chemical Equilibrium Problem as MCP MCP 127 Wall, T W
Econometrics
Name Description Type No. Author
least Nonlinear Regression Problem NLP 24 Bracken, J
like Maximum Likelihood Estimation NLP 25 Bracken, J
linear Linear Regression with Various Criteria DNLP 23 Bracken, J
Economic Development
Name Description Type No. Author
chenery Substitution and Structural Change NLP 33 Chenery, H B
dinam DINAMICO A Dynamic Multi-Sectoral Multi-Skill Model LP 35 Manne, A S
dmcmge Accounting for economic growth with new inputs MPSGE 143 Markusen, J R
pak Optimal Patterns of Growth and Aid LP 34 Chenery, H B
prolog Market Equilibrium and Activity Analysis NLP 41 Norton, R D
ramsey Savings Model by Ramsey NLP 63 Murtagh, B A
Energy Economics
Name Description Type No. Author
etamac ETA-MACRO Energy Model for the USA NLP 80 Manne, A S
etamge ETA-MACRO Energy Model for the USA - MPSGE Format MPSGE 144 Manne, A S
gtm International Gas Trade Model NLP 53 Manne, A S
korpet Investment Planning in the Korean Oil-Petro Industry MIP 48 Suh, J S
launch Launch Vehicle Design and Costing NLP 161 Bracken, J
markov Strategic Petroleum Reserve LP 82 Teisberg, T J
otpop OPEC Trade and Production NLP 47 Blitzer, C
pindyck Optimal Pricing and Extraction for OPEC NLP 28 Pindyck, R S
pollut Industrial Pollution Control NLP 96 Mangasarian, O L
shale Investment Planning in the Oil Shale Industry LP 46 Melton, J W
166 The GAMS Model Library
Engineering
Name Description Type No. Author
bearing Hydrostatic Thrust Bearing Design for a Turbogenerator NLP 202 Coello Coello, C A
camshape Shape optimization of a cam COPS 2.0 #4 NLP 232 Anitescu, M
chain Hanging Chain COPS 2.0 #3 NLP 231 Cesari, L
dispatch Economic Load Dispatch Including Transmission Losses NLP 166 Wood, A J
elec Distribution of electrons on a sphere COPS 2.0 #2 NLP 230 Dolan, E D
flowchan Flow in a channel COPS 2.0 #7 NLP 235 Ascher, U M
fuel Fuel Scheduling and Unit Commitment Problem MINLP 168 Wood, A J
gasnet Optimal Design of a Gas Transmission Network MINLP 223 Edgar, T F
gastrans Gas Transmission Problem - Belgium MINLP 217 de Wolf, D
gear Gear Train Design MINLP 203 Deb, K
glider Hang glider COPS 2.0 #11 NLP 239 Bulirsch, R
herves Herves (Transposable Element) Activity Calculations DNLP 305 Meeraus, W H
hydro Hydrothermal Scheduling Problem NLP 167 Wood, A J
jbearing Journal bearing COPS 2.0 #16 NLP 244 Averick, B M
lnts Particle steering COPS 2.0 #9 NLP 237 Betts, J
minsurf Minimal surface with obstacle COPS 2.0 #17 NLP 245 Dolan, E D
pump Pump Network Synthesis MINLP 205 Floudas, C A
robot Robot arm COPS 2.0 #8 NLP 236 Dolan, E D
rocket Goddard rocket COPS 2.0 #10 NLP 238 Bryson, A E
ship Structural Optimization NLP 22 Bracken, J
spring Coil Compression String Design MINLP 206 Sandgren, E
torsion Elastic-plastic torsion COPS 2.0 #15 NLP 243 Averick, B M
trimloss Trim Loss Minimization MINLP 204 Floudas, C A
water Design of a Water Distribution Network DNLP 68 Brooke, A
waterx Design of a Water Distribution Network (MINLP) MINLP 125 Brooke, A
windfac Winding Factor of Electrical Machines MINLP 224 Michna, M
Finance
Name Description Type No. Author
alan A Quadratic Programming Model for Portfolio Analysis MINLP 124 Manne, A S
cmo Financial Optimization: Financial Engineering MIP 114 Dahl, H
immun Financial Optimization: Risk Management NLP 110 Dahl, H
meanvar Financial Optimization: Risk Management NLP 112 Dahl, H
meanvarx Financial Optimization: Risk Management MINLP 113 Dahl, H
port Simple Portfolio Model LP 50 CDC
qalan A Quadratic Programming Model for Portfolio Analysis QCP MIQCP 282 Manne, A S
qcp1 Standard QP Model QCP QCP 283 Kalvelagen, E
qmeanvar Financial Optimization: Risk Management using MIQCP MIQCP 291 Dahl, H
qp1 Standard QP Model NLP 171 Kalvelagen, E
qp1x Standard QP Model with GDX data input NLP 246 Kalvelagen, E
qp2 Standard QP Model - symmetry exploitations NLP 172 Kalvelagen, E
qp3 Standard QP Model - intermediate variables NLP 173 Kalvelagen, E
qp4 Standard QP Model - no covariance matrix NLP 174 Kalvelagen, E
qp5 Standard QP Model - linear approximation LP 175 Kalvelagen, E
qp6 Standard QP Model - LCP formulation of QP4 MCP 184 Kalvelagen, E
qp7 Standard QP Model - conic programming formulation LP 271 Andersen, E
repay Repayment Factors for Loans GAMS 72 GAMS Development
Corporation
tba Financial Optimization: Financial Engineering MIP 115 Dahl, H
worst Financial Optimization: Risk Management NLP 111 Dahl, H
Forestry
Name Description Type No. Author
tfordy Antalya Forestry Model - Dynamic LP 62 Bergendorff, H
tforss Antalya Forestry Model - Steady State LP 61 Bergendorff, H
B.1 The Models 167
International Trade
Name Description Type No. Author
alum World Aluminum Model MIP 31 Brown, M
copper Modeling Investment in the World Copper Industry MIP 45 Dammert, A
co2mge Carbon-related trade model (static) MPSGE 142 Perroni, C
nsmge North-South trade and capital flows MPSGE 150 Manne, A S
Macro Economics
Name Description Type No. Author
abel Linear Quadratic Control Problem NLP 64 Kendrick, D
chakra Optimal Growth Model NLP 43 Chakravarty, S
hhfair Household Optimization Problem by Fair NLP 69 Fair, R C
orani A Miniature Version of Orani 78 LP 40 Kendrick, D
qabel Linear Quadratic Control Problem QCP 293 Kendrick, D
Mathematics
Name Description Type No. Author
absmip Discontinous functions abs() min() max() sign() as MIPs MIP 208 GAMS Dev. Corp.
circle Circle Enclosing Points - SNOPT Example NLP 201 Gill, P E
crazy Examples of Extended Arithmetic GAMS 67 GAMS Dev. Corp.
cube Three-Dimensional Noughts and Crosses MIP 42 Williams, H P
dice Non-transitive Dice Design MIP 176 Gardner, M
dicex Non-transitive Dice Design - Enhanced MIP 272 Bosh, R A
fct LGO Interface Example DNLP 265 Pinter, J D
gauss Matrix Inversion with Full Pivoting GAMS 71 GAMS Dev. Corp.
great Great Circle Distances GAMS 73 Brooke, A
himmel11 Himmelblau Test Problem Number 11 NLP 95 Himmelblau, D M
himmel16 Area of Hexagon Test Problem NLP 36 Himmelblau, D M
house House Plan Design NLP 99 Borland
hs62 Hock - Schittkowski Problem 62 NLP 264 Hock, W
icut Integer Cut Example MIP 160 GAMS Dev. Corp.
imsl Piecewise Linear Approximation LP 59 IMSL Inc
knights Maximum Knights Problem MIP 158 Dudeney, H E
latin The Orthogonal Latin-Square Problem MIP 159 Dantzig, G B
lmp1 Linear Multiplicative Model - Type 1 NLP 251 Konno, H
lmp2 Linear Multiplicative Model - Type 2 NLP 252 Tawarmalani, M
lmp3 Linear Multiplicative Model - Type 3 NLP 253 Benson, H P
mathopt1 MathOptimizer Example 1 NLP 255 Mathematica
mathopt2 MathOptimizer Example 2 NLP 256 Mathematica
mathopt3 MathOptimizer Example 3 NLP 257 Mathematica
mathopt4 MathOptimizer Example 4 NLP 258 Mathematica
mathopt5 MathOptimizer Example 5 DNLP 259 Mathematica
mathopt6 MathOptimizer Example 6 DNLP 260 Mathematica
maxmin Max Min Location of Points in Unit Square DNLP 263 Pinter, J D
mhw4d Nonlinear Test Problem NLP 84 Wright, M H
mhw4dx MHW4D with additional Tests NLP 267 Wright, M H
mhw4dxx MHW4DX with multiple solutions NLP 295 Wright, M H
pentium Pentium Error Test GAMS 162 House, B
pivot Simple Gaussian Elimination GAMS 70 GAMS Dev. Corp.
polygon Largest small polygon COPS 2.0 #1 NLP 229 Dolan, E D
popdynm Marine population dynamics COPS 2.0 #6 NLP 234 Dolan, E D
prime Prime Number Generation GAMS 157 GAMS Dev. Corp.
qp6 Standard QP Model - LCP formulation of QP4 MCP 184 Kalvelagen, E
rbrock Rosenbrock Test Function NLP 83 Rosenbrock, H H
t1000 LGO Test t1000 10 NLP 266 Pinter, J D
trig Simple Trigonometric Example NLP 261 Pinter, J D
Micro Economics
Name Description Type No. Author
andean Andean Fertilizer Model MIP 44 Mennes, L B M
bid Bid Evaluation MIP 19 Bracken, J
bidsos Bid Evaluation with SOS2 Sets MIP 163 Bracken, J
cesam Cross Entropy SAM Estimation MCP 228 Golan, A
diet Stiglers Nutrition Model LP 7 Dantzig, G B
fertd Egypt - Dynamic Fertilizer Model MIP 14 Choksi, A M
ferts Egypt - Static Fertilizer Model LP 13 Choksi, A M
harker Models of Spatial Competition NLP 85 Harker, P T
marco Mini Oil Refining Model LP 32 Aronofsky, J S
mexls Mexico Steel - Large Static LP 17 Kendrick, D
mexsd Mexico Steel - Small Dynamic MIP 16 Kendrick, D
mexss Mexico Steel - Small Static LP 15 Kendrick, D
msm Morocco Fertilizer Distribution - Mode Selection LP 66 Meeraus, A
phosdis Sea Distances for World Phosphate Model LP 78 Meeraus, A
B.1 The Models 169
Stochastic Programming
Name Description Type No. Author
airsp Aircraft Allocation OSLSE 189 Dantzig, G B
airsp2 Aircraft Allocation - stochastic optimization with DECIS DECIS 196 Dantzig, G B
apl1p Stochastic Programming Example for DECIS DECIS 197 Infanger, G
apl1pca Stochastic Programming Example for DECIS DECIS 198 Infanger, G
clearlak Scenario Reduction: ClearLake exercise LP 249 Birge, J R
farm The Farmers Problem formulated for DECIS DECIS 199 Birge, J R
kand Stochastic Programming OSLSE 187 Kall, P
lands Optimal Investment OSLSE 188 Louveaux, F V
prodsp Stochastic Programming Example OSLSE 186 King, A J
prodsp2 Stochastic Programming Example - reformulated for DECIS DECIS 200 King, A J
srkandw Stochastic Programming Scenario Reduction LP 248 Kall, P
Statistics
Name Description Type No. Author
dea Data Envelopment Analysis - DEA LP 192 Dyson
mingamma Minimal y of GAMMA(x) DNLP 299 Sloane, N J A
mlbeta Fitting of beta distribution through maximum likelihood NLP 302 Kotz, J
mlgamma Max. Likelihood est. of parameters of the gamma distribution NLP 303 Cox, D R
procmean Optimal Process Mean NLP 301 Chen, C H
quantum An application from quantum mechanics DNLP 300 Ogura, A
qsambal Social Accounting Matrix Balancing Problem QCP QCP 285 Zenios, S A
sambal Social Accounting Matrix Balancing Problem NLP 77 Zenios, S A
sample Stratified Sample Design NLP 27 Bracken, J
Recreational Models
Name Description Type No. Author
alphamet Alphametics - a Mathematical Puzzle MIP 170 de Wetering, A V
coex Peacefully Coexisting Armies of Queens MIP 219 Bosch, R
coexx Peacefully Coexisting Armies of Queens - tight MIP 218 Bosch, R
cross Alcuins River Crossing MIP 191 Borndoerfer, R
cube Three-Dimensional Noughts and Crosses MIP 42 Williams, H P
dice Non-transitive Dice Design MIP 176 Gardner, M
dicex Non-transitive Dice Design - Enhanced MIP 272 Bosh, R A
knights Maximum Knights Problem MIP 158 Dudeney, H E
latin The Orthogonal Latin-Square Problem MIP 159 Dantzig, G B
marilyn Numerical Puzzle MIP 193 Savant, M
queens Maximum Queens Chess Problem MIP 103 Beauvais, J
tsp1 Traveling Salesman Problem - One MIP 177 Kalvelagen, E
tsp2 Traveling Salesman Problem - Two MIP 178 Kalvelagen, E
tsp3 Traveling Salesman Problem - Three MIP 179 Kalvelagen, E
tsp4 Traveling Salesman Problem - Four MIP 180 Kalvelagen, E
tsp42 TSP solution with subtour elimination MIP 213 Dantzig, G B
GAMS Tools
Name Description Type No. Author
awkqap Input file generation with AWK for the QAP MIQCP 296 Burkard, R E
awktsp Traveling Salesman Problem Instance prepared with AWK MIP 298 Juenger, M
seders ERS Data Manipulations with SED GAMS 297 Somwaru, A
C
The entire GAMS system appears to the user as a single call that reads the input file and produces an output
file. Several options are available at this level to define the overall layout of the output page, and when to save
and restore the entire environment. Although details will vary with the type of computer and operating system
used, the general operating principles are the same on all machines.
>gams myfile
from the system prompt and GAMS will compile and execute the GAMS statements in the file myfile. If a file
with this name cannot be found, GAMS will look for a file with the extended name myfile.gms. The output will
be written by default on the file myfile.lst on PC and Unix systems, and myfile.lis on OpenVMS systems.
For example, the following statement compiles and executes the example problem [TRNSPORT] from the GAMS
model library,
>gams trnsport
where key1 is the name of the option that is being set on the command line, and value1 is the value to which the
option is set. Depending on the option, value1 could be a character string, or an integer number. For example,
consider the following commands to run [TRNSPORT] from the GAMS model library,
All the four commands above are equivalent, and each directs the output listing to the file myrun.lst. o is the
name of the option, and it is set to myfile.lst. In addition, in each case, the log of the run is redirected to the
file myrun.log.
172 The GAMS Call
Table C.1 briefly describes the various options in each of the categories. Section C.3 contains a reference list of
all options available through the command line with detailed description for each.
Model Generation During this stage, the variables and equations involved in the model being solved are
generated.
The a=e setting can only be used during restart on files that have previously been compiled, since models first
need to be compiled before they can be executed.
To illustrate the use of the dumpopt option, [TRNSPORT] has been split into two files. The first file (say
trans1.gms) contains most of the original file except for the solve statement, and looks as follows,
sets
i canning plants / seattle, san-diego /
j markets / new-york, chicago, topeka / ;
parameters
a(i) capacity of plant i in cases
/ seattle 350
san-diego 600 /
variables
x(i,j) shipment quantities in cases
z total transportation costs in thousands of dollars ;
positive variable x ;
equations
cost define objective function
supply(i) observe supply limit at plant i
demand(j) satisfy demand at market j ;
All comments have been removed from [TRNSPORT] for brevity. Running this model and saving the work files
through the save parameter leads to the generation of eight work files. The second file (say trans2.gms) generated
from [TRNSPORT] looks as follows,
One can then run trans2.gms restarting from the saved work files generated from running trans1.gms. The
result obtained is equivalent to running [TRNSPORT].
In order to use the dumpopt parameter effectively, it is required that the first line in the restart file
be the solve statement.
To illustrate the use of the dumpopt option, run the second model using the following command
where trans is the name of the saved files generated through the save parameter from trans1.gms. A new file
trans2.dmp is created as a result of this call, and looks as follows,
model transport;
Note that all the data entering the model in the solve statement has been regenerated. The parameter d has
not been regenerated since it does not appear in the model, but the parameter c is. Changing the value of the
parameter dumpopt will result in alternate names being used for the identifiers in the regenerated file.
Note that with dp=2, all the file operations are listed including the full path of each file on which any operation
is performed.
The listing file that results from running this model contains the following section,
Note that numbers ($170 and $120) flags the two errors as they occur, but the errors are explained only at the
end of the source listing. However, if the code is run using the option errmsg=1, the resulting listing file contains
the following,
Note that the explanation for each error is provided immediately following the error marker.
parameter a ; a = 0 ;
$include file2.inc
$include file2.inc
a = a+1 ;
display a ;
Running the model with the command line flag expand myfile.fil results in the creation of the file myfile.fil.
The content of this file is provided below,
1 INPUT 0 0 0 1 7 E:\TEMP\FILE1.GMS
2 INCLUDE 1 1 2 2 4 E:\TEMP\FILE2.INC
3 INCLUDE 1 1 3 5 7 E:\TEMP\FILE2.INC
The first column gives the sequence number of the input files encountered. The first row always refers the parent
file called by the GAMS call. The second column refers to the type of file being referenced. The various types of
files are
178 The GAMS Call
0 INPUT 1 INCLUDE
2 BATINCLUDE 3 LIBINCLUDE
4 SYSINCLUDE
The third column provides the sequence number of the parent file for the file being referenced. The fifth column
gives the local line number in the parent file where the $include appeared. The sixth column gives the global
(expanded) line number which contained the $include statement. The seventh column provides the total number
of lines in the file after it is processed. The eighth and last column provides the name of the file In the example
listed above, the include files file1.inc and file2.inc were included on lines 1 and 4 of the parent file test1.gms.
0 0 0 0 D:\GAMS\NEW.LST
1 1 170 31 D:\GAMS\NEW.GMS
2 2 120 14 D:\GAMS\NEW.GMS
The first column refers to the global row number of the error in the listing file. The second column refers to
the row number of the error in the individual file where the problem occurs. This will be different from the first
column only if the error occurs in an include file. In this case, the second column will contain the line number in
the include file where the error occurs, while the first number will contain the global line number (as reported in
the listing file) where the error occurs. The number in the third column refers to the error number of the error.
The fourth number refers to the column number of the error in the source file. The fifth column contains the
individual file in which the error occurred.
The option value 1 is mainly used by solvers that can be interrupted from the terminal.
that have since become keywords will cause errors when run with the latest version of GAMS. This option will
allow one to run such models.
Values: 0 latest syntax
1 syntax from Release 2.05 only
2 syntax from the first version of Release 2.25 only
For example, the word if is a key word in GAMS introduced with the first version of Release 2.25. Setting
the g205=1 option allows if to be used as an identifier since it was not a keyword in Release 2.05. As another
example, the word for is a key word in GAMS introduced with the later versions of Release 2.25. Setting the
g205=2 option allows for to be used as an identifier since it was not a keyword in the first version of Release 2.25.
Using values of 1 or 2 for g205 will not permit the use of enhancements to the language introduced
in the later versions.
The search order for the file myfile (or myfile.gms) and all included files in PC systems is as follows: current di-
rectory directories specified by inputdir (\mydir and \mydir2 directories) in order under Unix, the corresponding
command is
1. current directory
180 The GAMS Call
2. mydir1
3. mydir2
1. current directory
2. mydir2
3. mydir1
Note that it is not the order in which they are specified that matters but the number of the inputdir that they
have been assigned to.
Unlike idir, additional directories cannot be set with ldir. The string passed will be treated as one
directory. Passing additional directories will cause errors.
Note that if the ldir parameter is set, the default library include directory is not searched.
GAMS searches for any referenced $libinclude file in the directory mydir.
The log file that results from running [TRNSPORT] with the option ll=0 is shown below,
4 1.53675000E+02 nopt 0
SOLVER STATUS: 1 NORMAL COMPLETION
MODEL STATUS : 1 OPTIMAL
OBJECTIVE VALUE 153.67500
Comparing this output to the one shown in the example of option logfile, one can see that the line numbers
are absent from the log file.
To illustrate the use of the lo option, run [TRNSPORT] with the options lo=2. The resulting log file, trnsport.log,
looks exactly as shown in the example of option logfile.
By default, running a file containing just the two statements shown above results in the following listing file,
1 a(i) = b(i)*5 ;
**** $140$120$140
2 b(i) = c(j) ;
**** $140$120$149
None of the sets i, or j have been defined or initialized, and the identifiers a,b and c have not been defined.
Further, an assignment cannot be made without the right hand side of the assignment being known. In both the
assignments in the example above, there is no data available for the right hand side. Running the model with the
setting mp=1 results in the following listing file,
1 a(i) = b(i)*5 ;
2 b(i) = c(j) ;
**** $149
Error Messages
149 Uncontrolled set entered as constant
**** 1 ERROR(S) 0 WARNING(S)
Note that the statements in the example have now been processed independently of its context. They are now
checked only for consistency. GAMS now assumes that sets i and j, as well as the identifiers a, b, and c are
defined and, if necessary, initialized elsewhere. The only error that is reported is the inconsistency of indices in
the second statement.
A pool of n unique nonlinear constants is kept. Lookup for first n constants only.
Since there is no action called q, GAMS will complain and provide the following message:
model m /all/ ;
solve m using nlp maximizing dollars ;
The option file that is being used after this assignment is solvername.opt, where solvername is the name of the
solver that is specified. For CONOPT, the option file is called conopt.opt; for MINOS, it is minos.opt. The
names that you can use are listed in the Solver Manual.
Setting modelname.optfile in the GAMS input file overrides the value of the optfile parameter
passed through the command line.
To allow different option file names for the same solver, the optfile parameter can take other values as well.
Formally, the rule is optfile=n will use solvename.opt if n = 1, and solvername.opX, solvername.oXX or
solvername.XXX, where Xs are the characters representing the value of n, for n > 1 and will use no option file
at all for n = 0. This scheme implies an upper bound on n of 999. For example, the following optfile values
profile the option file names for the CONOPT solver
0 no option file used 1 conopt.opt
2 conopt.op2 26 conopt.o26
345 conopt.345
gams trnsport
gams trnsport o=trnsport.out
gams trnsport o=c:\test\trnsport.out
The first call will create an output file called trnsport.lst (for PC and Unix platforms) in the current directory.
The second call will create a file called trnsport.out in the current directory. The last call will create the file as
listed. If the directory c:\test does not exist, GAMS will exit with a parameter error.
Setting the profile option through the option statement in the GAMS input file overrides the value
of the profile parameter passed through the command line.
Values: 0 no profiling
1 minimum profiling
2 detailed profiling
A value of 0 does not cause an execution profile to be generated. A value of 1 reports execution times for each
statement and the number of set elements over which the particular statement is executed. A value of 2 reports
specific times for statements inside control structures like loops. Running [TRNSPORT] with profile=1 provides
the following additional information in the listing file,
---- 1 EXEC-INIT 0.010 0.010 SECONDS
---- 43 ASSIGNMENT C 0.000 0.010 SECONDS 6
---- 63 ASSIGNMENT TRANSPORT 0.000 0.010 SECONDS
---- 65 SOLVE INIT TRANSPORT 0.000 0.020 SECONDS
---- 56 EQUATION COST 0.000 0.020 SECONDS 1
---- 58 EQUATION SUPPLY 0.030 0.050 SECONDS 2
---- 60 EQUATION DEMAND 0.000 0.050 SECONDS 3
---- 65 SOLVE FINI TRANSPORT 0.040 0.090 SECONDS
---- 65 GAMS FINI 0.030 0.120 SECONDS
---- 1 EXEC-INIT 0.000 0.000 SECONDS
---- 65 SOLVE READ TRANSPORT 0.020 0.020 SECONDS
---- 67 DISPLAY 0.010 0.030 SECONDS
---- 67 GAMS FINI 0.000 0.030 SECONDS
The first column provides the line number in the input file of the statement being executed. The second column
provides the type of statement being executed.
EXEC-INIT denotes the beginning of the execution phase of the GAMS input file.
GAMS-FINI denotes the end of this phase.
Note that GAMS finishes processing of an input file as soon as a solve statement is processed, and passes control
to the solver being called. After the solver is done, GAMS restarts. This causes two EXEC-INIT GAMS-FINI pairs
to be generated for [TRNSPORT].
ASSIGNMENT C denotes an assignment statement involving the identifier c.
SOLVE INIT, SOLVE FINI are book ends enclosing the generation of the model [TRNSPORT].
Note that only equations are listed, and not variables. This happens because GAMS uses an equation based
scheme to generate a model. The third and fourth columns provide the individual time needed to execute the
statement, and the cumulative time taken by the GAMS system so far. The last column gives the number of
assignments generated in the specified line.
186 The GAMS Call
1 47 I SETS DECLARED 26 26 9 0 1
E:\WORK\TRNSPORT.GMS
2 47 I SETS DEFINED 26 26 29 0 1
E:\WORK\TRNSPORT.GMS
3 48 J SETS DECLARED 27 27 9 0 1
E:\WORK\TRNSPORT.GMS
4 48 J SETS DEFINED 27 27 29 0 1
E:\WORK\TRNSPORT.GMS
On Unix platforms the ? character is a special character and may require a backslash character (\)
in front of it in order to be interpreted correctly. The name myfile? should be written on this
platform as myfile\?.
To illustrate the use of the stepsum option, the default GAMS log file from running [TRNSPORT] with the option
stepsum=1 contains the following step summaries.
The first step summary occurs before the model is sent to the solver, and the second occurs after the solver
completes its task and returns control back to GAMS. The first column reports time for the individual section of
the run, while the second column reports accumulated times including previous sections.
The $offlisting and $onlisting dollar control options effect the listing file only if suppress is set
to 0. If suppress is set to 1, the input file(s) is not echoed to the listing file, and these dollar control
options have no effect on the listing file.
Unlike idir, additional directories cannot be set with sdir. The string passed will be treated as one
directory. Passing additional directories will cause errors.
Note that if the sdir parameter is set, the default system include directory is not searched.
GAMS searches for any referenced $sysinclude file in the directory mydir.
workdir (wdir=curdir) T
his option sets the working directory. This option is useful when GAMS is called from an external system like
Visual Basic. If not specified, it will be set to the curdir directory.
D.1 Introduction
The Dollar Control Options are used to indicated compiler directives and options. Dollar control options are not
part of the GAMS language and must be entered on separate lines recognized by a $ symbol in the first column.
A dollar control option line may be placed anywhere within a GAMS program and it is processed during the
compilation of the program. The $ symbol is followed by one or more options identified by spaced. Since the
dollar control options are not part of the GAMS language, they do not appear on the compiler listing unless an
error had been detected. Dollar control option lines are not case sensitive and a continued compilation uses the
previous settings.
D.1.1 Syntax
In general, the syntax in GAMS for dollar control options is as follows,
$option_name argument_list {option_name argument_list}
where option name is the name of the dollar control option, while argument list is the list of arguments for the
option. Depending on the particular option, the number of arguments required can vary from 0 to many.
No blank space is permitted between the $ character and the first option that follows.
In most cases, multiple dollar control options can be processed on a line. However, some dollar
control options require that they be the first option on a line.
The effect of the dollar control option is felt immediately after the option is processed.
Note that there is no blank space between the $ character and the option that follows. The first dollar control
option $title sets the title of the pages in the listing file to the text that follows the option name. In the second
line of the example above, two options are set - $onsymxref and $onsymlist. Both these options turn of the
echoing of the symbol reference table and listings to the listing file.
Table D.1 briefly describes the various options in each of the categories. Section D.3 contains a reference list of
all dollar control options in alphabetical order with detailed description for each.
Non-default settings are reported before the file summary at the end of a GAMS listing as a reminder for future
continued compilations. This is only relevant if a restart file has been requested with the GAMS call.
bracketed.
abort
This option will issue a compilation error and abort the compilation.
Consider the following example,
$if not %system.filesys% == UNIX
$abort We only do UNIX
This attempts to stop compilation if the operating system is not Unix. Running the above example on a
non-Unix platform results in the compilation being aborted, and the following listing file,
2 $abort We only do UNIX
**** $343
Error Messages
343 Abort triggered by above statement
batinclude
The $batinclude facility performs the same task as the $include facility in that it inserts the contents of
the specified text file at the location of the call. In addition, however, it also passes on arguments which
can be used inside the include file:
$batinclude file arg1 arg2 ...
The $batinclude option can appear in any place the $include option can appear. The name of the batch
include file may be quoted or unquoted, while arg1, arg2,.. are arguments that are passed on to the
batch include file. These arguments are treated as character strings that are substituted by number inside
the included file. These arguments can be single unbroken strings (quoted or unquoted) or quoted
multi-part strings.
The syntax has been modeled after the DOS batch facility. Inside the batch file, a parameter substitution
is indicated by using the character % followed immediately by an integer value corresponding to the order
of parameters on the list where %1 refers to the first argument, %2 to the second argument, and so on. If
the integer value is specified that does not correspond to a passed parameter, then the parameter flag is
substituted with a null string. The parameter flag %0 is a special case that will substitute a fully expanded
file name specification of the current batch included file. The flag %$ is the current $ symbol (see $dollar).
Parameters are substituted independent of context, and the entire line is processed before it is passed to
the compiler. The exception to this is that parameter flags appearing in comments are not substituted.
GAMS requires that processing the substitutions must result in a line of less than or equal to the
maximum input line length (currently 255 characters).
The case of the passed parameters is preserved for use in string comparisons.
5 display b ;
BATINCLUDE D:\GAMS\INC2.INC
7 b = sqr(c) - c ;
8 display b ;
BATINCLUDE D:\GAMS\INC2.INC
10 b = sqr(a+5) - a+5 ;
11 display b ;
Note that the three calls to $batinclude with the various arguments lead to GAMS interpreting the
contents of batch include file in turn as
b = sqr(a) - a ;
b = sqr(c) - c ;
b = sqr(a+5) - a+5 ;
Note that third call is not interpreted as sqr(a+5)-(a+5), but instead as sqr(a+5)-a+5. The results of
the display statement are shown at the end of the listing file,
---- 5 PARAMETER B = 0.000
---- 8 PARAMETER B = 2.000
---- 11 PARAMETER B = 40.000
The third call leads to b = sqr(6)-1+5 which results in b taking a value of 40. If the statement in the
batch include file was modified to read as follows,
%1 = sqr(%2) - %2 ;
The results of the display statement in the listing file would read,
---- 5 PARAMETER B = 0.000
---- 8 PARAMETER B = 2.000
---- 11 PARAMETER B = 30.000
The third call leads to b = sqr(6)-6 which results in b taking a value of 30.
call
Passes a followed string command to the current operating system command processor and interrupts
compilation until the command has been completed. If the command string is empty or omitted, a new
interactive command processor will be loaded.
Consider the following slice of code,
$call dir
This command makes a directory listing on a PC.
The command string can be passed to the system and executed directly without using a command
processor by prefixing the command with an = sign. Compilation errors are issued if the command or the
command processor cannot be loaded and executed properly.
$call gams trnsport
$call =gams trnsport
The first call runs [TRNSPORT] in a new command shell. The DOS command shell does not send any
return codes from the run back to GAMS. Therefore any errors in the run are not reported back. The
second call, however, sends the command directly to the system. The return codes from the system are
intercepted correctly and available to the GAMS system through the errorlevel DOS batch function.
Some commands (like copy on a PC and cd in Unix) are shell commands and cannot be spawned off
to the system. Using these in a system call will create a compilation error.
id1, id2, ... are the identifiers whose data is being reset. Note that this is carried out during compile
time, and not when the GAMS program executes. Not all data types can be cleared - only set,
parameter, equation and variable types can be reset.
Consider the following example,
set i /1*20/ ; scalar a /2/ ;
$clear i a
display i, a ;
The $clear option resets i and a to their default values. The result of the display statement in the listing
file shows that i is now an empty set, and a takes a value of 0.
---- 3 SET I
(EMPTY)
---- 3 PARAMETER A = 0.000
The two-pass processing of a GAMS file can lead to seemingly unexpected results. Both the dollar
control options and the data initialization is done in the first pass, and assignments in the second,
irrespective of their relative locations. This is an issue particularly with $clear since data can be
both initialized and assigned.
comment (*)
This option changes the start-of-line comment to a single character which follows immediately th command
keyword. This should be used with great care, and one should reset it quickly to the default symbol *.
The case of the start-of-line comment character does not matter when being used.
dollar ($)
This option changes the current $ symbol to a single character which follows immediately the $dollar
keyword. When the include file is inserted, all dollar control options are inherited, and the current $
symbol may not be known. The special %$ substitution symbol can be used to get the correct symbol (see
$batinclude).
Consider the following example,
$dollar #
#hidden now we can use # as the $ symbol
double
The lines following the $double statement will be echoed double spaced to the listing file.
Consider the following example,
198 Dollar Control Options
set i /1*2/ ;
scalar a /1/ ;
$double
set j /10*15/ ;
scalar b /2/ ;
The resulting listing file looks as follows,
1 set i /1*2/ ;
2 scalar a /1/ ;
4 set j /10*15/ ;
5 scalar b /2/ ;
Note that lines before the $double option are listed singly spaced, while the lines after the option are
listed with double space.
echo
The echo option allows to write text to a file:
$echo text > file
$echo text >> file
These options send the message text to the file file. Both the text and the file name can be quoted or
unquoted. The file name is expanded using the working directory. The $echo statement tries to minimize
file operations by keeping the file open in anticipation of another $echo to be appended to the same file.
The file will be closed at the end of the compilation or when a $call or any kind of $include statement is
encountered. The redirection symbols > and >> have the usual meaning of starting at the beginning or
appending to an existing file.
Consider the following example,
$echo > echo
$echo The message written goes from the first non blank >> echo
$echo to the first > or >> symbol unless the text is >> echo
$echo "is quoted. The Listing File is %gams.input%. The" >> echo
$echo file name "echo" will be completed with >> echo
$echo %gams.workdir%. >> echo
$echo >> echo
The contents of the resulting file echo are as follows,
The message written goes from the first non blank
to the first > or >> symbol unless the text is
is quoted. The Listing File is C:\PROGRAM FILES\GAMSIDE\CC.GMS. The
file name "echo" will be completed with
C:\PROGRAM FILES\GAMSIDE\.
eject
Advances the output to the next page.
Consider the following example,
$eject
Set i,j ;
Parameter data(i,j) ;
$eject
This will force the statements between the two $eject calls to be reported on a separated page in the
listing file.
eolcom (!!)
This option redefines the end-of-line comment symbol, which can be a one or two character sequence. By
default the system is initialized to !! but not active. The $oneolcom option is used to activate the
end-of-line comment. The $eolcom option sets $oneolcom automatically.
Consider the following example,
$eolcom ->
set i /1*2/ ; -> set declaration
parameter a(i) ; -> parameter declaration
The character set -> serves as the end-of-line-comment indicator.
GAMS requires that one not reset the $eolcom option to the existing symbol.
The following code is illegal since $eolcom is being reset to the same symbol as it is currently,
D.3 Detailed Description of Dollar Control Options 199
$eolcom ->
$eolcom ->
error
This option will issue a compilation error and will continue with the next line.
Consider the following example,
$if not exist myfile
$error File myfile not found - will continue anyway
This checks if the file myfile exists, and if not, it will generate an error with the comment File not
found - will continue anyway, and then compilation continues with the following line.
exit
This option will cause the compiler to exit (stop reading) from the current file. This is equivalent to
having reached the end of file.
Consider the following example,
scalar a ; a = 5 ;
display a ;
$exit end
a = a+5 ; display a ;
The lines following the $exit will not be compiled.
goto id
This option will cause GAMS to search for a line starting with $label id and then continue reading
from there. This option can be used to skip over or repeat sections of the input files. In batch include files,
the target labels or label arguments can be passed as parameters because of the manner in which
parameter substitution occurs in such files. In order to avoid infinite loops, one can only jump a maximum
of 100 times to the same label.
Consider the following example,
scalar a ; a = 5;
display a ;
$goto next
a = a+5 ; display a ;
$label next
a = a+10 ; display a ;
On reaching the $goto next option, GAMS continues from $label next. All lines in between are
ignored. On running the example, a finally takes a value of 15.
The $goto and $label have to be in the same file. If the target label is not found in the current file,
and error is issued.
hidden
This line will be ignored and will not be echoed to the listing file. This option is used to enter information
only relevant to the person manipulating the file.
Consider the following example,
$hidden You need to edit the following lines if you want to:
$hidden
$hidden 1. Change form a to b
$hidden 2. Expand the set
The lines above serve as comments to the person who wrote the file. However, these comments will not be
visible in the listing file, and is therefore hidden from view.
if
The $if dollar control option provides the greatest amount of control over conditional processing of the
input file(s). The syntax is similar to the IF statement of the DOS Batch language:
$if [not] {exist filename | string1 == string2} new_input_line
The syntax allows for negating the conditional with a not operator followed either of two types of
conditional expressions: a file operation or a string comparison. The result of the conditional test is used
to determine whether to read the remainder of the line, which can be any valid GAMS, input line.
200 Dollar Control Options
The exist file operator can be used to check for the existence of the given file name specification. The
string compare consists of two strings (quoted or unquoted) for which the comparison result is true only if
the strings match exactly. Null (empty) strings can be indicated by an empty quote: ""
The case of the strings provided either explicitly or, more likely, through a parameter substitution, is
preserved and therefore will effect the string comparison.
Quoted strings with leading and trailing blanks are not trimmed and the blanks are considered part
of the string.
If the string to be compared is a possibly empty parameter, the parameter operator must be quoted.
New input line is the remainder of the line containing the $if option, and could be any valid GAMS
input line.
The first non-blank character on the line following the conditional expression is considered to be the
first column position of the GAMS input line. Therefore, if the first character encountered is a
comment character the rest of the line is treated as a comment line. Likewise if the first character
encountered is the dollar control character, the line is treated as a dollar control line.
An alternative to placing new input line on the same line as the conditional is to leave the remainder of
the line blank and place new input line on the line immediately following the if line. If the conditional
is found to be false, either the remainder of the line (if any) is skipped or the next line is not read.
Consider the following example,
$if exist myfile.dat $include myfile.dat
The statement above includes the file myfile.dat if the file exists. Note that the $ character at the
beginning of the $include option is the first non-blank character after the conditional expression, if
exist myfile.dat and is therefore treated as the first column position. The above statement can also be
written as
$if exist myfile.dat
$include myfile.dat
Consider the following additional examples,
$if not "%1a" == a $goto labelname
$if not exist "%1" display "file %1 not found" ;
The first statement illustrates the use of the $if option inside a batch include file where parameters are
passed through the $batinclude call from the parent file. The $if condition checks if the parameter is
empty, and if not processes the $goto option. Note that the string comparison attempted, "%1a" == a,
can also be done using %1 == "".
The second statement illustrates using standard GAMS statements if the conditional is valid. If the file
name passed as a parameter through the $batinclude call does not exist, the GAMS display statement is
processed.
In line and end of line comments are stripped out of the input file before processing for
new input line. If either of these forms of comments appears, it will be treated as blanks.
Consider the following example,
parameter a ; a=10 ;
$eolcom ! inlinecom /* */
$if exist myfile.dat
/* in line comments */ ! end of line comments
a = 4 ;
display a;
The fourth line is ignored, and the fifth line involving an assignment setting a to 4 will be treated as the
result of the conditional. So the result of the display statement would be the listing of a with a value of 4
if the file myfile.dat exists, and a value of 10 if the file does not exist.
It is suggested that a $label not appear as part of the conditional input line. The result is that if
the $label appears on the $if line, a $goto to this label will re-scan the entire line thus causing a
reevaluation of the conditional expression. On the other hand, if the $label appears on the next
line, the condition will not be reevaluated on subsequent gotos to the label.
D.3 Detailed Description of Dollar Control Options 201
The following example illustrates how an unknown number of file specifications can be passed on to a
batch include file that will include each of them if they exist. The batch include file could look as follows,
/* Batch Include File - inclproc.bch */
/* Process and INCLUDE an unknown number of input files */
$label nextfile
$if "%1a" == a $goto end
$if exist "%1" $include "%1" /* name might have blanks */
$shift goto nextfile
$label end
The call to this file in the parent file could look like:
$batinclude inclproc.bch fil1.inc fil2.inc fil3.inc fil4.inc
include
The $include option inserts the contents of a specified text file at the location of the call. The name of
the file to be included which follows immediately the keyword include may be quoted or unquoted.
Include files can be nested.
The include file names are processed in the same way as the input file is handled. The names are
expanded using the working directory. If the file cannot be found and no extension is given, the standard
GAMS input extension is tried. However, if an incomplete path is given, the file name is completed using
the include directory. By default, the library include directory is set to the working directory. The default
directory can be reset with the idir command line parameter.
The start of the include file is marked in the compiler listing. This reference to the include file can be
omitted by using the $offinclude option.
The following example illustrates the use of an include statement,
$include myfile
$include "myfile"
Both statements above are equivalent, and the search order for the include file is as follows:
1. myfile in current working directory
2. myfile.gms in current working directory
3. myfile and myfile.gms (in that order) in directories specified by idir parameter.
The current settings of the dollar control options are passed on to the lower level include files.
However, the dollar control options set in the lower level include file are passed on to the parent file
only if the $onglobal option is set.
Compiler errors in include files have additional information about the name of the include file and the
local line number.
At the end of the compiler listing, an include file summary shows the context and type of include files.
The line number where an include file has been called is given. For example, in the Include File Summary
below we see that:
SEQ GLOBAL TYPE PARENT LOCAL FILENAME
1 1 INPUT 0 0 C:\TEST\TEST1.GMS
2 1 INCLUDE 1 1 .C:\TEST\FILE1.INC
3 6 INCLUDE 1 4 .C:\TEST\FILE2.INC
The first column named SEQ gives the sequence number of the input files encountered. The first row
always refers the parent file called by the GAMS call. The second column named GLOBAL gives the global
(expanded) line number which contained the $include statement. The third column named TYPE refers to
the type of file being referenced. The various types of files are INPUT, INCLUDE, BATINCLUDE, LIBINCLUDE,
and SYSINCLUDE. The fourth column named PARENT provides the sequence number of the parent file for
the file being referenced. The fifth column named LOCAL gives the local line number in the parent file
where the $include appeared. In the example listed above, the include files file1.inc and file2.inc
were included on lines 1 and 4 of the parent file test1.gms.
inlinecom (/* */)
This option redefines the in-line comment symbols, which are a pair of one or two character sequence. By
default, the system is initialized to /* and */, but is not active. The $oninline option is used to activate
the in-line comments. The $inlinecom option sets the $oninline automatically.
202 Dollar Control Options
GAMS requires that one not reset the $inlinecom option to an existing symbol.
The following code is illegal since $inlinecom is being reset to the same symbol as it is currently,
$inlinecom {{ }}
$inlinecom {{ }}
kill
Removes all data for an identifier and resets the identifier, only the type and dimension are retained. Note
that this is carried out during compile time, and not when the GAMS program executes. Not all data
types can be killed - only set, parameter, equation and variable types can be reset.
Consider the following example,
set i / 1*20 /; scalar a /2/
$kill i a
Note that this is different from $clear in that after setting $kill, i and a are treated as though they
have been only defined and have not been initialized or assigned. The result of the $kill statement above
is equivalent to i and a being defined as follows,
set i ; scalar a ;
Unlike the $clear, a display statement for i and a after they are killed will trigger an error.
label id
This option marks a line to be jumped to by a $goto statement. Any number of labels can be used in files
and not all of them need to be referenced. Re-declaration of a label identifier will not generate an error,
and only the first occurrence encountered by the GAMS compiler will be used for future $goto references.
Consider the following example,
scalar a ; a = 5 ;
display a ;
$goto next
a = a+5 ; display a ;
$label next
a = a+10 ; display a ;
On reaching the $goto next option, GAMS continues from $label next. All lines in between are
ignored. On running the example, a finally takes a value of 15.
The $label statement has to be the first dollar control option of multiple dollar control options that
appear on the same line.
libinclude
Equivalent to $batinclude:
$libinclude file arg1 arg2 ...
However, if an incomplete path is given, the file name is completed using the library include directory. By
default, the library include directory is set to the inclib directory in the GAMS system directory. The
default directory can be reset with the ldir command line parameter.
Consider the following example,
$libinclude abc x y
This call first looks for the include file [GAMS System Directory]/inclib/abc, and if this file does not
exist, GAMS looks for the file [GAMS System Directory]/inclib/abc.gms. The arguments x and y are
passed on to the include file to interpret as explained for the $batinclude option.
Consider the following example,
$libinclude c:\abc\myinc.inc x y
D.3 Detailed Description of Dollar Control Options 203
This call first looks specifically for the include file c:\abc\myfile.inc. The arguments x and y are passed
on to the include file to interpret as explained for the $batinclude option.
lines n
This option starts a new page in the listing file if less than n lines are available on the current page.
Consider the following example,
$hidden Never split the first few lines of the following table
$lines 5
table io(i,j) Transaction matrix
This will ensure that the if there are less than five lines available on the current page in the listing file
before the next statement (in this case, the table statement) is echoed to it, the contents of this statement
are echoed to a new page.
log
This option will send a message to the log file. By default, the log file is the console. The default log file
can be reset with the lo and lf command line parameters.
Leading blanks are ignored when the text is written out to the log file using the $log option.
All special % symbols will be substituted out before the text passed through the $log option is sent
to the log file.
maxcol n (255)
Sets the right margin for the input file. All valid data is before and including column n in the input file.
All text after column n is treated as comment and ignored.
Consider the following example,
$maxcol 30
set i /vienna, rome/ set definition
scalar a /2.3/ ; scalar definition
The text strings set definition and scalar definition are treated as comments and ignored since
they begin on or after column 31.
Any changes in the margins via maxcol or mincol will be reported in the listing file with the message that
gives the valid range of input columns. For example, the dollar control option $mincol 20 maxcol 110
will trigger the message:
NEW MARGIN = 20-110
204 Dollar Control Options
GAMS requires that the right margin set by $maxcol is greater than 15.
GAMS requires that the right margin set by $maxcol is greater than the left margin set by $mincol.
mincol n (1)
Sets the left margin for the input file. All valid data is after and including column n in the input file. All
text before column n is treated as comment and ignored.
Consider the following example,
$mincol 30
set definition set i /vienna, rome/
scalar definition scalar a /2.3/ ;
The text strings set definition and scalar definition are treated as comments and ignored since
they begin before column 30.
Any changes in the margins via maxcol or mincol will be reported in the listing file with the message that
gives the valid range of input columns. For example, the dollar control option $mincol 20 maxcol 110
will trigger the message:
NEW MARGIN = 20-110
GAMS requires that the left margin set by $mincol is smaller than the right margin set by $maxcol.
set i /1,2,3/ ;
set j(i) / / ;
parameter x(i) empty parameter / / ;
table y(i,i) headers only
1 2 3
;
$onempty
set k(i) / / ;
parameter xx(i) empty parameter / / ;
table yy(i,i)
1 2 3
;
The listing file that results looks like,
1 set i /1,2,3/ ;
2 set j(i) / / ;
**** $460
3 parameter x(i) empty parameter / / ;
**** $460
4 table y(i,i) headers only
5 1 2 3
6 ;
**** $462
8 set k(i) / / ;
9 parameter xx(i) empty parameter / / ;
10 table yy(i,i)
11 1 2 3
12 ;
Error Messages
460 Empty data statements not allowed. You may want to use $ON/OFFEMPTY
462 The row section in the previous table is missing
Note that empty data statements are not allowed for sets, parameters or tables. These are most likely
to occur when data is being entered into the GAMS model by an external program. Using the $onempty
dollar control option allows one to overcome this problem.
The empty data statement can only be used with symbols, which have a known dimension. If the
dimension is also derived from the data, the $phantom dollar control option should be used to
generate phantom set elements.
Note that the alternate syntax is more in line with syntax used in some of the popular programming
languages.
Both forms of the syntax will never be valid simultaneously. Setting the $onend option will make the
alternate syntax valid, but makes the standard syntax invalid.
INCLUDE D:\GAMS\INC.INC
2 $ondollar
3 $hidden text inside include file
INCLUDE D:\GAMS\INC.INC
7 $ondollar
8 $hidden text inside include file
9 $hidden after second call to include file
Note that the effect of the $ondollar dollar control option inside the include file does not affect the
parent file until $onglobal is turned on. The $hidden text is then echoed to the listing file.
[on|off ]include ($oninclude)
Controls the listing of the expanded include file name in the listing file.
Consider running the following slice of code,
$include inc.inc
$offinclude
$include inc.inc
where the file inc.inc contains the line,
$ondollar
$hidden text inside include file
The resulting listing file is as follows,
INCLUDE D:\GAMS\INC.INC
2 $ondollar
3 $hidden text inside include file
6 $ondollar
7 $hidden text inside include file
Note that the include file name is echoed the first time the include file is used. However, the include file
name is not echoed after $offinclude is set.
[on|off ]inline ($offinline)
Switch to control the use of in-line comments. By default, the in-line comments are set to the two
character pairs /* and */ but the processing is disabled. These comments can span lines till the
end-of-comment characters are encountered.
Consider running the following slice of code,
$oninline
/* the default comment symbols are now
active. These comments can continue
to additional lines till the closing
comments are found */
The two-pass processing of a GAMS file can lead to seemingly unexpected results. Both the dollar
control options and the data initialization is done in the first pass, and assignments in the second,
irrespective of their relative locations. This is an issue particularly with $onmulti since it allows
data initializations to be performed more than once.
Consider the following example,
scalar a /12/ ;
a=a+1;
$onmulti
scalar a /20/ ;
display a ;
The two scalar data initialization statements and the $onmulti option are processed before the
assignment statement a=a+1. In the order that it is processed, the example above is read by GAMS as,
* compilation step
scalar a /12/ ;
$onmulti
scalar a /20/ ;
* execution step
a=a+1;
display a ;
D.3 Detailed Description of Dollar Control Options 209
The example results in a taking a value of 21. The display statement in the resulting listing file is as
follows,
---- 5 PARAMETER A = 21.000
FUNCTIONS
**********
ABS
ARCTAN
CEIL
SETS
I canning plants
J markets
PARAMETERS
A capacity of plant i in cases
B demand at market j in cases
C transport cost in thousands of dollars per case
D distance in thousands of miles
F freight in dollars per case per thousand miles
VARIABLES
X shipment quantities in cases
Z total transportation costs in thousands of dollars
EQUATIONS
COST define objective function
DEMAND satisfy demand at market j
SUPPLY observe supply limit at plant i
MODELS
TRANSPORT
FILES
FILE Current file name for FILE.xxx use
PREDEFINED
DIAG
SAMEAS
This serves as a simple description of the symbols used in a model, and can be used in reports and other
documentation.
[on|off ]symxref ($onsymxref)
This option controls the following,
Collection of cross references for identifiers like sets, parameters, and variables.
Cross-reference report of all collected symbols in listing file
Listing of all referenced symbols and their explanatory text by symbol type in listing file. This is
also reported by using $onsymlist.
Consider the following slice of code,
210 Dollar Control Options
$offsymxref
set j(i) will not show / 1*3 /
display i;
$onsymxref
k(1) = yes;
The resulting listing file will contain the following symbol reference reports,
SYMBOL TYPE REFERENCES
I SET DECLARED 1 DEFINED 1 REF 1
K SET DECLARED 1 ASSIGNED 6
SETS
I this is set declaration
K some more
[on|off ]text
The $ontext - $offtext pair encloses comment lines. Line numbers in the compiler listing are suppressed
to mark skipped lines.
Consider the following,
* standard comment line
$ontext
Everything here is a comment
until we encounter the closing $offtext
like the one below
$offtext
* another standard comment line
The resulting listing file is as follows,
1 * standard comment line
Everything here is a comment
until we encounter the closing $offtext
like the one below
7 * another standard comment line
GAMS requires that every $ontext has a matching $offtext, and vice versa.
$onupper
* now we list everything in upper case
$offupper
* now we are back to list lines as entered
The resulting listing file is as follows,
2 * NOW WE LIST EVERYTHING IN UPPER CASE
4 * now we are back to list lines as entered
Note that all the characters in the lines between $onupper and $offupper are capitalized.
This switch affects three types of domain errors usually referred to as error numbers 116, 170 and
171.
This can have serious side affects and one has to exercise great care when using this feature.
Error Messages
116 Label is unknown
170 Domain violation for element
**** 0 ERROR(S) 6 WARNING(S)
E x e c u t i o n
---- 7 SET I
one , two , three
---- 7 SET J
four, five
$phantom null
set i / null/
j / a,b,null/ ;
display i,j ;
The resulting section of the listing file is shown below,
---- 4 SET I
(EMPTY)
---- 4 SET J
a, b
a 1.000
shift
The $shift option is similar to the DOS batch shift operator. It shifts the order of all parameters passed
once to the left. This effectively drops the lowest numbered parameter in the list.
Consider the following example,
scalar a, b, c ; a = 1 ;
$batinclude inc.inc a b c
display a, b, c ;
where the batch include file inc.inc is as follows,
%2 = %1 + 1 ;
$shift
%2 = %1 + 1 ;
The resulting listing file is,
1 scalar a, b, c ; a = 1 ;
BATINCLUDE C:\PROGRAM FILES\GAMSIDE\INC.INC
3 b = a + 1 ;
5 c = b + 1 ;
6 display a, b, c ;
In the first statement in the include file, %1 is the first argument in the $batinclude call and is
interpreted in this case as a. %2 is the second argument in the $batinclude call and is interpreted as b.
This leads to the overall assignment being interpreted as b=a+1.
The $shift option shifts the arguments to the left. So now, %1 is interpreted as b, and %2 is interpreted as
c. This leads to the second assignment being interpreted as c=b+1.
The result of the display statement in the input file is therefore,
---- 6 PARAMETER A = 1.000
PARAMETER B = 2.000
PARAMETER C = 3.000
single
The lines following a $single option will be echoed single spaced on the compiler listing. This option is
the default, and is only useful as a switch to turn off the $double option.
Consider the following example,
set i /1*2/ ;
scalar a /1/ ;
$double
set j /10*15/ ;
scalar b /2/ ;
$single
set k /5*10/ ;
scalar c /3/ ;
D.3 Detailed Description of Dollar Control Options 213
4 set j /10*15/ ;
5 scalar b /2/ ;
7 set k /5*10/ ;
8 scalar c /3/ ;
Note that lines between the $double and $single options are listed double spaced, while the lines after
the $single option revert back to being listed singly spaced.
stars (****)
This option is used to redefine the **** marker in the GAMS listing file. By default, important lines like
those denote errors, and the solver/model status are prefixed with ****.
Consider the following example,
$stars ####
garbage
The resulting listing file looks as follows,
2 garbage
#### $140
#### $36,299 UNEXPECTED END OF FILE (1)
Error Messages
36 = or .. or := or $= operator expected
rest of statement ignored
140 Unknown symbol
299 Unexpected end of file
stitle
This option sets the subtitle in the page header of the listing file to text which follows immediately the
keyword stitle. The next output line will appear on a new page in the listing file.
Consider the following example,
$stitle data tables for input/output
sysinclude
Equivalent to $batinclude:
$sysinclude file arg1 arg2 ...
However, if an incomplete path is given, the file name is completed using the system include directory. By
default, the system include directory is set to the GAMS system directory. The default directory can be
reset with the sdir command line parameter.
Consider the following example,
$sysinclude abc x y
This call first looks for the include file [GAMS System Directory]/abc, and if this file does not exist,
looks for [GAMS System Directory]/abc.gms. The arguments x and y are passed on to the include file to
interpret as explained for the $batinclude option.
Consider the following example,
$sysinclude c:\abc\myinc.inc x y
This call first looks specifically for the include file c:\abc\myfile.inc.
title
This option sets the title in the page header of the listing file to text which follows immediately the
keyword title. The next output line will appear on a new page in the listing file.
Consider the following example,
$title Production Planning Model
$stitle Set Definitions
use205
This option sets the GAMS syntax to that of Release 2.05. This is mainly used for backward
compatibility. New keywords have been introduced in the GAMS language since Release 2.05. Models
214 Dollar Control Options
developed earlier that use identifiers that have since become keywords will cause errors when run with the
latest version of GAMS. This option will allow one to run such models.
Consider the following example,
$use205
set if /1.2.3/; scalar x ;
The word if is a keyword in GAMS introduced with the first version of Release 2.25. The setting of the
$use205 option allows if to be used as an identifier since it was not a keyword in Release 2.05.
use225
This option sets the GAMS syntax to that of first version of Release 2.25. This is mainly used for
backward compatibility. New keywords have been introduced in the GAMS language since the first version
of Release 2.25. Models developed earlier that use identifiers that have since become keywords will cause
errors when run with the latest version of GAMS. This option will allow one to run such models.
Consider the following example,
$use225
set for /1.2.3/; scalar x ;
The word for is a keyword in GAMS introduced with the later versions of Release 2.25. The setting of the
$use225 option allows for to be used as an identifier since it was not a keyword in the first version of
Release 2.25.
use999
This option sets the GAMS syntax to that of the latest version of the compiler. This option is the default.
Consider the following example,
$use225
set for /1.2.3/; scalar x ;
$use999
for (x=1 to 3, display x) ;
The word for is used as a set identifier by setting the option $use225, and later the keyword for is used
as a looping construct by setting the language syntax to that of the latest version by setting $use999.
E
E.1 Introduction
The option statement is used to set various global system parameters that control output detail, solution process
and the layout of displays. They are processed at execution time unlike the dollar control options discussed in
Appendix D. They are provided to give flexibility to the user who would like to change the way GAMS would
normally do things. GAMS does provide default values that are adequate for the most purposes, but there are
always cases when the user would like to maintain control of aspects of the run.
where the keyword1 and keyword2 are recognized option names (but not reserved words) and the value1
and value2 are valid values for each of the respective options. Note that commas or end-of-line characters are
both legal separators between options.
Option names are not reserved words and therefore do not conflict with other uses of their name.
1. a display specifier.
2. a recognized name, number following = sign or value
3. a recognized name, number following an = sign, then an unsigned integer value.
4. a recognized name, number following an = sign, then an unsigned real number.
5. a recognized name, number following an = sign, then either of two recognized words.
The values associated with an option can be changed as often as necessary, with the new value
replacing the older one each time.
option profit:0:3:2;
option eject
iterlim = 100 , solprint = off ;
solve mymodel using lp maximizing profit ;
display profit.l ;
input("val1") = 5.3 ;
option iterlim = 50 ;
solve mymodel using lp maximizing profit ;
The option statement in the second line affects the display format of the identifier profit. More details on this
option can be found under the heading <identifier> in the following section. The option on the second line has
no value associated with it, and serves to advance the output in the listing file to the next page. The third line
contains two options - iterlim, and solprint. The values associated with the two options on the fourth line are
of different types - iterlim has an integer value while solprint requires a character string as a value. Note also
that the end of line and the comma serve as legal separators between two options.
The option iterlim serves to limit the number of iterations taken by the solver while attempting to solve the
lp model mymodel. After mymodel is solved for the first time, some of the input data is changed and the model
is solved again. However, before the second solve statement, the option iterlim is changed to 50. The effect of
the sequence above is to limit the first solve to less than 100 iterations and the second to less than 50.
Table E.1 briefly describes the various options in each of the categories. Section E.3 contains a reference list of
all options available through the option statement in alphabetical order with detailed description for each.
Options controlling output detail:
<identifier> controls print format profile lists program execution profile
decimals global control of print format profiletol sets tolerance for execution profile
eject advances output to next page solprint controls printing of solution
limcol number of columns listed solslack controls type of equation information
limrow number of rows listed sysout controls printing of solver status file
Options controlling solver specific parameters:
bratio use of advanced basis optca sets absolute optimality tolerance
domlim limits number of domain errors optcr sets relative optimality tolerance
iterlim limits number of solver iterations reslim limits amount of solver time
Options controlling choice of solver:
cns sets solver for cns model type mip sets solver for mip model type
dnlp sets solver for dnlp model type mpec sets solver for mpec model type
lp sets solver for lp model type nlp sets solver for nlp model type
mcp sets solver for mcp model type rminlp sets solver for rminlp model type
minlp sets solver for minlp model type rmip sets solver for rmip model type
Options affecting input program control:
seed resets seed for pseudo random number generator
solveopt controls return of solution values to GAMS
<identifier>
Display specifier: identifier:d, identifier:d:r:c Defines print formats for identifier when used in a
display statement. d is the number of decimal places, r is the number of index positions printed as row
labels, c is the number of index positions printed as column labels; the remaining index positions (if any)
will be used to index the planes (index order: plane, row, column); if r is zero list format will be used.
The default setting is described in Section 14.4.
bratio (0.25)
Certain solution procedures can restart from an advanced basis that is constructed automatically. This
option is used to specify whether or not basis information (probably from an earlier solve) is used. The
use of this basis is rejected if the number of basic variables is smaller than bratio times the size of the
basis. Setting bratio to 1 will cause all existing basis information to be discarded, which is sometimes
needed with nonlinear problems. A bratio of 0 accepts any basis, and a bratio of 1 always rejects the
basis. Setting bratio to 0 forces GAMS to construct a basis using whatever information is available. If
bratio has been set to 0 and there was no previous solve, an all slack (sometimes called all logical )
basis will be provided. This option is not useful for MIP solvers.
Range:
0-1
cns (default)
The default cns solver is set during installation. The user can change this default by setting this option to
the required solver. The list of cns solvers available with your system can be obtained by reading the
gamscomp.txt file that is present in the GAMS system directory. A value of default will change the cns
solver back to the default one as specified in gamscomp.txt.
decimals (3)
Number of decimals printed for symbols not having a specific print format attached.
Range:
0-8
dnlp (default)
This option controls the solver used to solve dnlp models. For details cf. option cns.
domlim (0)
This controls the number of domain errors (undefined operations like division by zero) a nonlinear solver
will perform, while calculating function and derivative values, before it terminates the run. Nonlinear
solvers have difficulty recovering after attempting an undefined operation.
eject
Advances output in the listing file to the next page.
iterlim (1000)
This option will cause the solver to interrupt the solution process after iterlim iterations and return the
current solution values to GAMS.
limcol (3)
This controls the number of columns that are listed for each variable in the COLUMN LISTING section of the
listing file. Specify zero to suppress the COLUMN LISTING altogether.
limrow (3)
This controls the number of rows that are listed for each equation in the EQUATION LISTING section of the
listing file. Specify zero to suppress the EQUATION LISTING altogether.
lp (default)
This option controls the solver used to solve lp models. For details cf. option cns.
mcp (default)
This option controls the solver used to solve mcp models. For details cf. option cns.
218 The Option Statement
minlp (default)
This option controls the solver used to solve minlp models. For details cf. option cns.
mip (default)
This option controls the solver used to solve mip models. For details cf. option cns.
miqcp (default)
This option controls the solver used to solve miqcp models. For details cf. option cns.
nlp (default)
This option controls the solver used to solve nlp models. For details cf. option cns.
optca (0.0)
This option is only used with problems containing discrete variables (i.e. the GAMS model type mip).
General mixed integer problems are often extremely difficult to solve, and proving that a solution found is
the best possible can use enormous amounts of resources. This option sets an absolute termination
tolerance, which means that the solver will stop and report on the first solution found whose objective
value is within optca of the best possible solution.
optcr (0.1)
This option sets a relative termination tolerance for problems containing discrete variables, which means
that the solver will stop and report on the first solution found whose objective value is within 100*optcr
of the best possible solution.
profile (0)
This option is used to generate more information on program execution profiles. This option is equivalent
in function to the profile command line parameter.
0 No execution profile is generated in listing file
1 The listing file reports execution times for each statement and the number of set elements
over which the particular statement is executed.
2 Specific times for statements inside control structures like loops.
profiletol (0.0)
This option sets profile tolerance in seconds. All statements that take less time to execute than this
tolerance are not reported in the listing file.
qcp (default)
This option controls the solver used to solve qcp models. For details cf. option cns.
reslim (1000)
This option causes the solver to terminate the solution process after reslim units of processor time have
been used, and the current solution values are returned to GAMS. The units are seconds on the wall
clock for PCs, or CPU seconds for larger machines. The SOLUTION SUMMARY part of the listing file shows
the limit was used.
rmip (default)
This option controls the solver used to solve rmip models. For details cf. option cns.
rminlp (default)
This option controls the solver used to solve rminlp models. For details cf. option cns.
seed (3141)
This option resets the seed for the pseudo random number generator.
solprint (on)
This option controls the printing of the model solution in the listing file. Using this specification
suppresses the list of the solution following a solve.
on The solution is printed one line per row and column in the listing file.
off Solution details are not printed. Although this saves paper, we do not recommend it unless
you understand your model very well and solve it often.
E.3 Detailed Description of Options 219
solslack (0)
This option causes the equation output in the listing file to contain slack variable values instead of level
values.
0 Equation output in listing file contains level values between lower and upper bound values
1 Equation output in listing file contains slack values between lower and upper bound values
solveopt (merge)
Controls the way solution values from a solve are returned to GAMS.
merge Old and new values merged together, and new values over-write old ones
replace All old values associated with a variable or equation are reset to default values before new
solution values are returned
sysout (off)
This option controls the printing of the solver status file as part of the listing file. The contents of the
solver status file are useful if you are interested in the behavior of the solver. If the solver crashes or
encounters any difficulty, the contents of the solver status file will be automatically sent to the listing file.
on Prints the system output file of the solver
off No subsystem output appears on output file unless a subsystem error has occurred.
220 The Option Statement
F
F.1 Introduction
GAMS saves the information provided in the input files in intermediate, mostly binary, files. These files are
referred to as work files or scratch files. Some of these files are used to exchange information between GAMS and
the various solvers. Just before a GAMS run is complete, these files are usually deleted.
Input files can be processed in parts through the use of these intermediate files. This is an extremely powerful
feature that can reduce the time needed when several runs of the same model are being made with different data.
It may be clearer if the process is described in a different way. Imagine taking a large GAMS program and running
it, producing one output file. Then think of splitting the program into two pieces. The first piece is run and the
resulting work file is saved along with the resulting listing file. Then the second piece is run after reading in the
data from the work file saved previously. A new listing file is generated for the second piece. The content of the
output that results is the same, though slightly rearranged, as the case when the large file was run. Splitting the
files allows one to interrupt a GAMS task and restart it later without loss of information. Further, changes could
be made, or errors corrected, to the later parts.
Sets
i "canning plants" / seattle, san-diego /
j "markets" / new-york, chicago, topeka / ;
Parameters
a(i) "capacity of plant i in cases"
/ seattle 350
san-diego 600 /
Variables
x(i,j) "shipment quantities in cases"
z "total transportation costs in 1000$" ;
Positive Variable x ;
Equations
cost "define objective function"
supply(i) "observe supply limit at plant i"
demand(j) "satisfy demand at market j" ;
Note that [TRNSPORT] results from appending file2.gms at the end of file1.gms.
The Work file preserves all information (including declarations, values, option settings and compiler
dollar directives) known to GAMS at the end of the run that created them.
The work file is not machine specific, it is portable between platforms. For example, a work file
generated on a PC running Windows can be re-used on a Sun machine running Solaris.
GAMS reads the work file named trans.g00 and regenerates the information stored in file1.gms. Then
file2.gms is run and the result is as if the two files were concatenated.
A restarted run also requires a continuation input file. The restart does not alter work files. They can be used
repeatedly to continue a particular run many times, possibly with many different continuation input files.
It is the responsibility of the modeler to ensure that the contents of the input file matches that of the
work file, although the compiler will issue errors if it detects any inconsistencies, such as references
to symbols not previously declared.
F.3 Ways in which a Work File is Useful 223
A Work file can be used only by GAMS tasks requesting a restarted run.
A File can be saved following a restarted run, thus producing another work file that reflects the state
of the job following completion of the statements in the continuation file.
Positive Variable x ;
Note that this representation does not contain any data, and is a purely algebraic representation of the trans-
portation problem. Running this model and saving the resulting work file will allow the model to be used with
the data stored in a separate file file2.gms.
Table d(i,j)
new-york chicago topeka
seattle 2.5 1.7 1.8
san-diego 2.5 1.8 1.4 ;
Scalar f / 90 / ;
224 The Save and Restart Feature
This file contains the data for the model and the solve statement.
G.1 Introduction
This document describes the GDX (GAMS Data Exchange) facilities available in GAMS. In addition to these
facilities, there are a few utilities to work with GDX files.
A GDX file is a file that stores the values of one or more GAMS symbols such as sets, parameters variables and
equations. GDX files can be used to prepare data for a GAMS model, present results of a GAMS model, store
results of the same model using different parameters etc. A GDX file does not store a model formulation or
executable statements.
GDX files are binary files that are portable between different platforms. They are written using the byte ordering
native to the hardware platform they are created on, but can be read on a platform using a different byte ordering.
Users can also write their own programs using GDX files by using the gdxio.dll (gdxio.so shared objects). The
interface and usage for these libraries is described in a separate document.
Compression: Starting with version 22.3 of GAMS, gdx files can be written in a compressed format. Compression
is controlled by the environment variable GDXCOMPRESS. A non-zero value indicates compression.
GDX files can be converted to a compressed format too.
During compilation, we can use dollar control options to specify the GDX file and the symbols to read or write.
Reading during the compilation phase also allows us to define the elements of a set and the subsequent use of
such a set as a domain.
The following directives are available for reading data from a GDX file into GAMS during compilation of a GAMS
model:
226 GDX Facilities and Tools
When reading data, the symbol to be read has to be defined in GAMS already
Compile Phase Example 1 The trnsport.gms model has been modified to use the demand data from an
external source. Only the relevant declarations are shown.
The parameter B is read from the GDX file using the name demand, and only those elements that are in the
domain J will be used. Values for parameter B that are outside the domain J will be ignored without generating
any error messages.
*Example 1
Set
j markets / new-york, chicago, topeka / ;
Parameter
B(j) demand at market j in cases ;
$GDXIN demanddata.gdx
$LOAD b=demand
$GDXIN
Compile Phase Example 2 In this example, the set J is also read from the GDX file, and is used as the
domain for parameter B. All elements read for the set J will be used. Values for the parameter B that are outside
the domain J will be ignored. Note that the dimension of set J is set to one by specifying its domain.
*Example 2
$GDXIN demanddata.gdx
Set
J(*) markets;
$LOAD j=markets
Parameter
B(j) demand at market j in cases ;
$LOAD b=demand
$GDXIN
*Example 3
$GDXIN trnsport.gdx
$LOAD
When writing data, an existing GDX file will be overwritten with the new data; there is no merge or
append option.
Read a solution:
228 GDX Facilities and Tools
execute\_loadpoint filename;
execute_loadpoint filename,id1,id2=gdxid2,..;
The execute loadpoint allows you to merge solution points into any GAMS database. Loading the data acts
like an assignment statement and it merges/replaces data with current data. For variables and equations, only
level values and marginal values (.L and .M) are used. If no symbols are specified, all symbols that match in type
and dimensionality will be loaded into the GAMS database.
The gdx file that can be used is not limited to files created with SAVE POINT; any gdx file can be used.
Example 4:
This example again uses the trnsport.gms model. After solving the model, we write the sets I and J and the
variables Z and X to the GDX file:
*Example 4
Set I /. . ./,
J / . . . /;
Variable X(I,J),
Z;
. . .
Solve trnsport using LP minimizing Z;
Execute_Unload results.gdx,I,J,Z,X;
Using the GDX option in the GAMS call, will cause all sets, parameters, variables and equations to be written
to the GDX file.
For example:
Gams trnsport gdx=trnsport
or
Gams trnsport a=c gdx=trnsport
Using the GDX parameter when running the model using the GAMSIDE, the process window will show the GDX
filename in blue indicating that the file can be opened using a double-click with the mouse.
GDXXRW allows reading and writing of an Excel spreadsheet. This utility requires the presence of
Microsoft Excel and therefore can only be used on a PC running the Windows operating system with
Microsoft Excel installed.
GDXDUMP writes the contents of a GDX file as a GAMS formatted text file.
GDXDIFF compares the data of symbols with the same name, type and dimension in two GDX files
and writes the differences to a third GDX file.
GDXMERGE combines multiple GDX files into one file. Symbols with the same name, dimension and
type are combined into a single symbol of a higher dimension. The added dimension has the file name
of the combined file as its unique element.
GDXRANK reads one or more one dimensional parameters from a GDX file, sorts each parameter and
writes the sorted indices as a one dimensional parameters to the output GDX file.
GDXCOPY Copies/Converts one or more GDX files to a different format.
G.5 GDXXRW
GDXXRW is a utility to read and write Excel spreadsheet data. GDXXRW can read multiple ranges in a
spreadsheet and write the data to a GDX file, or read from a GDX file, and write the data to different ranges in
a spreadsheet.
Usage:
gdxxrw Inputfile {Outputfile} {Options} [Symbols]
Parameters can also be read from a text file; the use of a file for parameters is indicated by preceding the file
name with a @ (At sign). When reading parameters from a text file, lines starting with an * (asterisk) will be
ignored and act as a comment.
Parameters can also be read from an area in a spreadsheet.
Files without a full path name are assumed to be in the current directory when using a DOS command prompt.
When using the GAMSIDE, these files are assumed to be in the current project directory. The use of file names
with embedded blanks is allowed as long as the file name is enclosed in doublequotes ().
Immediate options are command line options recognized independant of their position on the command line.
These options are global and they can only be specified once
The file extension of the input file is required and determines the action taken by the program. The extension
.gdx for the input file will read data from a GDX file and write data to a spreadsheet. The extension .xls
for the input file will read a spreadsheet and write the data to a .gdx file. In addition to the .xls input file
extension, the following file extensions are also valid for spreadsheet input files: .wk1, .wk2, .wk3 and .dbf.
A file sharing conflict will arise when writing to a spreadsheet with the target file open in Excel. Either close the
file in Excel before executing GDXXRW, or mark the spreadsheet as a shared workbook in Excel. To change the
shared status of a workbook, use the Excel commands available under: Tools Share Workbook.
Writing to a shared workbook can be painfully slow; simply closing the file and reopen the file after GDXXRW
is finished is often a better option.
GDXXRW Symbols
To write data to a spreadsheet or to a GDX file, one or more symbols and their associated Excel range need to
be specified. See also Excel Ranges
The general syntax for a Symbol specification is:
DataType = SymbolName {DataRange}{Dimensions}{SymbolOptions}
DataType
Par = GAMS Param
Specify a GAMS parameter to be read from a GDX file and written to spreadsheet, or to be read from a
spreadsheet and written to a GDX file.
When writing to a spreadsheet, special values such as Eps, NA and Inf will be written in ASCII. When
reading data from a spreadsheet, the ASCII strings will be used to write the corresponding special values
to the GDX file.
Equ = GAMS Equation
The following options apply to the symbol preceding the option, and only affect that symbol:
DataRange (Optional)
Rng = Excel Range
The Excel Range for the data for the symbol. Note that an empty range is equivalent to the first cell of
the first sheet.
Dimensions (Optional)
Dim = integer
The total dimension for the symbol
Cdim = Integer
Column dimension: the number of rows in the data range that will be used to define the labels for the
columns. The first Cdim rows of the data range will be used for labels.
Rdim = Integer
Row dimension: the number of columns in the data range that will be used to define the labels for the
rows. The first Rdim columns of the data range will be used for the labels.
IndexRange (Optional)
Index = Excel Range
The Index option is used to obtain the parameters by reading from the spreadsheet directly. The
parameters are read using the specified range, and treated as if they were specified directly on the
command line. The first three columns of the range have a fixed interpretation: DataType, Symbol
identifier and Data range. The fourth and following columns can be used for additional parameters. The
column header contains the keyword when necessary, and the Cell content is used as the option value.
G.5 GDXXRW 233
GDXXRW Options
The following options affect the symbols that follow the option. They remain in effect unless they are used again
for another symbol.
EpsOut = string (default = Eps)
String to be used when writing the value for Epsilon.
IncRC = flag (default = N)
Valid only when reading a spreadsheet. Include Excel row and column indices when a symbol is written to
the GDX file. For example, when we write a parameter P with indices I and J, without this option it will
be written a P(I, J). When IncRC is enabled, the parameter will be written as P(Excel Rows, I,
Excel Columns, J). Note that the sets Excel Rows and Excel Columns will be added the GDX file
automatically.
MinfOut = string (default = -Inf )
String to be used when writing the value for Negative infinity.
NameConv = flag or NC = flag
q The naming convention parameter is used to change the interpretation of an Excel range that does not
contain an ! (exclamation mark). For details see Ranges below. The default value is false.
NaOut = string (default = NA)
String to be used when writing the value for Not available.
PinfOut = string (default = +Inf )
String to be used when writing the value for Positive infinity.
ResetOut
Reset the output strings for special values to their defaults.
Squeeze = flag (default = Y) or SQ = flag
Writing to a spreadsheet: The squeeze option affects the writing of sub-fields of variables and equations. A
value for the field that is the default value for that type of variable or equation will not be written to the
spreadsheet. For example, the default for .L (Level value) is 0.0, and therefore zero will not be written to
the spreadsheet. When we set SQ=N, all values will be written to the spreadsheet.
Reading a spreadsheet: When the squeeze option is enabled, zero values for parameters will not be written
to the GDX file. When the squeeze option is disabled, zero values will be written to the GDX file. In
either case, empty cells, or cells containing blanks only, will never be written to the GDX file.
SkipEmpty = integer (default = 1) or SE = integer
The SkipEmpty option can be used when reading a spreadsheet, and the range is specified using the top
left corner instead of a block range. The value defines the number of empty row or column cells signal the
end of a block. Valid values are 0..n.
UndfOut = string (default = Undf )
String to be used when writing the value for Undefined.
UseAllUELs = flag (default = Y)
Valid only when reading a spreadsheet. When enabled, all unique elements found in a range will be
entered in the GDX file. When disabled, only those unique elements that are used in conjunction with a
value will be entered in the GDX file.
ZeroOut = string (default = 0)
String to be used when writing the value for Zero; by default this is 0.
Syntax elements
GDXXRW Ranges
Warning: When executing gdxxrw.exe twice and redirecting output to the same log file may result in a fatal error.
For example:
Gdxxrw step1 parameters > logfile
Gdxxrw step2 parameters > logfile
The execution of step2 may fail, because Excel will close the logfile in step1 in a delayed fashion, but return
control to gdxxrw.exe immediately. Using the Log or LogAppend parameter will avoid this problem.
Assuming we want to read parameter Data1 from the file test1.xls and write the data to test1.gdx. The sheet
name in a range can be omitted when it refers to the first sheet.
A B C D E
1 a1 a2 a3
2 i1 1 2 3
3 i2 4 5 6
4
The same data as in the previous example, but organized differently. We use the Dset option to read set I (in
column A) and set A (in column B).
G.5 GDXXRW 235
A B C D
1 i1 a1 1
2 i1 a2 2
3 i1 a3 3
4 i2 a1 4
5 i2 a2 5
6 i2 a3 6
7
GDXXRW test1.xls par=Data2 rng=EX2!A1 Rdim=2 Dset=I rng=EX2!A1 Rdim=1 Dset=A rng=EX2!B1 Rdim=1
When using a few symbols, the command line can become too long to be practical. In such case, use a text file to
hold the parameters. A parameter file can contain multiple lines to increase readability and a line starting with
a * will be ignored.
*file example6.txt
par =Data2 rng=EX2!A1 RDim=2
Dset=I rng=EX2!A1 Rdim=1
Dset=A rng=EX2!B1 Rdim=1
A parameter file can also be written during the execution of a GAMS model using the GAMS PUT
A B C D E F G
1 q1 q1 q2 q2
2 r1 r2 r1 r2
3 a1 b1 1 2 3 4
4 a1 b2 5 6 7 8
5 a2 b1 9 10 11 12
6 a2 b2 13 14 15 16
7
Special values can be read and written; the division by zero error in the spreadsheet will be written as Undefined.
A B C D E F G
1 v1 v2 v3 v4 v5 v6
2 Eps NA Eps +inf -inf #DIV/0!
3
Example of reading a set; the result will only contain the element s1.
A B C
1
2 s1 Y
3 s2 N
4 s3 No
5 s4 0
6 s5
7
The Index option is used to read a number of parameters and sets based on information stored in the spreadsheet
itself. The first row of the range is used for column headings indicating additional parameters.
A B C D E F G
1 Dim Rdim
2 par D1 Ex1a1 2
3 par D2 Ex2a1 3
4 par D2 Ex3a1 4 2
5 set S1 Ex5a1 1
6 set S2 Ex5a1:a6 1 Values=NoData
7
First, we create a GDX file using the GDX parameter in the GAMS call:
*file makedata.gms
set i /i1*i4/
j /j1*j4/
k /k1*k4/;
parameter v(i,j,k);
v(i,j,k)$(uniform(0,1) < 0.30) = uniform(0,1);
When we run this GAMS model, the file test2.gdx will be created at the end of the run.
GAMS makedata gdx=test2
Using the file test2.gdx, we can write to a spreadsheet:
Write parameter V to the first cell in the first sheet; because we only specify the top left corner of the sheet, the
complete sheet can be used to store the data. We do not specify the row and column dimension, so they will
default to rdim=2 and cdim=1.
Before executing this example, open the Excel file (test2.xls) and use the Excel Tools menu to make this a shared
notebook. After writing to the spreadsheet, use the Excel File Save command to verify the changes made.
GDXXRW test2.gdx par=V rng=a1
The steps above can be combined in a single GAMS model using the Execute Unload and Execute statements as
follows:
G.5 GDXXRW 237
set i /i1*i4/
j /j1*j4/
k /k1*k4/;
parameter v(i,j,k);
v(i,j,k)$(uniform(0,1) < 0.30) = uniform(0,1);
Execute_Unload "test2.gdx",I,J,K,V;
Execute GDXXRW.EXE test2.gdx par=V rng=a1;
The second sheet of this spreadsheet contains a number of labels to illustrate the use of the merge option. Note
that the values written are no longer in the same cells because they have been matched with the column and row
labels. Cells that were not changed by the merge option still contain xxx.
GDXXRW test2.gdx par=V rng=sheet2!a1 merge
In the previous example, the cells that were not changed still contained xxx. We can clear the data range before
a merge by using the Clear option.
GDXXRW test2.gdx par=V rng=sheet2!a1 clear
In the following example, we read data from a spreadsheet and save the data in a GDX file. Using the $GDXIN
and $LOAD GAMS directives, we read data from the GDX file into GAMS. The GAMS program modifies the data
and at the end of the run the data is saved in a new GDX file (tmp.gdx). The last step updates the spreadsheet
with the modified parameter.
The data in spreadsheet test1.xls:
A B C D E
1 a1 a2 a3
2 i1 1 2 3
3 i2 4 5 6
4
$CALL GDXXRW test1.xls Set=I rng=A2:A3 Rdim=1 Set=A rng=B1:D1 Cdim=1 Par=X
rng=A1:D3 Rdim=1 Cdim=1
$GDXIN test1.gdx
Set I(*),A(*);
$LOAD I A
Parameter X(I,A);
$LOAD X
Display I,A,X;
$GDXIN
X(I,A) = - X(I,A);
Execute_Unload tmp.gdx,I,A,X;
Execute GDXXRW.EXE tmp.gdx O=test1.xls par=X rng=EX6!A1:D3 rdim=1 cdim=1;
In this example we use a modified version of the trnsport.gms model from the GAMS model library. This example
illustrates:
Compilation phase
A B C D E
1 Distance values
2 new-york chicago topeka
3 seattle 2.5 1.7 1.8
4 san diego 2.5 1.8 1.4
5
6
7 plant capacities
8 seattle 350
9 san diego 600
10
11
12 demand at markets
13 new york 325
14 chicago 300
15 topeka 275
16
17
18 freight factor
19 90
20
Data for trnsport model
$gdxin trnsportdata.gdx
sets
i(*) canning plants
j(*) markets;
$load i j
display i,j;
Parameters
a(i) capacity of plant i in cases
b(j) demand at market j in cases
d(i,j) distance in thousands of miles
scalar f freight in dollars per case per thousand miles
$load d a b f
$gdxin
VARIABLES
x(i,j) shipment quantities in cases
z total transportation costs in thousands of dollars ;
POSITIVE VARIABLE x ;
EQUATIONS
cost define objective function
supply(i) observe supply limit at plant i
demand(j) satisfy demand at market j ;
execute_unload trnsportdata.gdx, x;
execute gdxxrw.exe trnsportdata.gdx var=x.l rng=sheet2!a1 ;
A B C D E
1 new-york chicago topeka
2 seattle 50 300
3 san diego 275 275
4
5
Solution
G.6 GDXDUMP
The program gdxdump will write scalars, sets and parameters (tables) to standard output formatted as a GAMS
program with data statements. To write to a file, use the output redirection provided by the operating system.
Usage
gdxdump filename {options}
Options
Symb = identifier
Selects a single identifier to be written.
UelTable = identifier
Write all unique elements found in the GDX file to a set using identifier as the name for the set.
Delim = [period, comma, tab, blank]
Selects a different delimiter to separate unique elements; period is the default.
NoHeader
Suppress the header information for a single symbol; only the data for the symbol will be written, not its
declaration.
Symbols
Generate an alphabetical list of all symbols in the GDX file.
G.7 GDXDIFF
The gdxdiff program compares the data of symbols with the same name, type and dimension in two GDX files
and writes the differences to a third GDX file. A summary report will be written to standard output.
Usage:
gdxdiff file1 file2 {diffile} {Eps = value} {RelEps = value} {Field = FieldName}
GDXDIFF requires two parameters, the file names of two GDX files. An optional third parameter is the name
of the GDX difference file. Without the third parameter, the difference file will be diffile.gdx in the current
directory.
Eps = value
A value used to test if two values are equal
RelEps = value
A value used to modify the test if two values are equal
The use of Eps and RelEps is best described by the code fragment below.
AbsDiff := Abs(V1 - V2);
Result := AbsDiff <= EpsAbsolute;
if not Result
then
begin
if EpsRelative = 0.0
then
Result := AbsDiff / (1.0 + Min(Abs(V1), Abs(V2))) <= EpsAbsolute
else
Result := AbsDiff / Max(Abs(V1), Abs(V2)) <= EpsRelative;
end;
Field = FieldName
A field can be specified to limit the comparison of variables and equations to that field. Field names follow
the GAMS conventions (Lo, L, Up, M, Scale and Prior)
G.8 GDXMERGE 241
Only symbols with the same name, type and dimension will be compared. Tuples with different values are written
to the GDX difference file, and a dimension is added to describe the difference using the following labels:
ins1 indicates that the tuple only occurs in the first file.
ins2 indicates that the tuple only occurs in the second file.
dif1 indicates that the tuple occurs in both files; contains the value from the first file.
dif2 indicates that the tuple occurs in both files; contains the value from the second file.
Example 16:
In the following example, the trnsport model is solved twice with different capacity data. GDX files are saved for
each run, and compared using gdxdiff. The shipments variable is loaded into a new variable used for a display
statement. We introduce four new unique elements that are used in the difference file.
G.8 GDXMERGE
Combines multiple GDX files into one file. Symbols with the same name, dimension and type are combined into
a single symbol of a higher dimension. The added dimension has the file name of the combined file as its unique
element.
Usage:
gdxmerge filepattern1 filepattern2...filepatternn [id=<ident1>,<ident2>...][big=<integer>]
Each file pattern represents a file name or a wildcard representation using ? and *. A parameter of the form
@filename, will process the commands from the text file specified.
The result of the merge will be written to a file called merged.gdx.
All symbols with matching type and dimension will be merged. By specifying the parameter id=ident1 the merge
process will only be performed for the identifier(s) specified.
By default, the program reads all gdx once and stores all data in memory before writing the merged.gdx file. The
big parameter is used to specify a cutoff for symbols that will be written one at a time. Each symbol that exceeds
the size will be processed by reading each gdx file and only process the data for that symbol. This can lead to
reading the same gdx file many times, but it allows the merging of large data sets. The formula used to calculate
the cutoff is: Dimension * TotalNumberOfElements. The number id doubled for variables and equations.
In addition to the symbols written, a set is added to the GDX file representing all the files processed during
the merge operation. The name of the set is Merged set 1, and is made unique by changing the number. The
explanatory text for each set element contains the date and time of the GDX file processed.
The file merged.gdx will never be used in a merge operation even if the name matches a file pattern.
Symbols with dimension 10 cannot be merged, because the resulting symbol will have dimension 11
which exceeds the maximum dimension allowed by GAMS.
242 GDX Facilities and Tools
Instead of using the display statement, we can also use the GAMSIDE to view the merged.gdx file. After selecting
the level field to be displayed and arranging the display:
seattle san-diego
new-york chicago new-york topeka
bdmlp 50 300 275 275
conopt 300 325 275
cplex 50 300 275 275
minos 50 300 275 275
snopt 50 300 275 275
xpress 300 325 275
G.9 GDXRANK
Reads one or more one dimensional parameters from a GDX file, sorts each parameter and writes the sorted
indices as a one dimensional parameter to the output GDX file.
Usage:
gdxrank inputfile outputfile
Each one dimensional parameter is read from the input file, sorted and the corresponding integer permutation
index written to the output file using the same name for the symbol. GAMS special values such as Eps, +Inf and
-Inf are recognized.
G.10 GDXCOPY
Copy GDX files to a different directory converting the file format used.
Usage:
gdxcopy option infile outdir
Output option:
Version 6 formatted files were introduced with version 22.3 of GAMS. Prior versions used version 5.
dir
1,219 t1.gdx
1,740 t0.gdx
889,973 i.gdx
1,740 pv.gdx
894,672 bytes
H.1 Introduction
When models are distributed to users other than the original developers or embedded in applications to be
deployed by other developers, issues of privacy, security, data integrity and ownership arise. We may have to
hide, protect or purge some parts of the model before it can be released. The information to be protected can be
of numeric or symbolic nature. For example:
Privacy
A Social Accounting Matrix supplied by a Statistical Office is required in a general equilibrium model to
be used by the Ministry of Finance. The data from the statistical office needs to be protected for obvious
privacy reasons and the model experiments are used to evaluate policy options that are highly
confidential. Most of the model structure is public, most of the data however is private and model results
need to be transformed in such a way as to prohibit the discovery of the original data.
Security
Components of a model contain proprietary information that describes mathematically a chemical
reaction. The associated algebra and some of the data are considered of strategic importance and need to
be hidden completely. The final model however, will be used at different locations around the world.
Integrity
Data integrity safeguards are needed to assure the proper functioning of a model. Certain data and
symbolic information needs to be protected from accidental changes that would compromise the operation
of the model.
To address these issues, access control at a symbol level and secure restart files have been added to the GAMS
system.
Access Control
The access to GAMS symbols like sets, variables, parameters and equations can be changed once with the
compile time commands $purge, $hide, $protect and $expose. $Purge will remove any information
associated with this symbol. $Hide will make the symbol and all its information invisible. $Protect
prevents changes to information. $Expose will revert the symbol to its original state.
Secure Restart Files
The GAMS licensing mechanism can be used to save a secure model in a secure work file. A secure work
file1 behaves like any other work file but is locked to a specific users license file. A privacy license, the
license file of the target users, is required to create a secure work file. The content of a secure work file is
disguised and protected against unauthorized access via the GAMS license mechanism.
A special license is required to set the access controls and to create a corresponding secure work file. Reporting
features have been added to allow audits and traces during generation and use of secure work files.
1 Work files are used to save and restart the state of a GAMS program. Depending on the context, we refer to those files as work
We continue to enter access control commands in a file called t2.gms and create a secure work file with the name
t2.g00:
$eolcom //
$protect all // make all symbol read only
$purge d f // remove items d and f
$hide cost supply a // make objective invisible
$expose transport b // allow changes to b
The access control commands are activated by the use of the privacy GAMS license option PLICENSE. This
option specifies the name of the target user license file. The save/restart file t2.g00 can only be read with the
target license file.
The three lines starting with *** are a recap of the content of the target license file. From now on, the source
and the target licensees are burned into this file and all its descendants. We are ready to send the restart file to
the target user or system.
The target user can now run the model with new data, add new GAMS statements, and make new save/restart
files. The only restrictions are that some of the symbols are hidden and that this model can only be executed
using the target license file. For example, the target user may want to half the demand and compare the original
solution with the new one. We will call this program t3.gms and it will be executed on the target system:
2 When using the GAMS/IDE interface, the GAMS parameters are entered and mainened in the text window just right to the Run
Note that the originator/owner of the secure work file is mentioned by name on the log file. A similar message is
contained in the listing file:
A more detailed inspection of the listing file will show that the hidden variables and equations do not appear in
the usual equation/variable listings and the solution print. The hidden items can only be accessed via a public
(exposed) model and a solve statement.
In the following two sections we will describe secure work files and the access control commands in more detail.
The keyword ALL applies the ACC to all identifiers defined up to this point in the GAMS source code. ACCs
can be changed and redefined within the same GAMS program. Identifiers inherited from a restart file cannot be
changed, however.
$include trnsport.gms
variable newc(i,j) new tansport data;
equation defnewc(i,j) definition of new transport data;
model getc compute new transport data / defnewc /;
defnewc(i,j).. newc(i,j) =e= f*d(i,j)/1000;
solve getc using cns;
Next, we change the objective function of the original model to a more complicated nonlinear function. Fur-
thermore, we will compute a base case value to be used later in the reporting model. Note the reference to
newc.l(i,j), since nexc is a variable we have to specify that we only want the level value:
We will save the above GAMS code under the name p1.gms, execute and make a save/restart file with the name
p1.g00 as follows:
> gams p1 s=p1
Now we are ready to make some test runs similar to those we expect to be defined by the target user. We will
define three scenarios to be solved in a loop and name the file u1.gms:
set s / one,two,three /;
parameter sbeta(s) / one 1.25, two 1.5, three 2.0 /
sf(s) / one 85, two 75, three 50 /;
H.6 Limitations and Future Requirements 249
When executing the above GAMS code together with the original transport model from the GAMS model library
we will get the following results.
Note that all symbols are still completely exposed. We need to add access controls to the model p1.gms before
we can ship it to the target client. The information to be protected is the original distance matrix and derived
information. We start out by hiding everything and then give access to selected parts of the model. We collect the
access control information in the file s1.gms shown below and save the secure work file under the name s1.g00.
Since we are still testing, we use our own license as target user. This will allows us to test the system the same
way the target user will use it:
$hide all
$expose getc newtrans rep
$expose i j z delta
$expose f beta a b
To test the new secure file, we run again the problem u1.gms. When doing so you will observe that equation,
variable and solution listings related to the hidden variables are not shown any more. Any attempt to reference
a hidden variable will case a compilation error.
> gams u1 r=s1
Before we can ship a secure work file we need a copy of the target user license file. We then will restart again
from p1.gms, zip the resulting secure files and we are ready to distribute the model:
implies that, in principle, certain information could be recovered, if one has knowledge of GAMS solvers internals
and is willing to expand considerable programming effort. In this section we will discuss current limitations and
possible extension to the security features.
The following limitations exist:
Solvers are not security aware and it would be possible to write a special GAMS solver that extracts
information about a specific model instance. Primal and duals values as well as first partial derivatives
could be extracted and displayed.
The names and explanatory text of all GAMS symbols are retained in the work file and could be
accessed by a special GAMS solver.
The source and target license files locked to the secure work file cannot be changed. If the target user
upgrades the GAMS system and receives a new license file, the secure work file cannot be read any
more.
I
I.1 Introduction
When models are distributed to users other than the original developers, issues of privacy, security, data integrity
and ownership arise. Besides using secure work files, one can compress and encrypt GAMS input files. The
compression and decompression of files is available to any GAMS user. The encryption follows the work file
security model and requires special licensing. Three new Dollar Control Options have been introduced:
$Encrypt <source> <target> Encrypts into a GAMS system file
$Compress <source> <target> Compresses into a GAMS system file
$Decompress <source> <target> Decompresses a GAMS system file
Encryption is only available if a system is licensed for secure work files and usually requires a target license file
which will contain the user or target encryption key. Once a file has been encrypted it cannot be decrypted any
more.
The use of a PLICENSE parameter will specify the target or privacy license to be used as a user key for encrypting.
Decompression and encrypting is done on the fly into memory when reading the GAMS system files. GAMS
will recognize if a file is just plain text or compressed and/or encrypted and will validate and process the files
accordingly.
Finally, all compressed and encrypted files are, of course, platform independent as any other GAMS input file.
...
--- Compress Source: C:\support\28Dec\trnsport.gms
--- Compress Target: C:\support\28Dec\t1.gms.gms
--- Compress Time : 0msec
...
Now we can treat the compressed input files like any other GAMS input file and the listing files will be identical
because the decompressed input is echoed just like any normal input line:
252 Compressed and Encrypted Input Files
Finally, we may want to encrypt the input file in a way to hide the equation definitions. To do this we just insert
$Offlisting and $Onlisting around the blocks of GAMS code we want to hide. We now can encrypt the modified
model file by using a privacy or target license file to lock the new encrypted file to this license key.
This new version of t1.gms can only be used with the target license file. To simulate the target environment we
can force the target license to be used instead of the installed one.
Note the use of dumopt=19 which is sometimes used for debugging and maintenance, writes a clean copy of the
input to the file t1.dmp (all include files and macros are expanded). In the file t1.dmp you will see that the input
text between $offlisting and $onlisting is suppressed.
An attempt to $Decompress and decrypt the file will fail as well. Once a file has been encrypted, it cannot be
decrypted any more. For example, trying to decompress as in the example before will fail:
> gams t4
...
--- Decompress Source: C:\support\28Dec\t1.gms
--- Decompress Target: C:\support\28Dec\t3.org
--- Decompress Error : Security violation
$ontext
This model demonstrates the use of compressed input files.
Remember, if the file names contain spaces you need
to use single or double quotes around the file names.
$offtext
$ontext
Input files can be encrypted and use the save/privacy license
file mechanism for managing the user password. Similar to
compression, we offer an $encrypt utility to lock any file to a
specific target license file. Once a file has been encrypted it
can only be read by a gams program that has the matching license
file. There is no inverse operation possible: you cannot recover
the original GAMS file from the encrypted version.
$eolcom //
$if NOT errorfree $abort pending errors
$decompress t1.gms t1.org // this has to fail
$if errorfree $abort decompress did not fail
$clearerror
$if errorlevel 1 $abort results for trnsport and t1 are not equal
J.1 Introduction
As systems with multiple CPUs and High Performance Computing Grids are becoming available more widely, the
GAMS language has been extended to take advantage of these new environments. New language features facilitate
the management of asynchronous submission and collection of model solution tasks in a platform independent
fashion. A simple architecture, relying on existing operating system functionality allows for rapid introduction of
new environments and provides for an open research architecture.
A typical application uses a coarse grain approach involving hundreds or thousands of model solutions tasks which
can be carried out in parallel. For example:
Scenario Analysis
Monte Carlo Simulations
Lagrangian Relaxation
Decomposition Algorithms
Advanced Solution Approaches
The grid features work on all GAMS platforms and have been tailored to many different environments, like the
Condor Resource Manager, a system for high throughput computing from the University of Wisconsin or the Sun
Grid Engine. Researchers using Condor reported a delivery of 5000 CPU hours in 20 hours wall clock time.
Disclaimer. The use of the term grid computing may be offensive to some purists in the computer science world.
We use it very loosely to refer to a collection of computing components that allow us to provide high throughput
to certain applications. One may also think of it as a resurrection of the commercial service bureau concept of
some 30 years ago.
Caution. Although these features have been tested on all platforms and are part of our standard release we may
change the approach and introduce alternative mechanisms in the future.
Acknowledgments. Prof. Monique Guignard-Spielberg from the Wharton School at U Penn introduced us
to parallel Lagrangian Relaxation on the SUN Grid Environment. Prof. Michael Ferris from the University of
Wisconsin at Madison adopted our original GAMS grid approach to the high throughput system Condor and
helped to make this approach a practical proposition.
1. Generation. The symbolic equations of the model are used to instantiate the model using the current state
of the GAMS data base. This instance contains all information and services needed by a solution method to
attempt a solution. This representation is independent of the solution subsystem and computing platform.
2. Solution. The model instance is handed over to a solution subsystem and GAMS will wait until the solver
subsystem terminates.
3. Update. The detailed solution and statistics are used to update the GAMS data base.
In most cases, the time taken to generate the model and update the data base with the solution will be much
smaller than the actual time spent in a specific solution subsystem. The model generation takes a few seconds,
whereas the time to obtain an optimal solution may take a few minutes to several hours. If sequential model
solutions do not depend on each other, we can solve in parallel and update the data base in random order. All we
need is a facility to generate models, submit them for solution and continue. At a convenient point in our GAMS
program we will then look for the completed solution and update the data base accordingly. We will term this
first phase the submission loop and the subsequent phase the collection loop:
Submission Loop. In this phase we will generate and submit models for solutions that can be solved indepen-
dently.
Collection Loop. The solutions of the previously submitted models are collected as soon a solution is available.
It may be necessary to wait for some solutions to complete by putting the GAMS program to sleep.
Note that we have assumed that there will be no errors in any of those steps. This, of course, will not always be
the case and elaborate mechanisms are in place to make the operation fail-safe.
Loop(p(pp),
ret.fx = rmin + (rmax-rmin)/(card(pp)+1)*ord(pp) ;
Solve minvar min var using miqcp ;
xres(i,p) = x.l(i);
report(p,i,inc) = xi.l(i);
report(p,i,dec) = xd.l(i) );
This loop will save the solutions to the model MINVAR for different returns RET. Since the solutions do not
depend on the order in which they are carried out, we can rewrite this loop to operate in parallel. The first step
is to write the submit loop:
The model attribute .solvelink controls the behavior of the solve statement. A value of 3 tells GAMS to
generate and submit the model for solution and continue without waiting for the completion of the solution step.
There is a new model attribute .handle which provides a unique identification of the submitted solution request.
We need to store those handle values, in this case in the parameter h, to be used later to collect the solutions
once completed. This is then done with a collection loop:
J.3 A First Example 257
loop(pp$handlecollect(h(pp)),
xres(i,pp) = x.l(i);
report(pp,i,inc) = xi.l(i);
report(pp,i,dec) = xd.l(i) );
The function handlecollect interrogates the solution process. If the solution process has been completed the
results will be retrieved and the function returns a value of 1. If the solution is not ready to be retrieved the
value zero will be returned.
The above collection loop has one big flaw. If a solution was not ready it will not be retrieved. We need to
call this loop several times until all solutions have been retrieved or we get tired of it and quit. We will use a
repeat-until construct and the handle parameter h to control the loop to look only for the not yet loaded solutions
as shown below:
Repeat
loop(pp$handlecollect(h(pp)),
xres(i,pp) = x.l(i);
report(pp,i,inc) = xi.l(i);
report(pp,i,dec) = xd.l(i);
display$handledelete(h(pp)) trouble deleting handles ;
h(pp) = 0 ) ;
display$sleep(card(h)*0.2) sleep some time;
until card(h) = 0 or timeelapsed > 100;
xres(i,pp)$h(pp) = na;
Once we have extracted a solution we will set the handle parameter to zero. In addition, we want to remove the
instance from the system by calling the function handledelete which returns zero if successful (see definition).
No harm is done if it fails but we want to be notified via the conditional display statement. Before running the
collection loop again, we may want to wait a while to give the system time to complete more solution steps. This
is done with the sleep command that sleeps some time. The final wrinkle is to terminate after 100 seconds elapsed
time, even if we did not get all solutions. This is important, because if one of the solution steps fails our program
would never terminate. The last statement sets the results of the missed solves to NA to signal the failed solve.
The parameter h will now contain the handles of the failed solvers for later analysis.
Alternatively, we could have used the function handlestatus and collect the solution which is stored in a GDX
file. For example we could write:
loop(pp$(handlestatus(h(pp))=2),
minvar.handle = h(pp);
execute_loadhandle minvar;
xres(i,pp) = x.l(i);
report(pp,i,inc) = xi.l(i);
report(pp,i,dec) = xd.l(i) );
The function handlestatus interrogates the solution process and returns 2 if the solution process has been
completed and the results can be retrieved. The solution is stored in a GDX file which can be loaded in a
way similar to other gdx solution points. First we need to tell GAMS what solution to retrieve by setting
the minvar.handle to the appropriate value. Then we can use execute loadhandle to load the solution for
model minvar back into the GAMS data base. Using handlestatus and loadhandle instead of the simpler
handlecollect adds one more layer of control to the final collection loop. We now need one additional if statement
inside the above collection loop:
Repeat
loop(pp$h(pp),
if(handlestatus(h(pp))=2,
minvar.handle = h(pp);
execute_loadhandle minvar;
xres(i,pp) = x.l(i);
report(pp,i,inc) = xi.l(i);
report(pp,i,dec) = xd.l(i);
display$handledelete(h(pp)) trouble deleting handles ;
h(pp) = 0 ) ) ;
258 The GAMS Grid Computing Facility
Now we are ready to run the modified model. The execution log will contain some new information that may be
useful on more advanced applications:
--- LOOPS pp = p1
--- 46 rows 37 columns 119 non-zeroes
--- 311 nl-code 7 nl-non-zeroes
--- 14 discrete-columns
--- Submitting model minvar with handle grid137000002
--- Executing after solve
...
--- GDXin=C:\answerv5\gams_srcdev\225j\grid137000003\gmsgrid.gdx
--- Removing handle grid137000003
The log will now contain some additional information about the submission, retrieval and removal of the solution
instance. In the following sections we will make use of this additional information. You can find a complete
example of a grid enabled transport model in the GAMS model library.
At a final note, we have made no assumptions about what kind of solvers and what kind of computing environment
we will operate. The above example is completely platform and solver independent and it runs on your Windows
laptop or on a massive grid network like the Condor system without any changes in the GAMS source code.
status(pp,modelstat) = minvar.modelstat;
status(pp,seconds) = minvar.resusd;
else
handle := Waiting );
status(pp,status) = handle );
display status;
To run the above program we will restart from the previous save file by using the restart or r parameter.
> gams qcheck r=submit gdir=c:\test\grid
The output may then look like:
You may want to do some more detailed analysis on one of the solved model instances. Then we may have
a qanalyze.gms program that may look like and be called using the double dash option, which sets a GAMS
environment variable:
loop(pp$handlestatus(h(pp)),
minvar.handle = h(pp);
execute_loadhandle minvar;
xres(i,pp) = x.l(i);
report(pp,i,inc) = xi.l(i);
report(pp,i,dec) = xd.l(i);
display$handledelete(h(pp)) trouble deleting handles ;
h(pp) = 0 );
xres(i,pp)$h(pp) = na;
. . .
We would restart the above program from the save file that was created by the submitting job like:
> gams qreport r=submit gdir=c:\test\grid
Note that it would not be necessary to run the job out of the same directory we did the initial submission. We
dont even have to run the same operating system.
HandleStatus(handle) returns the status of the solve identified by handle. An execution error is triggered if
GAMS cannot retrieve the status of the handle.
0 the model instance is not known to the system
1 the model instance exists but no solution process is incomplete
2 the solution process has terminated and the solution is ready for retrieval
3 the solution process signaled completion but the solution cannot be retrieved
HandleDelete(handle) returns the status of the deletion of the handle model instance. In case of a nonzero
return an execution error is triggered.
0 the model instance has been removed
1 the argument is not a legal handle
2 the model instance is not known to the system
3 the deletion of the model instance encountered errors
HandleSubmit(handle) resubmits a previously created instance for solution. In case of a nonzero return an
execution error is triggered.
0 the model instance has been resubmitted for solution
1 the argument is not a legal handle
2 the model instance is not known to the system
3 the completion signal could not be removed
4 the resubmit procedure could not be found
5 the resubmit process could not be started
In addition, GAMS will issue execution errors which will give additional information that may help to identify
the source of problems. The property execerror can be used to get and set the number of execution errors.
1. Create a subdirectory in the GridDir with the name gridnnn. Where nnn stands for the numeric value of
the handle. The handle value is the internal symbol ID number x 1e6 + the model instance number. For
example, in the QMEANVAR example the first grid subdirectory was grid137000002.
2. Remove the completion signal in case the file already exists. Currently the signal is a file called finished.
For example, grid137000002/finished.
3. Create or replace a gdx file called gmsgrid.gdx which will contain a dummy solution with failed model
and solver status. This file will be overwritten by the final step of the solution process and will be read
when calling execute loadhandle.
262 The GAMS Grid Computing Facility
4. Place all standard GAMS solver interface files into the above instance directory.
5. Execute the submission wrapper called gmsgrid.cmd under Windows or gmsgrid.run under Unix. These
submission scripts are usually located in the GAMS system directory but are found via the current path if
not found in the GAMS system directory.
The grid submission script gmsgrid.cmd or gmsgrid.run is called with three arguments needed to make a
standard GAMS solver call:
The submission script then does the final submission to the operating system. This final script will perform the
following steps:
If we want to use the function handlesubmit() we also have to create the gmsrerun.cmd or gmsrerun.run
script which could later be used to resubmit the job.
For example, the default submission script for Windows is shown below:
@echo off
: gams grid submission script
:
: arg1 solver executable
: 2 control file
: 3 scratch directory
:
: gmscr_nx.exe processes the solution and produces gmsgrid.gdx
:
: note: %3 will be the short name because START cannot
: handle spaces and/or "...". We could use the original
: and use %~s3 which will strip ".." and makes the name :
: short
: gmsrerun.cmd will resubmit runit.cmd
echo @echo off > %3runit.cmd
echo %1 %2 >> %3runit.cmd
echo gmscr_nx.exe %2 >> %3runit.cmd
echo echo OK ^> %3finished >> %3runit.cmd
echo exit >> %3runit.cmd
echo @start /b %3runit.cmd ^> nul > %3gmsrerun.cmd
start /b %3runit.cmd > nul
exit
2. Copy the license file: If no license file is found in the GAMS directory, you will be prompted for one during
the installation. If you are not sure if you have a license file, or do not have it yet, choose No when asked if
you wish to copy a license file. You can always do this later. If no valid license file is found, GAMS will still
function in demonstration mode and will only solve small problems. For example, all demonstration and
student systems are shipped without a license. If you have a license file that you wish to copy to the GAMS
directory at this time, answer Yes to the license file prompt. You will now be given the opportunity to
browse the file system and find the license file you wish to copy. For example, if your system came with a
license file on a diskette, browse to the A: drive and select gamslice.txt. When you have found the correct
file, choose open to perform the copy.
3. Create project files: If this is the first installation of GAMS with the GAMS IDE on your system, setup
will create a default GAMS project in a subdirectory of the WINDOWS directory. If this is not the first
such installation, your existing GAMS projects will be preserved.
4. Choose default solvers: Run the GAMS IDE by double clicking gamside.exe from the GAMS directory.
To view or edit the default solvers, choose File Options Solvers from the IDE. You can accept the
existing defaults if you wish, but most users will want to pick default solvers for each model type. It is a
good idea to review the solver defaults when installing a new GAMS system or when updating a license file.
5. Run a few models to test the GAMS system: The on-line help for the IDE (Help Help Topics
Guided Tour) describes how to copy a model from the GAMS model library, run it, and view the solution.
To test your installation, run the following models from the GAMS model library:
If there are any problems during these test runs, read the section on Troubleshooting.
266 Installation and System Notes
1. Run the program gamsinst: gamsinst is a command line program used to install and configure GAMS.
It prompts the user for default solvers to be used for each model type. If possible choose solvers you have
licensed since unlicensed solvers will only run in demonstration mode.These solver defaults can be changed
by:
(a) rerunning gamsinst and resetting the default values
(b) setting a command line default, e.g. gams trnsport lp=bdmlp
(c) by an option statement in the GAMS model, e.g: option lp=bdmlp
If you have to support different operating systems from the same installation, please use gamsinst -sys
all A complete log of the installation is stored in gamsinst.log.
N.B. The system wide solver defaults are shared by the command line and GUI versions of GAMS, so you
can also choose these defaults using GAMS IDE.
2. Add the GAMS directory to your path. To avoid having to type in an absolute path name each time you
run GAMS, we recommend adding the GAMS directory to your PATH when using the console mode (not
the GAMS IDE) version of GAMS. In case more than one GAMS system is installed on the machine, then
separate paths have to be set before invoking each version. Under Windows 95/98, this change requires
editing the autoexec.bat file and adding the GAMS directory to the path, as the following example
illustrates:
path=your\current\path\setting
path=%PATH%;c:\gams
Under Windows 2000/XP the following procedure must be applied to add the GAMS directory to your
path,
Open the System Properties under the Control Panel.
On the Advanced tab click on the Environment Variables button and select the existing variable PATH.
Click edit.
In the Value box, add the GAMS directory to the path, as the following example illustrates:
c:\your\current\path\setting;C:\gams
Click OK.
Network Support
GAMS 2.50 supports the Uniform Naming Convention (UNC) both for the system and for the solvers. GAMS
can be installed and accessed over a network of PCs. For example: GAMS has been installed on a machine
called bach in the directory gams. Now you are sitting in front of another machine and want to run a model
(trnsport.gms), which is located in the current directory. To run this model, you can use the command
>\\bach\gams\gams trnsport
or, if the path to the remote machines GAMS system is in your PATH statement, simply
>gams trnsport
Troubleshooting
GAMS rarely causes problems, and most of the problems that occur can be corrected quite easily by following
the notes in this section.
267
GAMS itself is running out of memory. If the memory installed on your machine is not sufficient for your problems,
try to run the model on a machine with more memory. Note that using virtual memory will allow you to solve
larger problems than possible using real memory alone. However, using virtual memory will slow down the GAMS
system.
Environment space When running the console mode version of GAMS, it may under some circumstances
be necessary to adjust the settings for the initial environment to 4096. In order to change this under Windows
95/98, click the right mouse button on the title bar of MS DOS Prompt Window and select Properties. Open
the Memory tab and change the value of the Initial environment box to 4096.
Other technical problems: For other technical questions like insufficient disk space, licensing problems or
solver failures please check the GAMS web sites or contact our support staff directly.
1. Choose a location for the GAMS system directory (the directory where the GAMS system files should reside).
We recommend that you choose a name that indicates what distribution of GAMS you are installing, and
that anticipates the installation of updated GAMS systems in the future. For example, if you are installing
the 21.0 distribution (from the 21.0 CD, say), a good choice for the GAMS system directory would be
/usr/gams/21.0. The total GAMS system uses 40-100 MB. The command df will tell you how much disk
space is available on each file system. If the directory where you want to install GAMS is not below your
home directory, you may need to have super user privileges on the machine.
2. Create the GAMS system directory, for instance /usr/gams/21.0. Go to this directory. Make sure pwd
returns the name of this directory correctly.
3. Transfer the distribution files into the GAMS system directory. These files are available, whether from the
GAMS CD or via ftp, in two forms: in one large self extracting zip archive with a sfx.exe file extension,
or as a zip archive split into diskette sized pieces, along with some routines to put the pieces together. If you
use the self extracting archive (e.g. lnxgams sfx.exe on a Linux system), you can run the archive directly
from the CD while in the GAMS system directory to extract the necessary files to the system directory. If
you use the split archive, copy or download the files into the GAMS system directory. For example, you
might execute the following commands:
mkdir /usr/gams/21.0
cd /usr/gams/21.0
/dev/cdrom/systems/lnx/lnx2.2/lnxgams_sfx.exe
4. To mount the GAMS CD, you will usually need to be logged in as root. We assume you want to mount the
CD over the directory /cdrom. If the directory you want to mount over does not exist, you must create it
now. Once this directory is created, mount the CD, using the appropriate command. The correct arguments
for the mount command vary from machine to machine. Some examples are given below:
After mounting the CD, view the README.TXT file on it to find the subdirectory containing the GAMS
system for your machine.
5. If you transferred the split files via ftp, check that you have execute permission set on the files gamsinst,
gamsunpack, and gmsunzip. If you are not sure how to do this, just type in the command
6. Check if the file gamslice.txt exists in the GAMS system directory. If you have a licensed (professional
or evaluation) version of GAMS, a license file is required to solve large models. Copy the license file from
the diskette in MS-DOS format to the GAMS system directory. Some license files are sent via email. As a
last resort, the license may have been sent on paper and must be typed in by hand.
If no license file is present, GAMS will still function in demonstration mode and can only solve small
problems. For example, student and demonstration systems are sent without a license file. A license file
can easily be added later and ./gamsinst rerun, so if you cannot find a license file, you can safely proceed
without one.
7. Run the program ./gamsinst. This will unpack files if necessary. It will also prompt you for default solvers
to be used for each class of models. If possible, choose solvers you have licensed since unlicensed solvers will
only run in demonstration mode. These solver defaults can be changed or overridden by:
(a) rerunning ./gamsinst and resetting the default values
(b) setting a command line default, e.g. gams trnsport lp=bdmlp
(c) an option statement in the GAMS model, e.g: option lp=bdmlp
8. Add the GAMS system directory to your path (see ACCESS TO GAMS below).
9. To test the installation, login as a normal user and run a few models from a scratch directory, not the GAMS
system directory.
10. If you move the GAMS system to another directory, remember to rerun ./gamsinst. It is also good practice
to rerun ./gamsinst when you add or change your license file if this has changed the set of solvers you are
licensed to run.
Access to GAMS
To run GAMS you must be able to execute the GAMS programs located in the GAMS system directory. There are
several ways to do this. Remember that the GAMS system directory in the examples below may not correspond
to the directory where you have installed your GAMS system.
1. Those using the C shell (csh) and its variants can modify their .cshrc file by adding the second of the two
lines given below:
Those using the Bourne (sh) or Korn (ksh) shells and their variants can modify their .profile file by
adding the second of the three lines below:
PATH=/your/previous/path/setting
PATH=$PATH:/usr/gams/21.0 # new
export PATH
You should logout and login again after you make any changes to your path, unless you make the same
changes locally in a shell.
2. You may prefer to alias the names of the programs instead of adding to the path as described above. C
shell users can use the following commands on the command line or in their .cshrc file:
The correct Bourne or Korn shell syntax (either command line or .profile) is:
alias gams=/usr/gams/21.0/gams
alias gamslib=/usr/gams/21.0/gamslib
alias gamsbatch=/usr/gams/21.0/gamsbatch
Again, you should logout and login again in order for the alias settings in .cshrc or .profile to take effect.
3. Casual users can always type the absolute path names of the GAMS programs. For example,
/usr/gams/21.0/gams trnsport
Technical problems: For any technical questions like insufficient disk space, licensing problems or solver
failures please check the GAMS web sites or contact our support staff directly.
270 Installation and System Notes
Index
text, 31 infeasible, 88
GAMSLIB, 163 solution, 88
GDX facilities and tools, 225 variable, 61, 72
GDXCOPY, 229, 243 intermediate
GDXDIFF, 228, 240 infeasible, 88
GDXDUMP, 228, 239 noninteger, 88
GDXMERGE, 229, 241 nonoptimal, 88
GDXRANK, 229, 242 intersection, set operation, 111
GDXXRW, 228, 229 iteration
ge, relational operator, 30, 97 default limit, 88
goto interrupt, 89
dollar control option, 199 iterlim
grid computing, 255 GAMS option, 217
gt, relational operator, 30 model attribute, 73
GTM, example from GAMSLIB, 104, 105 option, 88, 89
iterusd, model attribute, 73
hidden
dollar control option, 199 keep
holdfixed, model attribute, 73 GAMS call parameter, 180
kill
ide dollar control option, 202
GAMS call parameter, 179 KORPET, example from GAMSLIB, 46
identifiers, 159
if label, 29, 31, 40
dollar control option, 199 dollar control option, 202
if-elseif-else order, 125
a statement, 144 order on displays, 122
example, 145 quoted, 31, 36
syntax, 144 row and column, 47
impl-asn, reference type, 82 unquoted, 31
include using in equations, 67
dollar control option, 201 le, relational operator, 30, 97
index position(s), 159 leftmargin
indices, controlling, 54 GAMS call parameter, 180
INDUS, example from GAMSLIB, 48 legal characters, 30
inequality constraint, 159 level, 62
inf libincdir
an extended range value, 95 GAMS call parameter, 180
as a variable bound, 70, 122 libinclude
extended range value, 32, 56 dollar control option, 202
variable bound, 60 license
infeasible, 86, 88, 92, 159 GAMS call parameter, 180
infes, solution marker, 92 limcol
initial values, 28, 62 GAMS call parameter, 180
initialization, 82, 159 GAMS option, 217
of data, 43 model attribute, 73
of parameters, 44 limrow
inlinecom GAMS call parameter, 181
dollar control option, 201 GAMS option, 217
input model attribute, 73
GAMS call parameter, 179 lines
inputdir dollar control option, 203
GAMS call parameter, 179 list, 159
inputdir1 to inputdir18 of labels using Asterisks, 37
GAMS call parameter, 179 list format, 159
integer locally
276 INDEX
variable
binary, 60
free, 60
integer, 60
negative, 60
positive, 60
statement, 59, 60
styles for declaration, 60
suffix, 61
syntax of declaration, 59
types, 60
variable attributes
activity level (.l), 61
branching priority value (.prior), 61
fixed value (.fx), 61
lower bound (.lo), 61
marginal or dual value (.m), 61
scale value (.scale), 61
upper bound (.up), 61
variable bounds
activity level, 62
fixing, 61
variable type, 161
vector, 161
while
example, 146
statement, 146
syntax, 146
Windows Installation Notes, 265
workdir
GAMS call parameter, 190
workspace, model attribute, 73