/* -*- C++ -*-

 * lemon/lemon_reader.h - Part of LEMON, a generic C++ optimization library

 *

 * Copyright (C) 2005 Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

 * (Egervary Research Group on Combinatorial Optimization, EGRES).

 *

 * Permission to use, modify and distribute this software is granted

 * provided that this copyright notice appears in all copies. For

 * precise terms see the accompanying LICENSE file.

 *

 * This software is provided "AS IS" with no warranty of any kind,

 * express or implied, and with no claim as to its suitability for any

 * purpose.

 *

 */

///\ingroup io_group

///\file

///\brief Lemon Format reader.

#ifndef LEMON_LEMON_READER_H

#define LEMON_LEMON_READER_H

#include <iostream>

#include <fstream>

#include <string>

#include <vector>

#include <algorithm>

#include <map>

#include <memory>

#include <lemon/error.h>

#include <lemon/graph_utils.h>

#include <lemon/utility.h>

#include <lemon/bits/item_reader.h>

#include <lemon/concept_check.h>

#include <lemon/concept/maps.h>

namespace lemon {

  namespace _reader_bits {

    template <typename Item>

    class ItemIdReader {

    public:

      bool isIdReader() { return true; }

      void readId(std::istream&, Item&) {}

      template <class _ItemIdReader>

      struct Constraints {

	void constraints() {

	  bool b = reader.isIdReader();

	  ignore_unused_variable_warning(b);

	  Item item;

	  reader.readId(is, item);

	}

	_ItemIdReader& reader;

	std::istream& is;

      };

    };

    template <typename Item>

    class ItemReader {

    public:

      void read(std::istream&, Item&) {}

      template <class _ItemReader>

      struct Constraints {

	void constraints() {

	  Item item;

	  reader.read(is, item);

	}

	_ItemReader& reader;

	std::istream& is;

      };

    };

    template <typename T>

    bool operator<(T, T) {

      throw DataFormatError("Id is not comparable");

    }

    template <typename T>

    struct Less {

      bool operator()(const T& p, const T& q) const {

	return p < q;

      }

    };

    template <typename M1, typename M2>

    class WriteComposeMap {

    public:

      typedef True NeedCopy;

      typedef typename M2::Key Key;

      typedef typename M1::Value Value;

      WriteComposeMap(typename SmartParameter<M1>::Type _m1, const M2& _m2) 

	: m1(_m1), m2(_m2) {}

      void set(const Key& key, const Value& value) {

	m1.set(m2[key], value);

      }

    private:

      typename SmartReference<M1>::Type m1;

      typename SmartConstReference<M2>::Type m2;

    };

    template <typename M1, typename M2>

    WriteComposeMap<M1, M2> writeComposeMap(M1& m1, const M2& m2) {

      return WriteComposeMap<M1, M2>(m1, m2);

    }

    template <typename M1, typename M2>

    WriteComposeMap<M1, M2> writeComposeMap(const M1& m1, const M2& m2) {

      return WriteComposeMap<M1, M2>(m1, m2);

    }

  }

  /// \ingroup io_group

  /// \brief Lemon Format reader class.

  /// 

  /// The Lemon Format contains several sections. We do not want to

  /// determine what sections are in a lemon file we give only a framework

  /// to read a section oriented format.

  ///

  /// In the Lemon Format each section starts with a line contains a \c \@

  /// character on the first not white space position. This line is the

  /// header line of the section. Each next lines belong to this section

  /// while it does not starts with \c \@ character. This line can start a 

  /// new section or if it can close the file with the \c \@end line.

  /// The file format ignore the empty and comment lines. The line is

  /// comment line if it starts with a \c # character. 

  ///

  /// The framework provides an abstract LemonReader::SectionReader class

  /// what defines the interface of a SectionReader. The SectionReader

  /// has the \c header() member function what get a header line string and

  /// decides if it want to process the next section. Several SectionReaders

  /// can be attached to an LemonReader and the first attached what can

  /// process the section will be used. Its \c read() member will called

  /// with a stream contains the section. From this stream the empty and

  /// comment lines are filtered out.

  ///

