| 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_ */ | 
|---|