1 | /** \file bond.cpp
|
---|
2 | *
|
---|
3 | * Function implementations for the classes BondLeaf, BondTree and bond.
|
---|
4 | *
|
---|
5 | */
|
---|
6 |
|
---|
7 | #include "molecules.hpp"
|
---|
8 |
|
---|
9 |
|
---|
10 | /***************************************** Functions for class bond ********************************/
|
---|
11 |
|
---|
12 | /** Empty Constructor for class bond.
|
---|
13 | */
|
---|
14 | bond::bond()
|
---|
15 | {
|
---|
16 | leftatom = NULL;
|
---|
17 | rightatom = NULL;
|
---|
18 | previous = NULL;
|
---|
19 | next = NULL;
|
---|
20 | nr = -1;
|
---|
21 | HydrogenBond = 0;
|
---|
22 | BondDegree = 0;
|
---|
23 | Used = white;
|
---|
24 | Cyclic = false;
|
---|
25 | Type = Undetermined;
|
---|
26 | };
|
---|
27 |
|
---|
28 | /** Constructor for class bond, taking right and left bond partner
|
---|
29 | * \param *left left atom
|
---|
30 | * \param *right right atom
|
---|
31 | * \param degree bond degree
|
---|
32 | * \param number increasing index
|
---|
33 | */
|
---|
34 | bond::bond(atom *left, atom *right, int degree=1, int number=0)
|
---|
35 | {
|
---|
36 | leftatom = left;
|
---|
37 | rightatom = right;
|
---|
38 | previous = NULL;
|
---|
39 | next = NULL;
|
---|
40 | HydrogenBond = 0;
|
---|
41 | if ((left != NULL) && (right != NULL)) {
|
---|
42 | if ((left->type != NULL) && (left->type->Z == 1))
|
---|
43 | HydrogenBond++;
|
---|
44 | if ((right->type != NULL) && (right->type->Z == 1))
|
---|
45 | HydrogenBond++;
|
---|
46 | }
|
---|
47 | BondDegree = degree;
|
---|
48 | nr = number;
|
---|
49 | Used = white;
|
---|
50 | Cyclic = false;
|
---|
51 | };
|
---|
52 | bond::bond(atom *left, atom *right)
|
---|
53 | {
|
---|
54 | leftatom = left;
|
---|
55 | rightatom = right;
|
---|
56 | previous = NULL;
|
---|
57 | next = NULL;
|
---|
58 | HydrogenBond = 0;
|
---|
59 | if ((left != NULL) && (right != NULL)) {
|
---|
60 | if ((left->type != NULL) && (left->type->Z == 1))
|
---|
61 | HydrogenBond++;
|
---|
62 | if ((right->type != NULL) && (right->type->Z == 1))
|
---|
63 | HydrogenBond++;
|
---|
64 | }
|
---|
65 | BondDegree = 1;
|
---|
66 | nr = 0;
|
---|
67 | Used = white;
|
---|
68 | Cyclic = false;
|
---|
69 | };
|
---|
70 |
|
---|
71 | /** Empty Destructor for class bond.
|
---|
72 | */
|
---|
73 | bond::~bond()
|
---|
74 | {
|
---|
75 | // remove this node from the list structure
|
---|
76 | if (previous != NULL) {
|
---|
77 | previous->next = next;
|
---|
78 | }
|
---|
79 | if (next != NULL) {
|
---|
80 | next->previous = previous;
|
---|
81 | }
|
---|
82 | };
|
---|
83 |
|
---|
84 | ostream & operator << (ostream &ost, bond &b)
|
---|
85 | {
|
---|
86 | ost << "[" << b.leftatom->Name << " <" << b.BondDegree << "(H" << b.HydrogenBond << ")>" << b.rightatom->Name << "]";
|
---|
87 | return ost;
|
---|
88 | };
|
---|
89 |
|
---|
90 | /** Get the other atom in a bond if one is specified.
|
---|
91 | * \param *Atom the pointer to the one atom
|
---|
92 | * \return pointer to the other atom in the bond, NULL if no match (indicates something's wrong with the bond)
|
---|
93 | */
|
---|
94 | atom * bond::GetOtherAtom(atom *Atom) const
|
---|
95 | {
|
---|
96 | if(leftatom == Atom)
|
---|
97 | return rightatom;
|
---|
98 | if(rightatom == Atom)
|
---|
99 | return leftatom;
|
---|
100 | return NULL;
|
---|
101 | };
|
---|
102 |
|
---|
103 | /** Get the other atom in a bond if one is specified.
|
---|
104 | * \param *Atom the pointer to the one atom
|
---|
105 | * \return pointer to the other atom in the bond, NULL if no match (indicates something's wrong with the bond)
|
---|
106 | */
|
---|
107 | bond * bond::GetFirstBond()
|
---|
108 | {
|
---|
109 | return GetFirst(this);
|
---|
110 | };
|
---|
111 |
|
---|
112 | /** Get the other atom in a bond if one is specified.
|
---|
113 | * \param *Atom the pointer to the one atom
|
---|
114 | * \return pointer to the other atom in the bond, NULL if no match (indicates something's wrong with the bond)
|
---|
115 | */
|
---|
116 | bond * bond::GetLastBond()
|
---|
117 | {
|
---|
118 | return GetLast(this);
|
---|
119 | };
|
---|
120 |
|
---|
121 | /** Returns whether vertex was used in DFS.
|
---|
122 | * \return bond::Used
|
---|
123 | */
|
---|
124 | enum Shading bond::IsUsed()
|
---|
125 | {
|
---|
126 | return Used;
|
---|
127 | };
|
---|
128 |
|
---|
129 | /** Checks if an atom exists in a bond.
|
---|
130 | * \param *ptr pointer to atom
|
---|
131 | * \return true if it is either bond::leftatom or bond::rightatom, false otherwise
|
---|
132 | */
|
---|
133 | bool bond::Contains(const atom *ptr)
|
---|
134 | {
|
---|
135 | return ((leftatom == ptr) || (rightatom == ptr));
|
---|
136 | };
|
---|
137 |
|
---|
138 | /** Checks if an atom exists in a bond.
|
---|
139 | * \param nr index of atom
|
---|
140 | * \return true if it is either bond::leftatom or bond::rightatom, false otherwise
|
---|
141 | */
|
---|
142 | bool bond::Contains(const int number)
|
---|
143 | {
|
---|
144 | return ((leftatom->nr == number) || (rightatom->nr == number));
|
---|
145 | };
|
---|
146 |
|
---|
147 | /** Masks vertex as used in DFS.
|
---|
148 | * \return bond::Used, false if bond was already marked used
|
---|
149 | */
|
---|
150 | bool bond::MarkUsed(enum Shading color) {
|
---|
151 | if (Used == black) {
|
---|
152 | cerr << "ERROR: Bond " << this << " was already marked black!." << endl;
|
---|
153 | return false;
|
---|
154 | } else {
|
---|
155 | Used = color;
|
---|
156 | return true;
|
---|
157 | }
|
---|
158 | };
|
---|
159 |
|
---|
160 | /** Resets used flag in DFS.
|
---|
161 | * \return bond::Used
|
---|
162 | */
|
---|
163 | void bond::ResetUsed() {
|
---|
164 | Used = white;
|
---|
165 | };
|
---|