  /// \relates GraphReader

  /// \relates NodeSetReader

  /// \relates EdgeSetReader

  /// \relates NodesReader

  /// \relates EdgesReader

  /// \relates AttributeReader

  class LemonReader {

  private:

    class FilterStreamBuf : public std::streambuf {

    public:

      typedef std::streambuf Parent;

      typedef Parent::char_type char_type;

      FilterStreamBuf(std::istream& is, int& num) 

	: _is(is), _base(0), _eptr(0), 

	  _num(num), skip_state(after_endl) {}

    protected:

      enum skip_state_type {

	no_skip,

	after_endl,

	comment_line

      };

      char_type small_buf[1];

      std::istream& _is;

      char_type* _base;

      char_type* _eptr;

      int& _num;

      skip_state_type skip_state;

      char_type* base() { return _base; }

      char_type* eptr() { return _eptr; }

      int blen() { return _eptr - _base; }

      void setb(char_type* buf, int len) {

	_base = buf;

	_eptr = buf + len;

      }

      virtual std::streambuf* setbuf(char *buf, int len) {

	if (base()) return 0;

	if (buf != 0 && len >= (int)sizeof(small_buf)) {

	  setb(buf, len);

	} else {

	  setb(small_buf, sizeof(small_buf));

	}

	setg(0, 0, 0);

	return this;

      }

      bool put_char(char c) {

	switch (skip_state) {

	case no_skip:

	  switch (c) {

	  case '\n': 

	    skip_state = after_endl;

	    return true;

	  default:

	    return true;

	  }

	case after_endl:

	  switch (c) {

	  case '@':

	    return false;

	  case '\n': 

	    return false;

	  case '#':

	    skip_state = comment_line;

	    return false;

	  default:

	    if (!isspace(c)) {

	      skip_state = no_skip;

	      return true;

	    } else {

	      return false;

	    }

	  }

	  break;

	case comment_line:

	  switch (c) {

	  case '\n': 

	    skip_state = after_endl;

	    return false;

	  default:

	    return false;

	  }

	}

	return false;

      }

      virtual int underflow() {

	char c;

	if (_is.read(&c, 1)) {

	  _is.putback(c);

	  if (c == '@') {

	    return EOF;

	  }

	} else {

	  return EOF;

	}

	char_type *ptr;

	for (ptr = base(); ptr != eptr(); ++ptr) {

	  if (_is.read(&c, 1)) {

	    if (c == '\n') ++_num;

	    if (put_char(c)) {

	      *ptr = c;

	    } else {

	      if (skip_state == after_endl && c == '@') {

		_is.putback('@');

		break;

	      }

	      --ptr;

	    }

	  } else {

	    break;

	  }

	}

	setg(base(), base(), ptr);

	return *base();

      }

      virtual int sync() {

	return EOF;

      }

    };

  public:

    /// \brief Abstract base class for reading a section.

    ///

    /// This class has an \c header() member function what get a 

    /// header line string and decides if it want to process the next 

    /// section. Several SectionReaders can be attached to an LemonReader 

    /// and the first attached what can process the section will be used. 

    /// Its \c read() member will called with a stream contains the section. 

    /// From this stream the empty lines and comments are filtered out.

    class SectionReader {

      friend class LemonReader;

    protected:

      /// \brief Constructor for SectionReader.

      ///

      /// Constructor for SectionReader. It attach this reader to

      /// the given LemonReader.

      SectionReader(LemonReader& reader) {

	reader.attach(*this);

      }

      /// \brief Gives back true when the SectionReader can process 

      /// the section with the given header line.

      ///

      /// It gives back true when the SectionReader can process

      /// the section with the given header line.

      virtual bool header(const std::string& line) = 0;

      /// \brief Reader function of the section.

      ///

      /// It reads the content of the section.

      virtual void read(std::istream& is) = 0;

    };

    /// \brief Constructor for LemonReader.

    ///

    /// Constructor for LemonReader which reads from the given stream.

    LemonReader(std::istream& _is) 

      : is(&_is), own_is(false) {}

    /// \brief Constructor for LemonReader.

    ///

