| 1 | /*
 | 
|---|
| 2 |  * Cluster.hpp
 | 
|---|
| 3 |  *
 | 
|---|
| 4 |  *  Created on: Jan 16, 2012
 | 
|---|
| 5 |  *      Author: heber
 | 
|---|
| 6 |  */
 | 
|---|
| 7 | 
 | 
|---|
| 8 | #ifndef CLUSTER_HPP_
 | 
|---|
| 9 | #define CLUSTER_HPP_
 | 
|---|
| 10 | 
 | 
|---|
| 11 | 
 | 
|---|
| 12 | // include config.h
 | 
|---|
| 13 | #ifdef HAVE_CONFIG_H
 | 
|---|
| 14 | #include <config.h>
 | 
|---|
| 15 | #endif
 | 
|---|
| 16 | 
 | 
|---|
| 17 | #include "types.hpp"
 | 
|---|
| 18 | 
 | 
|---|
| 19 | #include <set>
 | 
|---|
| 20 | 
 | 
|---|
| 21 | #include "ClusterInterface.hpp"
 | 
|---|
| 22 | 
 | 
|---|
| 23 | #include "Atom/CopyAtoms/CopyAtomsInterface.hpp"
 | 
|---|
| 24 | #include "LinearAlgebra/Vector.hpp"
 | 
|---|
| 25 | #include "Shapes/Shape.hpp"
 | 
|---|
| 26 | 
 | 
|---|
| 27 | class atom;
 | 
|---|
| 28 | class RealSpaceMatrix;
 | 
|---|
| 29 | 
 | 
|---|
| 30 | /** A Cluster is a set of atoms within a given Shape that contains them.
 | 
|---|
| 31 |  *
 | 
|---|
| 32 |  */
 | 
|---|
| 33 | class Cluster : public ClusterInterface
 | 
|---|
| 34 | {
 | 
|---|
| 35 |   //!> grant unit test access
 | 
|---|
| 36 |   friend class ClusterTest;
 | 
|---|
| 37 | public:
 | 
|---|
| 38 |   typedef std::set<atomId_t> atomIdSet;
 | 
|---|
| 39 |   typedef std::vector<atom *> AtomVector;
 | 
|---|
| 40 | 
 | 
|---|
| 41 |   Cluster(const Cluster & _cluster);
 | 
|---|
| 42 |   Cluster(const Shape &s);
 | 
|---|
| 43 |   Cluster(const atomIdSet &atoms, const Shape &s);
 | 
|---|
| 44 |   ~Cluster();
 | 
|---|
| 45 | 
 | 
|---|
| 46 |   void insert(const atomId_t id);
 | 
|---|
| 47 |   void erase(const atomId_t id);
 | 
|---|
| 48 | 
 | 
|---|
| 49 |   /** Getter for the internal set of atomic ids.
 | 
|---|
| 50 |    *
 | 
|---|
| 51 |    * @return atomIdSet with \a atoms
 | 
|---|
| 52 |    */
 | 
|---|
| 53 |   const atomIdSet & getAtomIds() const { return atoms; }
 | 
|---|
| 54 | 
 | 
|---|
| 55 |   Cluster_impl clone(CopyAtomsInterface& copyMethod, const Vector &offset = zeroVec) const;
 | 
|---|
| 56 | 
 | 
|---|
| 57 |   void translate(const Vector &offset);
 | 
|---|
| 58 | 
 | 
|---|
| 59 |   void transform(const RealSpaceMatrix &M);
 | 
|---|
| 60 | 
 | 
|---|
| 61 |   bool isInside(const atomId_t id) const {
 | 
|---|
| 62 |     return IsInShape(id);
 | 
|---|
| 63 |   }
 | 
|---|
| 64 | 
 | 
|---|
| 65 |   /** Getter for the internal Shape.
 | 
|---|
| 66 |    *
 | 
|---|
| 67 |    * @return Shape \a s
 | 
|---|
| 68 |    */
 | 
|---|
| 69 |   const Shape& getShape() const { return s; }
 | 
|---|
| 70 | 
 | 
|---|
| 71 | private:
 | 
|---|
| 72 |   /** Getter for the contained atoms.
 | 
|---|
| 73 |    *
 | 
|---|
| 74 |    * @return vector of atomic refs
 | 
|---|
| 75 |    */
 | 
|---|
| 76 |   AtomVector getAtomRefs() const;
 | 
|---|
| 77 | 
 | 
|---|
| 78 |   //!> do not allow Cluster without Shape
 | 
|---|
| 79 |   Cluster();
 | 
|---|
| 80 | 
 | 
|---|
| 81 |   bool IsInShape(const atomId_t id) const;
 | 
|---|
| 82 |   atom * const getAtomById(const atomId_t id) const;
 | 
|---|
| 83 | 
 | 
|---|
| 84 |   atomIdSet atoms;
 | 
|---|
| 85 |   Shape s;
 | 
|---|
| 86 | };
 | 
|---|
| 87 | 
 | 
|---|
| 88 | 
 | 
|---|
| 89 | #endif /* CLUSTER_HPP_ */
 | 
|---|