| 1 | /*
 | 
|---|
| 2 |  * AdjacencyList.hpp
 | 
|---|
| 3 |  *
 | 
|---|
| 4 |  *  Created on: Mar 3, 2011
 | 
|---|
| 5 |  *      Author: heber
 | 
|---|
| 6 |  */
 | 
|---|
| 7 | 
 | 
|---|
| 8 | #ifndef ADJACENCYLIST_HPP_
 | 
|---|
| 9 | #define ADJACENCYLIST_HPP_
 | 
|---|
| 10 | 
 | 
|---|
| 11 | // include config.h
 | 
|---|
| 12 | #ifdef HAVE_CONFIG_H
 | 
|---|
| 13 | #include <config.h>
 | 
|---|
| 14 | #endif
 | 
|---|
| 15 | 
 | 
|---|
| 16 | #include <iosfwd>
 | 
|---|
| 17 | #include <map>
 | 
|---|
| 18 | #include <set>
 | 
|---|
| 19 | 
 | 
|---|
| 20 | #include "types.hpp"
 | 
|---|
| 21 | #include "World.hpp"
 | 
|---|
| 22 | 
 | 
|---|
| 23 | class atom;
 | 
|---|
| 24 | 
 | 
|---|
| 25 | /** This class contains the adjacency structure inside an internal map of atoms.
 | 
|---|
| 26 |  *
 | 
|---|
| 27 |  * The adjacency structure is either from a file or from a given set of atoms.
 | 
|---|
| 28 |  *
 | 
|---|
| 29 |  * We may compare a subset of atoms against this internal bond structure. It is
 | 
|---|
| 30 |  * true it is a true subset of the bond structure.
 | 
|---|
| 31 |  */
 | 
|---|
| 32 | class AdjacencyList
 | 
|---|
| 33 | {
 | 
|---|
| 34 |   //!> Unit test is granted access to internal data
 | 
|---|
| 35 |   friend class AdjacencyListTest;
 | 
|---|
| 36 | public:
 | 
|---|
| 37 |   typedef std::vector<atomId_t> atomids_t;
 | 
|---|
| 38 | 
 | 
|---|
| 39 |   /** Default constructor for class AdjacencyList.
 | 
|---|
| 40 |    *
 | 
|---|
| 41 |    * We simply have an empty adjacency list here.
 | 
|---|
| 42 |    */
 | 
|---|
| 43 |   AdjacencyList() {}
 | 
|---|
| 44 |   AdjacencyList(std::istream &File);
 | 
|---|
| 45 |   AdjacencyList(const atomids_t &atoms);
 | 
|---|
| 46 |   ~AdjacencyList();
 | 
|---|
| 47 | 
 | 
|---|
| 48 |   bool operator<(const AdjacencyList &other) const;
 | 
|---|
| 49 | 
 | 
|---|
| 50 |   /** Comparison operator whether this adjacency list is a subset of \a other.
 | 
|---|
| 51 |    *
 | 
|---|
| 52 |    * @return true - is subset, false - is not subset
 | 
|---|
| 53 |    */
 | 
|---|
| 54 |   bool operator>(const AdjacencyList &other) const {
 | 
|---|
| 55 |     return other < *this;
 | 
|---|
| 56 |   }
 | 
|---|
| 57 |   /** Equality operator, determines whether both adjacencies are the same.
 | 
|---|
| 58 |    *
 | 
|---|
| 59 |    * @return true - both are the same, false - at least one is not a subset of the other
 | 
|---|
| 60 |    */
 | 
|---|
| 61 |   bool operator==(const AdjacencyList &other) const {
 | 
|---|
| 62 |     return (other < *this) && (*this < other);
 | 
|---|
| 63 |   }
 | 
|---|
| 64 |   /** Inquality operator, determines whether both adjacencies are not equal.
 | 
|---|
| 65 |    *
 | 
|---|
| 66 |    * @return true - both are not equal, false - both are the subset of one another
 | 
|---|
| 67 |    */
 | 
|---|
| 68 |   bool operator!=(const AdjacencyList &other) const {
 | 
|---|
| 69 |     return !(*this == other);
 | 
|---|
| 70 |   }
 | 
|---|
| 71 | 
 | 
|---|
| 72 |   /** Stores the adjacency contained in this instance to file.
 | 
|---|
| 73 |    *
 | 
|---|
| 74 |    * @param File stream to write to
 | 
|---|
| 75 |    * @return true - File is good, false - else
 | 
|---|
| 76 |    */
 | 
|---|
| 77 |   bool StoreToFile(std::ostream &File) const;
 | 
|---|
| 78 | 
 | 
|---|
| 79 | private:
 | 
|---|
| 80 |   typedef std::set<atomId_t> KeysSet;
 | 
|---|
| 81 |   typedef std::set<atomId_t> ValuesSet;
 | 
|---|
| 82 |   typedef std::pair<atomId_t, atomId_t> AtomBondPair;
 | 
|---|
| 83 |   typedef std::multimap< atomId_t, atomId_t > AtomBondMap;
 | 
|---|
| 84 |   typedef std::pair<AtomBondMap::const_iterator, AtomBondMap::const_iterator> AtomBondRange;
 | 
|---|
| 85 |   AtomBondMap atombondmap;
 | 
|---|
| 86 | 
 | 
|---|
| 87 |   KeysSet getKeys(const AtomBondRange &_range) const;
 | 
|---|
| 88 |   ValuesSet getValues(const AtomBondRange&_range) const;
 | 
|---|
| 89 | 
 | 
|---|
| 90 |   void CreateMap(atomids_t atoms);
 | 
|---|
| 91 |   bool ParseFromFile(std::istream &File);
 | 
|---|
| 92 | };
 | 
|---|
| 93 | 
 | 
|---|
| 94 | #endif /* ADJACENCYLIST_HPP_ */
 | 
|---|