    /// Constructor for LemonReader which reads from the given file.

    LemonReader(const std::string& filename) 

      : is(0), own_is(true) {

      is = new std::ifstream(filename.c_str());

    }

    /// \brief Desctructor for LemonReader.

    ///

    /// Desctructor for LemonReader.

    ~LemonReader() {

      if (own_is) {

	delete is;

      }

    }

  private:

    LemonReader(const LemonReader&);

    void operator=(const LemonReader&);

    void attach(SectionReader& reader) {

      readers.push_back(&reader);

    }

  public:

    /// \brief Executes the LemonReader.

    /// 

    /// It executes the LemonReader.

    void run() {

      int line_num = 0;

      std::string line;

      try {

	while ((++line_num, getline(*is, line)) && line.find("@end") != 0) {

	  SectionReaders::iterator it;

	  for (it = readers.begin(); it != readers.end(); ++it) {

	    if ((*it)->header(line)) {

	      char buf[2048];

	      FilterStreamBuf buffer(*is, line_num);

	      buffer.pubsetbuf(buf, sizeof(buf));

	      std::istream is(&buffer);

	      (*it)->read(is);

	      break;

	    }

	  }

	}

      } catch (DataFormatError& error) {

	error.line(line_num);

	throw error;

      }	

    }

  private:

    std::istream* is;

    bool own_is;

    typedef std::vector<SectionReader*> SectionReaders;

    SectionReaders readers;

  };

  /// \brief Helper class for implementing the common SectionReaders.

  ///

  /// Helper class for implementing the common SectionReaders.

  class CommonSectionReaderBase : public LemonReader::SectionReader {

    typedef LemonReader::SectionReader Parent;

  protected:

    /// \brief Constructor for CommonSectionReaderBase.

    ///

    /// Constructor for CommonSectionReaderBase. It attach this reader to

    /// the given LemonReader.

    CommonSectionReaderBase(LemonReader& _reader) 

      : Parent(_reader) {}

    template <typename _Item>

    class ReaderBase;

    template <typename _Item>

    class InverterBase : public ReaderBase<_Item> {

    public:

      typedef _Item Item;

      virtual void read(std::istream&, const Item&) = 0;

      virtual Item read(std::istream&) const = 0;

      virtual InverterBase<_Item>* getInverter() {

	return this;

      }

    };

    template <typename _Item, typename _Map, typename _Reader>

    class MapReaderInverter : public InverterBase<_Item> {

    public:

      typedef _Item Item;

      typedef _Reader Reader;

      typedef typename Reader::Value Value;

      typedef _Map Map;

      typedef std::map<Value, Item, _reader_bits::Less<Value> > Inverse;

      typename SmartReference<Map>::Type map;

      Reader reader;

      Inverse inverse;

      MapReaderInverter(typename SmartParameter<Map>::Type _map,

			const Reader& _reader) 

	: map(_map), reader(_reader) {}

      virtual ~MapReaderInverter() {}

      virtual void read(std::istream& is, const Item& item) {

	Value value;

	reader.read(is, value);

	map.set(item, value);

	typename Inverse::iterator it = inverse.find(value);

	if (it == inverse.end()) {

	  inverse.insert(std::make_pair(value, item));

	} else {

	  throw DataFormatError("Multiple ID occurence");

	}

      }

      virtual Item read(std::istream& is) const {

	Value value;

	reader.read(is, value);	

	typename Inverse::const_iterator it = inverse.find(value);

	if (it != inverse.end()) {

	  return it->second;

	} else {

	  throw DataFormatError("Invalid ID error");

	}

      }      

    };

    template <typename _Item, typename _Reader>

    class SkipReaderInverter : public InverterBase<_Item> {

    public:

      typedef _Item Item;

      typedef _Reader Reader;

      typedef typename Reader::Value Value;

      typedef std::map<Value, Item, _reader_bits::Less<Value> > Inverse;

      Reader reader;

      SkipReaderInverter(const Reader& _reader) 

	: reader(_reader) {}

      virtual ~SkipReaderInverter() {}

