doc/quicktour.dox
author zsuzska
Wed, 20 Jul 2005 22:36:37 +0000
changeset 1578 1d3a1bcbc874
parent 1541 305ce06287c9
child 1580 a9e4208cf4e3
permissions -rw-r--r--
kruskal_demo corrected, quicktour filled with kruskal
athos@1169
     1
/**
athos@1169
     2
alpar@1170
     3
\page quicktour Quick Tour to LEMON
alpar@1170
     4
athos@1175
     5
Let us first answer the question <b>"What do I want to use LEMON for?"
athos@1175
     6
</b>. 
athos@1175
     7
LEMON is a C++ library, so you can use it if you want to write C++ 
athos@1175
     8
programs. What kind of tasks does the library LEMON help to solve? 
athos@1175
     9
It helps to write programs that solve optimization problems that arise
athos@1175
    10
frequently when <b>designing and testing certain networks</b>, for example
athos@1175
    11
in telecommunication, computer networks, and other areas that I cannot
athos@1175
    12
think of now. A very natural way of modelling these networks is by means
athos@1183
    13
of a <b> graph</b> (we will always mean a directed graph by that and say
athos@1183
    14
<b> undirected graph </b> otherwise). 
athos@1175
    15
So if you want to write a program that works with 
athos@1183
    16
graphs then you might find it useful to use our library LEMON. LEMON 
athos@1183
    17
defines various graph concepts depending on what you want to do with the 
athos@1183
    18
graph: a very good description can be found in the page
athos@1183
    19
about \ref graphs "graphs".
athos@1175
    20
athos@1514
    21
You will also want to assign data to the edges or nodes of the graph, for
athos@1514
    22
example a length or capacity function defined on the edges. You can do this in
athos@1514
    23
LEMON using so called \b maps. You can define a map on the nodes or on the edges of the graph and the value of the map (the range of the function) can be practically almost of any type. Read more about maps \ref maps-page "here".
athos@1175
    24
athos@1528
    25
In this quick tour we want to show you some facilities LEMON library can provide through examples (simple demo programs). The examples will only show part of the functionality, but links will always be given to reach complete details. 
athos@1528
    26
You will find links next to the code fragments that help to download full demo programs: save them on your computer and compile them according to the description in the page about \ref getstart "How to start using LEMON". 
athos@1528
    27
athos@1528
    28
Have fun!
athos@1175
    29
athos@1522
    30
<ul> <li> The first thing to discuss is the way one can create data structures
athos@1522
    31
like graphs and maps in a program using LEMON. 
athos@1522
    32
//There are more graph types
athos@1522
    33
//implemented in LEMON and you can implement your own graph type just as well:
athos@1522
    34
//read more about this in the already mentioned page on \ref graphs "graphs".
athos@1522
    35
athos@1522
    36
First we show how to add nodes and edges to a graph manually. We will also
athos@1522
    37
define a map on the edges of the graph. After this we show the way one can
athos@1522
    38
read a graph (and perhaps maps on it) from a stream (e.g. a file). Of course
athos@1522
    39
we also have routines that write a graph (and perhaps maps) to a stream
athos@1522
    40
(file): this will also be shown. LEMON supports the DIMACS file formats to
athos@1534
    41
read network optimization problems, but more importantly we also have our own
athos@1522
    42
file format that gives a more flexible way to store data related to network
athos@1522
    43
optimization.
athos@1522
    44
athos@1530
    45
<ol> <li>The following code shows how to build a graph from scratch
athos@1530
    46
and iterate on its nodes and edges.  This example also shows how to
athos@1530
    47
give a map on the edges of the graph.  The type Listgraph is one of
athos@1530
    48
the LEMON graph types: the typedefs in the beginning are for
athos@1530
    49
convenience and we will assume them later as well.
athos@1522
    50
athos@1530
    51
\include hello_lemon.cc
athos@1522
    52
athos@1530
    53
See the whole program in file \ref hello_lemon.cc in the \c demo subdir of
athos@1526
    54
LEMON package.
athos@1175
    55
athos@1526
    56
    If you want to read more on the LEMON graph structures and
athos@1526
    57
concepts, read the page about \ref graphs "graphs".
athos@1522
    58
athos@1530
    59
athos@1530
    60
<li>LEMON has an own file format for storing graphs, maps on edges/nodes and some other things. Instead of any explanation let us give a
athos@1530
    61
short example file in this format: read the detailed description of the LEMON
athos@1530
    62
graph file format and input-output routines here: \ref graph-io-page.
athos@1530
    63
athos@1530
    64
So here is a file describing a graph of 6 nodes (0 to 5), two nodemaps
athos@1530
    65
(called \c coordinates_x and \c coordinates_y), several edges, an edge map
athos@1530
    66
called \c capacity and two designated nodes (called \c source and \c target).
athos@1530
    67
athos@1541
    68
\verbatim
athos@1541
    69
@nodeset
athos@1541
    70
id      coordinates_x   coordinates_y
athos@1541
    71
5       796.398 208.035
athos@1541
    72
4       573.002 63.002
athos@1541
    73
3       568.549 401.748
athos@1541
    74
2       277.889 68.476
athos@1541
    75
1       288.248 397.327
athos@1541
    76
0       102.239 257.532
athos@1541
    77
@edgeset
athos@1541
    78
                id      capacity
athos@1541
    79
4       5       6       8
athos@1541
    80
3       5       5       8
athos@1541
    81
2       4       4       5
athos@1541
    82
1       4       3       8
athos@1541
    83
1       3       2       5
athos@1541
    84
0       2       1       10
athos@1541
    85
0       1       0       10
athos@1541
    86
#This is a comment here
athos@1541
    87
@nodes
athos@1541
    88
source 0
athos@1541
    89
target 5
athos@1541
    90
@edges 
athos@1541
    91
@attributes 
athos@1541
    92
author "Attila BERNATH"
athos@1541
    93
@end
athos@1541
    94
\endverbatim
athos@1530
    95
athos@1530
    96
Finally let us give a simple example that reads a graph from a file and writes
athos@1530
    97
it to the standard output.
athos@1530
    98
athos@1530
    99
\include reader_writer_demo.cc
athos@1530
   100
athos@1530
   101
See the whole program in file \ref reader_writer_demo.cc.
athos@1530
   102
athos@1526
   103
<li> The following code shows how to read a graph from a stream
athos@1526
   104
(e.g. a file) in the DIMACS file format (find the documentation of the
athos@1526
   105
DIMACS file formats on the web).
athos@1522
   106
athos@1181
   107
\code
athos@1181
   108
Graph g;
athos@1181
   109
std::ifstream f("graph.dim");
athos@1181
   110
readDimacs(f, g);
athos@1181
   111
\endcode
athos@1522
   112
athos@1526
   113
One can also store network (graph+capacity on the edges) instances and
athos@1526
   114
other things (minimum cost flow instances etc.) in DIMACS format and
athos@1534
   115
read these in LEMON: to see the details read the documentation of the
athos@1534
   116
\ref dimacs.h "Dimacs file format reader". 
athos@1522
   117
athos@1514
   118
</ol>
athos@1514
   119
<li> If you want to solve some transportation problems in a network then 
athos@1175
   120
you will want to find shortest paths between nodes of a graph. This is 
athos@1175
   121
usually solved using Dijkstra's algorithm. A utility
athos@1175
   122
that solves this is  the \ref lemon::Dijkstra "LEMON Dijkstra class".
athos@1522
   123
The following code is a simple program using the 
athos@1530
   124
\ref lemon::Dijkstra "LEMON Dijkstra class": it calculates the shortest path between node \c s and \c t in a graph \c g.
athos@1530
   125
We omit the part reading the graph  \c g and the length map \c len.
athos@1175
   126
athos@1528
   127
\dontinclude dijkstra_demo.cc
athos@1528
   128
\skip ListGraph
athos@1530
   129
\until Graph g
athos@1530
   130
...
athos@1530
   131
\skip Dijkstra algorithm
athos@1528
   132
\until std::cout << g.id(s)
athos@1175
   133
alpar@1287
   134
See the whole program in \ref dijkstra_demo.cc.
athos@1183
   135
athos@1530
   136
Some explanation: after instantiating a member of the Dijkstra class
athos@1530
   137
we run the Dijkstra algorithm from node \c s. After this we read some
athos@1530
   138
of the results.  You can do much more with the Dijkstra class, for
athos@1530
   139
example you can run it step by step and gain full control of the
athos@1530
   140
execution. For a detailed description, see the documentation of the
athos@1530
   141
\ref lemon::Dijkstra "LEMON Dijkstra class".
athos@1183
   142
athos@1183
   143
zsuzska@1578
   144
<li> If you want to design a network and want to minimize the total
zsuzska@1578
   145
length of wires then you might be looking for a <b>minimum spanning
zsuzska@1578
   146
tree</b> in an undirected graph. This can be found using the Kruskal
zsuzska@1578
   147
algorithm: the function \ref lemon::kruskal "LEMON Kruskal " does
zsuzska@1578
   148
this job for you.  After we had a graph \c g and a cost map \c
zsuzska@1578
   149
edge_cost_map , the following code fragment shows an example how to get weight of the minmum spanning tree, if the costs are uniform:
athos@1175
   150
zsuzska@1578
   151
\dontinclude kruskal_demo.cc
zsuzska@1578
   152
\skip std::cout 
zsuzska@1578
   153
\until kruskal
zsuzska@1578
   154
zsuzska@1578
   155
It gives back a edge bool map, which contains the edges of the tree.
zsuzska@1578
   156
If the costs are non-uniform, for example  the cost is given by \c
zsuzska@1578
   157
edge_cost_map_2 , or the edges of the tree are have to be given in a
zsuzska@1578
   158
vector, then we can give to the kruskal a vector \c tree_edge_vec , instead of
zsuzska@1578
   159
an edge bool map:
zsuzska@1578
   160
zsuzska@1578
   161
\skip edge_cost_map_2 
zsuzska@1578
   162
\until edge_cost_map_2, std::back_inserter
zsuzska@1578
   163
zsuzska@1578
   164
And finally the next fragment shows how to use the functions \c makeKruskalMapInput and \c makeKruskalSequenceOutPut:
zsuzska@1578
   165
zsuzska@1578
   166
\skip makeKruskalSequenceOutput
zsuzska@1578
   167
\until tree_edge_vec
zsuzska@1578
   168
zsuzska@1578
   169
See the whole program in \ref kruskal_demo.cc.
zsuzska@1578
   170
zsuzska@1578
   171
athos@1511
   172
athos@1517
   173
<li>Many problems in network optimization can be formalized by means
athos@1517
   174
of a linear programming problem (LP problem, for short). In our
athos@1517
   175
library we decided not to write an LP solver, since such packages are
athos@1517
   176
available in the commercial world just as well as in the open source
athos@1517
   177
world, and it is also a difficult task to compete these. Instead we
athos@1517
   178
decided to develop an interface that makes it easier to use these
athos@1517
   179
solvers together with LEMON. The advantage of this approach is
athos@1517
   180
twofold. Firstly our C++ interface is more comfortable than the
athos@1517
   181
solvers' native interface. Secondly, changing the underlying solver in
athos@1517
   182
a certain software using LEMON's LP interface needs zero effort. So,
athos@1517
   183
for example, one may try his idea using a free solver, demonstrate its
athos@1517
   184
usability for a customer and if it works well, but the performance
athos@1517
   185
should be improved, then one may decide to purchase and use a better
athos@1517
   186
commercial solver.
athos@1517
   187
athos@1517
   188
So far we have an
athos@1526
   189
interface for the commercial LP solver software \b CPLEX (developed by ILOG)
athos@1514
   190
and for the open source solver \b GLPK (a shorthand for Gnu Linear Programming
athos@1517
   191
Toolkit).
athos@1514
   192
athos@1514
   193
We will show two examples, the first one shows how simple it is to formalize
athos@1514
   194
and solve an LP problem in LEMON, while the second one shows how LEMON
athos@1514
   195
facilitates solving network optimization problems using LP solvers.
athos@1514
   196
athos@1514
   197
<ol>
athos@1514
   198
<li>The following code shows how to solve an LP problem using the LEMON lp
athos@1517
   199
interface. The code together with the comments is self-explanatory.
athos@1511
   200
athos@1530
   201
\dontinclude lp_demo.cc
athos@1530
   202
\skip A default solver is taken
athos@1530
   203
\until End of LEMON style code
athos@1175
   204
athos@1514
   205
See the whole code in \ref lp_demo.cc.
athos@1514
   206
athos@1517
   207
<li>The second example shows how easy it is to formalize a max-flow
athos@1517
   208
problem as an LP problem using the LEMON LP interface: we are looking
athos@1517
   209
for a real valued function defined on the edges of the digraph
athos@1517
   210
satisfying the nonnegativity-, the capacity constraints and the
athos@1517
   211
flow-conservation constraints and giving the largest flow value
athos@1517
   212
between to designated nodes.
athos@1517
   213
athos@1517
   214
In the following code we suppose that we already have the graph \c g,
athos@1517
   215
the capacity map \c cap, the source node \c s and the target node \c t
athos@1517
   216
in the memory. We will also omit the typedefs.
athos@1517
   217
athos@1530
   218
\dontinclude lp_maxflow_demo.cc
athos@1530
   219
\skip Define a map on the edges for the variables of the LP problem
athos@1530
   220
\until lp.max();
athos@1530
   221
\skip Solve with the underlying solver
athos@1530
   222
\until lp.solve();
athos@1517
   223
athos@1517
   224
athos@1517
   225
The complete program can be found in file \ref lp_maxflow_demo.cc. After compiling run it in the form:
athos@1517
   226
athos@1528
   227
<tt>./lp_maxflow_demo < sample.lgf</tt>
athos@1517
   228
athos@1528
   229
where sample.lgf is a file in the lemon format containing a maxflow instance (designated "source", "target" nodes and "capacity" map on the edges).
athos@1517
   230
athos@1517
   231
athos@1514
   232
athos@1514
   233
</ol>
athos@1514
   234
</ul>
athos@1175
   235
athos@1175
   236
*/