alpar@209: /* -*- mode: C++; indent-tabs-mode: nil; -*-
deba@127:  *
alpar@209:  * This file is a part of LEMON, a generic C++ optimization library.
deba@127:  *
alpar@949:  * Copyright (C) 2003-2011
deba@127:  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
deba@127:  * (Egervary Research Group on Combinatorial Optimization, EGRES).
deba@127:  *
deba@127:  * Permission to use, modify and distribute this software is granted
deba@127:  * provided that this copyright notice appears in all copies. For
deba@127:  * precise terms see the accompanying LICENSE file.
deba@127:  *
deba@127:  * This software is provided "AS IS" with no warranty of any kind,
deba@127:  * express or implied, and with no claim as to its suitability for any
deba@127:  * purpose.
deba@127:  *
deba@127:  */
deba@127: 
deba@127: ///\ingroup lemon_io
deba@127: ///\file
ladanyi@236: ///\brief \ref lgf-format "LEMON Graph Format" reader.
deba@127: 
deba@127: 
deba@127: #ifndef LEMON_LGF_READER_H
deba@127: #define LEMON_LGF_READER_H
deba@127: 
deba@127: #include <iostream>
deba@127: #include <fstream>
deba@127: #include <sstream>
deba@127: 
deba@127: #include <set>
deba@127: #include <map>
deba@127: 
deba@220: #include <lemon/core.h>
deba@127: 
deba@127: #include <lemon/lgf_writer.h>
deba@127: 
deba@127: #include <lemon/concept_check.h>
deba@127: #include <lemon/concepts/maps.h>
deba@127: 
deba@127: namespace lemon {
deba@127: 
deba@127:   namespace _reader_bits {
deba@127: 
deba@127:     template <typename Value>
deba@127:     struct DefaultConverter {
deba@127:       Value operator()(const std::string& str) {
alpar@209:         std::istringstream is(str);
alpar@209:         Value value;
deba@290:         if (!(is >> value)) {
deba@290:           throw FormatError("Cannot read token");
deba@290:         }
alpar@209: 
alpar@209:         char c;
alpar@209:         if (is >> std::ws >> c) {
deba@290:           throw FormatError("Remaining characters in token");
alpar@209:         }
alpar@209:         return value;
deba@127:       }
deba@127:     };
deba@127: 
deba@127:     template <>
deba@127:     struct DefaultConverter<std::string> {
deba@127:       std::string operator()(const std::string& str) {
alpar@209:         return str;
deba@127:       }
deba@127:     };
deba@127: 
alpar@209:     template <typename _Item>
deba@127:     class MapStorageBase {
deba@127:     public:
deba@127:       typedef _Item Item;
deba@127: 
deba@127:     public:
deba@127:       MapStorageBase() {}
deba@127:       virtual ~MapStorageBase() {}
deba@127: 
deba@127:       virtual void set(const Item& item, const std::string& value) = 0;
deba@127: 
deba@127:     };
deba@127: 
alpar@209:     template <typename _Item, typename _Map,
alpar@209:               typename _Converter = DefaultConverter<typename _Map::Value> >
deba@127:     class MapStorage : public MapStorageBase<_Item> {
deba@127:     public:
deba@127:       typedef _Map Map;
deba@127:       typedef _Converter Converter;
deba@127:       typedef _Item Item;
alpar@209: 
deba@127:     private:
deba@127:       Map& _map;
deba@127:       Converter _converter;
deba@127: 
deba@127:     public:
alpar@209:       MapStorage(Map& map, const Converter& converter = Converter())
alpar@209:         : _map(map), _converter(converter) {}
deba@127:       virtual ~MapStorage() {}
deba@127: 
deba@127:       virtual void set(const Item& item ,const std::string& value) {
alpar@209:         _map.set(item, _converter(value));
deba@127:       }
deba@127:     };
deba@127: 
deba@598:     template <typename _GR, bool _dir, typename _Map,
alpar@209:               typename _Converter = DefaultConverter<typename _Map::Value> >
deba@598:     class GraphArcMapStorage : public MapStorageBase<typename _GR::Edge> {
deba@165:     public:
deba@165:       typedef _Map Map;
deba@165:       typedef _Converter Converter;
deba@598:       typedef _GR GR;
deba@598:       typedef typename GR::Edge Item;
deba@165:       static const bool dir = _dir;
alpar@209: 
deba@165:     private:
deba@598:       const GR& _graph;
deba@165:       Map& _map;
deba@165:       Converter _converter;
deba@165: 
deba@165:     public:
deba@598:       GraphArcMapStorage(const GR& graph, Map& map,
alpar@209:                          const Converter& converter = Converter())
alpar@209:         : _graph(graph), _map(map), _converter(converter) {}
deba@165:       virtual ~GraphArcMapStorage() {}
deba@165: 
deba@165:       virtual void set(const Item& item ,const std::string& value) {
alpar@209:         _map.set(_graph.direct(item, dir), _converter(value));
deba@165:       }
deba@165:     };
deba@165: 
deba@127:     class ValueStorageBase {
deba@127:     public:
deba@127:       ValueStorageBase() {}
deba@127:       virtual ~ValueStorageBase() {}
deba@127: 
deba@127:       virtual void set(const std::string&) = 0;
deba@127:     };
deba@127: 
deba@127:     template <typename _Value, typename _Converter = DefaultConverter<_Value> >
deba@127:     class ValueStorage : public ValueStorageBase {
deba@127:     public:
deba@127:       typedef _Value Value;
deba@127:       typedef _Converter Converter;
deba@127: 
deba@127:     private:
deba@127:       Value& _value;
deba@127:       Converter _converter;
deba@127: 
deba@127:     public:
deba@127:       ValueStorage(Value& value, const Converter& converter = Converter())
kpeter@212:         : _value(value), _converter(converter) {}
deba@127: 
deba@127:       virtual void set(const std::string& value) {
alpar@209:         _value = _converter(value);
deba@127:       }
deba@127:     };
deba@127: 
deba@1030:     template <typename Value,
deba@1030:               typename Map = std::map<std::string, Value> >
deba@127:     struct MapLookUpConverter {
deba@1030:       const Map& _map;
deba@1030: 
deba@1030:       MapLookUpConverter(const Map& map)
deba@127:         : _map(map) {}
deba@127: 
deba@127:       Value operator()(const std::string& str) {
deba@1030:         typename Map::const_iterator it = _map.find(str);
deba@127:         if (it == _map.end()) {
deba@127:           std::ostringstream msg;
deba@127:           msg << "Item not found: " << str;
deba@290:           throw FormatError(msg.str());
deba@127:         }
deba@127:         return it->second;
deba@127:       }
deba@127:     };
deba@127: 
deba@1030:     template <typename Value,
deba@1030:               typename Map1 = std::map<std::string, Value>,
deba@1030:               typename Map2 = std::map<std::string, Value> >
deba@1030:     struct DoubleMapLookUpConverter {
deba@1030:       const Map1& _map1;
deba@1030:       const Map2& _map2;
deba@1030: 
deba@1030:       DoubleMapLookUpConverter(const Map1& map1, const Map2& map2)
deba@1030:         : _map1(map1), _map2(map2) {}
deba@1030: 
deba@1030:       Value operator()(const std::string& str) {
deba@1030:         typename Map1::const_iterator it1 = _map1.find(str);
deba@1030:         typename Map2::const_iterator it2 = _map2.find(str);
deba@1030:         if (it1 == _map1.end()) {
deba@1030:           if (it2 == _map2.end()) {
deba@1030:             std::ostringstream msg;
deba@1030:             msg << "Item not found: " << str;
deba@1030:             throw FormatError(msg.str());
deba@1030:           } else {
deba@1030:             return it2->second;
deba@1030:           }
deba@1030:         } else {
deba@1030:           if (it2 == _map2.end()) {
deba@1030:             return it1->second;
deba@1030:           } else {
deba@1030:             std::ostringstream msg;
deba@1030:             msg << "Item is ambigous: " << str;
deba@1030:             throw FormatError(msg.str());
deba@1030:           }
deba@1030:         }
deba@1030:       }
deba@1030:     };
deba@1030: 
deba@598:     template <typename GR>
deba@165:     struct GraphArcLookUpConverter {
deba@598:       const GR& _graph;
deba@598:       const std::map<std::string, typename GR::Edge>& _map;
deba@598: 
deba@598:       GraphArcLookUpConverter(const GR& graph,
alpar@209:                               const std::map<std::string,
deba@598:                                              typename GR::Edge>& map)
alpar@209:         : _graph(graph), _map(map) {}
alpar@209: 
deba@598:       typename GR::Arc operator()(const std::string& str) {
alpar@209:         if (str.empty() || (str[0] != '+' && str[0] != '-')) {
deba@290:           throw FormatError("Item must start with '+' or '-'");
alpar@209:         }
deba@598:         typename std::map<std::string, typename GR::Edge>
alpar@209:           ::const_iterator it = _map.find(str.substr(1));
alpar@209:         if (it == _map.end()) {
deba@290:           throw FormatError("Item not found");
alpar@209:         }
alpar@209:         return _graph.direct(it->second, str[0] == '+');
deba@165:       }
deba@165:     };
deba@165: 
deba@197:     inline bool isWhiteSpace(char c) {
alpar@209:       return c == ' ' || c == '\t' || c == '\v' ||
alpar@209:         c == '\n' || c == '\r' || c == '\f';
deba@127:     }
alpar@209: 
deba@197:     inline bool isOct(char c) {
alpar@209:       return '0' <= c && c <='7';
deba@127:     }
alpar@209: 
deba@197:     inline int valueOct(char c) {
deba@127:       LEMON_ASSERT(isOct(c), "The character is not octal.");
deba@127:       return c - '0';
deba@127:     }
deba@127: 
deba@197:     inline bool isHex(char c) {
alpar@209:       return ('0' <= c && c <= '9') ||
alpar@209:         ('a' <= c && c <= 'z') ||
alpar@209:         ('A' <= c && c <= 'Z');
deba@127:     }
alpar@209: 
deba@197:     inline int valueHex(char c) {
deba@127:       LEMON_ASSERT(isHex(c), "The character is not hexadecimal.");
deba@127:       if ('0' <= c && c <= '9') return c - '0';
deba@127:       if ('a' <= c && c <= 'z') return c - 'a' + 10;
deba@127:       return c - 'A' + 10;
deba@127:     }
deba@127: 
deba@197:     inline bool isIdentifierFirstChar(char c) {
deba@127:       return ('a' <= c && c <= 'z') ||
alpar@209:         ('A' <= c && c <= 'Z') || c == '_';
deba@127:     }
deba@127: 
deba@197:     inline bool isIdentifierChar(char c) {
deba@127:       return isIdentifierFirstChar(c) ||
alpar@209:         ('0' <= c && c <= '9');
deba@127:     }
deba@127: 
deba@197:     inline char readEscape(std::istream& is) {
deba@127:       char c;
deba@127:       if (!is.get(c))
deba@290:         throw FormatError("Escape format error");
deba@127: 
deba@127:       switch (c) {
deba@127:       case '\\':
alpar@209:         return '\\';
deba@127:       case '\"':
alpar@209:         return '\"';
deba@127:       case '\'':
alpar@209:         return '\'';
deba@127:       case '\?':
alpar@209:         return '\?';
deba@127:       case 'a':
alpar@209:         return '\a';
deba@127:       case 'b':
alpar@209:         return '\b';
deba@127:       case 'f':
alpar@209:         return '\f';
deba@127:       case 'n':
alpar@209:         return '\n';
deba@127:       case 'r':
alpar@209:         return '\r';
deba@127:       case 't':
alpar@209:         return '\t';
deba@127:       case 'v':
alpar@209:         return '\v';
deba@127:       case 'x':
alpar@209:         {
alpar@209:           int code;
alpar@209:           if (!is.get(c) || !isHex(c))
deba@290:             throw FormatError("Escape format error");
alpar@209:           else if (code = valueHex(c), !is.get(c) || !isHex(c)) is.putback(c);
alpar@209:           else code = code * 16 + valueHex(c);
alpar@209:           return code;
alpar@209:         }
deba@127:       default:
alpar@209:         {
alpar@209:           int code;
alpar@209:           if (!isOct(c))
deba@290:             throw FormatError("Escape format error");
alpar@209:           else if (code = valueOct(c), !is.get(c) || !isOct(c))
alpar@209:             is.putback(c);
alpar@209:           else if (code = code * 8 + valueOct(c), !is.get(c) || !isOct(c))
alpar@209:             is.putback(c);
alpar@209:           else code = code * 8 + valueOct(c);
alpar@209:           return code;
alpar@209:         }
alpar@209:       }
deba@127:     }
alpar@209: 
deba@197:     inline std::istream& readToken(std::istream& is, std::string& str) {
deba@127:       std::ostringstream os;
deba@127: 
deba@127:       char c;
deba@127:       is >> std::ws;
alpar@209: 
alpar@209:       if (!is.get(c))
alpar@209:         return is;
deba@127: 
deba@127:       if (c == '\"') {
alpar@209:         while (is.get(c) && c != '\"') {
alpar@209:           if (c == '\\')
alpar@209:             c = readEscape(is);
alpar@209:           os << c;
alpar@209:         }
alpar@209:         if (!is)
deba@290:           throw FormatError("Quoted format error");
deba@127:       } else {
alpar@209:         is.putback(c);
alpar@209:         while (is.get(c) && !isWhiteSpace(c)) {
alpar@209:           if (c == '\\')
alpar@209:             c = readEscape(is);
alpar@209:           os << c;
alpar@209:         }
alpar@209:         if (!is) {
alpar@209:           is.clear();
alpar@209:         } else {
alpar@209:           is.putback(c);
alpar@209:         }
deba@127:       }
deba@127:       str = os.str();
deba@127:       return is;
deba@127:     }
deba@162: 
deba@162:     class Section {
deba@162:     public:
deba@162:       virtual ~Section() {}
deba@162:       virtual void process(std::istream& is, int& line_num) = 0;
deba@162:     };
deba@162: 
deba@162:     template <typename Functor>
deba@162:     class LineSection : public Section {
deba@162:     private:
deba@162: 
deba@162:       Functor _functor;
deba@162: 
deba@162:     public:
alpar@209: 
deba@162:       LineSection(const Functor& functor) : _functor(functor) {}
deba@162:       virtual ~LineSection() {}
deba@162: 
deba@162:       virtual void process(std::istream& is, int& line_num) {
alpar@209:         char c;
alpar@209:         std::string line;
alpar@209:         while (is.get(c) && c != '@') {
alpar@209:           if (c == '\n') {
alpar@209:             ++line_num;
alpar@209:           } else if (c == '#') {
alpar@209:             getline(is, line);
alpar@209:             ++line_num;
alpar@209:           } else if (!isWhiteSpace(c)) {
alpar@209:             is.putback(c);
alpar@209:             getline(is, line);
alpar@209:             _functor(line);
alpar@209:             ++line_num;
alpar@209:           }
alpar@209:         }
alpar@209:         if (is) is.putback(c);
alpar@209:         else if (is.eof()) is.clear();
deba@162:       }
deba@162:     };
deba@162: 
deba@162:     template <typename Functor>
deba@162:     class StreamSection : public Section {
deba@162:     private:
deba@162: 
deba@162:       Functor _functor;
deba@162: 
deba@162:     public:
alpar@209: 
deba@162:       StreamSection(const Functor& functor) : _functor(functor) {}
alpar@209:       virtual ~StreamSection() {}
deba@162: 
deba@162:       virtual void process(std::istream& is, int& line_num) {
alpar@209:         _functor(is, line_num);
alpar@209:         char c;
alpar@209:         std::string line;
alpar@209:         while (is.get(c) && c != '@') {
alpar@209:           if (c == '\n') {
alpar@209:             ++line_num;
alpar@209:           } else if (!isWhiteSpace(c)) {
alpar@209:             getline(is, line);
alpar@209:             ++line_num;
alpar@209:           }
alpar@209:         }
alpar@209:         if (is) is.putback(c);
alpar@209:         else if (is.eof()) is.clear();
deba@162:       }
deba@162:     };
alpar@209: 
deba@127:   }
alpar@156: 
deba@598:   template <typename DGR>
deba@190:   class DigraphReader;
deba@190: 
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, std::istream& is = std::cin);
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, const std::string& fn);
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, const char *fn);
deba@190: 
alpar@156:   /// \ingroup lemon_io
alpar@209:   ///
kpeter@192:   /// \brief \ref lgf-format "LGF" reader for directed graphs
alpar@156:   ///
alpar@156:   /// This utility reads an \ref lgf-format "LGF" file.
alpar@156:   ///
alpar@156:   /// The reading method does a batch processing. The user creates a
alpar@156:   /// reader object, then various reading rules can be added to the
alpar@156:   /// reader, and eventually the reading is executed with the \c run()
alpar@156:   /// member function. A map reading rule can be added to the reader
alpar@156:   /// with the \c nodeMap() or \c arcMap() members. An optional
deba@162:   /// converter parameter can also be added as a standard functor
kpeter@192:   /// converting from \c std::string to the value type of the map. If it
deba@162:   /// is set, it will determine how the tokens in the file should be
kpeter@192:   /// converted to the value type of the map. If the functor is not set,
deba@162:   /// then a default conversion will be used. One map can be read into
deba@162:   /// multiple map objects at the same time. The \c attribute(), \c
deba@162:   /// node() and \c arc() functions are used to add attribute reading
deba@162:   /// rules.
alpar@156:   ///
alpar@156:   ///\code
deba@598:   /// DigraphReader<DGR>(digraph, std::cin).
kpeter@192:   ///   nodeMap("coordinates", coord_map).
kpeter@192:   ///   arcMap("capacity", cap_map).
kpeter@192:   ///   node("source", src).
kpeter@192:   ///   node("target", trg).
kpeter@192:   ///   attribute("caption", caption).
kpeter@192:   ///   run();
alpar@156:   ///\endcode
alpar@156:   ///
kpeter@786:   /// By default, the reader uses the first section in the file of the
alpar@156:   /// proper type. If a section has an optional name, then it can be
deba@162:   /// selected for reading by giving an optional name parameter to the
deba@189:   /// \c nodes(), \c arcs() or \c attributes() functions.
alpar@156:   ///
alpar@156:   /// The \c useNodes() and \c useArcs() functions are used to tell the reader
alpar@156:   /// that the nodes or arcs should not be constructed (added to the
alpar@156:   /// graph) during the reading, but instead the label map of the items
alpar@156:   /// are given as a parameter of these functions. An
kpeter@192:   /// application of these functions is multipass reading, which is
kpeter@192:   /// important if two \c \@arcs sections must be read from the
kpeter@192:   /// file. In this case the first phase would read the node set and one
alpar@156:   /// of the arc sets, while the second phase would read the second arc
alpar@156:   /// set into an \e ArcSet class (\c SmartArcSet or \c ListArcSet).
alpar@156:   /// The previously read label node map should be passed to the \c
alpar@156:   /// useNodes() functions. Another application of multipass reading when
alpar@210:   /// paths are given as a node map or an arc map.
alpar@210:   /// It is impossible to read this in
alpar@156:   /// a single pass, because the arcs are not constructed when the node
alpar@156:   /// maps are read.
deba@598:   template <typename DGR>
deba@127:   class DigraphReader {
deba@127:   public:
deba@127: 
deba@598:     typedef DGR Digraph;
kpeter@559: 
kpeter@559:   private:
kpeter@559: 
deba@598:     TEMPLATE_DIGRAPH_TYPEDEFS(DGR);
alpar@209: 
deba@127:     std::istream* _is;
deba@127:     bool local_is;
deba@290:     std::string _filename;
deba@127: 
deba@598:     DGR& _digraph;
deba@127: 
deba@127:     std::string _nodes_caption;
deba@127:     std::string _arcs_caption;
deba@127:     std::string _attributes_caption;
deba@127: 
deba@127:     typedef std::map<std::string, Node> NodeIndex;
deba@127:     NodeIndex _node_index;
deba@127:     typedef std::map<std::string, Arc> ArcIndex;
deba@127:     ArcIndex _arc_index;
alpar@209: 
alpar@209:     typedef std::vector<std::pair<std::string,
alpar@209:       _reader_bits::MapStorageBase<Node>*> > NodeMaps;
alpar@209:     NodeMaps _node_maps;
deba@127: 
deba@127:     typedef std::vector<std::pair<std::string,
deba@127:       _reader_bits::MapStorageBase<Arc>*> >ArcMaps;
deba@127:     ArcMaps _arc_maps;
deba@127: 
alpar@209:     typedef std::multimap<std::string, _reader_bits::ValueStorageBase*>
deba@127:       Attributes;
deba@127:     Attributes _attributes;
deba@127: 
deba@127:     bool _use_nodes;
deba@127:     bool _use_arcs;
deba@127: 
deba@188:     bool _skip_nodes;
deba@188:     bool _skip_arcs;
deba@188: 
deba@127:     int line_num;
deba@127:     std::istringstream line;
deba@127: 
deba@127:   public:
deba@127: 
alpar@156:     /// \brief Constructor
alpar@156:     ///
alpar@156:     /// Construct a directed graph reader, which reads from the given
alpar@156:     /// input stream.
deba@598:     DigraphReader(DGR& digraph, std::istream& is = std::cin)
deba@127:       : _is(&is), local_is(false), _digraph(digraph),
alpar@209:         _use_nodes(false), _use_arcs(false),
alpar@209:         _skip_nodes(false), _skip_arcs(false) {}
deba@127: 
alpar@156:     /// \brief Constructor
alpar@156:     ///
alpar@156:     /// Construct a directed graph reader, which reads from the given
alpar@156:     /// file.
deba@598:     DigraphReader(DGR& digraph, const std::string& fn)
deba@290:       : _is(new std::ifstream(fn.c_str())), local_is(true),
deba@290:         _filename(fn), _digraph(digraph),
kpeter@212:         _use_nodes(false), _use_arcs(false),
deba@290:         _skip_nodes(false), _skip_arcs(false) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
alpar@209: 
alpar@156:     /// \brief Constructor
alpar@156:     ///
alpar@156:     /// Construct a directed graph reader, which reads from the given
alpar@156:     /// file.
deba@598:     DigraphReader(DGR& digraph, const char* fn)
deba@290:       : _is(new std::ifstream(fn)), local_is(true),
deba@290:         _filename(fn), _digraph(digraph),
kpeter@212:         _use_nodes(false), _use_arcs(false),
deba@290:         _skip_nodes(false), _skip_arcs(false) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
deba@127: 
alpar@156:     /// \brief Destructor
deba@127:     ~DigraphReader() {
alpar@209:       for (typename NodeMaps::iterator it = _node_maps.begin();
alpar@209:            it != _node_maps.end(); ++it) {
alpar@209:         delete it->second;
deba@127:       }
deba@127: 
alpar@209:       for (typename ArcMaps::iterator it = _arc_maps.begin();
alpar@209:            it != _arc_maps.end(); ++it) {
alpar@209:         delete it->second;
deba@127:       }
deba@127: 
alpar@209:       for (typename Attributes::iterator it = _attributes.begin();
alpar@209:            it != _attributes.end(); ++it) {
alpar@209:         delete it->second;
deba@127:       }
deba@127: 
deba@127:       if (local_is) {
alpar@209:         delete _is;
deba@127:       }
deba@127: 
deba@127:     }
deba@127: 
deba@127:   private:
deba@190: 
deba@598:     template <typename TDGR>
deba@598:     friend DigraphReader<TDGR> digraphReader(TDGR& digraph, std::istream& is);
deba@598:     template <typename TDGR>
alpar@877:     friend DigraphReader<TDGR> digraphReader(TDGR& digraph,
deba@598:                                              const std::string& fn);
deba@598:     template <typename TDGR>
deba@598:     friend DigraphReader<TDGR> digraphReader(TDGR& digraph, const char *fn);
alpar@209: 
alpar@209:     DigraphReader(DigraphReader& other)
deba@190:       : _is(other._is), local_is(other.local_is), _digraph(other._digraph),
alpar@209:         _use_nodes(other._use_nodes), _use_arcs(other._use_arcs),
alpar@209:         _skip_nodes(other._skip_nodes), _skip_arcs(other._skip_arcs) {
deba@190: 
deba@190:       other._is = 0;
deba@190:       other.local_is = false;
alpar@209: 
deba@190:       _node_index.swap(other._node_index);
deba@190:       _arc_index.swap(other._arc_index);
deba@190: 
deba@190:       _node_maps.swap(other._node_maps);
deba@190:       _arc_maps.swap(other._arc_maps);
deba@190:       _attributes.swap(other._attributes);
deba@190: 
deba@190:       _nodes_caption = other._nodes_caption;
deba@190:       _arcs_caption = other._arcs_caption;
deba@190:       _attributes_caption = other._attributes_caption;
deba@190: 
deba@190:     }
deba@190: 
deba@127:     DigraphReader& operator=(const DigraphReader&);
deba@127: 
deba@127:   public:
deba@127: 
kpeter@584:     /// \name Reading Rules
alpar@156:     /// @{
alpar@209: 
alpar@156:     /// \brief Node map reading rule
alpar@156:     ///
alpar@156:     /// Add a node map reading rule to the reader.
deba@127:     template <typename Map>
deba@127:     DigraphReader& nodeMap(const std::string& caption, Map& map) {
deba@127:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Node>* storage =
alpar@209:         new _reader_bits::MapStorage<Node, Map>(map);
deba@127:       _node_maps.push_back(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Node map reading rule
alpar@156:     ///
alpar@156:     /// Add a node map reading rule with specialized converter to the
alpar@156:     /// reader.
deba@127:     template <typename Map, typename Converter>
alpar@209:     DigraphReader& nodeMap(const std::string& caption, Map& map,
alpar@209:                            const Converter& converter = Converter()) {
deba@127:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Node>* storage =
alpar@209:         new _reader_bits::MapStorage<Node, Map, Converter>(map, converter);
deba@127:       _node_maps.push_back(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Arc map reading rule
alpar@156:     ///
alpar@156:     /// Add an arc map reading rule to the reader.
deba@127:     template <typename Map>
deba@127:     DigraphReader& arcMap(const std::string& caption, Map& map) {
deba@127:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Arc>* storage =
alpar@209:         new _reader_bits::MapStorage<Arc, Map>(map);
deba@127:       _arc_maps.push_back(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Arc map reading rule
alpar@156:     ///
alpar@156:     /// Add an arc map reading rule with specialized converter to the
alpar@156:     /// reader.
deba@127:     template <typename Map, typename Converter>
alpar@209:     DigraphReader& arcMap(const std::string& caption, Map& map,
alpar@209:                           const Converter& converter = Converter()) {
deba@127:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Arc>* storage =
alpar@209:         new _reader_bits::MapStorage<Arc, Map, Converter>(map, converter);
deba@127:       _arc_maps.push_back(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Attribute reading rule
alpar@156:     ///
alpar@156:     /// Add an attribute reading rule to the reader.
deba@127:     template <typename Value>
deba@127:     DigraphReader& attribute(const std::string& caption, Value& value) {
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Value>(value);
deba@127:       _attributes.insert(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Attribute reading rule
alpar@156:     ///
alpar@156:     /// Add an attribute reading rule with specialized converter to the
alpar@156:     /// reader.
deba@127:     template <typename Value, typename Converter>
alpar@209:     DigraphReader& attribute(const std::string& caption, Value& value,
alpar@209:                              const Converter& converter = Converter()) {
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Value, Converter>(value, converter);
deba@127:       _attributes.insert(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Node reading rule
alpar@156:     ///
alpar@156:     /// Add a node reading rule to reader.
deba@127:     DigraphReader& node(const std::string& caption, Node& node) {
deba@127:       typedef _reader_bits::MapLookUpConverter<Node> Converter;
deba@127:       Converter converter(_node_index);
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Node, Converter>(node, converter);
deba@127:       _attributes.insert(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Arc reading rule
alpar@156:     ///
alpar@156:     /// Add an arc reading rule to reader.
deba@127:     DigraphReader& arc(const std::string& caption, Arc& arc) {
deba@127:       typedef _reader_bits::MapLookUpConverter<Arc> Converter;
deba@127:       Converter converter(_arc_index);
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Arc, Converter>(arc, converter);
deba@127:       _attributes.insert(std::make_pair(caption, storage));
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// @}
alpar@156: 
kpeter@584:     /// \name Select Section by Name
alpar@156:     /// @{
alpar@156: 
alpar@156:     /// \brief Set \c \@nodes section to be read
alpar@156:     ///
alpar@156:     /// Set \c \@nodes section to be read
deba@127:     DigraphReader& nodes(const std::string& caption) {
deba@127:       _nodes_caption = caption;
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Set \c \@arcs section to be read
alpar@156:     ///
alpar@156:     /// Set \c \@arcs section to be read
deba@127:     DigraphReader& arcs(const std::string& caption) {
deba@127:       _arcs_caption = caption;
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Set \c \@attributes section to be read
alpar@156:     ///
alpar@156:     /// Set \c \@attributes section to be read
deba@127:     DigraphReader& attributes(const std::string& caption) {
deba@127:       _attributes_caption = caption;
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// @}
alpar@156: 
kpeter@584:     /// \name Using Previously Constructed Node or Arc Set
alpar@156:     /// @{
alpar@156: 
alpar@156:     /// \brief Use previously constructed node set
alpar@156:     ///
alpar@156:     /// Use previously constructed node set, and specify the node
alpar@156:     /// label map.
deba@127:     template <typename Map>
deba@127:     DigraphReader& useNodes(const Map& map) {
deba@127:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@127:       _use_nodes = true;
deba@127:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@127:       for (NodeIt n(_digraph); n != INVALID; ++n) {
alpar@209:         _node_index.insert(std::make_pair(converter(map[n]), n));
deba@127:       }
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Use previously constructed node set
alpar@156:     ///
alpar@156:     /// Use previously constructed node set, and specify the node
alpar@156:     /// label map and a functor which converts the label map values to
kpeter@192:     /// \c std::string.
deba@127:     template <typename Map, typename Converter>
alpar@209:     DigraphReader& useNodes(const Map& map,
alpar@209:                             const Converter& converter = Converter()) {
deba@127:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@127:       _use_nodes = true;
deba@127:       for (NodeIt n(_digraph); n != INVALID; ++n) {
alpar@209:         _node_index.insert(std::make_pair(converter(map[n]), n));
deba@127:       }
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Use previously constructed arc set
alpar@156:     ///
alpar@156:     /// Use previously constructed arc set, and specify the arc
alpar@156:     /// label map.
deba@127:     template <typename Map>
deba@127:     DigraphReader& useArcs(const Map& map) {
deba@127:       checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
deba@127:       LEMON_ASSERT(!_use_arcs, "Multiple usage of useArcs() member");
deba@127:       _use_arcs = true;
deba@127:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@127:       for (ArcIt a(_digraph); a != INVALID; ++a) {
alpar@209:         _arc_index.insert(std::make_pair(converter(map[a]), a));
deba@127:       }
deba@127:       return *this;
deba@127:     }
deba@127: 
alpar@156:     /// \brief Use previously constructed arc set
alpar@156:     ///
alpar@156:     /// Use previously constructed arc set, and specify the arc
alpar@156:     /// label map and a functor which converts the label map values to
kpeter@192:     /// \c std::string.
deba@127:     template <typename Map, typename Converter>
alpar@209:     DigraphReader& useArcs(const Map& map,
alpar@209:                            const Converter& converter = Converter()) {
deba@127:       checkConcept<concepts::ReadMap<Arc, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_arcs, "Multiple usage of useArcs() member");
deba@127:       _use_arcs = true;
deba@127:       for (ArcIt a(_digraph); a != INVALID; ++a) {
alpar@209:         _arc_index.insert(std::make_pair(converter(map[a]), a));
deba@127:       }
deba@127:       return *this;
deba@127:     }
deba@127: 
deba@188:     /// \brief Skips the reading of node section
deba@188:     ///
deba@188:     /// Omit the reading of the node section. This implies that each node
kpeter@192:     /// map reading rule will be abandoned, and the nodes of the graph
deba@188:     /// will not be constructed, which usually cause that the arc set
kpeter@192:     /// could not be read due to lack of node name resolving.
kpeter@192:     /// Therefore \c skipArcs() function should also be used, or
kpeter@192:     /// \c useNodes() should be used to specify the label of the nodes.
deba@188:     DigraphReader& skipNodes() {
alpar@209:       LEMON_ASSERT(!_skip_nodes, "Skip nodes already set");
deba@188:       _skip_nodes = true;
deba@188:       return *this;
deba@188:     }
deba@188: 
deba@188:     /// \brief Skips the reading of arc section
deba@188:     ///
deba@188:     /// Omit the reading of the arc section. This implies that each arc
kpeter@192:     /// map reading rule will be abandoned, and the arcs of the graph
deba@188:     /// will not be constructed.
deba@188:     DigraphReader& skipArcs() {
alpar@209:       LEMON_ASSERT(!_skip_arcs, "Skip arcs already set");
deba@188:       _skip_arcs = true;
deba@188:       return *this;
deba@188:     }
deba@188: 
alpar@156:     /// @}
alpar@156: 
deba@127:   private:
deba@127: 
deba@127:     bool readLine() {
deba@127:       std::string str;
deba@127:       while(++line_num, std::getline(*_is, str)) {
alpar@209:         line.clear(); line.str(str);
alpar@209:         char c;
alpar@209:         if (line >> std::ws >> c && c != '#') {
alpar@209:           line.putback(c);
alpar@209:           return true;
alpar@209:         }
deba@127:       }
deba@127:       return false;
deba@127:     }
deba@127: 
deba@127:     bool readSuccess() {
deba@127:       return static_cast<bool>(*_is);
deba@127:     }
alpar@209: 
deba@127:     void skipSection() {
deba@127:       char c;
deba@127:       while (readSuccess() && line >> c && c != '@') {
alpar@209:         readLine();
deba@127:       }
deba@427:       if (readSuccess()) {
deba@427:         line.putback(c);
deba@427:       }
deba@127:     }
deba@127: 
deba@127:     void readNodes() {
deba@127: 
deba@127:       std::vector<int> map_index(_node_maps.size());
deba@127:       int map_num, label_index;
deba@127: 
deba@186:       char c;
deba@186:       if (!readLine() || !(line >> c) || c == '@') {
alpar@209:         if (readSuccess() && line) line.putback(c);
alpar@209:         if (!_node_maps.empty())
deba@290:           throw FormatError("Cannot find map names");
alpar@209:         return;
deba@186:       }
deba@186:       line.putback(c);
deba@186: 
deba@127:       {
alpar@209:         std::map<std::string, int> maps;
alpar@209: 
alpar@209:         std::string map;
alpar@209:         int index = 0;
alpar@209:         while (_reader_bits::readToken(line, map)) {
alpar@209:           if (maps.find(map) != maps.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Multiple occurence of node map: " << map;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           maps.insert(std::make_pair(map, index));
alpar@209:           ++index;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {
alpar@209:           std::map<std::string, int>::iterator jt =
alpar@209:             maps.find(_node_maps[i].first);
alpar@209:           if (jt == maps.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Map not found: " << _node_maps[i].first;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           map_index[i] = jt->second;
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           std::map<std::string, int>::iterator jt = maps.find("label");
alpar@209:           if (jt != maps.end()) {
alpar@209:             label_index = jt->second;
alpar@209:           } else {
alpar@209:             label_index = -1;
alpar@209:           }
alpar@209:         }
alpar@209:         map_num = maps.size();
deba@127:       }
deba@127: 
deba@127:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::vector<std::string> tokens(map_num);
alpar@209:         for (int i = 0; i < map_num; ++i) {
alpar@209:           if (!_reader_bits::readToken(line, tokens[i])) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Column not found (" << i + 1 << ")";
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:         }
alpar@209:         if (line >> std::ws >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         Node n;
alpar@209:         if (!_use_nodes) {
alpar@209:           n = _digraph.addNode();
alpar@209:           if (label_index != -1)
alpar@209:             _node_index.insert(std::make_pair(tokens[label_index], n));
alpar@209:         } else {
alpar@209:           if (label_index == -1)
kpeter@291:             throw FormatError("Label map not found");
alpar@209:           typename std::map<std::string, Node>::iterator it =
alpar@209:             _node_index.find(tokens[label_index]);
alpar@209:           if (it == _node_index.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Node with label not found: " << tokens[label_index];
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           n = it->second;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {
alpar@209:           _node_maps[i].second->set(n, tokens[map_index[i]]);
alpar@209:         }
deba@127: 
deba@127:       }
deba@127:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@127:       }
deba@127:     }
deba@127: 
deba@127:     void readArcs() {
deba@127: 
deba@127:       std::vector<int> map_index(_arc_maps.size());
deba@127:       int map_num, label_index;
deba@127: 
deba@186:       char c;
deba@186:       if (!readLine() || !(line >> c) || c == '@') {
alpar@209:         if (readSuccess() && line) line.putback(c);
alpar@209:         if (!_arc_maps.empty())
deba@290:           throw FormatError("Cannot find map names");
alpar@209:         return;
deba@186:       }
deba@186:       line.putback(c);
alpar@209: 
deba@127:       {
alpar@209:         std::map<std::string, int> maps;
alpar@209: 
alpar@209:         std::string map;
alpar@209:         int index = 0;
alpar@209:         while (_reader_bits::readToken(line, map)) {
alpar@949:           if(map == "-") {
alpar@949:               if(index!=0)
alpar@949:                 throw FormatError("'-' is not allowed as a map name");
alpar@949:               else if (line >> std::ws >> c)
alpar@949:                 throw FormatError("Extra character at the end of line");
alpar@949:               else break;
alpar@949:             }
alpar@209:           if (maps.find(map) != maps.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Multiple occurence of arc map: " << map;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           maps.insert(std::make_pair(map, index));
alpar@209:           ++index;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_arc_maps.size()); ++i) {
alpar@209:           std::map<std::string, int>::iterator jt =
alpar@209:             maps.find(_arc_maps[i].first);
alpar@209:           if (jt == maps.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Map not found: " << _arc_maps[i].first;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           map_index[i] = jt->second;
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           std::map<std::string, int>::iterator jt = maps.find("label");
alpar@209:           if (jt != maps.end()) {
alpar@209:             label_index = jt->second;
alpar@209:           } else {
alpar@209:             label_index = -1;
alpar@209:           }
alpar@209:         }
alpar@209:         map_num = maps.size();
deba@127:       }
deba@127: 
deba@127:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::string source_token;
alpar@209:         std::string target_token;
alpar@209: 
alpar@209:         if (!_reader_bits::readToken(line, source_token))
deba@290:           throw FormatError("Source not found");
alpar@209: 
alpar@209:         if (!_reader_bits::readToken(line, target_token))
deba@290:           throw FormatError("Target not found");
alpar@209: 
alpar@209:         std::vector<std::string> tokens(map_num);
alpar@209:         for (int i = 0; i < map_num; ++i) {
alpar@209:           if (!_reader_bits::readToken(line, tokens[i])) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Column not found (" << i + 1 << ")";
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:         }
alpar@209:         if (line >> std::ws >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         Arc a;
alpar@209:         if (!_use_arcs) {
deba@127: 
deba@127:           typename NodeIndex::iterator it;
alpar@209: 
deba@127:           it = _node_index.find(source_token);
deba@127:           if (it == _node_index.end()) {
deba@127:             std::ostringstream msg;
deba@127:             msg << "Item not found: " << source_token;
deba@290:             throw FormatError(msg.str());
deba@127:           }
deba@127:           Node source = it->second;
deba@127: 
deba@127:           it = _node_index.find(target_token);
alpar@209:           if (it == _node_index.end()) {
alpar@209:             std::ostringstream msg;
deba@127:             msg << "Item not found: " << target_token;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           Node target = it->second;
alpar@209: 
alpar@209:           a = _digraph.addArc(source, target);
alpar@209:           if (label_index != -1)
alpar@209:             _arc_index.insert(std::make_pair(tokens[label_index], a));
alpar@209:         } else {
alpar@209:           if (label_index == -1)
kpeter@291:             throw FormatError("Label map not found");
alpar@209:           typename std::map<std::string, Arc>::iterator it =
alpar@209:             _arc_index.find(tokens[label_index]);
alpar@209:           if (it == _arc_index.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Arc with label not found: " << tokens[label_index];
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           a = it->second;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_arc_maps.size()); ++i) {
alpar@209:           _arc_maps[i].second->set(a, tokens[map_index[i]]);
alpar@209:         }
deba@127: 
deba@127:       }
deba@127:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@127:       }
deba@127:     }
deba@127: 
deba@127:     void readAttributes() {
deba@127: 
deba@127:       std::set<std::string> read_attr;
deba@127: 
deba@127:       char c;
deba@127:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::string attr, token;
alpar@209:         if (!_reader_bits::readToken(line, attr))
deba@290:           throw FormatError("Attribute name not found");
alpar@209:         if (!_reader_bits::readToken(line, token))
deba@290:           throw FormatError("Attribute value not found");
alpar@209:         if (line >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         {
alpar@209:           std::set<std::string>::iterator it = read_attr.find(attr);
alpar@209:           if (it != read_attr.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Multiple occurence of attribute: " << attr;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           read_attr.insert(attr);
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           typename Attributes::iterator it = _attributes.lower_bound(attr);
alpar@209:           while (it != _attributes.end() && it->first == attr) {
alpar@209:             it->second->set(token);
alpar@209:             ++it;
alpar@209:           }
alpar@209:         }
deba@127: 
deba@127:       }
deba@127:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@127:       }
deba@127:       for (typename Attributes::iterator it = _attributes.begin();
alpar@209:            it != _attributes.end(); ++it) {
alpar@209:         if (read_attr.find(it->first) == read_attr.end()) {
alpar@209:           std::ostringstream msg;
kpeter@291:           msg << "Attribute not found: " << it->first;
deba@290:           throw FormatError(msg.str());
alpar@209:         }
deba@127:       }
deba@127:     }
deba@127: 
deba@127:   public:
alpar@156: 
kpeter@584:     /// \name Execution of the Reader
alpar@156:     /// @{
alpar@156: 
alpar@156:     /// \brief Start the batch processing
alpar@156:     ///
alpar@156:     /// This function starts the batch processing
deba@127:     void run() {
deba@127:       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");
alpar@209: 
deba@188:       bool nodes_done = _skip_nodes;
deba@188:       bool arcs_done = _skip_arcs;
deba@127:       bool attributes_done = false;
deba@127: 
alpar@209:       line_num = 0;
deba@127:       readLine();
deba@172:       skipSection();
deba@127: 
deba@127:       while (readSuccess()) {
alpar@209:         try {
alpar@209:           char c;
alpar@209:           std::string section, caption;
alpar@209:           line >> c;
alpar@209:           _reader_bits::readToken(line, section);
alpar@209:           _reader_bits::readToken(line, caption);
alpar@209: 
alpar@209:           if (line >> c)
kpeter@291:             throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:           if (section == "nodes" && !nodes_done) {
alpar@209:             if (_nodes_caption.empty() || _nodes_caption == caption) {
alpar@209:               readNodes();
alpar@209:               nodes_done = true;
alpar@209:             }
alpar@209:           } else if ((section == "arcs" || section == "edges") &&
alpar@209:                      !arcs_done) {
alpar@209:             if (_arcs_caption.empty() || _arcs_caption == caption) {
alpar@209:               readArcs();
alpar@209:               arcs_done = true;
alpar@209:             }
alpar@209:           } else if (section == "attributes" && !attributes_done) {
alpar@209:             if (_attributes_caption.empty() || _attributes_caption == caption) {
alpar@209:               readAttributes();
alpar@209:               attributes_done = true;
alpar@209:             }
alpar@209:           } else {
alpar@209:             readLine();
alpar@209:             skipSection();
alpar@209:           }
deba@290:         } catch (FormatError& error) {
alpar@209:           error.line(line_num);
deba@290:           error.file(_filename);
alpar@209:           throw;
alpar@209:         }
deba@127:       }
deba@127: 
deba@127:       if (!nodes_done) {
deba@290:         throw FormatError("Section @nodes not found");
deba@127:       }
deba@127: 
deba@127:       if (!arcs_done) {
deba@290:         throw FormatError("Section @arcs not found");
deba@127:       }
deba@127: 
deba@127:       if (!attributes_done && !_attributes.empty()) {
deba@290:         throw FormatError("Section @attributes not found");
deba@127:       }
deba@127: 
deba@127:     }
alpar@156: 
alpar@156:     /// @}
alpar@209: 
deba@127:   };
alpar@877: 
deba@598:   /// \ingroup lemon_io
deba@598:   ///
deba@598:   /// \brief Return a \ref DigraphReader class
deba@598:   ///
deba@598:   /// This function just returns a \ref DigraphReader class.
deba@598:   ///
alpar@877:   /// With this function a digraph can be read from an
deba@598:   /// \ref lgf-format "LGF" file or input stream with several maps and
deba@598:   /// attributes. For example, there is network flow problem on a
deba@598:   /// digraph, i.e. a digraph with a \e capacity map on the arcs and
deba@598:   /// \e source and \e target nodes. This digraph can be read with the
deba@598:   /// following code:
deba@598:   ///
deba@598:   ///\code
deba@598:   ///ListDigraph digraph;
deba@598:   ///ListDigraph::ArcMap<int> cm(digraph);
deba@598:   ///ListDigraph::Node src, trg;
deba@598:   ///digraphReader(digraph, std::cin).
deba@598:   ///  arcMap("capacity", cap).
deba@598:   ///  node("source", src).
deba@598:   ///  node("target", trg).
deba@598:   ///  run();
deba@598:   ///\endcode
deba@598:   ///
deba@598:   /// For a complete documentation, please see the \ref DigraphReader
deba@598:   /// class documentation.
deba@598:   /// \warning Don't forget to put the \ref DigraphReader::run() "run()"
deba@598:   /// to the end of the parameter list.
deba@598:   /// \relates DigraphReader
deba@598:   /// \sa digraphReader(TDGR& digraph, const std::string& fn)
deba@598:   /// \sa digraphReader(TDGR& digraph, const char* fn)
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, std::istream& is) {
deba@598:     DigraphReader<TDGR> tmp(digraph, is);
deba@598:     return tmp;
deba@598:   }
deba@127: 
kpeter@498:   /// \brief Return a \ref DigraphReader class
kpeter@498:   ///
kpeter@498:   /// This function just returns a \ref DigraphReader class.
kpeter@498:   /// \relates DigraphReader
deba@598:   /// \sa digraphReader(TDGR& digraph, std::istream& is)
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, const std::string& fn) {
deba@598:     DigraphReader<TDGR> tmp(digraph, fn);
kpeter@498:     return tmp;
kpeter@498:   }
kpeter@498: 
kpeter@498:   /// \brief Return a \ref DigraphReader class
kpeter@498:   ///
kpeter@498:   /// This function just returns a \ref DigraphReader class.
kpeter@498:   /// \relates DigraphReader
deba@598:   /// \sa digraphReader(TDGR& digraph, std::istream& is)
deba@598:   template <typename TDGR>
deba@598:   DigraphReader<TDGR> digraphReader(TDGR& digraph, const char* fn) {
deba@598:     DigraphReader<TDGR> tmp(digraph, fn);
kpeter@498:     return tmp;
kpeter@498:   }
kpeter@498: 
deba@598:   template <typename GR>
ladanyi@303:   class GraphReader;
alpar@877: 
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, std::istream& is = std::cin);
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, const std::string& fn);
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, const char *fn);
deba@165: 
deba@165:   /// \ingroup lemon_io
alpar@209:   ///
kpeter@192:   /// \brief \ref lgf-format "LGF" reader for undirected graphs
deba@165:   ///
deba@165:   /// This utility reads an \ref lgf-format "LGF" file.
kpeter@192:   ///
kpeter@192:   /// It can be used almost the same way as \c DigraphReader.
kpeter@192:   /// The only difference is that this class can handle edges and
kpeter@192:   /// edge maps as well as arcs and arc maps.
deba@201:   ///
deba@201:   /// The columns in the \c \@edges (or \c \@arcs) section are the
deba@201:   /// edge maps. However, if there are two maps with the same name
deba@201:   /// prefixed with \c '+' and \c '-', then these can be read into an
deba@201:   /// arc map.  Similarly, an attribute can be read into an arc, if
deba@201:   /// it's value is an edge label prefixed with \c '+' or \c '-'.
kpeter@559:   template <typename GR>
deba@165:   class GraphReader {
deba@165:   public:
deba@165: 
kpeter@559:     typedef GR Graph;
kpeter@559: 
kpeter@559:   private:
kpeter@559: 
deba@598:     TEMPLATE_GRAPH_TYPEDEFS(GR);
alpar@209: 
deba@165:     std::istream* _is;
deba@165:     bool local_is;
deba@290:     std::string _filename;
deba@165: 
deba@598:     GR& _graph;
deba@165: 
deba@165:     std::string _nodes_caption;
deba@165:     std::string _edges_caption;
deba@165:     std::string _attributes_caption;
deba@165: 
deba@165:     typedef std::map<std::string, Node> NodeIndex;
deba@165:     NodeIndex _node_index;
deba@165:     typedef std::map<std::string, Edge> EdgeIndex;
deba@165:     EdgeIndex _edge_index;
alpar@209: 
alpar@209:     typedef std::vector<std::pair<std::string,
alpar@209:       _reader_bits::MapStorageBase<Node>*> > NodeMaps;
alpar@209:     NodeMaps _node_maps;
deba@165: 
deba@165:     typedef std::vector<std::pair<std::string,
deba@165:       _reader_bits::MapStorageBase<Edge>*> > EdgeMaps;
deba@165:     EdgeMaps _edge_maps;
deba@165: 
alpar@209:     typedef std::multimap<std::string, _reader_bits::ValueStorageBase*>
deba@165:       Attributes;
deba@165:     Attributes _attributes;
deba@165: 
deba@165:     bool _use_nodes;
deba@165:     bool _use_edges;
deba@165: 
deba@188:     bool _skip_nodes;
deba@188:     bool _skip_edges;
deba@188: 
deba@165:     int line_num;
deba@165:     std::istringstream line;
deba@165: 
deba@165:   public:
deba@165: 
deba@165:     /// \brief Constructor
deba@165:     ///
kpeter@192:     /// Construct an undirected graph reader, which reads from the given
deba@165:     /// input stream.
deba@598:     GraphReader(GR& graph, std::istream& is = std::cin)
deba@165:       : _is(&is), local_is(false), _graph(graph),
alpar@209:         _use_nodes(false), _use_edges(false),
alpar@209:         _skip_nodes(false), _skip_edges(false) {}
deba@165: 
deba@165:     /// \brief Constructor
deba@165:     ///
kpeter@192:     /// Construct an undirected graph reader, which reads from the given
deba@165:     /// file.
deba@598:     GraphReader(GR& graph, const std::string& fn)
deba@290:       : _is(new std::ifstream(fn.c_str())), local_is(true),
deba@290:         _filename(fn), _graph(graph),
kpeter@212:         _use_nodes(false), _use_edges(false),
deba@290:         _skip_nodes(false), _skip_edges(false) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
alpar@209: 
deba@165:     /// \brief Constructor
deba@165:     ///
kpeter@192:     /// Construct an undirected graph reader, which reads from the given
deba@165:     /// file.
deba@598:     GraphReader(GR& graph, const char* fn)
deba@290:       : _is(new std::ifstream(fn)), local_is(true),
deba@290:         _filename(fn), _graph(graph),
kpeter@212:         _use_nodes(false), _use_edges(false),
deba@290:         _skip_nodes(false), _skip_edges(false) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
deba@165: 
deba@165:     /// \brief Destructor
deba@165:     ~GraphReader() {
alpar@209:       for (typename NodeMaps::iterator it = _node_maps.begin();
alpar@209:            it != _node_maps.end(); ++it) {
alpar@209:         delete it->second;
deba@165:       }
deba@165: 
alpar@209:       for (typename EdgeMaps::iterator it = _edge_maps.begin();
alpar@209:            it != _edge_maps.end(); ++it) {
alpar@209:         delete it->second;
deba@165:       }
deba@165: 
alpar@209:       for (typename Attributes::iterator it = _attributes.begin();
alpar@209:            it != _attributes.end(); ++it) {
alpar@209:         delete it->second;
deba@165:       }
deba@165: 
deba@165:       if (local_is) {
alpar@209:         delete _is;
deba@165:       }
deba@165: 
deba@165:     }
deba@165: 
deba@165:   private:
deba@598:     template <typename TGR>
deba@598:     friend GraphReader<TGR> graphReader(TGR& graph, std::istream& is);
deba@598:     template <typename TGR>
alpar@877:     friend GraphReader<TGR> graphReader(TGR& graph, const std::string& fn);
deba@598:     template <typename TGR>
deba@598:     friend GraphReader<TGR> graphReader(TGR& graph, const char *fn);
alpar@209: 
alpar@209:     GraphReader(GraphReader& other)
deba@190:       : _is(other._is), local_is(other.local_is), _graph(other._graph),
alpar@209:         _use_nodes(other._use_nodes), _use_edges(other._use_edges),
alpar@209:         _skip_nodes(other._skip_nodes), _skip_edges(other._skip_edges) {
deba@190: 
deba@190:       other._is = 0;
deba@190:       other.local_is = false;
alpar@209: 
deba@190:       _node_index.swap(other._node_index);
deba@190:       _edge_index.swap(other._edge_index);
deba@190: 
deba@190:       _node_maps.swap(other._node_maps);
deba@190:       _edge_maps.swap(other._edge_maps);
deba@190:       _attributes.swap(other._attributes);
deba@190: 
deba@190:       _nodes_caption = other._nodes_caption;
deba@190:       _edges_caption = other._edges_caption;
deba@190:       _attributes_caption = other._attributes_caption;
deba@190: 
deba@190:     }
deba@190: 
deba@165:     GraphReader& operator=(const GraphReader&);
deba@165: 
deba@165:   public:
deba@165: 
kpeter@584:     /// \name Reading Rules
deba@165:     /// @{
alpar@209: 
deba@165:     /// \brief Node map reading rule
deba@165:     ///
deba@165:     /// Add a node map reading rule to the reader.
deba@165:     template <typename Map>
deba@165:     GraphReader& nodeMap(const std::string& caption, Map& map) {
deba@165:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Node>* storage =
alpar@209:         new _reader_bits::MapStorage<Node, Map>(map);
deba@165:       _node_maps.push_back(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Node map reading rule
deba@165:     ///
deba@165:     /// Add a node map reading rule with specialized converter to the
deba@165:     /// reader.
deba@165:     template <typename Map, typename Converter>
alpar@209:     GraphReader& nodeMap(const std::string& caption, Map& map,
alpar@209:                            const Converter& converter = Converter()) {
deba@165:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Node>* storage =
alpar@209:         new _reader_bits::MapStorage<Node, Map, Converter>(map, converter);
deba@165:       _node_maps.push_back(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Edge map reading rule
deba@165:     ///
deba@165:     /// Add an edge map reading rule to the reader.
deba@165:     template <typename Map>
deba@165:     GraphReader& edgeMap(const std::string& caption, Map& map) {
deba@165:       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Edge>* storage =
alpar@209:         new _reader_bits::MapStorage<Edge, Map>(map);
deba@165:       _edge_maps.push_back(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Edge map reading rule
deba@165:     ///
deba@165:     /// Add an edge map reading rule with specialized converter to the
deba@165:     /// reader.
deba@165:     template <typename Map, typename Converter>
alpar@209:     GraphReader& edgeMap(const std::string& caption, Map& map,
alpar@209:                           const Converter& converter = Converter()) {
deba@165:       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Edge>* storage =
alpar@209:         new _reader_bits::MapStorage<Edge, Map, Converter>(map, converter);
deba@165:       _edge_maps.push_back(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Arc map reading rule
deba@165:     ///
deba@165:     /// Add an arc map reading rule to the reader.
deba@165:     template <typename Map>
deba@165:     GraphReader& arcMap(const std::string& caption, Map& map) {
deba@165:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Edge>* forward_storage =
alpar@209:         new _reader_bits::GraphArcMapStorage<Graph, true, Map>(_graph, map);
deba@165:       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
alpar@209:       _reader_bits::MapStorageBase<Edge>* backward_storage =
deba@598:         new _reader_bits::GraphArcMapStorage<GR, false, Map>(_graph, map);
deba@165:       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Arc map reading rule
deba@165:     ///
deba@165:     /// Add an arc map reading rule with specialized converter to the
deba@165:     /// reader.
deba@165:     template <typename Map, typename Converter>
alpar@209:     GraphReader& arcMap(const std::string& caption, Map& map,
alpar@209:                           const Converter& converter = Converter()) {
deba@165:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
alpar@209:       _reader_bits::MapStorageBase<Edge>* forward_storage =
deba@598:         new _reader_bits::GraphArcMapStorage<GR, true, Map, Converter>
alpar@209:         (_graph, map, converter);
deba@165:       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
alpar@209:       _reader_bits::MapStorageBase<Edge>* backward_storage =
deba@598:         new _reader_bits::GraphArcMapStorage<GR, false, Map, Converter>
alpar@209:         (_graph, map, converter);
deba@165:       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Attribute reading rule
deba@165:     ///
deba@165:     /// Add an attribute reading rule to the reader.
deba@165:     template <typename Value>
deba@165:     GraphReader& attribute(const std::string& caption, Value& value) {
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Value>(value);
deba@165:       _attributes.insert(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Attribute reading rule
deba@165:     ///
deba@165:     /// Add an attribute reading rule with specialized converter to the
deba@165:     /// reader.
deba@165:     template <typename Value, typename Converter>
alpar@209:     GraphReader& attribute(const std::string& caption, Value& value,
alpar@209:                              const Converter& converter = Converter()) {
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Value, Converter>(value, converter);
deba@165:       _attributes.insert(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Node reading rule
deba@165:     ///
deba@165:     /// Add a node reading rule to reader.
deba@165:     GraphReader& node(const std::string& caption, Node& node) {
deba@165:       typedef _reader_bits::MapLookUpConverter<Node> Converter;
deba@165:       Converter converter(_node_index);
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Node, Converter>(node, converter);
deba@165:       _attributes.insert(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Edge reading rule
deba@165:     ///
deba@165:     /// Add an edge reading rule to reader.
deba@165:     GraphReader& edge(const std::string& caption, Edge& edge) {
deba@165:       typedef _reader_bits::MapLookUpConverter<Edge> Converter;
deba@165:       Converter converter(_edge_index);
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Edge, Converter>(edge, converter);
deba@165:       _attributes.insert(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Arc reading rule
deba@165:     ///
deba@165:     /// Add an arc reading rule to reader.
deba@165:     GraphReader& arc(const std::string& caption, Arc& arc) {
deba@598:       typedef _reader_bits::GraphArcLookUpConverter<GR> Converter;
deba@165:       Converter converter(_graph, _edge_index);
alpar@209:       _reader_bits::ValueStorageBase* storage =
alpar@209:         new _reader_bits::ValueStorage<Arc, Converter>(arc, converter);
deba@165:       _attributes.insert(std::make_pair(caption, storage));
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// @}
deba@165: 
kpeter@584:     /// \name Select Section by Name
deba@165:     /// @{
deba@165: 
deba@165:     /// \brief Set \c \@nodes section to be read
deba@165:     ///
kpeter@192:     /// Set \c \@nodes section to be read.
deba@165:     GraphReader& nodes(const std::string& caption) {
deba@165:       _nodes_caption = caption;
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Set \c \@edges section to be read
deba@165:     ///
kpeter@192:     /// Set \c \@edges section to be read.
deba@165:     GraphReader& edges(const std::string& caption) {
deba@165:       _edges_caption = caption;
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Set \c \@attributes section to be read
deba@165:     ///
kpeter@192:     /// Set \c \@attributes section to be read.
deba@165:     GraphReader& attributes(const std::string& caption) {
deba@165:       _attributes_caption = caption;
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// @}
deba@165: 
kpeter@584:     /// \name Using Previously Constructed Node or Edge Set
deba@165:     /// @{
deba@165: 
deba@165:     /// \brief Use previously constructed node set
deba@165:     ///
deba@165:     /// Use previously constructed node set, and specify the node
deba@165:     /// label map.
deba@165:     template <typename Map>
deba@165:     GraphReader& useNodes(const Map& map) {
deba@165:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@165:       _use_nodes = true;
deba@165:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@165:       for (NodeIt n(_graph); n != INVALID; ++n) {
alpar@209:         _node_index.insert(std::make_pair(converter(map[n]), n));
deba@165:       }
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Use previously constructed node set
deba@165:     ///
deba@165:     /// Use previously constructed node set, and specify the node
deba@165:     /// label map and a functor which converts the label map values to
kpeter@192:     /// \c std::string.
deba@165:     template <typename Map, typename Converter>
alpar@209:     GraphReader& useNodes(const Map& map,
alpar@209:                             const Converter& converter = Converter()) {
deba@165:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@165:       _use_nodes = true;
deba@165:       for (NodeIt n(_graph); n != INVALID; ++n) {
alpar@209:         _node_index.insert(std::make_pair(converter(map[n]), n));
deba@165:       }
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Use previously constructed edge set
deba@165:     ///
deba@165:     /// Use previously constructed edge set, and specify the edge
deba@165:     /// label map.
deba@165:     template <typename Map>
deba@165:     GraphReader& useEdges(const Map& map) {
deba@165:       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
deba@165:       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");
deba@165:       _use_edges = true;
deba@165:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@165:       for (EdgeIt a(_graph); a != INVALID; ++a) {
alpar@209:         _edge_index.insert(std::make_pair(converter(map[a]), a));
deba@165:       }
deba@165:       return *this;
deba@165:     }
deba@165: 
deba@165:     /// \brief Use previously constructed edge set
deba@165:     ///
deba@165:     /// Use previously constructed edge set, and specify the edge
deba@165:     /// label map and a functor which converts the label map values to
kpeter@192:     /// \c std::string.
deba@165:     template <typename Map, typename Converter>
alpar@209:     GraphReader& useEdges(const Map& map,
alpar@209:                             const Converter& converter = Converter()) {
deba@165:       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
alpar@209:       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");
deba@165:       _use_edges = true;
deba@165:       for (EdgeIt a(_graph); a != INVALID; ++a) {
alpar@209:         _edge_index.insert(std::make_pair(converter(map[a]), a));
deba@165:       }
deba@165:       return *this;
deba@165:     }
deba@165: 
kpeter@192:     /// \brief Skip the reading of node section
deba@188:     ///
deba@188:     /// Omit the reading of the node section. This implies that each node
kpeter@192:     /// map reading rule will be abandoned, and the nodes of the graph
deba@188:     /// will not be constructed, which usually cause that the edge set
deba@188:     /// could not be read due to lack of node name
kpeter@192:     /// could not be read due to lack of node name resolving.
kpeter@192:     /// Therefore \c skipEdges() function should also be used, or
kpeter@192:     /// \c useNodes() should be used to specify the label of the nodes.
deba@188:     GraphReader& skipNodes() {
alpar@209:       LEMON_ASSERT(!_skip_nodes, "Skip nodes already set");
deba@188:       _skip_nodes = true;
deba@188:       return *this;
deba@188:     }
deba@188: 
kpeter@192:     /// \brief Skip the reading of edge section
deba@188:     ///
deba@188:     /// Omit the reading of the edge section. This implies that each edge
kpeter@192:     /// map reading rule will be abandoned, and the edges of the graph
deba@188:     /// will not be constructed.
deba@188:     GraphReader& skipEdges() {
alpar@209:       LEMON_ASSERT(!_skip_edges, "Skip edges already set");
deba@188:       _skip_edges = true;
deba@188:       return *this;
deba@188:     }
deba@188: 
deba@165:     /// @}
deba@165: 
deba@165:   private:
deba@165: 
deba@165:     bool readLine() {
deba@165:       std::string str;
deba@165:       while(++line_num, std::getline(*_is, str)) {
alpar@209:         line.clear(); line.str(str);
alpar@209:         char c;
alpar@209:         if (line >> std::ws >> c && c != '#') {
alpar@209:           line.putback(c);
alpar@209:           return true;
alpar@209:         }
deba@165:       }
deba@165:       return false;
deba@165:     }
deba@165: 
deba@165:     bool readSuccess() {
deba@165:       return static_cast<bool>(*_is);
deba@165:     }
alpar@209: 
deba@165:     void skipSection() {
deba@165:       char c;
deba@165:       while (readSuccess() && line >> c && c != '@') {
alpar@209:         readLine();
deba@165:       }
deba@427:       if (readSuccess()) {
deba@427:         line.putback(c);
deba@427:       }
deba@165:     }
deba@165: 
deba@165:     void readNodes() {
deba@165: 
deba@165:       std::vector<int> map_index(_node_maps.size());
deba@165:       int map_num, label_index;
deba@165: 
deba@186:       char c;
deba@186:       if (!readLine() || !(line >> c) || c == '@') {
alpar@209:         if (readSuccess() && line) line.putback(c);
alpar@209:         if (!_node_maps.empty())
deba@290:           throw FormatError("Cannot find map names");
alpar@209:         return;
deba@186:       }
deba@186:       line.putback(c);
alpar@209: 
deba@165:       {
alpar@209:         std::map<std::string, int> maps;
alpar@209: 
alpar@209:         std::string map;
alpar@209:         int index = 0;
alpar@209:         while (_reader_bits::readToken(line, map)) {
alpar@209:           if (maps.find(map) != maps.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Multiple occurence of node map: " << map;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           maps.insert(std::make_pair(map, index));
alpar@209:           ++index;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {
alpar@209:           std::map<std::string, int>::iterator jt =
alpar@209:             maps.find(_node_maps[i].first);
alpar@209:           if (jt == maps.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Map not found: " << _node_maps[i].first;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           map_index[i] = jt->second;
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           std::map<std::string, int>::iterator jt = maps.find("label");
alpar@209:           if (jt != maps.end()) {
alpar@209:             label_index = jt->second;
alpar@209:           } else {
alpar@209:             label_index = -1;
alpar@209:           }
alpar@209:         }
alpar@209:         map_num = maps.size();
deba@165:       }
deba@165: 
deba@165:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::vector<std::string> tokens(map_num);
alpar@209:         for (int i = 0; i < map_num; ++i) {
alpar@209:           if (!_reader_bits::readToken(line, tokens[i])) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Column not found (" << i + 1 << ")";
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:         }
alpar@209:         if (line >> std::ws >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         Node n;
alpar@209:         if (!_use_nodes) {
alpar@209:           n = _graph.addNode();
alpar@209:           if (label_index != -1)
alpar@209:             _node_index.insert(std::make_pair(tokens[label_index], n));
alpar@209:         } else {
alpar@209:           if (label_index == -1)
kpeter@291:             throw FormatError("Label map not found");
alpar@209:           typename std::map<std::string, Node>::iterator it =
alpar@209:             _node_index.find(tokens[label_index]);
alpar@209:           if (it == _node_index.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Node with label not found: " << tokens[label_index];
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           n = it->second;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_node_maps.size()); ++i) {
alpar@209:           _node_maps[i].second->set(n, tokens[map_index[i]]);
alpar@209:         }
deba@165: 
deba@165:       }
deba@165:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@165:       }
deba@165:     }
deba@165: 
deba@165:     void readEdges() {
deba@165: 
deba@165:       std::vector<int> map_index(_edge_maps.size());
deba@165:       int map_num, label_index;
deba@165: 
deba@186:       char c;
deba@186:       if (!readLine() || !(line >> c) || c == '@') {
alpar@209:         if (readSuccess() && line) line.putback(c);
alpar@209:         if (!_edge_maps.empty())
deba@290:           throw FormatError("Cannot find map names");
alpar@209:         return;
deba@186:       }
deba@186:       line.putback(c);
alpar@209: 
deba@165:       {
alpar@209:         std::map<std::string, int> maps;
alpar@209: 
alpar@209:         std::string map;
alpar@209:         int index = 0;
alpar@209:         while (_reader_bits::readToken(line, map)) {
alpar@949:           if(map == "-") {
alpar@949:               if(index!=0)
alpar@949:                 throw FormatError("'-' is not allowed as a map name");
alpar@949:               else if (line >> std::ws >> c)
alpar@949:                 throw FormatError("Extra character at the end of line");
alpar@949:               else break;
alpar@949:             }
alpar@209:           if (maps.find(map) != maps.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Multiple occurence of edge map: " << map;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           maps.insert(std::make_pair(map, index));
alpar@209:           ++index;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {
alpar@209:           std::map<std::string, int>::iterator jt =
alpar@209:             maps.find(_edge_maps[i].first);
alpar@209:           if (jt == maps.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Map not found: " << _edge_maps[i].first;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           map_index[i] = jt->second;
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           std::map<std::string, int>::iterator jt = maps.find("label");
alpar@209:           if (jt != maps.end()) {
alpar@209:             label_index = jt->second;
alpar@209:           } else {
alpar@209:             label_index = -1;
alpar@209:           }
alpar@209:         }
alpar@209:         map_num = maps.size();
deba@165:       }
deba@165: 
deba@165:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::string source_token;
alpar@209:         std::string target_token;
alpar@209: 
alpar@209:         if (!_reader_bits::readToken(line, source_token))
deba@290:           throw FormatError("Node u not found");
alpar@209: 
alpar@209:         if (!_reader_bits::readToken(line, target_token))
deba@290:           throw FormatError("Node v not found");
alpar@209: 
alpar@209:         std::vector<std::string> tokens(map_num);
alpar@209:         for (int i = 0; i < map_num; ++i) {
alpar@209:           if (!_reader_bits::readToken(line, tokens[i])) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Column not found (" << i + 1 << ")";
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:         }
alpar@209:         if (line >> std::ws >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         Edge e;
alpar@209:         if (!_use_edges) {
deba@165: 
deba@165:           typename NodeIndex::iterator it;
alpar@209: 
deba@165:           it = _node_index.find(source_token);
deba@165:           if (it == _node_index.end()) {
deba@165:             std::ostringstream msg;
deba@165:             msg << "Item not found: " << source_token;
deba@290:             throw FormatError(msg.str());
deba@165:           }
deba@165:           Node source = it->second;
deba@165: 
deba@165:           it = _node_index.find(target_token);
alpar@209:           if (it == _node_index.end()) {
alpar@209:             std::ostringstream msg;
deba@165:             msg << "Item not found: " << target_token;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           Node target = it->second;
alpar@209: 
alpar@209:           e = _graph.addEdge(source, target);
alpar@209:           if (label_index != -1)
alpar@209:             _edge_index.insert(std::make_pair(tokens[label_index], e));
alpar@209:         } else {
alpar@209:           if (label_index == -1)
kpeter@291:             throw FormatError("Label map not found");
alpar@209:           typename std::map<std::string, Edge>::iterator it =
alpar@209:             _edge_index.find(tokens[label_index]);
alpar@209:           if (it == _edge_index.end()) {
alpar@209:             std::ostringstream msg;
alpar@209:             msg << "Edge with label not found: " << tokens[label_index];
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           e = it->second;
alpar@209:         }
alpar@209: 
alpar@209:         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {
alpar@209:           _edge_maps[i].second->set(e, tokens[map_index[i]]);
alpar@209:         }
deba@165: 
deba@165:       }
deba@165:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@165:       }
deba@165:     }
deba@165: 
deba@165:     void readAttributes() {
deba@165: 
deba@165:       std::set<std::string> read_attr;
deba@165: 
deba@165:       char c;
deba@165:       while (readLine() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209: 
alpar@209:         std::string attr, token;
alpar@209:         if (!_reader_bits::readToken(line, attr))
deba@290:           throw FormatError("Attribute name not found");
alpar@209:         if (!_reader_bits::readToken(line, token))
deba@290:           throw FormatError("Attribute value not found");
alpar@209:         if (line >> c)
kpeter@291:           throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:         {
alpar@209:           std::set<std::string>::iterator it = read_attr.find(attr);
alpar@209:           if (it != read_attr.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Multiple occurence of attribute: " << attr;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           read_attr.insert(attr);
alpar@209:         }
alpar@209: 
alpar@209:         {
alpar@209:           typename Attributes::iterator it = _attributes.lower_bound(attr);
alpar@209:           while (it != _attributes.end() && it->first == attr) {
alpar@209:             it->second->set(token);
alpar@209:             ++it;
alpar@209:           }
alpar@209:         }
deba@165: 
deba@165:       }
deba@165:       if (readSuccess()) {
alpar@209:         line.putback(c);
deba@165:       }
deba@165:       for (typename Attributes::iterator it = _attributes.begin();
alpar@209:            it != _attributes.end(); ++it) {
alpar@209:         if (read_attr.find(it->first) == read_attr.end()) {
alpar@209:           std::ostringstream msg;
kpeter@291:           msg << "Attribute not found: " << it->first;
deba@290:           throw FormatError(msg.str());
alpar@209:         }
deba@165:       }
deba@165:     }
deba@165: 
deba@165:   public:
deba@165: 
kpeter@584:     /// \name Execution of the Reader
deba@165:     /// @{
deba@165: 
deba@165:     /// \brief Start the batch processing
deba@165:     ///
deba@165:     /// This function starts the batch processing
deba@165:     void run() {
alpar@209: 
deba@165:       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");
alpar@209: 
deba@188:       bool nodes_done = _skip_nodes;
deba@188:       bool edges_done = _skip_edges;
deba@165:       bool attributes_done = false;
deba@165: 
alpar@209:       line_num = 0;
deba@165:       readLine();
deba@172:       skipSection();
deba@165: 
deba@165:       while (readSuccess()) {
alpar@209:         try {
alpar@209:           char c;
alpar@209:           std::string section, caption;
alpar@209:           line >> c;
alpar@209:           _reader_bits::readToken(line, section);
alpar@209:           _reader_bits::readToken(line, caption);
alpar@209: 
alpar@209:           if (line >> c)
kpeter@291:             throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:           if (section == "nodes" && !nodes_done) {
alpar@209:             if (_nodes_caption.empty() || _nodes_caption == caption) {
alpar@209:               readNodes();
alpar@209:               nodes_done = true;
alpar@209:             }
alpar@209:           } else if ((section == "edges" || section == "arcs") &&
alpar@209:                      !edges_done) {
alpar@209:             if (_edges_caption.empty() || _edges_caption == caption) {
alpar@209:               readEdges();
alpar@209:               edges_done = true;
alpar@209:             }
alpar@209:           } else if (section == "attributes" && !attributes_done) {
alpar@209:             if (_attributes_caption.empty() || _attributes_caption == caption) {
alpar@209:               readAttributes();
alpar@209:               attributes_done = true;
alpar@209:             }
alpar@209:           } else {
alpar@209:             readLine();
alpar@209:             skipSection();
alpar@209:           }
deba@290:         } catch (FormatError& error) {
alpar@209:           error.line(line_num);
deba@290:           error.file(_filename);
alpar@209:           throw;
alpar@209:         }
deba@165:       }
deba@165: 
deba@165:       if (!nodes_done) {
deba@290:         throw FormatError("Section @nodes not found");
deba@165:       }
deba@165: 
deba@165:       if (!edges_done) {
deba@290:         throw FormatError("Section @edges not found");
deba@165:       }
deba@165: 
deba@165:       if (!attributes_done && !_attributes.empty()) {
deba@290:         throw FormatError("Section @attributes not found");
deba@165:       }
deba@165: 
deba@165:     }
deba@165: 
deba@165:     /// @}
alpar@209: 
deba@165:   };
deba@165: 
deba@598:   /// \ingroup lemon_io
deba@598:   ///
deba@598:   /// \brief Return a \ref GraphReader class
deba@598:   ///
alpar@877:   /// This function just returns a \ref GraphReader class.
deba@598:   ///
alpar@877:   /// With this function a graph can be read from an
deba@598:   /// \ref lgf-format "LGF" file or input stream with several maps and
deba@598:   /// attributes. For example, there is weighted matching problem on a
deba@598:   /// graph, i.e. a graph with a \e weight map on the edges. This
deba@598:   /// graph can be read with the following code:
deba@598:   ///
deba@598:   ///\code
deba@598:   ///ListGraph graph;
deba@598:   ///ListGraph::EdgeMap<int> weight(graph);
deba@598:   ///graphReader(graph, std::cin).
deba@598:   ///  edgeMap("weight", weight).
deba@598:   ///  run();
deba@598:   ///\endcode
deba@598:   ///
deba@598:   /// For a complete documentation, please see the \ref GraphReader
deba@598:   /// class documentation.
deba@598:   /// \warning Don't forget to put the \ref GraphReader::run() "run()"
deba@598:   /// to the end of the parameter list.
deba@598:   /// \relates GraphReader
deba@598:   /// \sa graphReader(TGR& graph, const std::string& fn)
deba@598:   /// \sa graphReader(TGR& graph, const char* fn)
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, std::istream& is) {
deba@598:     GraphReader<TGR> tmp(graph, is);
deba@598:     return tmp;
deba@598:   }
deba@598: 
kpeter@498:   /// \brief Return a \ref GraphReader class
kpeter@498:   ///
kpeter@498:   /// This function just returns a \ref GraphReader class.
kpeter@498:   /// \relates GraphReader
deba@598:   /// \sa graphReader(TGR& graph, std::istream& is)
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, const std::string& fn) {
deba@598:     GraphReader<TGR> tmp(graph, fn);
kpeter@498:     return tmp;
kpeter@498:   }
kpeter@498: 
kpeter@498:   /// \brief Return a \ref GraphReader class
kpeter@498:   ///
kpeter@498:   /// This function just returns a \ref GraphReader class.
kpeter@498:   /// \relates GraphReader
deba@598:   /// \sa graphReader(TGR& graph, std::istream& is)
deba@598:   template <typename TGR>
deba@598:   GraphReader<TGR> graphReader(TGR& graph, const char* fn) {
deba@598:     GraphReader<TGR> tmp(graph, fn);
kpeter@498:     return tmp;
kpeter@498:   }
kpeter@498: 
deba@1024:   template <typename BGR>
deba@1024:   class BpGraphReader;
deba@1024: 
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, std::istream& is = std::cin);
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, const std::string& fn);
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, const char *fn);
deba@1024: 
deba@1024:   /// \ingroup lemon_io
deba@1024:   ///
deba@1024:   /// \brief \ref lgf-format "LGF" reader for bipartite graphs
deba@1024:   ///
deba@1024:   /// This utility reads an \ref lgf-format "LGF" file.
deba@1024:   ///
deba@1024:   /// It can be used almost the same way as \c GraphReader, but it
deba@1024:   /// reads the red and blue nodes from separate sections, and these
deba@1024:   /// sections can contain different set of maps.
deba@1024:   ///
deba@1026:   /// The red and blue node maps are read from the corresponding
deba@1024:   /// sections. If a map is defined with the same name in both of
deba@1024:   /// these sections, then it can be read as a node map.
deba@1024:   template <typename BGR>
deba@1024:   class BpGraphReader {
deba@1024:   public:
deba@1024: 
deba@1024:     typedef BGR Graph;
deba@1024: 
deba@1024:   private:
deba@1024: 
deba@1024:     TEMPLATE_BPGRAPH_TYPEDEFS(BGR);
deba@1024: 
deba@1024:     std::istream* _is;
deba@1024:     bool local_is;
deba@1024:     std::string _filename;
deba@1024: 
deba@1024:     BGR& _graph;
deba@1024: 
deba@1024:     std::string _nodes_caption;
deba@1024:     std::string _edges_caption;
deba@1024:     std::string _attributes_caption;
deba@1024: 
deba@1030:     typedef std::map<std::string, RedNode> RedNodeIndex;
deba@1030:     RedNodeIndex _red_node_index;
deba@1030:     typedef std::map<std::string, BlueNode> BlueNodeIndex;
deba@1030:     BlueNodeIndex _blue_node_index;
deba@1024:     typedef std::map<std::string, Edge> EdgeIndex;
deba@1024:     EdgeIndex _edge_index;
deba@1024: 
deba@1024:     typedef std::vector<std::pair<std::string,
deba@1030:       _reader_bits::MapStorageBase<RedNode>*> > RedNodeMaps;
deba@1030:     RedNodeMaps _red_node_maps;
deba@1030:     typedef std::vector<std::pair<std::string,
deba@1030:       _reader_bits::MapStorageBase<BlueNode>*> > BlueNodeMaps;
deba@1030:     BlueNodeMaps _blue_node_maps;
deba@1024: 
deba@1024:     typedef std::vector<std::pair<std::string,
deba@1024:       _reader_bits::MapStorageBase<Edge>*> > EdgeMaps;
deba@1024:     EdgeMaps _edge_maps;
deba@1024: 
deba@1024:     typedef std::multimap<std::string, _reader_bits::ValueStorageBase*>
deba@1024:       Attributes;
deba@1024:     Attributes _attributes;
deba@1024: 
deba@1024:     bool _use_nodes;
deba@1024:     bool _use_edges;
deba@1024: 
deba@1024:     bool _skip_nodes;
deba@1024:     bool _skip_edges;
deba@1024: 
deba@1024:     int line_num;
deba@1024:     std::istringstream line;
deba@1024: 
deba@1024:   public:
deba@1024: 
deba@1024:     /// \brief Constructor
deba@1024:     ///
deba@1024:     /// Construct an undirected graph reader, which reads from the given
deba@1024:     /// input stream.
deba@1024:     BpGraphReader(BGR& graph, std::istream& is = std::cin)
deba@1024:       : _is(&is), local_is(false), _graph(graph),
deba@1024:         _use_nodes(false), _use_edges(false),
deba@1024:         _skip_nodes(false), _skip_edges(false) {}
deba@1024: 
deba@1024:     /// \brief Constructor
deba@1024:     ///
deba@1024:     /// Construct an undirected graph reader, which reads from the given
deba@1024:     /// file.
deba@1024:     BpGraphReader(BGR& graph, const std::string& fn)
deba@1024:       : _is(new std::ifstream(fn.c_str())), local_is(true),
deba@1024:         _filename(fn), _graph(graph),
deba@1024:         _use_nodes(false), _use_edges(false),
deba@1024:         _skip_nodes(false), _skip_edges(false) {
deba@1024:       if (!(*_is)) {
deba@1024:         delete _is;
deba@1024:         throw IoError("Cannot open file", fn);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Constructor
deba@1024:     ///
deba@1024:     /// Construct an undirected graph reader, which reads from the given
deba@1024:     /// file.
deba@1024:     BpGraphReader(BGR& graph, const char* fn)
deba@1024:       : _is(new std::ifstream(fn)), local_is(true),
deba@1024:         _filename(fn), _graph(graph),
deba@1024:         _use_nodes(false), _use_edges(false),
deba@1024:         _skip_nodes(false), _skip_edges(false) {
deba@1024:       if (!(*_is)) {
deba@1024:         delete _is;
deba@1024:         throw IoError("Cannot open file", fn);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Destructor
deba@1024:     ~BpGraphReader() {
deba@1030:       for (typename RedNodeMaps::iterator it = _red_node_maps.begin();
deba@1030:            it != _red_node_maps.end(); ++it) {
deba@1024:         delete it->second;
deba@1024:       }
deba@1024: 
deba@1030:       for (typename BlueNodeMaps::iterator it = _blue_node_maps.begin();
deba@1030:            it != _blue_node_maps.end(); ++it) {
deba@1024:         delete it->second;
deba@1024:       }
deba@1024: 
deba@1024:       for (typename EdgeMaps::iterator it = _edge_maps.begin();
deba@1024:            it != _edge_maps.end(); ++it) {
deba@1024:         delete it->second;
deba@1024:       }
deba@1024: 
deba@1024:       for (typename Attributes::iterator it = _attributes.begin();
deba@1024:            it != _attributes.end(); ++it) {
deba@1024:         delete it->second;
deba@1024:       }
deba@1024: 
deba@1024:       if (local_is) {
deba@1024:         delete _is;
deba@1024:       }
deba@1024: 
deba@1024:     }
deba@1024: 
deba@1024:   private:
deba@1024:     template <typename TBGR>
deba@1024:     friend BpGraphReader<TBGR> bpGraphReader(TBGR& graph, std::istream& is);
deba@1024:     template <typename TBGR>
deba@1024:     friend BpGraphReader<TBGR> bpGraphReader(TBGR& graph,
deba@1024:                                              const std::string& fn);
deba@1024:     template <typename TBGR>
deba@1024:     friend BpGraphReader<TBGR> bpGraphReader(TBGR& graph, const char *fn);
deba@1024: 
deba@1024:     BpGraphReader(BpGraphReader& other)
deba@1024:       : _is(other._is), local_is(other.local_is), _graph(other._graph),
deba@1024:         _use_nodes(other._use_nodes), _use_edges(other._use_edges),
deba@1024:         _skip_nodes(other._skip_nodes), _skip_edges(other._skip_edges) {
deba@1024: 
deba@1024:       other._is = 0;
deba@1024:       other.local_is = false;
deba@1024: 
deba@1030:       _red_node_index.swap(other._red_node_index);
deba@1030:       _blue_node_index.swap(other._blue_node_index);
deba@1024:       _edge_index.swap(other._edge_index);
deba@1024: 
deba@1030:       _red_node_maps.swap(other._red_node_maps);
deba@1030:       _blue_node_maps.swap(other._blue_node_maps);
deba@1024:       _edge_maps.swap(other._edge_maps);
deba@1024:       _attributes.swap(other._attributes);
deba@1024: 
deba@1024:       _nodes_caption = other._nodes_caption;
deba@1024:       _edges_caption = other._edges_caption;
deba@1024:       _attributes_caption = other._attributes_caption;
deba@1024: 
deba@1024:     }
deba@1024: 
deba@1024:     BpGraphReader& operator=(const BpGraphReader&);
deba@1024: 
deba@1024:   public:
deba@1024: 
deba@1024:     /// \name Reading Rules
deba@1024:     /// @{
deba@1024: 
deba@1024:     /// \brief Node map reading rule
deba@1024:     ///
deba@1024:     /// Add a node map reading rule to the reader.
deba@1024:     template <typename Map>
deba@1024:     BpGraphReader& nodeMap(const std::string& caption, Map& map) {
deba@1024:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<RedNode>* red_storage =
deba@1030:         new _reader_bits::MapStorage<RedNode, Map>(map);
deba@1030:       _red_node_maps.push_back(std::make_pair(caption, red_storage));
deba@1030:       _reader_bits::MapStorageBase<BlueNode>* blue_storage =
deba@1030:         new _reader_bits::MapStorage<BlueNode, Map>(map);
deba@1030:       _blue_node_maps.push_back(std::make_pair(caption, blue_storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Node map reading rule
deba@1024:     ///
deba@1024:     /// Add a node map reading rule with specialized converter to the
deba@1024:     /// reader.
deba@1024:     template <typename Map, typename Converter>
deba@1024:     BpGraphReader& nodeMap(const std::string& caption, Map& map,
deba@1024:                            const Converter& converter = Converter()) {
deba@1024:       checkConcept<concepts::WriteMap<Node, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<RedNode>* red_storage =
deba@1030:         new _reader_bits::MapStorage<RedNode, Map, Converter>(map, converter);
deba@1030:       _red_node_maps.push_back(std::make_pair(caption, red_storage));
deba@1030:       _reader_bits::MapStorageBase<BlueNode>* blue_storage =
deba@1030:         new _reader_bits::MapStorage<BlueNode, Map, Converter>(map, converter);
deba@1030:       _blue_node_maps.push_back(std::make_pair(caption, blue_storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1026:     /// Add a red node map reading rule to the reader.
deba@1024:     template <typename Map>
deba@1026:     BpGraphReader& redNodeMap(const std::string& caption, Map& map) {
deba@1030:       checkConcept<concepts::WriteMap<RedNode, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<RedNode>* storage =
deba@1030:         new _reader_bits::MapStorage<RedNode, Map>(map);
deba@1030:       _red_node_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1026:     /// \brief Red node map reading rule
deba@1024:     ///
deba@1026:     /// Add a red node map node reading rule with specialized converter to
deba@1026:     /// the reader.
deba@1024:     template <typename Map, typename Converter>
deba@1026:     BpGraphReader& redNodeMap(const std::string& caption, Map& map,
deba@1026:                               const Converter& converter = Converter()) {
deba@1030:       checkConcept<concepts::WriteMap<RedNode, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<RedNode>* storage =
deba@1030:         new _reader_bits::MapStorage<RedNode, Map, Converter>(map, converter);
deba@1030:       _red_node_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1026:     /// Add a blue node map reading rule to the reader.
deba@1024:     template <typename Map>
deba@1026:     BpGraphReader& blueNodeMap(const std::string& caption, Map& map) {
deba@1030:       checkConcept<concepts::WriteMap<BlueNode, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<BlueNode>* storage =
deba@1030:         new _reader_bits::MapStorage<BlueNode, Map>(map);
deba@1030:       _blue_node_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1026:     /// \brief Blue node map reading rule
deba@1024:     ///
deba@1026:     /// Add a blue node map reading rule with specialized converter to
deba@1026:     /// the reader.
deba@1024:     template <typename Map, typename Converter>
deba@1026:     BpGraphReader& blueNodeMap(const std::string& caption, Map& map,
deba@1026:                                const Converter& converter = Converter()) {
deba@1030:       checkConcept<concepts::WriteMap<BlueNode, typename Map::Value>, Map>();
deba@1030:       _reader_bits::MapStorageBase<BlueNode>* storage =
deba@1030:         new _reader_bits::MapStorage<BlueNode, Map, Converter>(map, converter);
deba@1030:       _blue_node_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Edge map reading rule
deba@1024:     ///
deba@1024:     /// Add an edge map reading rule to the reader.
deba@1024:     template <typename Map>
deba@1024:     BpGraphReader& edgeMap(const std::string& caption, Map& map) {
deba@1024:       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();
deba@1024:       _reader_bits::MapStorageBase<Edge>* storage =
deba@1024:         new _reader_bits::MapStorage<Edge, Map>(map);
deba@1024:       _edge_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Edge map reading rule
deba@1024:     ///
deba@1024:     /// Add an edge map reading rule with specialized converter to the
deba@1024:     /// reader.
deba@1024:     template <typename Map, typename Converter>
deba@1024:     BpGraphReader& edgeMap(const std::string& caption, Map& map,
deba@1024:                           const Converter& converter = Converter()) {
deba@1024:       checkConcept<concepts::WriteMap<Edge, typename Map::Value>, Map>();
deba@1024:       _reader_bits::MapStorageBase<Edge>* storage =
deba@1024:         new _reader_bits::MapStorage<Edge, Map, Converter>(map, converter);
deba@1024:       _edge_maps.push_back(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Arc map reading rule
deba@1024:     ///
deba@1024:     /// Add an arc map reading rule to the reader.
deba@1024:     template <typename Map>
deba@1024:     BpGraphReader& arcMap(const std::string& caption, Map& map) {
deba@1024:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
deba@1024:       _reader_bits::MapStorageBase<Edge>* forward_storage =
deba@1024:         new _reader_bits::GraphArcMapStorage<Graph, true, Map>(_graph, map);
deba@1024:       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
deba@1024:       _reader_bits::MapStorageBase<Edge>* backward_storage =
deba@1024:         new _reader_bits::GraphArcMapStorage<BGR, false, Map>(_graph, map);
deba@1024:       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Arc map reading rule
deba@1024:     ///
deba@1024:     /// Add an arc map reading rule with specialized converter to the
deba@1024:     /// reader.
deba@1024:     template <typename Map, typename Converter>
deba@1024:     BpGraphReader& arcMap(const std::string& caption, Map& map,
deba@1024:                           const Converter& converter = Converter()) {
deba@1024:       checkConcept<concepts::WriteMap<Arc, typename Map::Value>, Map>();
deba@1024:       _reader_bits::MapStorageBase<Edge>* forward_storage =
deba@1024:         new _reader_bits::GraphArcMapStorage<BGR, true, Map, Converter>
deba@1024:         (_graph, map, converter);
deba@1024:       _edge_maps.push_back(std::make_pair('+' + caption, forward_storage));
deba@1024:       _reader_bits::MapStorageBase<Edge>* backward_storage =
deba@1024:         new _reader_bits::GraphArcMapStorage<BGR, false, Map, Converter>
deba@1024:         (_graph, map, converter);
deba@1024:       _edge_maps.push_back(std::make_pair('-' + caption, backward_storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Attribute reading rule
deba@1024:     ///
deba@1024:     /// Add an attribute reading rule to the reader.
deba@1024:     template <typename Value>
deba@1024:     BpGraphReader& attribute(const std::string& caption, Value& value) {
deba@1024:       _reader_bits::ValueStorageBase* storage =
deba@1024:         new _reader_bits::ValueStorage<Value>(value);
deba@1024:       _attributes.insert(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Attribute reading rule
deba@1024:     ///
deba@1024:     /// Add an attribute reading rule with specialized converter to the
deba@1024:     /// reader.
deba@1024:     template <typename Value, typename Converter>
deba@1024:     BpGraphReader& attribute(const std::string& caption, Value& value,
deba@1024:                              const Converter& converter = Converter()) {
deba@1024:       _reader_bits::ValueStorageBase* storage =
deba@1024:         new _reader_bits::ValueStorage<Value, Converter>(value, converter);
deba@1024:       _attributes.insert(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Node reading rule
deba@1024:     ///
deba@1024:     /// Add a node reading rule to reader.
deba@1024:     BpGraphReader& node(const std::string& caption, Node& node) {
deba@1030:       typedef _reader_bits::DoubleMapLookUpConverter<
deba@1030:         Node, RedNodeIndex, BlueNodeIndex> Converter;
deba@1030:       Converter converter(_red_node_index, _blue_node_index);
deba@1024:       _reader_bits::ValueStorageBase* storage =
deba@1024:         new _reader_bits::ValueStorage<Node, Converter>(node, converter);
deba@1024:       _attributes.insert(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1030:     /// \brief Red node reading rule
deba@1030:     ///
deba@1030:     /// Add a red node reading rule to reader.
deba@1030:     BpGraphReader& redNode(const std::string& caption, RedNode& node) {
deba@1030:       typedef _reader_bits::MapLookUpConverter<RedNode> Converter;
deba@1030:       Converter converter(_red_node_index);
deba@1030:       _reader_bits::ValueStorageBase* storage =
deba@1030:         new _reader_bits::ValueStorage<RedNode, Converter>(node, converter);
deba@1030:       _attributes.insert(std::make_pair(caption, storage));
deba@1030:       return *this;
deba@1030:     }
deba@1030: 
deba@1030:     /// \brief Blue node reading rule
deba@1030:     ///
deba@1030:     /// Add a blue node reading rule to reader.
deba@1030:     BpGraphReader& blueNode(const std::string& caption, BlueNode& node) {
deba@1030:       typedef _reader_bits::MapLookUpConverter<BlueNode> Converter;
deba@1030:       Converter converter(_blue_node_index);
deba@1030:       _reader_bits::ValueStorageBase* storage =
deba@1030:         new _reader_bits::ValueStorage<BlueNode, Converter>(node, converter);
deba@1030:       _attributes.insert(std::make_pair(caption, storage));
deba@1030:       return *this;
deba@1030:     }
deba@1030: 
deba@1024:     /// \brief Edge reading rule
deba@1024:     ///
deba@1024:     /// Add an edge reading rule to reader.
deba@1024:     BpGraphReader& edge(const std::string& caption, Edge& edge) {
deba@1024:       typedef _reader_bits::MapLookUpConverter<Edge> Converter;
deba@1024:       Converter converter(_edge_index);
deba@1024:       _reader_bits::ValueStorageBase* storage =
deba@1024:         new _reader_bits::ValueStorage<Edge, Converter>(edge, converter);
deba@1024:       _attributes.insert(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Arc reading rule
deba@1024:     ///
deba@1024:     /// Add an arc reading rule to reader.
deba@1024:     BpGraphReader& arc(const std::string& caption, Arc& arc) {
deba@1024:       typedef _reader_bits::GraphArcLookUpConverter<BGR> Converter;
deba@1024:       Converter converter(_graph, _edge_index);
deba@1024:       _reader_bits::ValueStorageBase* storage =
deba@1024:         new _reader_bits::ValueStorage<Arc, Converter>(arc, converter);
deba@1024:       _attributes.insert(std::make_pair(caption, storage));
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// @}
deba@1024: 
deba@1024:     /// \name Select Section by Name
deba@1024:     /// @{
deba@1024: 
deba@1024:     /// \brief Set \c \@nodes section to be read
deba@1024:     ///
deba@1024:     /// Set \c \@nodes section to be read.
deba@1024:     BpGraphReader& nodes(const std::string& caption) {
deba@1024:       _nodes_caption = caption;
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Set \c \@edges section to be read
deba@1024:     ///
deba@1024:     /// Set \c \@edges section to be read.
deba@1024:     BpGraphReader& edges(const std::string& caption) {
deba@1024:       _edges_caption = caption;
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Set \c \@attributes section to be read
deba@1024:     ///
deba@1024:     /// Set \c \@attributes section to be read.
deba@1024:     BpGraphReader& attributes(const std::string& caption) {
deba@1024:       _attributes_caption = caption;
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// @}
deba@1024: 
deba@1024:     /// \name Using Previously Constructed Node or Edge Set
deba@1024:     /// @{
deba@1024: 
deba@1024:     /// \brief Use previously constructed node set
deba@1024:     ///
deba@1024:     /// Use previously constructed node set, and specify the node
deba@1024:     /// label map.
deba@1024:     template <typename Map>
deba@1024:     BpGraphReader& useNodes(const Map& map) {
deba@1024:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@1024:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@1024:       _use_nodes = true;
deba@1024:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@1030:       for (RedNodeIt n(_graph); n != INVALID; ++n) {
deba@1030:         _red_node_index.insert(std::make_pair(converter(map[n]), n));
deba@1030:       }
deba@1030:       for (BlueNodeIt n(_graph); n != INVALID; ++n) {
deba@1030:         _blue_node_index.insert(std::make_pair(converter(map[n]), n));
deba@1024:       }
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Use previously constructed node set
deba@1024:     ///
deba@1024:     /// Use previously constructed node set, and specify the node
deba@1024:     /// label map and a functor which converts the label map values to
deba@1024:     /// \c std::string.
deba@1024:     template <typename Map, typename Converter>
deba@1024:     BpGraphReader& useNodes(const Map& map,
deba@1024:                             const Converter& converter = Converter()) {
deba@1024:       checkConcept<concepts::ReadMap<Node, typename Map::Value>, Map>();
deba@1024:       LEMON_ASSERT(!_use_nodes, "Multiple usage of useNodes() member");
deba@1024:       _use_nodes = true;
deba@1030:       for (RedNodeIt n(_graph); n != INVALID; ++n) {
deba@1030:         _red_node_index.insert(std::make_pair(converter(map[n]), n));
deba@1030:       }
deba@1030:       for (BlueNodeIt n(_graph); n != INVALID; ++n) {      
deba@1030:         _blue_node_index.insert(std::make_pair(converter(map[n]), n));
deba@1024:       }
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Use previously constructed edge set
deba@1024:     ///
deba@1024:     /// Use previously constructed edge set, and specify the edge
deba@1024:     /// label map.
deba@1024:     template <typename Map>
deba@1024:     BpGraphReader& useEdges(const Map& map) {
deba@1024:       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
deba@1024:       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");
deba@1024:       _use_edges = true;
deba@1024:       _writer_bits::DefaultConverter<typename Map::Value> converter;
deba@1024:       for (EdgeIt a(_graph); a != INVALID; ++a) {
deba@1024:         _edge_index.insert(std::make_pair(converter(map[a]), a));
deba@1024:       }
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Use previously constructed edge set
deba@1024:     ///
deba@1024:     /// Use previously constructed edge set, and specify the edge
deba@1024:     /// label map and a functor which converts the label map values to
deba@1024:     /// \c std::string.
deba@1024:     template <typename Map, typename Converter>
deba@1024:     BpGraphReader& useEdges(const Map& map,
deba@1024:                             const Converter& converter = Converter()) {
deba@1024:       checkConcept<concepts::ReadMap<Edge, typename Map::Value>, Map>();
deba@1024:       LEMON_ASSERT(!_use_edges, "Multiple usage of useEdges() member");
deba@1024:       _use_edges = true;
deba@1024:       for (EdgeIt a(_graph); a != INVALID; ++a) {
deba@1024:         _edge_index.insert(std::make_pair(converter(map[a]), a));
deba@1024:       }
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Skip the reading of node section
deba@1024:     ///
deba@1024:     /// Omit the reading of the node section. This implies that each node
deba@1024:     /// map reading rule will be abandoned, and the nodes of the graph
deba@1024:     /// will not be constructed, which usually cause that the edge set
deba@1024:     /// could not be read due to lack of node name
deba@1024:     /// could not be read due to lack of node name resolving.
deba@1024:     /// Therefore \c skipEdges() function should also be used, or
deba@1024:     /// \c useNodes() should be used to specify the label of the nodes.
deba@1024:     BpGraphReader& skipNodes() {
deba@1024:       LEMON_ASSERT(!_skip_nodes, "Skip nodes already set");
deba@1024:       _skip_nodes = true;
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// \brief Skip the reading of edge section
deba@1024:     ///
deba@1024:     /// Omit the reading of the edge section. This implies that each edge
deba@1024:     /// map reading rule will be abandoned, and the edges of the graph
deba@1024:     /// will not be constructed.
deba@1024:     BpGraphReader& skipEdges() {
deba@1024:       LEMON_ASSERT(!_skip_edges, "Skip edges already set");
deba@1024:       _skip_edges = true;
deba@1024:       return *this;
deba@1024:     }
deba@1024: 
deba@1024:     /// @}
deba@1024: 
deba@1024:   private:
deba@1024: 
deba@1024:     bool readLine() {
deba@1024:       std::string str;
deba@1024:       while(++line_num, std::getline(*_is, str)) {
deba@1024:         line.clear(); line.str(str);
deba@1024:         char c;
deba@1024:         if (line >> std::ws >> c && c != '#') {
deba@1024:           line.putback(c);
deba@1024:           return true;
deba@1024:         }
deba@1024:       }
deba@1024:       return false;
deba@1024:     }
deba@1024: 
deba@1024:     bool readSuccess() {
deba@1024:       return static_cast<bool>(*_is);
deba@1024:     }
deba@1024: 
deba@1024:     void skipSection() {
deba@1024:       char c;
deba@1024:       while (readSuccess() && line >> c && c != '@') {
deba@1024:         readLine();
deba@1024:       }
deba@1024:       if (readSuccess()) {
deba@1024:         line.putback(c);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     void readRedNodes() {
deba@1024: 
deba@1030:       std::vector<int> map_index(_red_node_maps.size());
deba@1024:       int map_num, label_index;
deba@1024: 
deba@1024:       char c;
deba@1024:       if (!readLine() || !(line >> c) || c == '@') {
deba@1024:         if (readSuccess() && line) line.putback(c);
deba@1030:         if (!_red_node_maps.empty())
deba@1024:           throw FormatError("Cannot find map names");
deba@1024:         return;
deba@1024:       }
deba@1024:       line.putback(c);
deba@1024: 
deba@1024:       {
deba@1024:         std::map<std::string, int> maps;
deba@1024: 
deba@1024:         std::string map;
deba@1024:         int index = 0;
deba@1024:         while (_reader_bits::readToken(line, map)) {
deba@1024:           if (maps.find(map) != maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1026:             msg << "Multiple occurence of red node map: " << map;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           maps.insert(std::make_pair(map, index));
deba@1024:           ++index;
deba@1024:         }
deba@1024: 
deba@1030:         for (int i = 0; i < static_cast<int>(_red_node_maps.size()); ++i) {
deba@1024:           std::map<std::string, int>::iterator jt =
deba@1030:             maps.find(_red_node_maps[i].first);
deba@1024:           if (jt == maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1030:             msg << "Map not found: " << _red_node_maps[i].first;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           map_index[i] = jt->second;
deba@1024:         }
deba@1024: 
deba@1024:         {
deba@1024:           std::map<std::string, int>::iterator jt = maps.find("label");
deba@1024:           if (jt != maps.end()) {
deba@1024:             label_index = jt->second;
deba@1024:           } else {
deba@1024:             label_index = -1;
deba@1024:           }
deba@1024:         }
deba@1024:         map_num = maps.size();
deba@1024:       }
deba@1024: 
deba@1024:       while (readLine() && line >> c && c != '@') {
deba@1024:         line.putback(c);
deba@1024: 
deba@1024:         std::vector<std::string> tokens(map_num);
deba@1024:         for (int i = 0; i < map_num; ++i) {
deba@1024:           if (!_reader_bits::readToken(line, tokens[i])) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Column not found (" << i + 1 << ")";
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:         }
deba@1024:         if (line >> std::ws >> c)
deba@1024:           throw FormatError("Extra character at the end of line");
deba@1024: 
deba@1030:         RedNode n;
deba@1024:         if (!_use_nodes) {
deba@1024:           n = _graph.addRedNode();
deba@1024:           if (label_index != -1)
deba@1030:             _red_node_index.insert(std::make_pair(tokens[label_index], n));
deba@1024:         } else {
deba@1024:           if (label_index == -1)
deba@1024:             throw FormatError("Label map not found");
deba@1030:           typename std::map<std::string, RedNode>::iterator it =
deba@1030:             _red_node_index.find(tokens[label_index]);
deba@1030:           if (it == _red_node_index.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Node with label not found: " << tokens[label_index];
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           n = it->second;
deba@1024:         }
deba@1024: 
deba@1030:         for (int i = 0; i < static_cast<int>(_red_node_maps.size()); ++i) {
deba@1030:           _red_node_maps[i].second->set(n, tokens[map_index[i]]);
deba@1024:         }
deba@1024: 
deba@1024:       }
deba@1024:       if (readSuccess()) {
deba@1024:         line.putback(c);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     void readBlueNodes() {
deba@1024: 
deba@1030:       std::vector<int> map_index(_blue_node_maps.size());
deba@1024:       int map_num, label_index;
deba@1024: 
deba@1024:       char c;
deba@1024:       if (!readLine() || !(line >> c) || c == '@') {
deba@1024:         if (readSuccess() && line) line.putback(c);
deba@1030:         if (!_blue_node_maps.empty())
deba@1024:           throw FormatError("Cannot find map names");
deba@1024:         return;
deba@1024:       }
deba@1024:       line.putback(c);
deba@1024: 
deba@1024:       {
deba@1024:         std::map<std::string, int> maps;
deba@1024: 
deba@1024:         std::string map;
deba@1024:         int index = 0;
deba@1024:         while (_reader_bits::readToken(line, map)) {
deba@1024:           if (maps.find(map) != maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1026:             msg << "Multiple occurence of blue node map: " << map;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           maps.insert(std::make_pair(map, index));
deba@1024:           ++index;
deba@1024:         }
deba@1024: 
deba@1030:         for (int i = 0; i < static_cast<int>(_blue_node_maps.size()); ++i) {
deba@1024:           std::map<std::string, int>::iterator jt =
deba@1030:             maps.find(_blue_node_maps[i].first);
deba@1024:           if (jt == maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1030:             msg << "Map not found: " << _blue_node_maps[i].first;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           map_index[i] = jt->second;
deba@1024:         }
deba@1024: 
deba@1024:         {
deba@1024:           std::map<std::string, int>::iterator jt = maps.find("label");
deba@1024:           if (jt != maps.end()) {
deba@1024:             label_index = jt->second;
deba@1024:           } else {
deba@1024:             label_index = -1;
deba@1024:           }
deba@1024:         }
deba@1024:         map_num = maps.size();
deba@1024:       }
deba@1024: 
deba@1024:       while (readLine() && line >> c && c != '@') {
deba@1024:         line.putback(c);
deba@1024: 
deba@1024:         std::vector<std::string> tokens(map_num);
deba@1024:         for (int i = 0; i < map_num; ++i) {
deba@1024:           if (!_reader_bits::readToken(line, tokens[i])) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Column not found (" << i + 1 << ")";
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:         }
deba@1024:         if (line >> std::ws >> c)
deba@1024:           throw FormatError("Extra character at the end of line");
deba@1024: 
deba@1030:         BlueNode n;
deba@1024:         if (!_use_nodes) {
deba@1024:           n = _graph.addBlueNode();
deba@1024:           if (label_index != -1)
deba@1030:             _blue_node_index.insert(std::make_pair(tokens[label_index], n));
deba@1024:         } else {
deba@1024:           if (label_index == -1)
deba@1024:             throw FormatError("Label map not found");
deba@1030:           typename std::map<std::string, BlueNode>::iterator it =
deba@1030:             _blue_node_index.find(tokens[label_index]);
deba@1030:           if (it == _blue_node_index.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Node with label not found: " << tokens[label_index];
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           n = it->second;
deba@1024:         }
deba@1024: 
deba@1030:         for (int i = 0; i < static_cast<int>(_blue_node_maps.size()); ++i) {
deba@1030:           _blue_node_maps[i].second->set(n, tokens[map_index[i]]);
deba@1024:         }
deba@1024: 
deba@1024:       }
deba@1024:       if (readSuccess()) {
deba@1024:         line.putback(c);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     void readEdges() {
deba@1024: 
deba@1024:       std::vector<int> map_index(_edge_maps.size());
deba@1024:       int map_num, label_index;
deba@1024: 
deba@1024:       char c;
deba@1024:       if (!readLine() || !(line >> c) || c == '@') {
deba@1024:         if (readSuccess() && line) line.putback(c);
deba@1024:         if (!_edge_maps.empty())
deba@1024:           throw FormatError("Cannot find map names");
deba@1024:         return;
deba@1024:       }
deba@1024:       line.putback(c);
deba@1024: 
deba@1024:       {
deba@1024:         std::map<std::string, int> maps;
deba@1024: 
deba@1024:         std::string map;
deba@1024:         int index = 0;
deba@1024:         while (_reader_bits::readToken(line, map)) {
deba@1024:           if (maps.find(map) != maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Multiple occurence of edge map: " << map;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           maps.insert(std::make_pair(map, index));
deba@1024:           ++index;
deba@1024:         }
deba@1024: 
deba@1024:         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {
deba@1024:           std::map<std::string, int>::iterator jt =
deba@1024:             maps.find(_edge_maps[i].first);
deba@1024:           if (jt == maps.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Map not found: " << _edge_maps[i].first;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           map_index[i] = jt->second;
deba@1024:         }
deba@1024: 
deba@1024:         {
deba@1024:           std::map<std::string, int>::iterator jt = maps.find("label");
deba@1024:           if (jt != maps.end()) {
deba@1024:             label_index = jt->second;
deba@1024:           } else {
deba@1024:             label_index = -1;
deba@1024:           }
deba@1024:         }
deba@1024:         map_num = maps.size();
deba@1024:       }
deba@1024: 
deba@1024:       while (readLine() && line >> c && c != '@') {
deba@1024:         line.putback(c);
deba@1024: 
deba@1024:         std::string source_token;
deba@1024:         std::string target_token;
deba@1024: 
deba@1024:         if (!_reader_bits::readToken(line, source_token))
deba@1024:           throw FormatError("Red node not found");
deba@1024: 
deba@1024:         if (!_reader_bits::readToken(line, target_token))
deba@1024:           throw FormatError("Blue node not found");
deba@1024: 
deba@1024:         std::vector<std::string> tokens(map_num);
deba@1024:         for (int i = 0; i < map_num; ++i) {
deba@1024:           if (!_reader_bits::readToken(line, tokens[i])) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Column not found (" << i + 1 << ")";
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:         }
deba@1024:         if (line >> std::ws >> c)
deba@1024:           throw FormatError("Extra character at the end of line");
deba@1024: 
deba@1024:         Edge e;
deba@1024:         if (!_use_edges) {
deba@1030:           typename RedNodeIndex::iterator rit =
deba@1030:             _red_node_index.find(source_token);
deba@1030:           if (rit == _red_node_index.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Item not found: " << source_token;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1030:           RedNode source = rit->second;
deba@1030:           typename BlueNodeIndex::iterator it =
deba@1030:             _blue_node_index.find(target_token);
deba@1030:           if (it == _blue_node_index.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Item not found: " << target_token;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1030:           BlueNode target = it->second;
deba@1024: 
poroszd@1029:           // It is checked that source is red and
poroszd@1029:           // target is blue, so this should be safe:
deba@1030:           e = _graph.addEdge(source, target);
deba@1024:           if (label_index != -1)
deba@1024:             _edge_index.insert(std::make_pair(tokens[label_index], e));
deba@1024:         } else {
deba@1024:           if (label_index == -1)
deba@1024:             throw FormatError("Label map not found");
deba@1024:           typename std::map<std::string, Edge>::iterator it =
deba@1024:             _edge_index.find(tokens[label_index]);
deba@1024:           if (it == _edge_index.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Edge with label not found: " << tokens[label_index];
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           e = it->second;
deba@1024:         }
deba@1024: 
deba@1024:         for (int i = 0; i < static_cast<int>(_edge_maps.size()); ++i) {
deba@1024:           _edge_maps[i].second->set(e, tokens[map_index[i]]);
deba@1024:         }
deba@1024: 
deba@1024:       }
deba@1024:       if (readSuccess()) {
deba@1024:         line.putback(c);
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:     void readAttributes() {
deba@1024: 
deba@1024:       std::set<std::string> read_attr;
deba@1024: 
deba@1024:       char c;
deba@1024:       while (readLine() && line >> c && c != '@') {
deba@1024:         line.putback(c);
deba@1024: 
deba@1024:         std::string attr, token;
deba@1024:         if (!_reader_bits::readToken(line, attr))
deba@1024:           throw FormatError("Attribute name not found");
deba@1024:         if (!_reader_bits::readToken(line, token))
deba@1024:           throw FormatError("Attribute value not found");
deba@1024:         if (line >> c)
deba@1024:           throw FormatError("Extra character at the end of line");
deba@1024: 
deba@1024:         {
deba@1024:           std::set<std::string>::iterator it = read_attr.find(attr);
deba@1024:           if (it != read_attr.end()) {
deba@1024:             std::ostringstream msg;
deba@1024:             msg << "Multiple occurence of attribute: " << attr;
deba@1024:             throw FormatError(msg.str());
deba@1024:           }
deba@1024:           read_attr.insert(attr);
deba@1024:         }
deba@1024: 
deba@1024:         {
deba@1024:           typename Attributes::iterator it = _attributes.lower_bound(attr);
deba@1024:           while (it != _attributes.end() && it->first == attr) {
deba@1024:             it->second->set(token);
deba@1024:             ++it;
deba@1024:           }
deba@1024:         }
deba@1024: 
deba@1024:       }
deba@1024:       if (readSuccess()) {
deba@1024:         line.putback(c);
deba@1024:       }
deba@1024:       for (typename Attributes::iterator it = _attributes.begin();
deba@1024:            it != _attributes.end(); ++it) {
deba@1024:         if (read_attr.find(it->first) == read_attr.end()) {
deba@1024:           std::ostringstream msg;
deba@1024:           msg << "Attribute not found: " << it->first;
deba@1024:           throw FormatError(msg.str());
deba@1024:         }
deba@1024:       }
deba@1024:     }
deba@1024: 
deba@1024:   public:
deba@1024: 
deba@1024:     /// \name Execution of the Reader
deba@1024:     /// @{
deba@1024: 
deba@1024:     /// \brief Start the batch processing
deba@1024:     ///
deba@1024:     /// This function starts the batch processing
deba@1024:     void run() {
deba@1024: 
deba@1024:       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");
deba@1024: 
deba@1024:       bool red_nodes_done = _skip_nodes;
deba@1024:       bool blue_nodes_done = _skip_nodes;
deba@1024:       bool edges_done = _skip_edges;
deba@1024:       bool attributes_done = false;
deba@1024: 
deba@1024:       line_num = 0;
deba@1024:       readLine();
deba@1024:       skipSection();
deba@1024: 
deba@1024:       while (readSuccess()) {
deba@1024:         try {
deba@1024:           char c;
deba@1024:           std::string section, caption;
deba@1024:           line >> c;
deba@1024:           _reader_bits::readToken(line, section);
deba@1024:           _reader_bits::readToken(line, caption);
deba@1024: 
deba@1024:           if (line >> c)
deba@1024:             throw FormatError("Extra character at the end of line");
deba@1024: 
deba@1024:           if (section == "red_nodes" && !red_nodes_done) {
deba@1024:             if (_nodes_caption.empty() || _nodes_caption == caption) {
deba@1024:               readRedNodes();
deba@1024:               red_nodes_done = true;
deba@1024:             }
deba@1024:           } else if (section == "blue_nodes" && !blue_nodes_done) {
deba@1024:             if (_nodes_caption.empty() || _nodes_caption == caption) {
deba@1024:               readBlueNodes();
deba@1024:               blue_nodes_done = true;
deba@1024:             }
deba@1024:           } else if ((section == "edges" || section == "arcs") &&
deba@1024:                      !edges_done) {
deba@1024:             if (_edges_caption.empty() || _edges_caption == caption) {
deba@1024:               readEdges();
deba@1024:               edges_done = true;
deba@1024:             }
deba@1024:           } else if (section == "attributes" && !attributes_done) {
deba@1024:             if (_attributes_caption.empty() || _attributes_caption == caption) {
deba@1024:               readAttributes();
deba@1024:               attributes_done = true;
deba@1024:             }
deba@1024:           } else {
deba@1024:             readLine();
deba@1024:             skipSection();
deba@1024:           }
deba@1024:         } catch (FormatError& error) {
deba@1024:           error.line(line_num);
deba@1024:           error.file(_filename);
deba@1024:           throw;
deba@1024:         }
deba@1024:       }
deba@1024: 
deba@1024:       if (!red_nodes_done) {
deba@1024:         throw FormatError("Section @red_nodes not found");
deba@1024:       }
deba@1024: 
deba@1024:       if (!blue_nodes_done) {
deba@1024:         throw FormatError("Section @blue_nodes not found");
deba@1024:       }
deba@1024: 
deba@1024:       if (!edges_done) {
deba@1024:         throw FormatError("Section @edges not found");
deba@1024:       }
deba@1024: 
deba@1024:       if (!attributes_done && !_attributes.empty()) {
deba@1024:         throw FormatError("Section @attributes not found");
deba@1024:       }
deba@1024: 
deba@1024:     }
deba@1024: 
deba@1024:     /// @}
deba@1024: 
deba@1024:   };
deba@1024: 
deba@1024:   /// \ingroup lemon_io
deba@1024:   ///
deba@1024:   /// \brief Return a \ref BpGraphReader class
deba@1024:   ///
deba@1024:   /// This function just returns a \ref BpGraphReader class.
deba@1024:   ///
deba@1024:   /// With this function a graph can be read from an
deba@1024:   /// \ref lgf-format "LGF" file or input stream with several maps and
deba@1024:   /// attributes. For example, there is bipartite weighted matching problem
deba@1024:   /// on a graph, i.e. a graph with a \e weight map on the edges. This
deba@1024:   /// graph can be read with the following code:
deba@1024:   ///
deba@1024:   ///\code
deba@1024:   ///ListBpGraph graph;
deba@1024:   ///ListBpGraph::EdgeMap<int> weight(graph);
deba@1024:   ///bpGraphReader(graph, std::cin).
deba@1024:   ///  edgeMap("weight", weight).
deba@1024:   ///  run();
deba@1024:   ///\endcode
deba@1024:   ///
deba@1024:   /// For a complete documentation, please see the \ref BpGraphReader
deba@1024:   /// class documentation.
deba@1024:   /// \warning Don't forget to put the \ref BpGraphReader::run() "run()"
deba@1024:   /// to the end of the parameter list.
deba@1024:   /// \relates BpGraphReader
deba@1024:   /// \sa bpGraphReader(TBGR& graph, const std::string& fn)
deba@1024:   /// \sa bpGraphReader(TBGR& graph, const char* fn)
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, std::istream& is) {
deba@1024:     BpGraphReader<TBGR> tmp(graph, is);
deba@1024:     return tmp;
deba@1024:   }
deba@1024: 
deba@1024:   /// \brief Return a \ref BpGraphReader class
deba@1024:   ///
deba@1024:   /// This function just returns a \ref BpGraphReader class.
deba@1024:   /// \relates BpGraphReader
deba@1024:   /// \sa bpGraphReader(TBGR& graph, std::istream& is)
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, const std::string& fn) {
deba@1024:     BpGraphReader<TBGR> tmp(graph, fn);
deba@1024:     return tmp;
deba@1024:   }
deba@1024: 
deba@1024:   /// \brief Return a \ref BpGraphReader class
deba@1024:   ///
deba@1024:   /// This function just returns a \ref BpGraphReader class.
deba@1024:   /// \relates BpGraphReader
deba@1024:   /// \sa bpGraphReader(TBGR& graph, std::istream& is)
deba@1024:   template <typename TBGR>
deba@1024:   BpGraphReader<TBGR> bpGraphReader(TBGR& graph, const char* fn) {
deba@1024:     BpGraphReader<TBGR> tmp(graph, fn);
deba@1024:     return tmp;
deba@1024:   }
deba@1024: 
deba@190:   class SectionReader;
deba@190: 
deba@190:   SectionReader sectionReader(std::istream& is);
deba@190:   SectionReader sectionReader(const std::string& fn);
deba@190:   SectionReader sectionReader(const char* fn);
alpar@209: 
kpeter@192:   /// \ingroup lemon_io
kpeter@192:   ///
deba@189:   /// \brief Section reader class
deba@189:   ///
alpar@209:   /// In the \ref lgf-format "LGF" file extra sections can be placed,
kpeter@192:   /// which contain any data in arbitrary format. Such sections can be
alpar@209:   /// read with this class. A reading rule can be added to the class
kpeter@192:   /// with two different functions. With the \c sectionLines() function a
kpeter@192:   /// functor can process the section line-by-line, while with the \c
deba@189:   /// sectionStream() member the section can be read from an input
deba@189:   /// stream.
deba@189:   class SectionReader {
deba@189:   private:
alpar@209: 
deba@189:     std::istream* _is;
deba@189:     bool local_is;
deba@290:     std::string _filename;
deba@189: 
deba@189:     typedef std::map<std::string, _reader_bits::Section*> Sections;
deba@189:     Sections _sections;
deba@189: 
deba@189:     int line_num;
deba@189:     std::istringstream line;
deba@189: 
deba@189:   public:
deba@189: 
deba@189:     /// \brief Constructor
deba@189:     ///
deba@189:     /// Construct a section reader, which reads from the given input
deba@189:     /// stream.
alpar@209:     SectionReader(std::istream& is)
deba@189:       : _is(&is), local_is(false) {}
deba@189: 
deba@189:     /// \brief Constructor
deba@189:     ///
deba@189:     /// Construct a section reader, which reads from the given file.
alpar@209:     SectionReader(const std::string& fn)
deba@290:       : _is(new std::ifstream(fn.c_str())), local_is(true),
deba@290:         _filename(fn) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
alpar@209: 
deba@189:     /// \brief Constructor
deba@189:     ///
deba@189:     /// Construct a section reader, which reads from the given file.
alpar@209:     SectionReader(const char* fn)
deba@290:       : _is(new std::ifstream(fn)), local_is(true),
deba@290:         _filename(fn) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
deba@189: 
deba@189:     /// \brief Destructor
deba@189:     ~SectionReader() {
alpar@209:       for (Sections::iterator it = _sections.begin();
alpar@209:            it != _sections.end(); ++it) {
alpar@209:         delete it->second;
deba@189:       }
deba@189: 
deba@189:       if (local_is) {
alpar@209:         delete _is;
deba@189:       }
deba@189: 
deba@189:     }
deba@189: 
deba@189:   private:
deba@190: 
deba@190:     friend SectionReader sectionReader(std::istream& is);
deba@190:     friend SectionReader sectionReader(const std::string& fn);
deba@190:     friend SectionReader sectionReader(const char* fn);
deba@190: 
alpar@209:     SectionReader(SectionReader& other)
deba@190:       : _is(other._is), local_is(other.local_is) {
deba@190: 
deba@190:       other._is = 0;
deba@190:       other.local_is = false;
alpar@209: 
deba@190:       _sections.swap(other._sections);
deba@190:     }
alpar@209: 
deba@189:     SectionReader& operator=(const SectionReader&);
deba@189: 
deba@189:   public:
deba@189: 
kpeter@584:     /// \name Section Readers
deba@189:     /// @{
deba@189: 
deba@189:     /// \brief Add a section processor with line oriented reading
deba@189:     ///
deba@189:     /// The first parameter is the type descriptor of the section, the
deba@189:     /// second is a functor, which takes just one \c std::string
deba@189:     /// parameter. At the reading process, each line of the section
deba@189:     /// will be given to the functor object. However, the empty lines
deba@189:     /// and the comment lines are filtered out, and the leading
deba@189:     /// whitespaces are trimmed from each processed string.
deba@189:     ///
kpeter@786:     /// For example, let's see a section, which contain several
deba@189:     /// integers, which should be inserted into a vector.
deba@189:     ///\code
deba@189:     ///  @numbers
deba@189:     ///  12 45 23
deba@189:     ///  4
deba@189:     ///  23 6
deba@189:     ///\endcode
deba@189:     ///
kpeter@192:     /// The functor is implemented as a struct:
deba@189:     ///\code
deba@189:     ///  struct NumberSection {
deba@189:     ///    std::vector<int>& _data;
deba@189:     ///    NumberSection(std::vector<int>& data) : _data(data) {}
deba@189:     ///    void operator()(const std::string& line) {
deba@189:     ///      std::istringstream ls(line);
deba@189:     ///      int value;
deba@189:     ///      while (ls >> value) _data.push_back(value);
deba@189:     ///    }
deba@189:     ///  };
deba@189:     ///
deba@189:     ///  // ...
deba@189:     ///
alpar@209:     ///  reader.sectionLines("numbers", NumberSection(vec));
deba@189:     ///\endcode
deba@189:     template <typename Functor>
deba@189:     SectionReader& sectionLines(const std::string& type, Functor functor) {
kpeter@192:       LEMON_ASSERT(!type.empty(), "Type is empty.");
alpar@209:       LEMON_ASSERT(_sections.find(type) == _sections.end(),
alpar@209:                    "Multiple reading of section.");
alpar@209:       _sections.insert(std::make_pair(type,
deba@189:         new _reader_bits::LineSection<Functor>(functor)));
deba@189:       return *this;
deba@189:     }
deba@189: 
deba@189: 
deba@189:     /// \brief Add a section processor with stream oriented reading
deba@189:     ///
deba@189:     /// The first parameter is the type of the section, the second is
kpeter@192:     /// a functor, which takes an \c std::istream& and an \c int&
deba@189:     /// parameter, the latter regard to the line number of stream. The
deba@189:     /// functor can read the input while the section go on, and the
deba@189:     /// line number should be modified accordingly.
deba@189:     template <typename Functor>
deba@189:     SectionReader& sectionStream(const std::string& type, Functor functor) {
kpeter@192:       LEMON_ASSERT(!type.empty(), "Type is empty.");
alpar@209:       LEMON_ASSERT(_sections.find(type) == _sections.end(),
alpar@209:                    "Multiple reading of section.");
alpar@209:       _sections.insert(std::make_pair(type,
alpar@209:          new _reader_bits::StreamSection<Functor>(functor)));
deba@189:       return *this;
alpar@209:     }
alpar@209: 
deba@189:     /// @}
deba@189: 
deba@189:   private:
deba@189: 
deba@189:     bool readLine() {
deba@189:       std::string str;
deba@189:       while(++line_num, std::getline(*_is, str)) {
alpar@209:         line.clear(); line.str(str);
alpar@209:         char c;
alpar@209:         if (line >> std::ws >> c && c != '#') {
alpar@209:           line.putback(c);
alpar@209:           return true;
alpar@209:         }
deba@189:       }
deba@189:       return false;
deba@189:     }
deba@189: 
deba@189:     bool readSuccess() {
deba@189:       return static_cast<bool>(*_is);
deba@189:     }
alpar@209: 
deba@189:     void skipSection() {
deba@189:       char c;
deba@189:       while (readSuccess() && line >> c && c != '@') {
alpar@209:         readLine();
deba@189:       }
deba@427:       if (readSuccess()) {
deba@427:         line.putback(c);
deba@427:       }
deba@189:     }
deba@189: 
deba@189:   public:
deba@189: 
deba@189: 
kpeter@584:     /// \name Execution of the Reader
deba@189:     /// @{
deba@189: 
deba@189:     /// \brief Start the batch processing
deba@189:     ///
kpeter@192:     /// This function starts the batch processing.
deba@189:     void run() {
alpar@209: 
deba@189:       LEMON_ASSERT(_is != 0, "This reader assigned to an other reader");
alpar@209: 
deba@189:       std::set<std::string> extra_sections;
deba@189: 
alpar@209:       line_num = 0;
deba@189:       readLine();
deba@189:       skipSection();
deba@189: 
deba@189:       while (readSuccess()) {
alpar@209:         try {
alpar@209:           char c;
alpar@209:           std::string section, caption;
alpar@209:           line >> c;
alpar@209:           _reader_bits::readToken(line, section);
alpar@209:           _reader_bits::readToken(line, caption);
alpar@209: 
alpar@209:           if (line >> c)
kpeter@291:             throw FormatError("Extra character at the end of line");
alpar@209: 
alpar@209:           if (extra_sections.find(section) != extra_sections.end()) {
alpar@209:             std::ostringstream msg;
kpeter@291:             msg << "Multiple occurence of section: " << section;
deba@290:             throw FormatError(msg.str());
alpar@209:           }
alpar@209:           Sections::iterator it = _sections.find(section);
alpar@209:           if (it != _sections.end()) {
alpar@209:             extra_sections.insert(section);
alpar@209:             it->second->process(*_is, line_num);
alpar@209:           }
alpar@209:           readLine();
alpar@209:           skipSection();
deba@290:         } catch (FormatError& error) {
alpar@209:           error.line(line_num);
deba@290:           error.file(_filename);
alpar@209:           throw;
alpar@209:         }
deba@189:       }
deba@189:       for (Sections::iterator it = _sections.begin();
alpar@209:            it != _sections.end(); ++it) {
alpar@209:         if (extra_sections.find(it->first) == extra_sections.end()) {
alpar@209:           std::ostringstream os;
alpar@209:           os << "Cannot find section: " << it->first;
deba@290:           throw FormatError(os.str());
alpar@209:         }
deba@189:       }
deba@189:     }
deba@189: 
deba@189:     /// @}
alpar@209: 
deba@189:   };
deba@189: 
deba@598:   /// \ingroup lemon_io
deba@598:   ///
deba@598:   /// \brief Return a \ref SectionReader class
deba@598:   ///
deba@598:   /// This function just returns a \ref SectionReader class.
deba@598:   ///
deba@598:   /// Please see SectionReader documentation about the custom section
deba@598:   /// input.
deba@598:   ///
deba@598:   /// \relates SectionReader
deba@598:   /// \sa sectionReader(const std::string& fn)
deba@598:   /// \sa sectionReader(const char *fn)
deba@598:   inline SectionReader sectionReader(std::istream& is) {
deba@598:     SectionReader tmp(is);
deba@598:     return tmp;
deba@598:   }
deba@598: 
kpeter@192:   /// \brief Return a \ref SectionReader class
alpar@209:   ///
kpeter@192:   /// This function just returns a \ref SectionReader class.
deba@189:   /// \relates SectionReader
deba@598:   /// \sa sectionReader(std::istream& is)
deba@598:   inline SectionReader sectionReader(const std::string& fn) {
deba@598:     SectionReader tmp(fn);
deba@189:     return tmp;
deba@189:   }
deba@189: 
kpeter@192:   /// \brief Return a \ref SectionReader class
alpar@209:   ///
kpeter@192:   /// This function just returns a \ref SectionReader class.
deba@189:   /// \relates SectionReader
deba@598:   /// \sa sectionReader(std::istream& is)
deba@189:   inline SectionReader sectionReader(const char* fn) {
deba@189:     SectionReader tmp(fn);
deba@189:     return tmp;
deba@189:   }
deba@189: 
deba@173:   /// \ingroup lemon_io
deba@173:   ///
alpar@209:   /// \brief Reader for the contents of the \ref lgf-format "LGF" file
deba@173:   ///
deba@173:   /// This class can be used to read the sections, the map names and
ladanyi@236:   /// the attributes from a file. Usually, the LEMON programs know
deba@173:   /// that, which type of graph, which maps and which attributes
deba@173:   /// should be read from a file, but in general tools (like glemon)
alpar@179:   /// the contents of an LGF file should be guessed somehow. This class
deba@173:   /// reads the graph and stores the appropriate information for
deba@173:   /// reading the graph.
deba@173:   ///
alpar@209:   ///\code
alpar@209:   /// LgfContents contents("graph.lgf");
alpar@179:   /// contents.run();
deba@173:   ///
kpeter@192:   /// // Does it contain any node section and arc section?
alpar@179:   /// if (contents.nodeSectionNum() == 0 || contents.arcSectionNum()) {
kpeter@192:   ///   std::cerr << "Failure, cannot find graph." << std::endl;
deba@173:   ///   return -1;
deba@173:   /// }
alpar@209:   /// std::cout << "The name of the default node section: "
alpar@179:   ///           << contents.nodeSection(0) << std::endl;
alpar@209:   /// std::cout << "The number of the arc maps: "
alpar@179:   ///           << contents.arcMaps(0).size() << std::endl;
alpar@209:   /// std::cout << "The name of second arc map: "
alpar@179:   ///           << contents.arcMaps(0)[1] << std::endl;
deba@173:   ///\endcode
alpar@209:   class LgfContents {
deba@173:   private:
deba@173: 
deba@173:     std::istream* _is;
deba@173:     bool local_is;
deba@173: 
deba@173:     std::vector<std::string> _node_sections;
deba@173:     std::vector<std::string> _edge_sections;
deba@173:     std::vector<std::string> _attribute_sections;
deba@173:     std::vector<std::string> _extra_sections;
deba@173: 
deba@173:     std::vector<bool> _arc_sections;
deba@173: 
deba@173:     std::vector<std::vector<std::string> > _node_maps;
deba@173:     std::vector<std::vector<std::string> > _edge_maps;
deba@173: 
deba@173:     std::vector<std::vector<std::string> > _attributes;
deba@173: 
deba@173: 
deba@173:     int line_num;
deba@173:     std::istringstream line;
alpar@209: 
deba@173:   public:
deba@173: 
deba@173:     /// \brief Constructor
deba@173:     ///
alpar@179:     /// Construct an \e LGF contents reader, which reads from the given
deba@173:     /// input stream.
alpar@209:     LgfContents(std::istream& is)
deba@173:       : _is(&is), local_is(false) {}
deba@173: 
deba@173:     /// \brief Constructor
deba@173:     ///
alpar@179:     /// Construct an \e LGF contents reader, which reads from the given
deba@173:     /// file.
alpar@209:     LgfContents(const std::string& fn)
deba@290:       : _is(new std::ifstream(fn.c_str())), local_is(true) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
deba@173: 
deba@173:     /// \brief Constructor
deba@173:     ///
alpar@179:     /// Construct an \e LGF contents reader, which reads from the given
deba@173:     /// file.
alpar@179:     LgfContents(const char* fn)
deba@290:       : _is(new std::ifstream(fn)), local_is(true) {
deba@295:       if (!(*_is)) {
deba@295:         delete _is;
deba@295:         throw IoError("Cannot open file", fn);
deba@295:       }
deba@290:     }
alpar@209: 
deba@173:     /// \brief Destructor
alpar@179:     ~LgfContents() {
deba@173:       if (local_is) delete _is;
deba@173:     }
deba@173: 
deba@190:   private:
alpar@209: 
deba@190:     LgfContents(const LgfContents&);
deba@190:     LgfContents& operator=(const LgfContents&);
deba@190: 
deba@190:   public:
deba@190: 
deba@173: 
kpeter@584:     /// \name Node Sections
deba@173:     /// @{
deba@173: 
deba@173:     /// \brief Gives back the number of node sections in the file.
deba@173:     ///
deba@173:     /// Gives back the number of node sections in the file.
deba@173:     int nodeSectionNum() const {
deba@173:       return _node_sections.size();
deba@173:     }
deba@173: 
alpar@209:     /// \brief Returns the node section name at the given position.
deba@173:     ///
alpar@209:     /// Returns the node section name at the given position.
deba@173:     const std::string& nodeSection(int i) const {
deba@173:       return _node_sections[i];
deba@173:     }
deba@173: 
deba@173:     /// \brief Gives back the node maps for the given section.
deba@173:     ///
deba@173:     /// Gives back the node maps for the given section.
alpar@182:     const std::vector<std::string>& nodeMapNames(int i) const {
deba@173:       return _node_maps[i];
deba@173:     }
deba@173: 
deba@173:     /// @}
deba@173: 
kpeter@584:     /// \name Arc/Edge Sections
deba@173:     /// @{
deba@173: 
alpar@181:     /// \brief Gives back the number of arc/edge sections in the file.
deba@173:     ///
alpar@181:     /// Gives back the number of arc/edge sections in the file.
alpar@181:     /// \note It is synonym of \c edgeSectionNum().
deba@173:     int arcSectionNum() const {
deba@173:       return _edge_sections.size();
deba@173:     }
deba@173: 
alpar@209:     /// \brief Returns the arc/edge section name at the given position.
deba@173:     ///
alpar@209:     /// Returns the arc/edge section name at the given position.
alpar@181:     /// \note It is synonym of \c edgeSection().
deba@173:     const std::string& arcSection(int i) const {
deba@173:       return _edge_sections[i];
deba@173:     }
deba@173: 
alpar@181:     /// \brief Gives back the arc/edge maps for the given section.
deba@173:     ///
alpar@181:     /// Gives back the arc/edge maps for the given section.
alpar@182:     /// \note It is synonym of \c edgeMapNames().
alpar@182:     const std::vector<std::string>& arcMapNames(int i) const {
deba@173:       return _edge_maps[i];
deba@173:     }
deba@173: 
deba@173:     /// @}
deba@173: 
alpar@181:     /// \name Synonyms
deba@173:     /// @{
deba@173: 
alpar@181:     /// \brief Gives back the number of arc/edge sections in the file.
deba@173:     ///
alpar@181:     /// Gives back the number of arc/edge sections in the file.
alpar@181:     /// \note It is synonym of \c arcSectionNum().
deba@173:     int edgeSectionNum() const {
deba@173:       return _edge_sections.size();
deba@173:     }
deba@173: 
alpar@209:     /// \brief Returns the section name at the given position.
deba@173:     ///
alpar@209:     /// Returns the section name at the given position.
alpar@181:     /// \note It is synonym of \c arcSection().
deba@173:     const std::string& edgeSection(int i) const {
deba@173:       return _edge_sections[i];
deba@173:     }
deba@173: 
deba@173:     /// \brief Gives back the edge maps for the given section.
deba@173:     ///
deba@173:     /// Gives back the edge maps for the given section.
alpar@182:     /// \note It is synonym of \c arcMapNames().
alpar@182:     const std::vector<std::string>& edgeMapNames(int i) const {
deba@173:       return _edge_maps[i];
deba@173:     }
deba@173: 
deba@173:     /// @}
deba@173: 
kpeter@584:     /// \name Attribute Sections
deba@173:     /// @{
deba@173: 
deba@173:     /// \brief Gives back the number of attribute sections in the file.
deba@173:     ///
deba@173:     /// Gives back the number of attribute sections in the file.
deba@173:     int attributeSectionNum() const {
deba@173:       return _attribute_sections.size();
deba@173:     }
deba@173: 
alpar@209:     /// \brief Returns the attribute section name at the given position.
deba@173:     ///
alpar@209:     /// Returns the attribute section name at the given position.
alpar@182:     const std::string& attributeSectionNames(int i) const {
deba@173:       return _attribute_sections[i];
deba@173:     }
deba@173: 
deba@173:     /// \brief Gives back the attributes for the given section.
deba@173:     ///
deba@173:     /// Gives back the attributes for the given section.
deba@173:     const std::vector<std::string>& attributes(int i) const {
deba@173:       return _attributes[i];
deba@173:     }
deba@173: 
deba@173:     /// @}
deba@173: 
kpeter@584:     /// \name Extra Sections
deba@173:     /// @{
deba@173: 
deba@173:     /// \brief Gives back the number of extra sections in the file.
deba@173:     ///
deba@173:     /// Gives back the number of extra sections in the file.
deba@173:     int extraSectionNum() const {
deba@173:       return _extra_sections.size();
deba@173:     }
deba@173: 
alpar@209:     /// \brief Returns the extra section type at the given position.
deba@173:     ///
alpar@209:     /// Returns the section type at the given position.
deba@173:     const std::string& extraSection(int i) const {
deba@173:       return _extra_sections[i];
deba@173:     }
deba@173: 
deba@173:     /// @}
deba@173: 
deba@173:   private:
deba@173: 
deba@173:     bool readLine() {
deba@173:       std::string str;
deba@173:       while(++line_num, std::getline(*_is, str)) {
alpar@209:         line.clear(); line.str(str);
alpar@209:         char c;
alpar@209:         if (line >> std::ws >> c && c != '#') {
alpar@209:           line.putback(c);
alpar@209:           return true;
alpar@209:         }
deba@173:       }
deba@173:       return false;
deba@173:     }
deba@173: 
deba@173:     bool readSuccess() {
deba@173:       return static_cast<bool>(*_is);
deba@173:     }
deba@173: 
deba@173:     void skipSection() {
deba@173:       char c;
deba@173:       while (readSuccess() && line >> c && c != '@') {
alpar@209:         readLine();
deba@173:       }
deba@427:       if (readSuccess()) {
deba@427:         line.putback(c);
deba@427:       }
deba@173:     }
deba@173: 
deba@173:     void readMaps(std::vector<std::string>& maps) {
deba@186:       char c;
deba@186:       if (!readLine() || !(line >> c) || c == '@') {
alpar@209:         if (readSuccess() && line) line.putback(c);
alpar@209:         return;
deba@186:       }
deba@186:       line.putback(c);
deba@173:       std::string map;
deba@173:       while (_reader_bits::readToken(line, map)) {
alpar@209:         maps.push_back(map);
deba@173:       }
deba@173:     }
deba@173: 
deba@173:     void readAttributes(std::vector<std::string>& attrs) {
deba@173:       readLine();
deba@173:       char c;
deba@173:       while (readSuccess() && line >> c && c != '@') {
alpar@209:         line.putback(c);
alpar@209:         std::string attr;
alpar@209:         _reader_bits::readToken(line, attr);
alpar@209:         attrs.push_back(attr);
alpar@209:         readLine();
deba@173:       }
deba@173:       line.putback(c);
deba@173:     }
deba@173: 
deba@173:   public:
deba@173: 
kpeter@584:     /// \name Execution of the Contents Reader
deba@173:     /// @{
deba@173: 
kpeter@192:     /// \brief Starts the reading
deba@173:     ///
kpeter@192:     /// This function starts the reading.
deba@173:     void run() {
deba@173: 
deba@173:       readLine();
deba@173:       skipSection();
deba@173: 
deba@173:       while (readSuccess()) {
deba@173: 
alpar@209:         char c;
alpar@209:         line >> c;
alpar@209: 
alpar@209:         std::string section, caption;
alpar@209:         _reader_bits::readToken(line, section);
alpar@209:         _reader_bits::readToken(line, caption);
alpar@209: 
alpar@209:         if (section == "nodes") {
alpar@209:           _node_sections.push_back(caption);
alpar@209:           _node_maps.push_back(std::vector<std::string>());
alpar@209:           readMaps(_node_maps.back());
alpar@209:           readLine(); skipSection();
alpar@209:         } else if (section == "arcs" || section == "edges") {
alpar@209:           _edge_sections.push_back(caption);
alpar@209:           _arc_sections.push_back(section == "arcs");
alpar@209:           _edge_maps.push_back(std::vector<std::string>());
alpar@209:           readMaps(_edge_maps.back());
alpar@209:           readLine(); skipSection();
alpar@209:         } else if (section == "attributes") {
alpar@209:           _attribute_sections.push_back(caption);
alpar@209:           _attributes.push_back(std::vector<std::string>());
alpar@209:           readAttributes(_attributes.back());
alpar@209:         } else {
alpar@209:           _extra_sections.push_back(section);
alpar@209:           readLine(); skipSection();
alpar@209:         }
deba@173:       }
deba@173:     }
deba@173: 
deba@173:     /// @}
alpar@209: 
deba@173:   };
deba@127: }
deba@127: 
deba@127: #endif