      virtual void read(std::istream& is, const Item& item) {

	Value value;

	reader.read(is, value);

	typename Inverse::iterator it = inverse.find(value);

	if (it == inverse.end()) {

	  inverse.insert(std::make_pair(value, item));

	} else {

	  throw DataFormatError("Multiple ID occurence error");

	}

      }

      virtual Item read(std::istream& is) const {

	Value value;

	reader.read(is, value);	

	typename Inverse::const_iterator it = inverse.find(value);

	if (it != inverse.end()) {

	  return it->second;

	} else {

	  throw DataFormatError("Invalid ID error");

	}

      }

    private:

      Inverse inverse;

    };

    template <typename _Item>    

    class ReaderBase {

    public:

      typedef _Item Item;

      virtual ~ReaderBase() {}

      virtual void read(std::istream& is, const Item& item) = 0;

      virtual InverterBase<_Item>* getInverter() = 0;

    };

    template <typename _Item, typename _Map, typename _Reader>

    class MapReader : public ReaderBase<_Item> {

    public:

      typedef _Map Map;

      typedef _Reader Reader;

      typedef typename Reader::Value Value;

      typedef _Item Item;

      typename SmartReference<Map>::Type map;

      Reader reader;

      MapReader(typename SmartParameter<Map>::Type _map, 

		const Reader& _reader) 

	: map(_map), reader(_reader) {}

      virtual ~MapReader() {}

      virtual void read(std::istream& is, const Item& item) {

	Value value;

	reader.read(is, value);

	map.set(item, value);

      }

      virtual InverterBase<_Item>* getInverter() {

	return new MapReaderInverter<Item, Map, Reader>(map, reader);

      }

    };

    template <typename _Item, typename _Reader>

    class SkipReader : public ReaderBase<_Item> {

    public:

      typedef _Reader Reader;

      typedef typename Reader::Value Value;

      typedef _Item Item;

      Reader reader;

      SkipReader(const Reader& _reader) : reader(_reader) {}

      virtual ~SkipReader() {}

      virtual void read(std::istream& is, const Item&) {

	Value value;

	reader.read(is, value);

      }      

      virtual InverterBase<Item>* getInverter() {

	return new SkipReaderInverter<Item, Reader>(reader);

      }

    };

    template <typename _Item>

    class IdReaderBase {

    public:

      typedef _Item Item;

      virtual Item read(std::istream& is) const = 0;

      virtual bool isIdReader() const = 0;

    };

    template <typename _Item, typename _BoxedIdReader>

    class IdReader : public IdReaderBase<_Item> {

    public:

      typedef _Item Item;

      typedef _BoxedIdReader BoxedIdReader;

      const BoxedIdReader& boxedIdReader;

      IdReader(const BoxedIdReader& _boxedIdReader) 

	: boxedIdReader(_boxedIdReader) {}

      virtual Item read(std::istream& is) const {

	Item item;

	boxedIdReader.readId(is, item);

	return item;

      }

      virtual bool isIdReader() const {

	return boxedIdReader.isIdReader();

      }

    };

    class ValueReaderBase {

    public:

      virtual void read(std::istream&) {};

    };

    template <typename _Value, typename _Reader>

    class ValueReader : public ValueReaderBase {

    public:

      typedef _Value Value;

      typedef _Reader Reader;

      ValueReader(Value& _value, const Reader& _reader)

 	: value(_value), reader(_reader) {}

      virtual void read(std::istream& is) {

	reader.read(is, value);

      }

    private:

      Value& value;

      Reader reader;

    };

  };

  /// \ingroup io_group

  /// \brief SectionReader for reading a graph's nodeset.

  ///

  /// The lemon format can store multiple graph nodesets with several maps.

  /// The nodeset section's header line is \c \@nodeset \c nodeset_id, but the

  /// \c nodeset_id may be empty.

  ///

  /// The first line of the section contains the names of the maps separated

  /// with white spaces. Each next lines describes a node in the nodeset, and

  /// contains the mapped values for each map.

  ///

  /// If the nodeset contains an \c "id" named map then it will be regarded

  /// as id map. This map should contain only unique values and when the 

