lemon/lgf_writer.h
author Alpar Juttner <alpar@cs.elte.hu>
Wed, 04 Jun 2008 11:28:08 +0100
changeset 166 579979fad316
parent 163 c82fd9568d75
child 185 33e45a9b868c
permissions -rw-r--r--
lemon/concept_check.h was installed into a wrong directory
deba@127
     1
/* -*- C++ -*-
deba@127
     2
 *
deba@127
     3
 * This file is a part of LEMON, a generic C++ optimization library
deba@127
     4
 *
deba@127
     5
 * Copyright (C) 2003-2008
deba@127
     6
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
deba@127
     7
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
deba@127
     8
 *
deba@127
     9
 * Permission to use, modify and distribute this software is granted
deba@127
    10
 * provided that this copyright notice appears in all copies. For
deba@127
    11
 * precise terms see the accompanying LICENSE file.
deba@127
    12
 *
deba@127
    13
 * This software is provided "AS IS" with no warranty of any kind,
deba@127
    14
 * express or implied, and with no claim as to its suitability for any
deba@127
    15
 * purpose.
deba@127
    16
 *
deba@127
    17
 */
deba@127
    18
deba@127
    19
///\ingroup lemon_io
deba@127
    20
///\file
deba@127
    21
///\brief Lemon Graph Format writer.
deba@127
    22
deba@127
    23
deba@127
    24
#ifndef LEMON_LGF_WRITER_H
deba@127
    25
#define LEMON_LGF_WRITER_H
deba@127
    26
deba@127
    27
#include <iostream>
deba@127
    28
#include <fstream>
deba@127
    29
#include <sstream>
deba@127
    30
deba@127
    31
#include <algorithm>
deba@127
    32
deba@127
    33
#include <vector>
deba@127
    34
#include <functional>
deba@127
    35
deba@127
    36
#include <lemon/assert.h>
deba@127
    37
#include <lemon/graph_utils.h>
deba@127
    38
deba@127
    39
