| 1 | /* | 
|---|
| 2 | * Project: MoleCuilder | 
|---|
| 3 | * Description: creates and alters molecular systems | 
|---|
| 4 | * Copyright (C)  2010 University of Bonn. All rights reserved. | 
|---|
| 5 | * Please see the LICENSE file or "Copyright notice" in builder.cpp for details. | 
|---|
| 6 | */ | 
|---|
| 7 |  | 
|---|
| 8 | /* | 
|---|
| 9 | * molecule_pointcloud.cpp | 
|---|
| 10 | * | 
|---|
| 11 | *  Created on: Oct 5, 2009 | 
|---|
| 12 | *      Author: heber | 
|---|
| 13 | */ | 
|---|
| 14 |  | 
|---|
| 15 | // include config.h | 
|---|
| 16 | #ifdef HAVE_CONFIG_H | 
|---|
| 17 | #include <config.h> | 
|---|
| 18 | #endif | 
|---|
| 19 |  | 
|---|
| 20 | #include "Helpers/MemDebug.hpp" | 
|---|
| 21 |  | 
|---|
| 22 | #include "atom.hpp" | 
|---|
| 23 | #include "config.hpp" | 
|---|
| 24 | #include "Helpers/Info.hpp" | 
|---|
| 25 | #include "molecule.hpp" | 
|---|
| 26 |  | 
|---|
| 27 | /************************************* Functions for class molecule *********************************/ | 
|---|
| 28 |  | 
|---|
| 29 | TesselPoint * molecule::getValue(const_iterator &rhs) const | 
|---|
| 30 | { | 
|---|
| 31 | return *rhs; | 
|---|
| 32 | } | 
|---|
| 33 |  | 
|---|
| 34 | TesselPoint * molecule::getValue(iterator &rhs) const | 
|---|
| 35 | { | 
|---|
| 36 | return *rhs; | 
|---|
| 37 | } | 
|---|
| 38 |  | 
|---|
| 39 | /** Returns a name for this point cloud, here the molecule's name. | 
|---|
| 40 | * \return name of point cloud | 
|---|
| 41 | */ | 
|---|
| 42 | const char * const molecule::GetName() const | 
|---|
| 43 | { | 
|---|
| 44 | return name; | 
|---|
| 45 | }; | 
|---|
| 46 |  | 
|---|
| 47 | /** Determine center of all atoms. | 
|---|
| 48 | * \param *out output stream for debugging | 
|---|
| 49 | * \return pointer to allocated with central coordinates | 
|---|
| 50 | */ | 
|---|
| 51 | Vector *molecule::GetCenter() const | 
|---|
| 52 | { | 
|---|
| 53 | Vector *center = DetermineCenterOfAll(); | 
|---|
| 54 | return center; | 
|---|
| 55 | }; | 
|---|
| 56 |  | 
|---|
| 57 |  | 
|---|
| 58 | /** PointCloud implementation of GoPoint | 
|---|
| 59 | * Uses atoms and STL stuff. | 
|---|
| 60 | */ | 
|---|
| 61 | TesselPoint* molecule::GetPoint() const | 
|---|
| 62 | { | 
|---|
| 63 | return (*InternalPointer); | 
|---|
| 64 | }; | 
|---|
| 65 |  | 
|---|
| 66 | /** PointCloud implementation of GoToNext. | 
|---|
| 67 | * Uses atoms and STL stuff. | 
|---|
| 68 | */ | 
|---|
| 69 | void molecule::GoToNext() const | 
|---|
| 70 | { | 
|---|
| 71 | if (InternalPointer != atoms.end()) | 
|---|
| 72 | InternalPointer++; | 
|---|
| 73 | }; | 
|---|
| 74 |  | 
|---|
| 75 | /** PointCloud implementation of GoToFirst. | 
|---|
| 76 | * Uses atoms and STL stuff. | 
|---|
| 77 | */ | 
|---|
| 78 | void molecule::GoToFirst() const | 
|---|
| 79 | { | 
|---|
| 80 | // evil hack necessary because | 
|---|
| 81 | // -# although InternalPointer is mutable | 
|---|
| 82 | // -# only const_iterator begin() is called due to const in the function declaration above | 
|---|
| 83 | // -# and there is no cast from const_iterator to const iterator | 
|---|
| 84 | atomSet::const_iterator test = begin(); | 
|---|
| 85 | InternalPointer = *(reinterpret_cast<atomSet::iterator *>(&test)); | 
|---|
| 86 | }; | 
|---|
| 87 |  | 
|---|
| 88 | /** PointCloud implementation of IsEmpty. | 
|---|
| 89 | * Uses atoms and STL stuff. | 
|---|
| 90 | */ | 
|---|
| 91 | bool molecule::IsEmpty() const | 
|---|
| 92 | { | 
|---|
| 93 | return (empty()); | 
|---|
| 94 | }; | 
|---|
| 95 |  | 
|---|
| 96 | /** PointCloud implementation of IsLast. | 
|---|
| 97 | * Uses atoms and STL stuff. | 
|---|
| 98 | */ | 
|---|
| 99 | bool molecule::IsEnd() const | 
|---|
| 100 | { | 
|---|
| 101 | return (InternalPointer == atoms.end()); | 
|---|
| 102 | }; | 
|---|
| 103 |  | 
|---|
| 104 | int molecule::GetMaxId() const { | 
|---|
| 105 | return getAtomCount(); | 
|---|
| 106 | } | 
|---|