  /// \c readId() member will read a value from the given stream it will

  /// give back that node which is mapped to this value.

  ///

  /// \relates LemonReader

  template <typename _Graph, typename _Traits = DefaultReaderTraits>

  class NodeSetReader : public CommonSectionReaderBase {

    typedef CommonSectionReaderBase Parent;

  public:

    typedef _Graph Graph;

    typedef _Traits Traits;

    typedef typename Graph::Node Node;

    typedef typename Traits::Skipper DefaultSkipper;

    /// \brief Constructor.

    ///

    /// Constructor for NodeSetReader. It creates the NodeSetReader and

    /// attach it into the given LemonReader. The nodeset reader will

    /// add the readed nodes to the given Graph. The reader will read

    /// the section when the \c section_id and the \c _id are the same. 

    NodeSetReader(LemonReader& _reader, 

		  typename SmartParameter<Graph>::Type _graph, 

		  const std::string& _id = std::string(),

		  const DefaultSkipper& _skipper = DefaultSkipper()) 

      : Parent(_reader), graph(_graph), id(_id), skipper(_skipper) {} 

    /// \brief Destructor.

    ///

    /// Destructor for NodeSetReader.

    virtual ~NodeSetReader() {

      for (typename MapReaders::iterator it = readers.begin(); 

	   it != readers.end(); ++it) {

	delete it->second;

      }

    }

  private:

    NodeSetReader(const NodeSetReader&);

    void operator=(const NodeSetReader&);

  public:

    /// \brief Add a new node map reader command for the reader.

    ///

    /// Add a new node map reader command for the reader.

    template <typename Map>

    NodeSetReader& readNodeMap(std::string name, Map& map) {

      return _readMap<

	typename Traits::template Reader<typename Map::Value>, Map,

	typename SmartParameter<Map>::Type>(name, map);

    }

    template <typename Map>

    NodeSetReader& readNodeMap(std::string name, const Map& map) {

      return _readMap<

	typename Traits::template Reader<typename Map::Value>, Map,

	typename SmartParameter<Map>::Type>(name, map);

    }

    /// \brief Add a new node map reader command for the reader.

    ///

    /// Add a new node map reader command for the reader.

    template <typename Reader, typename Map>

    NodeSetReader& readNodeMap(std::string name, Map& map, 

			       const Reader& reader = Reader()) {

      return _readMap<Reader, Map, typename SmartParameter<Map>::Type>

	(name, map, reader);

    }

    template <typename Reader, typename Map>

    NodeSetReader& readNodeMap(std::string name, const Map& map, 

			       const Reader& reader = Reader()) {

      return _readMap<Reader, Map, typename SmartParameter<Map>::Type>

	(name, map, reader);

    }

  private:

    template <typename Reader, typename Map, typename MapParameter>

    NodeSetReader& _readMap(std::string name, MapParameter map, 

			    const Reader& reader = Reader()) {

      checkConcept<concept::WriteMap<Node, typename Map::Value>, Map>();

      checkConcept<_reader_bits::ItemReader<typename Map::Value>, Reader>();

      if (readers.find(name) != readers.end()) {

	ErrorMessage msg;

	msg << "Multiple read rule for node map: " << name;

	throw IOParameterError(msg.message());

      }

      readers.insert(

	make_pair(name, new MapReader<Node, Map, Reader>(map, reader)));

      return *this;

    }

  public:

    /// \brief Add a new node map skipper command for the reader.

    ///

    /// Add a new node map skipper command for the reader.

    template <typename Reader>

    NodeSetReader& skipNodeMap(std::string name, 

			   const Reader& reader = Reader()) {

      if (readers.find(name) != readers.end()) {

	ErrorMessage msg;

	msg << "Multiple read rule for node map: " << name;

	throw IOParameterError(msg.message());

      }

      readers.insert(make_pair(name, new SkipReader<Node, Reader>(reader)));

      return *this;

    }

  protected:

    /// \brief Gives back true when the SectionReader can process 

    /// the section with the given header line.

    ///

    /// It gives back true when the header line starts with \c \@nodeset,

    /// and the header line's id and the nodeset's id are the same.