namespace lemon {
deba@127
    40
deba@127
    41
  namespace _writer_bits {
deba@127
    42
deba@127
    43
    template <typename Value>
deba@127
    44
    struct DefaultConverter {
deba@127
    45
      std::string operator()(const Value& value) {
deba@127
    46
	std::ostringstream os;
deba@127
    47
	os << value;
deba@127
    48
	return os.str();
deba@127
    49
      }
deba@127
    50
    };
deba@127
    51
deba@127
    52
    template <typename T>
deba@127
    53
    bool operator<(const T&, const T&) {
deba@127
    54
      throw DataFormatError("Label map is not comparable");
deba@127
    55
    }
deba@127
    56
deba@127
    57
    template <typename _Map>
deba@127
    58
    class MapLess {
deba@127
    59
    public:
deba@127
    60
      typedef _Map Map;
deba@127
    61
      typedef typename Map::Key Item;
deba@127
    62
deba@127
    63
    private:
deba@127
    64
      const Map& _map;
deba@127
    65
      
deba@127
    66
    public:
deba@127
    67
      MapLess(const Map& map) : _map(map) {}
deba@127
    68
deba@127
    69
      bool operator()(const Item& left, const Item& right) {
deba@127
    70
	return _map[left] < _map[right];
deba@127
    71
      }
deba@127
    72
    };
deba@127
    73
deba@165
    74
    template <typename _Graph, bool _dir, typename _Map>
deba@165
    75
    class GraphArcMapLess {
deba@165
    76
    public:
deba@165
    77
      typedef _Map Map;
deba@165
    78
      typedef _Graph Graph;
deba@165
    79
      typedef typename Graph::Edge Item;
deba@165
    80
deba@165
    81
    private:
deba@165
    82
      const Graph& _graph;
deba@165
    83
      const Map& _map;
deba@165
    84
      
deba@165
    85
    public:
deba@165
    86
      GraphArcMapLess(const Graph& graph, const Map& map) 
deba@165
    87
	: _graph(graph), _map(map) {}
deba@165
    88
deba@165
    89
      bool operator()(const Item& left, const Item& right) {
deba@165
    90
	return _map[_graph.direct(left, _dir)] < 
deba@165
    91
	  _map[_graph.direct(right, _dir)];
deba@165
    92
      }
deba@165
    93
    };
deba@165
    94
deba@127
    95
    template <typename _Item>    
deba@127
    96
    class MapStorageBase {
deba@127
    97
    public:
deba@127
    98
      typedef _Item Item;
deba@127
    99
deba@127
   100
    public:
deba@127
   101
      MapStorageBase() {}
deba@127
   102
      virtual ~MapStorageBase() {}
deba@127
   103
deba@127
   104
      virtual std::string get(const Item& item) = 0;
deba@127
   105
      virtual void sort(std::vector<Item>&) = 0;
deba@127
   106
    };
deba@127
   107
deba@127
   108
    template <typename _Item, typename _Map, 
deba@127
   109
	      typename _Converter = DefaultConverter<typename _Map::Value> >
deba@127
   110
    class MapStorage : public MapStorageBase<_Item> {
deba@127
   111
    public:
deba@127
   112
      typedef _Map Map;
deba@127
   113
      typedef _Converter Converter;
deba@127
   114
      typedef _Item Item;
deba@127
   115
      
deba@127
   116
    private:
deba@127
   117
      const Map& _map;
deba@127
   118
      Converter _converter;
deba@127
   119
deba@127
   120
    public:
deba@127
   121
      MapStorage(const Map& map, const Converter& converter = Converter()) 
deba@127
   122
	: _map(map), _converter(converter) {}
deba@127
   123
      virtual ~MapStorage() {}
deba@127
   124
deba@127
   125
      virtual std::string get(const Item& item) {
deba@127
   126
	return _converter(_map[item]);
deba@127
   127
      }
deba@127
   128
      virtual void sort(std::vector<Item>& items) {
deba@127
   129
	MapLess<Map> less(_map);
deba@127
   130
	std::sort(items.begin(), items.end(), less);
deba@127
   131
      }
deba@127
   132
    };
deba@127
   133
deba@165
   134
    template <typename _Graph, bool _dir, typename _Map, 
deba@165
   135
	      typename _Converter = DefaultConverter<typename _Map::Value> >
deba@165
   136
    class GraphArcMapStorage : public MapStorageBase<typename _Graph::Edge> {
deba@165
   137
    public:
deba@165
   138
      typedef _Map Map;
deba@165
   139
      typedef _Converter Converter;
deba@165
   140
      typedef _Graph Graph;
deba@165
   141
      typedef typename Graph::Edge Item;
deba@165
   142
      static const bool dir = _dir;
deba@165
   143
      
deba@165
   144
    private:
deba@165
   145
      const Graph& _graph;
deba@165
   146
      const Map& _map;
deba@165
   147
      Converter _converter;
deba@165
   148
deba@165
   149
    public:
deba@165
   150
      GraphArcMapStorage(const Graph& graph, const Map& map,  
deba@165
   151
			 const Converter& converter = Converter()) 
deba@165
   152
	: _graph(graph), _map(map), _converter(converter) {}
deba@165
   153
      virtual ~GraphArcMapStorage() {}
deba@165
   154
deba@165
   155
      virtual std::string get(const Item& item) {
deba@165
   156
	return _converter(_map[_graph.direct(item, dir)]);
deba@165
   157
      }
deba@165
   158
      virtual void sort(std::vector<Item>& items) {
deba@165
   159
	GraphArcMapLess<Graph, dir, Map> less(_graph, _map);
deba@165
   160
	std::sort(items.begin(), items.end(), less);
deba@165
   161
      }
deba@165
   162
    };
deba@165
   163
deba@127
   164
    class ValueStorageBase {
deba@127
   165
    public:
deba@127
   166
      ValueStorageBase() {}
deba@127
   167
      virtual ~ValueStorageBase() {}
deba@127
   168
deba@127
   169
      virtual std::string get() = 0;      
deba@127
   170
    };
deba@127
   171
deba@127
   172
    template <typename _Value, typename _Converter = DefaultConverter<_Value> >
deba@127
   173
    class ValueStorage : public ValueStorageBase {
deba@127
   174
    public:
deba@127
   175
      typedef _Value Value;
deba@127
   176
      typedef _Converter Converter;
deba@127
   177
deba@127
   178
    private:
deba@127
   179
      const Value& _value;
deba@127
   180
      Converter _converter;
deba@127
   181
deba@127
   182
    public:
deba@127
   183
      ValueStorage(const Value& value, const Converter& converter = Converter())
deba@127
   184
 	: _value(value), _converter(converter) {}
deba@127
   185
deba@127
   186
      virtual std::string get() {
deba@127
   187
	return _converter(_value);
deba@127
   188
      }
deba@127
   189
    };
deba@127
   190
deba@127
   191
    template <typename Value>
deba@127
   192
    struct MapLookUpConverter {
deba@127
   193
      const std::map<Value, std::string>& _map;
deba@127
   194
      
deba@127
   195
      MapLookUpConverter(const std::map<Value, std::string>& map) 
deba@127
   196
	: _map(map) {}
deba@127
   197
      
deba@127
   198
      std::string operator()(const Value& str) {
deba@127
   199
	typename std::map<Value, std::string>::const_iterator it = 
deba@127
   200
	  _map.find(str);
deba@127
   201
	if (it == _map.end()) {
deba@127
   202
	  throw DataFormatError("Item not found");
deba@127
   203
	}
deba@127
   204
	return it->second;
deba@127
   205
      }
deba@127
   206
    };
deba@127
   207
deba@165
   208
    template <typename Graph>
deba@165
   209
    struct GraphArcLookUpConverter {
deba@165
   210
      const Graph& _graph;
deba@165
   211
      const std::map<typename Graph::Edge, std::string>& _map;
deba@165
   212
      
deba@165
   213
      GraphArcLookUpConverter(const Graph& graph, 
deba@165
   214
			      const std::map<typename Graph::Edge, 
deba@165
   215
			                     std::string>& map) 
deba@165
   216
	: _graph(graph), _map(map) {}
deba@165
   217
      
deba@165
   218
      std::string operator()(const typename Graph::Arc& val) {
deba@165
   219
	typename std::map<typename Graph::Edge, std::string>
deba@165
   220
	  ::const_iterator it = _map.find(val);
deba@165
   221
	if (it == _map.end()) {
deba@165
   222
	  throw DataFormatError("Item not found");
deba@165
   223
	}
deba@165
   224
	return (_graph.direction(val) ? '+' : '-') + it->second;
deba@165
   225
      }
deba@165
   226
    };
deba@165
   227
deba@127
   228
    bool isWhiteSpace(char c) {
deba@127
   229
      return c == ' ' || c == '\t' || c == '\v' || 
deba@127
   230
        c == '\n' || c == '\r' || c == '\f'; 
deba@127
   231
    }
deba@127
   232
deba@127
   233
    bool isEscaped(char c) {
deba@127
   234
      return c == '\\' || c == '\"' || c == '\'' || 
deba@127
   235
	c == '\a' || c == '\b';
deba@127
   236
    }
deba@127
   237
deba@127
   238
    static void writeEscape(std::ostream& os, char c) {
deba@127
   239
      switch (c) {
deba@127
   240
      case '\\':
deba@127
   241
	os << "\\\\";
deba@127
   242
	return;
deba@127
   243
      case '\"':
deba@127
   244
	os << "\\\"";
deba@127
   245
	return;
deba@127
   246
      case '\a':
deba@127
   247
	os << "\\a";
deba@127
   248
	return;
deba@127
   249
      case '\b':
deba@127
   250
	os << "\\b";
deba@127
   251
	return;
deba@127
   252
      case '\f':
deba@127
   253
	os << "\\f";
deba@127
   254
	return;
deba@127
   255
      case '\r':
deba@127
   256
	os << "\\r";
deba@127
   257
	return;
deba@127
   258
      case '\n':
deba@127
   259
	os << "\\n";
deba@127
   260
	return;
deba@127
   261
      case '\t':
deba@127
   262
	os << "\\t";
deba@127
   263
	return;
deba@127
   264
      case '\v':
deba@127
   265
	os << "\\v";
deba@127
   266
	return;
deba@127
   267
      default:
deba@127
   268
	if (c < 0x20) {
deba@163
   269
	  std::ios::fmtflags flags = os.flags();
deba@127
   270
	  os << '\\' << std::oct << static_cast<int>(c);
deba@163
   271
	  os.flags(flags);
deba@127
   272
	} else {
deba@127
   273
	  os << c;
deba@127
   274
	}
deba@127
   275
	return;
deba@127
   276
      }     
deba@127
   277
    }
deba@127
   278
deba@127
   279
    bool requireEscape(const std::string& str) {
alpar@156
   280
      if (str.empty() || str[0] == '@') return true;
deba@127
   281
      std::istringstream is(str);
deba@127
   282
      char c;
deba@127
   283
      while (is.get(c)) {
deba@127
   284
	if (isWhiteSpace(c) || isEscaped(c)) {
deba@127
   285
	  return true;
deba@127
   286
	}
deba@127
   287
      }
deba@127
   288
      return false;
deba@127
   289
    }
deba@127
   290
    
deba@127
   291
    std::ostream& writeToken(std::ostream& os, const std::string& str) {
deba@127
   292
deba@127
   293
      if (requireEscape(str)) {
deba@127
   294
	os << '\"';
deba@127
   295
	for (std::string::const_iterator it = str.begin(); 
deba@127
   296
	     it != str.end(); ++it) {
deba@127
   297
	  writeEscape(os, *it);
deba@127
   298
	}	
deba@127
   299
	os << '\"';
deba@127
   300
      } else {
deba@127
   301
	os << str;
deba@127
   302
      }
deba@127
   303
      return os;
deba@127
   304
    }
deba@127
   305
deba@127
   306
  }
deba@127
   307
  
alpar@156
   308
  /// \ingroup lemon_io
alpar@156
   309
  ///  
alpar@156
   310
  /// \brief LGF writer for directed graphs
alpar@156
   311
  ///
alpar@156
   312
  /// This utility writes an \ref lgf-format "LGF" file.
alpar@156
   313
  ///
alpar@156
   314
  /// The writing method does a batch processing. The user creates a
alpar@156
   315
  /// writer object, then various writing rules can be added to the
alpar@156
   316
  /// writer, and eventually the writing is executed with the \c run()
alpar@156
   317
  /// member function. A map writing rule can be added to the writer
alpar@156
   318
  /// with the \c nodeMap() or \c arcMap() members. An optional
deba@163
   319
  /// converter parameter can also be added as a standard functor
deba@163
   320
  /// converting from the value type of the map to std::string. If it
deba@163
   321
  /// is set, it will determine how the map's value type is written to
deba@163
   322
  /// the output stream. If the functor is not set, then a default
deba@163
   323
  /// conversion will be used. The \c attribute(), \c node() and \c
deba@163
   324
  /// arc() functions are used to add attribute writing rules.
alpar@156
   325
  ///
alpar@156
   326
  ///\code
alpar@156
   327
  ///     DigraphWriter<Digraph>(std::cout, digraph).
alpar@156
   328
  ///       nodeMap("coordinates", coord_map).
alpar@156
   329
  ///       nodeMap("size", size).
alpar@156
   330
  ///       nodeMap("title", title).
alpar@156
   331
  ///       arcMap("capacity", cap_map).
alpar@156
   332
  ///       node("source", src).
alpar@156
   333
  ///       node("target", trg).
alpar@156
   334
  ///       attribute("caption", caption).
alpar@156
   335
  ///       run();
alpar@156
   336
  ///\endcode
alpar@156
   337
  ///
alpar@156
   338
  ///
alpar@156
   339
  /// By default, the writer does not write additional captions to the
alpar@156
   340
  /// sections, but they can be give as an optional parameter of
alpar@156
   341
  /// the \c nodes(), \c arcs() or \c
alpar@156
   342
  /// attributes() functions.
alpar@156
   343
  ///
alpar@156
   344
  /// The \c skipNodes() and \c skipArcs() functions forbid the
deba@163
   345
  /// writing of the sections. If two arc sections should be written
deba@163
   346
  /// to the output, it can be done in two passes, the first pass
deba@163
   347
  /// writes the node section and the first arc section, then the
deba@163
   348
  /// second pass skips the node section and writes just the arc
deba@163
   349
  /// section to the stream. The output stream can be retrieved with
deba@163
   350
  /// the \c ostream() function, hence the second pass can append its
deba@163
   351
  /// output to the output of the first pass.
deba@127
   352
  template <typename _Digraph>
deba@127
   353
  class DigraphWriter {
deba@127
   354
  public:
deba@127
   355
deba@127
   356
    typedef _Digraph Digraph;
deba@148
   357
    TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
deba@127
   358
    
deba@127
   359
  private:
deba@127
   360
deba@127
   361
deba@127
   362
    std::ostream* _os;
deba@127
   363
    bool local_os;
deba@127
   364
deba@127
   365
    Digraph& _digraph;
deba@127
   366
deba@127
   367
    std::string _nodes_caption;
deba@127
   368
    std::string _arcs_caption;
deba@127
   369
    std::string _attributes_caption;
deba@127
   370
    
deba@127
   371
    typedef std::map<Node, std::string> NodeIndex;
deba@127
   372
    NodeIndex _node_index;
deba@127
   373
    typedef std::map<Arc, std::string> ArcIndex;
deba@127
   374
    ArcIndex _arc_index;
deba@127
   375
deba@127
   376
    typedef std::vector<std::pair<std::string, 
deba@127
   377
      _writer_bits::MapStorageBase<Node>* > > NodeMaps;    
deba@127
   378
    NodeMaps _node_maps; 
deba@127
   379
deba@127
   380
    typedef std::vector<std::pair<std::string, 
deba@127
   381
      _writer_bits::MapStorageBase<Arc>* > >ArcMaps;
deba@127
   382
    ArcMaps _arc_maps;
deba@127
   383
deba@127
   384
    typedef std::vector<std::pair<std::string, 
deba@127
   385
      _writer_bits::ValueStorageBase*> > Attributes;
deba@127
   386
    Attributes _attributes;
deba@127
   387
deba@127
   388
    bool _skip_nodes;
deba@127
   389
    bool _skip_arcs;
deba@127
   390
deba@127
   391
  public:
deba@127
   392
alpar@156
   393
    /// \brief Constructor
alpar@156
   394
    ///
alpar@156
   395
    /// Construct a directed graph writer, which writes to the given
alpar@156
   396
    /// output stream.
deba@127
   397
    DigraphWriter(std::ostream& is, Digraph& digraph) 
deba@127
   398
      : _os(&is), local_os(false), _digraph(digraph),
deba@127
   399
	_skip_nodes(false), _skip_arcs(false) {}
deba@127
   400
alpar@156
   401
    /// \brief Constructor
alpar@156
   402
    ///
alpar@156
   403
    /// Construct a directed graph writer, which writes to the given
alpar@156
   404
    /// output file.
deba@127
   405
    DigraphWriter(const std::string& fn, Digraph& digraph) 
deba@127
   406
      : _os(new std::ofstream(fn.c_str())), local_os(true), _digraph(digraph),
deba@127
   407
	_skip_nodes(false), _skip_arcs(false) {}
deba@127
   408
alpar@156
   409
    /// \brief Constructor
alpar@156
   410
    ///
alpar@156
   411
    /// Construct a directed graph writer, which writes to the given
alpar@156
   412
    /// output file.
deba@127
   413
    DigraphWriter(const char* fn, Digraph& digraph) 
deba@127
   414
      : _os(new std::ofstream(fn)), local_os(true), _digraph(digraph),
deba@127
   415
	_skip_nodes(false), _skip_arcs(false) {}
deba@127
   416
alpar@156
   417
    /// \brief Copy constructor
alpar@156
   418
    ///
alpar@156
   419
    /// The copy constructor transfers all data from the other writer,
alpar@156
   420
    /// therefore the copied writer will not be usable more. 
deba@127
   421
    DigraphWriter(DigraphWriter& other) 
deba@127
   422
      : _os(other._os), local_os(other.local_os), _digraph(other._digraph),
deba@127
   423
	_skip_nodes(other._skip_nodes), _skip_arcs(other._skip_arcs) {
deba@127
   424
deba@163
   425
      other._os = 0;
deba@127
   426
      other.local_os = false;
deba@127
   427
deba@127
   428
      _node_index.swap(other._node_index);
deba@127
   429
      _arc_index.swap(other._arc_index);
deba@127
   430
deba@127
   431
      _node_maps.swap(other._node_maps);
deba@127
   432
      _arc_maps.swap(other._arc_maps);
deba@127
   433
      _attributes.swap(other._attributes);
deba@127
   434
deba@127
   435
      _nodes_caption = other._nodes_caption;
deba@127
   436
      _arcs_caption = other._arcs_caption;
deba@127
   437
      _attributes_caption = other._attributes_caption;
deba@127
   438
    }
deba@127
   439
alpar@156
   440
    /// \brief Destructor
deba@127
   441
    ~DigraphWriter() {
deba@127
   442
      for (typename NodeMaps::iterator it = _node_maps.begin(); 
deba@127
   443
	   it != _node_maps.end(); ++it) {
deba@127
   444
	delete it->second;
deba@127
   445
      }
deba@127
   446
deba@127
   447
      for (typename ArcMaps::iterator it = _arc_maps.begin(); 
deba@127
   448
	   it != _arc_maps.end(); ++it) {
deba@127
   449
	delete it->second;
deba@127
   450
      }
deba@127
   451
deba@127
   452
      for (typename Attributes::iterator it = _attributes.begin(); 
deba@127
   453
	   it != _attributes.end(); ++it) {
deba@127
   454
	delete it->second;
deba@127
   455
      }
deba@127
   456
deba@127
   457
      if (local_os) {
deba@127
   458
	delete _os;
deba@127
   459
      }
deba@127
   460
    }
deba@127
   461
deba@127
   462
  private:
deba@127
   463
    
deba@127
   464
    DigraphWriter& operator=(const DigraphWriter&);
deba@127
   465
deba@127
   466
  public:
deba@127
   467
alpar@156
   468
    /// \name Writing rules
alpar@156
   469
    /// @{
alpar@156
   470
    
alpar@156
   471
    /// \brief Node map reading rule
alpar@156
   472
    ///
alpar@156
   473
    /// Add a node map reading rule to the writer.
deba@127
   474
    template <typename Map>
deba@127
   475
    DigraphWriter& nodeMap(const std::string& caption, const Map& map) {
deba@127
   476
      checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@127
   477
      _writer_bits::MapStorageBase<Node>* storage = 
deba@127
   478
	new _writer_bits::MapStorage<Node, Map>(map);
deba@127
   479
      _node_maps.push_back(std::make_pair(caption, storage));
deba@127
   480
      return *this;
deba@127
   481
    }
deba@127
   482
alpar@156
   483
    /// \brief Node map writing rule
alpar@156
   484
    ///
alpar@156
   485
    /// Add a node map writing rule with specialized converter to the
alpar@156
   486
    /// writer.
deba@127
   487
    template <typename Map, typename Converter>
deba@127
   488
    DigraphWriter& nodeMap(const std::string& caption, const Map& map, 
deba@127
   489
			   const Converter& converter = Converter()) {
deba@127
   490
      checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@127
   491
      _writer_bits::MapStorageBase<Node>* storage = 
deba@127
   492
	new _writer_bits::MapStorage<Node, Map, Converter>(map, converter);
deba@127
   493
      _node_maps.push_back(std::make_pair(caption, storage));
deba@127
   494
      return *this;
deba@127
   495
    }
deba@127
   496
alpar@156
   497
    /// \brief Arc map writing rule
alpar@156
   498
    ///
alpar@156
   499
    /// Add an arc map writing rule to the writer.
deba@127
   500
    template <typename Map>
deba@127
   501
    DigraphWriter& arcMap(const std::string& caption, const Map& map) {
deba@127
   502
      checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
deba@127
   503
      _writer_bits::MapStorageBase<Arc>* storage = 
deba@127
   504
	new _writer_bits::MapStorage<Arc, Map>(map);
deba@127
   505
      _arc_maps.push_back(std::make_pair(caption, storage));
deba@127
   506
      return *this;
deba@127
   507
    }
deba@127
   508
alpar@156
   509
    /// \brief Arc map writing rule
alpar@156
   510
    ///
alpar@156
   511
    /// Add an arc map writing rule with specialized converter to the
alpar@156
   512
    /// writer.
deba@127
   513
    template <typename Map, typename Converter>
deba@127
   514
    DigraphWriter& arcMap(const std::string& caption, const Map& map, 
deba@127
   515
			  const Converter& converter = Converter()) {
deba@127
   516
      checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
deba@127
   517
      _writer_bits::MapStorageBase<Arc>* storage = 
deba@127
   518
	new _writer_bits::MapStorage<Arc, Map, Converter>(map, converter);
deba@127
   519
      _arc_maps.push_back(std::make_pair(caption, storage));
deba@127
   520
      return *this;
deba@127
   521
    }
deba@127
   522
alpar@156
   523
    /// \brief Attribute writing rule
alpar@156
   524
    ///
alpar@156
   525
    /// Add an attribute writing rule to the writer.
deba@127
   526
    template <typename Value>
deba@127
   527
    DigraphWriter& attribute(const std::string& caption, const Value& value) {
deba@127
   528
      _writer_bits::ValueStorageBase* storage = 
deba@127
   529
	new _writer_bits::ValueStorage<Value>(value);
deba@127
   530
      _attributes.push_back(std::make_pair(caption, storage));
deba@127
   531
      return *this;
deba@127
   532
    }
deba@127
   533
alpar@156
   534
    /// \brief Attribute writing rule
alpar@156
   535
    ///
alpar@156
   536
    /// Add an attribute writing rule with specialized converter to the
alpar@156
   537
    /// writer.
deba@127
   538
    template <typename Value, typename Converter>
deba@127
   539
    DigraphWriter& attribute(const std::string& caption, const Value& value, 
deba@127
   540
			     const Converter& converter = Converter()) {
deba@127
   541
      _writer_bits::ValueStorageBase* storage = 
deba@127
   542
	new _writer_bits::ValueStorage<Value, Converter>(value, converter);
deba@127
   543
      _attributes.push_back(std::make_pair(caption, storage));
deba@127
   544
      return *this;
deba@127
   545
    }
deba@127
   546
alpar@156
   547
    /// \brief Node writing rule
alpar@156
   548
    ///
alpar@156
   549
    /// Add a node writing rule to the writer.
deba@127
   550
    DigraphWriter& node(const std::string& caption, const Node& node) {
deba@127
   551
      typedef _writer_bits::MapLookUpConverter<Node> Converter;
deba@127
   552
      Converter converter(_node_index);
deba@127
   553
      _writer_bits::ValueStorageBase* storage = 
deba@127
   554
	new _writer_bits::ValueStorage<Node, Converter>(node, converter);
deba@127
   555
      _attributes.push_back(std::make_pair(caption, storage));
deba@127
   556
      return *this;
deba@127
   557
    }
deba@127
   558
alpar@156
   559
    /// \brief Arc writing rule
alpar@156
   560
    ///
alpar@156
   561
    /// Add an arc writing rule to writer.
deba@127
   562
    DigraphWriter& arc(const std::string& caption, const Arc& arc) {
deba@127
   563
      typedef _writer_bits::MapLookUpConverter<Arc> Converter;
deba@127
   564
      Converter converter(_arc_index);
deba@127
   565
      _writer_bits::ValueStorageBase* storage = 
deba@127
   566
	new _writer_bits::ValueStorage<Arc, Converter>(arc, converter);
deba@127
   567
      _attributes.push_back(std::make_pair(caption, storage));
deba@127
   568
      return *this;
deba@127
   569
    }
deba@127
   570
alpar@156
   571
    /// \name Select section by name
alpar@156
   572
    /// @{
alpar@156
   573
alpar@156
   574
    /// \brief Set \c \@nodes section to be read
alpar@156
   575
    ///
alpar@156
   576
    /// Set \c \@nodes section to be read
deba@127
   577
    DigraphWriter& nodes(const std::string& caption) {
deba@127
   578
      _nodes_caption = caption;
deba@127
   579
      return *this;
deba@127
   580
    }
deba@127
   581
alpar@156
   582
    /// \brief Set \c \@arcs section to be read
alpar@156
   583
    ///
alpar@156
   584
    /// Set \c \@arcs section to be read
deba@127
   585
    DigraphWriter& arcs(const std::string& caption) {
deba@127
   586
      _arcs_caption = caption;
deba@127
   587
      return *this;
deba@127
   588
    }
deba@127
   589
alpar@156
   590
    /// \brief Set \c \@attributes section to be read
alpar@156
   591
    ///
alpar@156
   592
    /// Set \c \@attributes section to be read
deba@127
   593
    DigraphWriter& attributes(const std::string& caption) {
deba@127
   594
      _attributes_caption = caption;
deba@127
   595
      return *this;
deba@127
   596
    }
deba@127
   597
alpar@156
   598
    /// \name Skipping section
alpar@156
   599
    /// @{
alpar@156
   600
alpar@156
   601
    /// \brief Skip writing the node set
alpar@156
   602
    ///
alpar@156
   603
    /// The \c \@nodes section will be not written to the stream.
deba@127
   604
    DigraphWriter& skipNodes() {
deba@127
   605
      LEMON_ASSERT(!_skip_nodes, "Multiple usage of skipNodes() member");
deba@127
   606
      return *this;
deba@127
   607
    }
deba@127
   608
alpar@156
   609
    /// \brief Skip writing arc set
alpar@156
   610
    ///
alpar@156
   611
    /// The \c \@arcs section will be not written to the stream.
deba@127
   612
    DigraphWriter& skipArcs() {
deba@127
   613
      LEMON_ASSERT(!_skip_arcs, "Multiple usage of skipArcs() member");
deba@127
   614
      return *this;
deba@127
   615
    }
deba@127
   616
alpar@156
   617
    /// @}
alpar@156
   618
deba@127
   619
  private:
deba@127
   620
deba@127
   621
    void writeNodes() {
deba@127
   622
      _writer_bits::MapStorageBase<Node>* label = 0;
deba@127
   623
      for (typename NodeMaps::iterator it = _node_maps.begin();
deba@127
   624
	   it != _node_maps.end(); ++it) {
deba@127
   625
        if (it->first == "label") {
deba@127
   626
	  label = it->second;
deba@127
   627
	  break;
deba@127
   628
	}
deba@127
   629
      }
deba@127
   630
deba@127
   631
      *_os << "@nodes";
deba@127
   632
      if (!_nodes_caption.empty()) {
alpar@156
   633
	_writer_bits::writeToken(*_os << ' ', _nodes_caption);
deba@127
   634
      }
deba@127
   635
      *_os << std::endl;
deba@127
   636
deba@127
   637
      if (label == 0) {
deba@127
   638
	*_os << "label" << '\t';
deba@127
   639
      }
deba@127
   640
      for (typename NodeMaps::iterator it = _node_maps.begin();
deba@127
   641
	   it != _node_maps.end(); ++it) {
alpar@156
   642
	_writer_bits::writeToken(*_os, it->first) << '\t';
deba@127
   643
      }
deba@127
   644
      *_os << std::endl;
deba@127
   645
deba@127
   646
      std::vector<Node> nodes;
deba@127
   647
      for (NodeIt n(_digraph); n != INVALID; ++n) {
deba@127
   648
	nodes.push_back(n);
deba@127
   649
      }
deba@127
   650
      
deba@127
   651
      if (label == 0) {
deba@127
   652
	IdMap<Digraph, Node> id_map(_digraph);
deba@127
   653
	_writer_bits::MapLess<IdMap<Digraph, Node> > id_less(id_map);
deba@127
   654
	std::sort(nodes.begin(), nodes.end(), id_less);
deba@127
   655
      } else {
deba@127
   656
	label->sort(nodes);
deba@127
   657
      }
deba@127
   658
deba@127
   659
      for (int i = 0; i < static_cast<int>(nodes.size()); ++i) {
deba@127
   660
	Node n = nodes[i];
deba@127
   661
	if (label == 0) {
deba@127
   662
	  std::ostringstream os;
deba@127
   663
	  os << _digraph.id(n);
deba@127
   664
	  _writer_bits::writeToken(*_os, os.str());
deba@127
   665
	  *_os << '\t';
deba@127
   666
	  _node_index.insert(std::make_pair(n, os.str()));
deba@127
   667
	}
deba@127
   668
	for (typename NodeMaps::iterator it = _node_maps.begin();
deba@127
   669
	     it != _node_maps.end(); ++it) {
deba@127
   670
	  std::string value = it->second->get(n);
deba@127
   671
	  _writer_bits::writeToken(*_os, value);
deba@127
   672
	  if (it->first == "label") {
deba@127
   673
	    _node_index.insert(std::make_pair(n, value));
deba@127
   674
	  }
deba@127
   675
	  *_os << '\t';
deba@127
   676
	}
deba@127
   677
	*_os << std::endl;
deba@127
   678
      }
deba@127
   679
    }
deba@127
   680
deba@127
   681
    void writeArcs() {
deba@127
   682
      _writer_bits::MapStorageBase<Arc>* label = 0;
deba@127
   683
      for (typename ArcMaps::iterator it = _arc_maps.begin();
deba@127
   684
	   it != _arc_maps.end(); ++it) {
deba@127
   685
        if (it->first == "label") {
deba@127
   686
	  label = it->second;
deba@127
   687
	  break;
deba@127
   688
	}
deba@127
   689
      }
deba@127
   690
deba@127
   691
      *_os << "@arcs";
deba@127
   692
      if (!_arcs_caption.empty()) {
alpar@156
   693
	_writer_bits::writeToken(*_os << ' ', _arcs_caption);
deba@127
   694
      }
deba@127
   695
      *_os << std::endl;
deba@127
   696
deba@127
   697
      *_os << '\t' << '\t';
deba@127
   698
      if (label == 0) {
deba@127
   699
	*_os << "label" << '\t';
deba@127
   700
      }
deba@127
   701
      for (typename ArcMaps::iterator it = _arc_maps.begin();
deba@127
   702
	   it != _arc_maps.end(); ++it) {
alpar@156
   703
	_writer_bits::writeToken(*_os, it->first) << '\t';
deba@127
   704
      }
deba@127
   705
      *_os << std::endl;
deba@127
   706
deba@127
   707
      std::vector<Arc> arcs;
deba@127
   708
      for (ArcIt n(_digraph); n != INVALID; ++n) {
deba@127
   709
	arcs.push_back(n);
deba@127
   710
      }
deba@127
   711
      
deba@127
   712
      if (label == 0) {
deba@127
   713
	IdMap<Digraph, Arc> id_map(_digraph);
deba@127
   714
	_writer_bits::MapLess<IdMap<Digraph, Arc> > id_less(id_map);
deba@127
   715
	std::sort(arcs.begin(), arcs.end(), id_less);
deba@127
   716
      } else {
deba@127
   717
	label->sort(arcs);
deba@127
   718
      }
deba@127
   719
deba@127
   720
      for (int i = 0; i < static_cast<int>(arcs.size()); ++i) {
deba@127
   721
	Arc a = arcs[i];
deba@127
   722
	_writer_bits::writeToken(*_os, _node_index.
deba@127
   723
				 find(_digraph.source(a))->second);
deba@127
   724
	*_os << '\t';
deba@127
   725
	_writer_bits::writeToken(*_os, _node_index.
deba@127
   726
				 find(_digraph.target(a))->second);
deba@127
   727
	*_os << '\t';
deba@127
   728
	if (label == 0) {
deba@127
   729
	  std::ostringstream os;
deba@127
   730
	  os << _digraph.id(a);
deba@127
   731
	  _writer_bits::writeToken(*_os, os.str());
deba@127
   732
	  *_os << '\t';
deba@127
   733
	  _arc_index.insert(std::make_pair(a, os.str()));
deba@127
   734
	}
deba@127
   735
	for (typename ArcMaps::iterator it = _arc_maps.begin();
deba@127
   736
	     it != _arc_maps.end(); ++it) {
deba@127
   737
	  std::string value = it->second->get(a);
deba@127
   738
	  _writer_bits::writeToken(*_os, value);
deba@127
   739
	  if (it->first == "label") {
deba@127
   740
	    _arc_index.insert(std::make_pair(a, value));
deba@127
   741
	  }
deba@127
   742
	  *_os << '\t';
deba@127
   743
	}
deba@127
   744
	*_os << std::endl;
deba@127
   745
      }
deba@127
   746
    }
deba@127
   747
deba@127
   748
    void writeAttributes() {
deba@127
   749
      if (_attributes.empty()) return;
deba@127
   750
      *_os << "@attributes";
deba@127
   751
      if (!_attributes_caption.empty()) {
alpar@156
   752
	_writer_bits::writeToken(*_os << ' ', _attributes_caption);
deba@127
   753
      }
deba@127
   754
      *_os << std::endl;
deba@127
   755
      for (typename Attributes::iterator it = _attributes.begin();
deba@127
   756
	   it != _attributes.end(); ++it) {
alpar@156
   757
	_writer_bits::writeToken(*_os, it->first) << ' ';
deba@127
   758
	_writer_bits::writeToken(*_os, it->second->get());
deba@127
   759
	*_os << std::endl;
deba@127
   760
      }
deba@127
   761
    }
deba@127
   762
    
deba@127
   763
  public:
deba@127
   764
    
alpar@156
   765
    /// \name Execution of the writer    
alpar@156
   766
    /// @{
alpar@156
   767
alpar@156
   768
    /// \brief Start the batch processing
alpar@156
   769
    ///
alpar@156
   770
    /// This function starts the batch processing
deba@127
   771
    void run() {
deba@127
   772
      if (!_skip_nodes) {
deba@127
   773
	writeNodes();
deba@127
   774
      }
deba@127
   775
      if (!_skip_arcs) {      
deba@127
   776
	writeArcs();
deba@127
   777
      }
deba@127
   778
      writeAttributes();
deba@127
   779
    }
deba@127
   780
alpar@156
   781
    /// \brief Gives back the stream of the writer
alpar@156
   782
    ///
alpar@156
   783
    /// Gives back the stream of the writer
alpar@156
   784
    std::ostream& ostream() {
deba@127
   785
      return *_os;
deba@127
   786
    }
alpar@156
   787
alpar@156
   788
    /// @}
deba@127
   789
  };
deba@127
   790
alpar@156
   791
  /// \relates DigraphWriter
deba@127
   792
  template <typename Digraph>
deba@163
   793
  DigraphWriter<Digraph> digraphWriter(std::ostream& os, Digraph& digraph) {
deba@163
   794
    DigraphWriter<Digraph> tmp(os, digraph);
deba@163
   795
    return tmp;
deba@127
   796
  }
deba@127
   797
alpar@156
   798
  /// \relates DigraphWriter
deba@127
   799
  template <typename Digraph>
deba@127
   800
  DigraphWriter<Digraph> digraphWriter(const std::string& fn, 
deba@127
   801
				       Digraph& digraph) {
deba@163
   802
    DigraphWriter<Digraph> tmp(fn, digraph);
deba@163
   803
    return tmp;
deba@127
   804
  }
deba@127
   805
alpar@156
   806
  /// \relates DigraphWriter
deba@127
   807
  template <typename Digraph>
deba@127
   808
  DigraphWriter<Digraph> digraphWriter(const char* fn, Digraph& digraph) {
deba@163
   809
    DigraphWriter<Digraph> tmp(fn, digraph);
deba@163
   810
    return tmp;
deba@127
   811
  }
deba@165
   812
deba@165
   813
  /// \ingroup lemon_io
deba@165
   814
  ///  
deba@165
   815
  /// \brief LGF writer for directed graphs
deba@165
   816
  ///
deba@165
   817
  /// This utility writes an \ref lgf-format "LGF" file.
deba@165
   818
  template <typename _Graph>
deba@165
   819
  class GraphWriter {
deba@165
   820
  public:
deba@165
   821
deba@165
   822
    typedef _Graph Graph;
deba@165
   823
    TEMPLATE_GRAPH_TYPEDEFS(Graph);
deba@165
   824
    
deba@165
   825
  private:
deba@165
   826
deba@165
   827
deba@165
   828
    std::ostream* _os;
deba@165
   829
    bool local_os;
deba@165
   830
deba@165
   831
    Graph& _graph;
deba@165
   832
deba@165
   833
    std::string _nodes_caption;
deba@165
   834
    std::string _edges_caption;
deba@165
   835
    std::string _attributes_caption;
deba@165
   836
    
deba@165
   837
    typedef std::map<Node, std::string> NodeIndex;
deba@165
   838
    NodeIndex _node_index;
deba@165
   839
    typedef std::map<Edge, std::string> EdgeIndex;
deba@165
   840
    EdgeIndex _edge_index;
deba@165
   841
deba@165
   842
    typedef std::vector<std::pair<std::string, 
deba@165
   843
      _writer_bits::MapStorageBase<Node>* > > NodeMaps;    
deba@165
   844
    NodeMaps _node_maps; 
deba@165
   845
deba@165
   846
    typedef std::vector<std::pair<std::string, 
deba@165
   847
      _writer_bits::MapStorageBase<Edge>* > >EdgeMaps;
deba@165
   848
    EdgeMaps _edge_maps;
deba@165
   849
deba@165
   850
    typedef std::vector<std::pair<std::string, 
deba@165
   851
      _writer_bits::ValueStorageBase*> > Attributes;
deba@165
   852
    Attributes _attributes;
deba@165
   853
deba@165
   854
    bool _skip_nodes;
deba@165
   855
    bool _skip_edges;
deba@165
   856
deba@165
   857
  public:
deba@165
   858
deba@165
   859
    /// \brief Constructor
deba@165
   860
    ///
deba@165
   861
    /// Construct a directed graph writer, which writes to the given
deba@165
   862
    /// output stream.
deba@165
   863
    GraphWriter(std::ostream& is, Graph& graph) 
deba@165
   864
      : _os(&is), local_os(false), _graph(graph),
deba@165
   865
	_skip_nodes(false), _skip_edges(false) {}
deba@165
   866
deba@165
   867
    /// \brief Constructor
deba@165
   868
    ///
deba@165
   869
    /// Construct a directed graph writer, which writes to the given
deba@165
   870
    /// output file.
deba@165
   871
    GraphWriter(const std::string& fn, Graph& graph) 
deba@165
   872
      : _os(new std::ofstream(fn.c_str())), local_os(true), _graph(graph),
deba@165
   873
	_skip_nodes(false), _skip_edges(false) {}
deba@165
   874
deba@165
   875
    /// \brief Constructor
deba@165
   876
    ///
deba@165
   877
    /// Construct a directed graph writer, which writes to the given
deba@165
   878
    /// output file.
deba@165
   879
    GraphWriter(const char* fn, Graph& graph) 
deba@165
   880
      : _os(new std::ofstream(fn)), local_os(true), _graph(graph),
deba@165
   881
	_skip_nodes(false), _skip_edges(false) {}
deba@165
   882
deba@165
   883
    /// \brief Copy constructor
deba@165
   884
    ///
deba@165
   885
    /// The copy constructor transfers all data from the other writer,
deba@165
   886
    /// therefore the copied writer will not be usable more. 
deba@165
   887
    GraphWriter(GraphWriter& other) 
deba@165
   888
      : _os(other._os), local_os(other.local_os), _graph(other._graph),
deba@165
   889
	_skip_nodes(other._skip_nodes), _skip_edges(other._skip_edges) {
deba@165
   890
deba@165
   891
      other._os = 0;
deba@165
   892
      other.local_os = false;
deba@165
   893
deba@165
   894
      _node_index.swap(other._node_index);
deba@165
   895
      _edge_index.swap(other._edge_index);
deba@165
   896
deba@165
   897
      _node_maps.swap(other._node_maps);
deba@165
   898
      _edge_maps.swap(other._edge_maps);
deba@165
   899
      _attributes.swap(other._attributes);
deba@165
   900
deba@165
   901
      _nodes_caption = other._nodes_caption;
deba@165
   902
      _edges_caption = other._edges_caption;
deba@165
   903
      _attributes_caption = other._attributes_caption;
deba@165
   904
    }
deba@165
   905
deba@165
   906
    /// \brief Destructor
deba@165
   907
    ~GraphWriter() {
deba@165
   908
      for (typename NodeMaps::iterator it = _node_maps.begin(); 
deba@165
   909
	   it != _node_maps.end(); ++it) {
deba@165
   910
	delete it->second;
deba@165
   911
      }
deba@165
   912
deba@165
   913
      for (typename EdgeMaps::iterator it = _edge_maps.begin(); 
deba@165
   914
	   it != _edge_maps.end(); ++it) {
deba@165
   915
	delete it->second;
deba@165
   916
      }
deba@165
   917
deba@165
   918
      for (typename Attributes::iterator it = _attributes.begin(); 
deba@165
   919
	   it != _attributes.end(); ++it) {
deba@165
   920
	delete it->second;
deba@165
   921
      }
deba@165
   922
deba@165
   923
      if (local_os) {
deba@165
   924
	delete _os;
deba@165
   925
      }
deba@165
   926
    }
deba@165
   927
deba@165
   928
  private:
deba@165
   929
    
deba@165
   930
    GraphWriter& operator=(const GraphWriter&);
deba@165
   931
deba@165
   932
  public:
deba@165
   933
deba@165
   934
    /// \name Writing rules
deba@165
   935
    /// @{
deba@165
   936
    
deba@165
   937
    /// \brief Node map reading rule
deba@165
   938
    ///
deba@165
   939
    /// Add a node map reading rule to the writer.
deba@165
   940
    template <typename Map>
deba@165
   941
    GraphWriter& nodeMap(const std::string& caption, const Map& map) {
deba@165
   942
      checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@165
   943
      _writer_bits::MapStorageBase<Node>* storage = 
deba@165
   944
	new _writer_bits::MapStorage<Node, Map>(map);
deba@165
   945
      _node_maps.push_back(std::make_pair(caption, storage));
deba@165
   946
      return *this;
deba@165
   947
    }
deba@165
   948
deba@165
   949
    /// \brief Node map writing rule
deba@165
   950
    ///
deba@165
   951
    /// Add a node map writing rule with specialized converter to the
deba@165
   952
    /// writer.
deba@165
   953
    template <typename Map, typename Converter>
deba@165
   954
    GraphWriter& nodeMap(const std::string& caption, const Map& map, 
deba@165
   955
			   const Converter& converter = Converter()) {
deba@165
   956
      checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@165
   957
      _writer_bits::MapStorageBase<Node>* storage = 
deba@165
   958
	new _writer_bits::MapStorage<Node, Map, Converter>(map, converter);
deba@165
   959
      _node_maps.push_back(std::make_pair(caption, storage));
deba@165
   960
      return *this;
deba@165
   961
    }
deba@165
   962
deba@165
   963
    /// \brief Edge map writing rule
deba@165
   964
    ///
deba@165
   965
    /// Add an edge map writing rule to the writer.
deba@165
   966
    template <typename Map>
deba@165
   967
    GraphWriter& edgeMap(const std::string& caption, const Map& map) {
deba@165
   968
      checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
deba@165
   969
      _writer_bits::MapStorageBase<Edge>* storage = 
deba@165
   970
	new _writer_bits::MapStorage<Edge, Map>(map);
deba@165
   971
      _edge_maps.push_back(std::make_pair(caption, storage));
deba@165
   972
      return *this;
deba@165
   973
    }
deba@165
   974
deba@165
   975
    /// \brief Edge map writing rule
deba@165
   976
    ///
deba@165
   977
    /// Add an edge map writing rule with specialized converter to the
deba@165
   978
    /// writer.
deba@165
   979
    template <typename Map, typename Converter>
deba@165
   980
    GraphWriter& edgeMap(const std::string& caption, const Map& map, 
deba@165
   981
			  const Converter& converter = Converter()) {
deba@165
   982
      checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
deba@165
   983
      _writer_bits::MapStorageBase<Edge>* storage = 
deba@165
   984
	new _writer_bits::MapStorage<Edge, Map, Converter>(map, converter);
deba@165
   985
      _edge_maps.push_back(std::make_pair(caption, storage));
deba@165
   986
      return *this;
deba@165
   987
    }
deba@165
   988
deba@165
   989
    /// \brief Arc map writing rule
deba@165
   990
    ///
deba@165
   991
    /// Add an arc map writing rule to the writer.
deba@165
   992
    template <typename Map>
deba@165
   993
    GraphWriter& arcMap(const std::string& caption, const Map& map) {
deba@165
   994
      checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
deba@165
   995
      _writer_bits::MapStorageBase<Edge>* forward_storage = 
deba@165
   996
	new _writer_bits::GraphArcMapStorage<Graph, true, Map>(_graph, map);
deba@165
   997
      _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
deba@165
   998
      _writer_bits::MapStorageBase<Edge>* backward_storage = 
deba@165
   999
	new _writer_bits::GraphArcMapStorage<Graph, false, Map>(_graph, map);
deba@165
  1000
      _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@165
  1001
      return *this;
deba@165
  1002
    }
deba@165
  1003
deba@165
  1004
    /// \brief Arc map writing rule
deba@165
  1005
    ///
deba@165
  1006
    /// Add an arc map writing rule with specialized converter to the
deba@165
  1007
    /// writer.
deba@165
  1008
    template <typename Map, typename Converter>
deba@165
  1009
    GraphWriter& arcMap(const std::string& caption, const Map& map, 
deba@165
  1010
			  const Converter& converter = Converter()) {
deba@165
  1011
      checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
deba@165
  1012
      _writer_bits::MapStorageBase<Edge>* forward_storage = 
deba@165
  1013
	new _writer_bits::GraphArcMapStorage<Graph, true, Map, Converter>
deba@165
  1014
	(_graph, map, converter);
deba@165
  1015
      _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
deba@165
  1016
      _writer_bits::MapStorageBase<Edge>* backward_storage = 
deba@165
  1017
	new _writer_bits::GraphArcMapStorage<Graph, false, Map, Converter>
deba@165
  1018
	(_graph, map, converter);
deba@165
  1019
      _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@165
  1020
      return *this;
deba@165
  1021
    }
deba@165
  1022
deba@165
  1023
    /// \brief Attribute writing rule
deba@165
  1024
    ///
deba@165
  1025
    /// Add an attribute writing rule to the writer.
deba@165
  1026
    template <typename Value>
deba@165
  1027
    GraphWriter& attribute(const std::string& caption, const Value& value) {
deba@165
  1028
      _writer_bits::ValueStorageBase* storage = 
deba@165
  1029
	new _writer_bits::ValueStorage<Value>(value);
deba@165
  1030
      _attributes.push_back(std::make_pair(caption, storage));
deba@165
  1031
      return *this;
deba@165
  1032
    }
deba@165
  1033
deba@165
  1034
    /// \brief Attribute writing rule
deba@165
  1035
    ///
deba@165
  1036
    /// Add an attribute writing rule with specialized converter to the
deba@165
  1037
    /// writer.
deba@165
  1038
    template <typename Value, typename Converter>
deba@165
  1039
    GraphWriter& attribute(const std::string& caption, const Value& value, 
deba@165
  1040
			     const Converter& converter = Converter()) {
deba@165
  1041
      _writer_bits::ValueStorageBase* storage = 
deba@165
  1042
	new _writer_bits::ValueStorage<Value, Converter>(value, converter);
deba@165
  1043
      _attributes.push_back(std::make_pair(caption, storage));
deba@165
  1044
      return *this;
deba@165
  1045
    }
deba@165
  1046
deba@165
  1047
    /// \brief Node writing rule
deba@165
  1048
    ///
deba@165
  1049
    /// Add a node writing rule to the writer.
deba@165
  1050
    GraphWriter& node(const std::string& caption, const Node& node) {
deba@165
  1051
      typedef _writer_bits::MapLookUpConverter<Node> Converter;
deba@165
  1052
      Converter converter(_node_index);
deba@165
  1053
      _writer_bits::ValueStorageBase* storage = 
deba@165
  1054
	new _writer_bits::ValueStorage<Node, Converter>(node, converter);
deba@165
  1055
      _attributes.push_back(std::make_pair(caption, storage));
deba@165
  1056
      return *this;
deba@165
  1057
    }
deba@165
  1058
deba@165
  1059
    /// \brief Edge writing rule
deba@165
  1060
    ///
deba@165
  1061
    /// Add an edge writing rule to writer.
deba@165
  1062
    GraphWriter& edge(const std::string& caption, const Edge& edge) {
deba@165
  1063
      typedef _writer_bits::MapLookUpConverter<Edge> Converter;
deba@165
  1064
      Converter converter(_edge_index);
deba@165
  1065
      _writer_bits::ValueStorageBase* storage = 
deba@165
  1066
	new _writer_bits::ValueStorage<Edge, Converter>(edge, converter);
deba@165
  1067
      _attributes.push_back(std::make_pair(caption, storage));
deba@165
  1068
      return *this;
deba@165
  1069
    }
deba@165
  1070
deba@165
  1071
    /// \brief Arc writing rule
deba@165
  1072
    ///
deba@165
  1073
    /// Add an arc writing rule to writer.
deba@165
  1074
    GraphWriter& arc(const std::string& caption, const Arc& arc) {
deba@165
  1075
      typedef _writer_bits::GraphArcLookUpConverter<Graph> Converter;
deba@165
  1076
      Converter converter(_graph, _edge_index);
deba@165
  1077
      _writer_bits::ValueStorageBase* storage = 
deba@165
  1078
	new _writer_bits::ValueStorage<Arc, Converter>(arc, converter);
deba@165
  1079
      _attributes.push_back(std::make_pair(caption, storage));
deba@165
  1080
      return *this;
deba@165
  1081
    }
deba@165
  1082
deba@165
  1083
    /// \name Select section by name
deba@165
  1084
    /// @{
deba@165
  1085
deba@165
  1086
    /// \brief Set \c \@nodes section to be read
deba@165
  1087
    ///
deba@165
  1088
    /// Set \c \@nodes section to be read
deba@165
  1089
    GraphWriter& nodes(const std::string& caption) {
deba@165
  1090
      _nodes_caption = caption;
deba@165
  1091
      return *this;
deba@165
  1092
    }
deba@165
  1093
deba@165
  1094
    /// \brief Set \c \@edges section to be read
deba@165
  1095
    ///
deba@165
  1096
    /// Set \c \@edges section to be read
deba@165
  1097
    GraphWriter& edges(const std::string& caption) {
deba@165
  1098
      _edges_caption = caption;
deba@165
  1099
      return *this;
deba@165
  1100
    }
deba@165
  1101
deba@165
  1102
    /// \brief Set \c \@attributes section to be read
deba@165
  1103
    ///
deba@165
  1104
    /// Set \c \@attributes section to be read
deba@165
  1105
    GraphWriter& attributes(const std::string& caption) {
deba@165
  1106
      _attributes_caption = caption;
deba@165
  1107
      return *this;
deba@165
  1108
    }
deba@165
  1109
deba@165
  1110
    /// \name Skipping section
deba@165
  1111
    /// @{
deba@165
  1112
deba@165
  1113
    /// \brief Skip writing the node set
deba@165
  1114
    ///
deba@165
  1115
    /// The \c \@nodes section will be not written to the stream.
deba@165
  1116
    GraphWriter& skipNodes() {
deba@165
  1117
      LEMON_ASSERT(!_skip_nodes, "Multiple usage of skipNodes() member");
deba@165
  1118
      return *this;
deba@165
  1119
    }
deba@165
  1120
deba@165
  1121
    /// \brief Skip writing edge set
deba@165
  1122
    ///
deba@165
  1123
    /// The \c \@edges section will be not written to the stream.
deba@165
  1124
    GraphWriter& skipEdges() {
deba@165
  1125
      LEMON_ASSERT(!_skip_edges, "Multiple usage of skipEdges() member");
deba@165
  1126
      return *this;
deba@165
  1127
    }
deba@165
  1128
deba@165
  1129
    /// @}
deba@165
  1130
deba@165
  1131
  private:
deba@165
  1132
deba@165
  1133
    void writeNodes() {
deba@165
  1134
      _writer_bits::MapStorageBase<Node>* label = 0;
deba@165
  1135
      for (typename NodeMaps::iterator it = _node_maps.begin();
deba@165
  1136
	   it != _node_maps.end(); ++it) {
deba@165
  1137
        if (it->first == "label") {
deba@165
  1138
	  label = it->second;
deba@165
  1139
	  break;
deba@165
  1140
	}
deba@165
  1141
      }
deba@165
  1142
deba@165
  1143
      *_os << "@nodes";
deba@165
  1144
      if (!_nodes_caption.empty()) {
deba@165
  1145
	_writer_bits::writeToken(*_os << ' ', _nodes_caption);
deba@165
  1146
      }
deba@165
  1147
      *_os << std::endl;
deba@165
  1148
deba@165
  1149
      if (label == 0) {
deba@165
  1150
	*_os << "label" << '\t';
deba@165
  1151
      }
deba@165
  1152
      for (typename NodeMaps::iterator it = _node_maps.begin();
deba@165
  1153
	   it != _node_maps.end(); ++it) {
deba@165
  1154
	_writer_bits::writeToken(*_os, it->first) << '\t';
deba@165
  1155
      }
deba@165
  1156
      *_os << std::endl;
deba@165
  1157
deba@165
  1158
      std::vector<Node> nodes;
deba@165
  1159
      for (NodeIt n(_graph); n != INVALID; ++n) {
deba@165
  1160
	nodes.push_back(n);
deba@165
  1161
      }
deba@165
  1162
      
deba@165
  1163
      if (label == 0) {
deba@165
  1164
	IdMap<Graph, Node> id_map(_graph);
deba@165
  1165
	_writer_bits::MapLess<IdMap<Graph, Node> > id_less(id_map);
deba@165
  1166
	std::sort(nodes.begin(), nodes.end(), id_less);
deba@165
  1167
      } else {
deba@165
  1168
	label->sort(nodes);
deba@165
  1169
      }
deba@165
  1170
deba@165
  1171
      for (int i = 0; i < static_cast<int>(nodes.size()); ++i) {
deba@165
  1172
	Node n = nodes[i];
deba@165
  1173
	if (label == 0) {
deba@165
  1174
	  std::ostringstream os;
deba@165
  1175
	  os << _graph.id(n);
deba@165
  1176
	  _writer_bits::writeToken(*_os, os.str());
deba@165
  1177
	  *_os << '\t';
deba@165
  1178
	  _node_index.insert(std::make_pair(n, os.str()));
deba@165
  1179
	}
deba@165
  1180
	for (typename NodeMaps::iterator it = _node_maps.begin();
deba@165
  1181
	     it != _node_maps.end(); ++it) {
deba@165
  1182
	  std::string value = it->second->get(n);
deba@165
  1183
	  _writer_bits::writeToken(*_os, value);
deba@165
  1184
	  if (it->first == "label") {
deba@165
  1185
	    _node_index.insert(std::make_pair(n, value));
deba@165
  1186
	  }
deba@165
  1187
	  *_os << '\t';
deba@165
  1188
	}
deba@165
  1189
	*_os << std::endl;
deba@165
  1190
      }
deba@165
  1191
    }
deba@165
  1192
deba@165
  1193
    void writeEdges() {
deba@165
  1194
      _writer_bits::MapStorageBase<Edge>* label = 0;
deba@165
  1195
      for (typename EdgeMaps::iterator it = _edge_maps.begin();
deba@165
  1196
	   it != _edge_maps.end(); ++it) {
deba@165
  1197
        if (it->first == "label") {
deba@165
  1198
	  label = it->second;
deba@165
  1199
	  break;
deba@165
  1200
	}
deba@165
  1201
      }
deba@165
  1202
deba@165
  1203
      *_os << "@edges";
deba@165
  1204
      if (!_edges_caption.empty()) {
deba@165
  1205
	_writer_bits::writeToken(*_os << ' ', _edges_caption);
deba@165
  1206
      }
deba@165
  1207
      *_os << std::endl;
deba@165
  1208
deba@165
  1209
      *_os << '\t' << '\t';
deba@165
  1210
      if (label == 0) {
deba@165
  1211
	*_os << "label" << '\t';
deba@165
  1212
      }
deba@165
  1213
      for (typename EdgeMaps::iterator it = _edge_maps.begin();
deba@165
  1214
	   it != _edge_maps.end(); ++it) {
deba@165
  1215
	_writer_bits::writeToken(*_os, it->first) << '\t';
deba@165
  1216
      }
deba@165
  1217
      *_os << std::endl;
deba@165
  1218
deba@165
  1219
      std::vector<Edge> edges;
deba@165
  1220
      for (EdgeIt n(_graph); n != INVALID; ++n) {
deba@165
  1221
	edges.push_back(n);
deba@165
  1222
      }
deba@165
  1223
      
deba@165
  1224
      if (label == 0) {
deba@165
  1225
	IdMap<Graph, Edge> id_map(_graph);
deba@165
  1226
	_writer_bits::MapLess<IdMap<Graph, Edge> > id_less(id_map);
deba@165
  1227
	std::sort(edges.begin(), edges.end(), id_less);
deba@165
  1228
      } else {
deba@165
  1229
	label->sort(edges);
deba@165
  1230
      }
deba@165
  1231
deba@165
  1232
      for (int i = 0; i < static_cast<int>(edges.size()); ++i) {
deba@165
  1233
	Edge e = edges[i];
deba@165
  1234
	_writer_bits::writeToken(*_os, _node_index.
deba@165
  1235
				 find(_graph.u(e))->second);
deba@165
  1236
	*_os << '\t';
deba@165
  1237
	_writer_bits::writeToken(*_os, _node_index.
deba@165
  1238
				 find(_graph.v(e))->second);
deba@165
  1239
	*_os << '\t';
deba@165
  1240
	if (label == 0) {
deba@165
  1241
	  std::ostringstream os;
deba@165
  1242
	  os << _graph.id(e);
deba@165
  1243
	  _writer_bits::writeToken(*_os, os.str());
deba@165
  1244
	  *_os << '\t';
deba@165
  1245
	  _edge_index.insert(std::make_pair(e, os.str()));
deba@165
  1246
	}
deba@165
  1247
	for (typename EdgeMaps::iterator it = _edge_maps.begin();
deba@165
  1248
	     it != _edge_maps.end(); ++it) {
deba@165
  1249
	  std::string value = it->second->get(e);
deba@165
  1250
	  _writer_bits::writeToken(*_os, value);
deba@165
  1251
	  if (it->first == "label") {
deba@165
  1252
	    _edge_index.insert(std::make_pair(e, value));
deba@165
  1253
	  }
deba@165
  1254
	  *_os << '\t';
deba@165
  1255
	}
deba@165
  1256
	*_os << std::endl;
deba@165
  1257
      }
deba@165
  1258
    }
deba@165
  1259
deba@165
  1260
    void writeAttributes() {
deba@165
  1261
      if (_attributes.empty()) return;
deba@165
  1262
      *_os << "@attributes";
deba@165
  1263
      if (!_attributes_caption.empty()) {
deba@165
  1264
	_writer_bits::writeToken(*_os << ' ', _attributes_caption);
deba@165
  1265
      }
deba@165
  1266
      *_os << std::endl;
deba@165
  1267
      for (typename Attributes::iterator it = _attributes.begin();
deba@165
  1268
	   it != _attributes.end(); ++it) {
deba@165
  1269
	_writer_bits::writeToken(*_os, it->first) << ' ';
deba@165
  1270
	_writer_bits::writeToken(*_os, it->second->get());
deba@165
  1271
	*_os << std::endl;
deba@165
  1272
      }
deba@165
  1273
    }
deba@165
  1274
    
deba@165
  1275
  public:
deba@165
  1276
    
deba@165
  1277
    /// \name Execution of the writer    
deba@165
  1278
    /// @{
deba@165
  1279
deba@165
  1280
    /// \brief Start the batch processing
deba@165
  1281
    ///
deba@165
  1282
    /// This function starts the batch processing
deba@165
  1283
    void run() {
deba@165
  1284
      if (!_skip_nodes) {
deba@165
  1285
	writeNodes();
deba@165
  1286
      }
deba@165
  1287
      if (!_skip_edges) {      
deba@165
  1288
	writeEdges();
deba@165
  1289
      }
deba@165
  1290
      writeAttributes();
deba@165
  1291
    }
deba@165
  1292
deba@165
  1293
    /// \brief Gives back the stream of the writer
deba@165
  1294
    ///
deba@165
  1295
    /// Gives back the stream of the writer
deba@165
  1296
    std::ostream& ostream() {
deba@165
  1297
      return *_os;
deba@165
  1298
    }
deba@165
  1299
deba@165
  1300
    /// @}
deba@165
  1301
  };
deba@165
  1302
deba@165
  1303
  /// \relates GraphWriter
deba@165
  1304
  template <typename Graph>
deba@165
  1305
  GraphWriter<Graph> graphWriter(std::ostream& os, Graph& graph) {
deba@165
  1306
    GraphWriter<Graph> tmp(os, graph);
deba@165
  1307
    return tmp;
deba@165
  1308
  }
deba@165
  1309
deba@165
  1310
  /// \relates GraphWriter
deba@165
  1311
  template <typename Graph>
deba@165
  1312
  GraphWriter<Graph> graphWriter(const std::string& fn, Graph& graph) {
deba@165
  1313
    GraphWriter<Graph> tmp(fn, graph);
deba@165
  1314
    return tmp;
deba@165
  1315
  }
deba@165
  1316
deba@165
  1317
  /// \relates GraphWriter
deba@165
  1318
  template <typename Graph>
deba@165
  1319
  GraphWriter<Graph> graphWriter(const char* fn, Graph& graph) {
deba@165
  1320
    GraphWriter<Graph> tmp(fn, graph);
deba@165
  1321
    return tmp;
deba@165
  1322
  }
deba@127
  1323
}
deba@127
  1324
deba@127
  1325
#endif