| 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 "CodePatterns/MemDebug.hpp"
 | 
|---|
| 21 | 
 | 
|---|
| 22 | #include "atom.hpp"
 | 
|---|
| 23 | #include "config.hpp"
 | 
|---|
| 24 | #include "CodePatterns/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 | }
 | 
|---|