    virtual bool header(const std::string& line) {

      std::istringstream ls(line);

      std::string command;

      std::string name;

      ls >> command >> name;

      return command == "@nodeset" && name == id;

    }

    /// \brief Reader function of the section.

    ///

    /// It reads the content of the section.

    virtual void read(std::istream& is) {

      std::vector<ReaderBase<Node>* > index;

      std::string line;

      getline(is, line);

      std::istringstream ls(line);	

      while (ls >> id) {

	typename MapReaders::iterator it = readers.find(id);

	if (it != readers.end()) {

	  index.push_back(it->second);

	} else {

	  index.push_back(&skipper);

	}

	if (id == "id" && inverter.get() == 0) {

	  inverter.reset(index.back()->getInverter());

	  index.back() = inverter.get();

	}

      }

      while (getline(is, line)) {	

	Node node = graph.addNode();

	std::istringstream ls(line);

	for (int i = 0; i < (int)index.size(); ++i) {

	  index[i]->read(ls, node);

	}

      }

    }

  public:

    /// \brief Returns true if the nodeset can give back the node by its id.

    ///

    /// Returns true if the nodeset can give back the node by its id.

    /// It is possible only if an "id" named map was read.

    bool isIdReader() const {

      return inverter.get() != 0;

    }

    /// \brief Gives back the node by its id.

    ///

    /// It reads an id from the stream and gives back which node belongs to

    /// it. It is possible only if there was read an "id" named map.

    void readId(std::istream& is, Node& node) const {

      node = inverter->read(is);

    } 

  private:

    typedef std::map<std::string, ReaderBase<Node>*> MapReaders;

    MapReaders readers;

    typename SmartReference<Graph>::Type graph;   

    std::string id;

    SkipReader<Node, DefaultSkipper> skipper;

    std::auto_ptr<InverterBase<Node> > inverter;

  };

  /// \ingroup io_group

  /// \brief SectionReader for reading a graph's edgeset.

  ///

  /// The lemon format can store multiple graph edgesets with several maps.

  /// The edgeset section's header line is \c \@edgeset \c edgeset_id, but the

  /// \c edgeset_id may be empty.

  ///

  /// The first line of the section contains the names of the maps separated

  /// with white spaces. Each next lines describes an edge in the edgeset. The

  /// line contains the source and the target nodes' id and the mapped 

  /// values for each map.

  ///

  /// If the edgeset contains an \c "id" named map then it will be regarded

  /// as id map. This map should contain only unique values and when the 

  /// \c readId() member will read a value from the given stream it will

  /// give back that edge which is mapped to this value.

  ///

  /// The edgeset reader needs a node id reader to identify which nodes

  /// have to be connected. If a NodeSetReader reads an "id" named map,

  /// it will be able to resolve the nodes by ids.

  ///

  /// \relates LemonReader

  template <typename _Graph, typename _Traits = DefaultReaderTraits>

  class EdgeSetReader : public CommonSectionReaderBase {

    typedef CommonSectionReaderBase Parent;

  public:

    typedef _Graph Graph;

    typedef _Traits Traits;

    typedef typename Graph::Node Node;

    typedef typename Graph::Edge Edge;

    typedef typename Traits::Skipper DefaultSkipper;

    /// \brief Constructor.

    ///

    /// Constructor for EdgeSetReader. It creates the EdgeSetReader and

    /// attach it into the given LemonReader. The edgeset reader will

    /// add the readed edges to the given Graph. It will use the given

    /// node id reader to read the source and target nodes of the edges.

    /// The reader will read the section only if the \c _id and the 

    /// \c edgset_id are the same. 

    template <typename NodeIdReader>

    EdgeSetReader(LemonReader& _reader, 

		  typename SmartParameter<Graph>::Type _graph, 

		  const NodeIdReader& _nodeIdReader, 

		  const std::string& _id = std::string(),

		  const DefaultSkipper& _skipper = DefaultSkipper()) 

      : Parent(_reader), graph(_graph), id(_id), skipper(_skipper) {

      checkConcept<_reader_bits::ItemIdReader<Node>, NodeIdReader>();

      nodeIdReader.reset(new IdReader<Node, NodeIdReader>(_nodeIdReader));

    }

    /// \brief Destructor.

    ///

    /// Destructor for EdgeSetReader.

    virtual ~EdgeSetReader() {

      for (typename MapReaders::iterator it = readers.begin(); 

	   it != readers.end(); ++it) {

	delete it->second;

      }

    }

  private:

    EdgeSetReader(const EdgeSetReader&);

    void operator=(const EdgeSetReader&);

  public:

    /// \brief Add a new edge map reader command for the reader.

    ///

    /// Add a new edge map reader command for the reader.

    template <typename Map>

    EdgeSetReader& readEdgeMap(std::string name, Map& map) {

      return _readMap<

	typename Traits::template Reader<typename Map::Value>, Map,

	typename SmartParameter<Map>::Type>(name, map);

    }

    template <typename Map>

    EdgeSetReader& readEdgeMap(std::string name, const Map& map) {

      return _readMap<

	typename Traits::template Reader<typename Map::Value>, Map,

	typename SmartParameter<Map>::Type>(name, map);

    }

    /// \brief Add a new edge map reader command for the reader.

    ///

    /// Add a new edge map reader command for the reader.

    template <typename Reader, typename Map>

    EdgeSetReader& readEdgeMap(std::string name, Map& map, 

			   const Reader& reader = Reader()) {

      return _readMap<Reader, Map,

	typename SmartParameter<Map>::Type>(name, map, reader);

    }

    template <typename Reader, typename Map>

    EdgeSetReader& readEdgeMap(std::string name, const Map& map, 

			       const Reader& reader = Reader()) {

      return _readMap<Reader, Map,

	typename SmartParameter<Map>::Type>(name, map, reader);

    }

  private:

    template <typename Reader, typename Map, typename MapParameter>

    EdgeSetReader& _readMap(std::string name, MapParameter map, 

			    const Reader& reader = Reader()) {

      checkConcept<concept::WriteMap<Edge, typename Map::Value>, Map>();

      checkConcept<_reader_bits::ItemReader<typename Map::Value>, Reader>();

      if (readers.find(name) != readers.end()) {

	ErrorMessage msg;

	msg << "Multiple read rule for edge map: " << name;

	throw IOParameterError(msg.message());

      }

      readers.insert(

	make_pair(name, new MapReader<Edge, Map, Reader>(map, reader)));

      return *this;

    }

  public:

    /// \brief Add a new edge map skipper command for the reader.

    ///

    /// Add a new edge map skipper command for the reader.

    template <typename Reader>

    EdgeSetReader& skipEdgeMap(std::string name, 

			       const Reader& reader = Reader()) {

      if (readers.find(name) != readers.end()) {

	ErrorMessage msg;

	msg << "Multiple read rule for edge map: " << name;

	throw IOParameterError(msg.message());

      }

      readers.insert(make_pair(name, new SkipReader<Edge, Reader>(reader)));

      return *this;

    }

  protected:

    /// \brief Gives back true when the SectionReader can process 

    /// the section with the given header line.

    ///

    /// It gives back true when the header line starts with \c \@edgeset,

    /// and the header line's id and the edgeset's id are the same.

    virtual bool header(const std::string& line) {

      std::istringstream ls(line);

      std::string command;

      std::string name;

      ls >> command >> name;

      return command == "@edgeset" && name == id;

    }

    /// \brief Reader function of the section.

    ///

    /// It reads the content of the section.

    virtual void read(std::istream& is) {

      if (!nodeIdReader->isIdReader()) {

	throw DataFormatError("Cannot find nodeset or ID map");

      }

      std::vector<ReaderBase<Edge>* > index;

      std::string line;

      getline(is, line);

      std::istringstream ls(line);	

      while (ls >> id) {

	typename MapReaders::iterator it = readers.find(id);

	if (it != readers.end()) {

	  index.push_back(it->second);

	} else {

	  index.push_back(&skipper);

	}

	if (id == "id" && inverter.get() == 0) {

	  inverter.reset(index.back()->getInverter());