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 | /** \file element.cpp
|
---|
9 | *
|
---|
10 | * Function implementations for the class element.
|
---|
11 | *
|
---|
12 | */
|
---|
13 |
|
---|
14 | // include config.h
|
---|
15 | #ifdef HAVE_CONFIG_H
|
---|
16 | #include <config.h>
|
---|
17 | #endif
|
---|
18 |
|
---|
19 | #include "CodePatterns/MemDebug.hpp"
|
---|
20 |
|
---|
21 | #include <iomanip>
|
---|
22 | #include <fstream>
|
---|
23 |
|
---|
24 | #include "CodePatterns/Assert.hpp"
|
---|
25 | #include "element.hpp"
|
---|
26 |
|
---|
27 | using namespace std;
|
---|
28 |
|
---|
29 | /************************************* Functions for class element **********************************/
|
---|
30 |
|
---|
31 | /** Constructor of class element.
|
---|
32 | */
|
---|
33 | element::element() :
|
---|
34 | mass(0),
|
---|
35 | CovalentRadius(0),
|
---|
36 | VanDerWaalsRadius(0),
|
---|
37 | Z(-1),
|
---|
38 | Valence(0),
|
---|
39 | NoValenceOrbitals(0)
|
---|
40 | {
|
---|
41 | };
|
---|
42 |
|
---|
43 | element::element(const element &src) :
|
---|
44 | mass(src.mass),
|
---|
45 | CovalentRadius(src.CovalentRadius),
|
---|
46 | VanDerWaalsRadius(src.VanDerWaalsRadius),
|
---|
47 | Z(src.Z),
|
---|
48 | Valence(src.Valence),
|
---|
49 | NoValenceOrbitals(src.NoValenceOrbitals),
|
---|
50 | name(src.name),
|
---|
51 | symbol(src.symbol)
|
---|
52 | {
|
---|
53 | period = src.period;
|
---|
54 | group = src.group;
|
---|
55 | block = src.block;
|
---|
56 | }
|
---|
57 |
|
---|
58 | /** Destructor of class element.
|
---|
59 | */
|
---|
60 | element::~element() {};
|
---|
61 |
|
---|
62 | element &element::operator=(const element &src){
|
---|
63 | if(this!=&src){
|
---|
64 | mass=src.mass;
|
---|
65 | CovalentRadius=src.CovalentRadius;
|
---|
66 | VanDerWaalsRadius=src.VanDerWaalsRadius;
|
---|
67 | Z=src.Z;
|
---|
68 | Valence=src.Valence;
|
---|
69 | NoValenceOrbitals=src.NoValenceOrbitals;
|
---|
70 | name=src.name;
|
---|
71 | symbol=src.symbol;
|
---|
72 | period = src.period;
|
---|
73 | group = src.group;
|
---|
74 | block = src.block;
|
---|
75 | }
|
---|
76 | return *this;
|
---|
77 | }
|
---|
78 |
|
---|
79 | /** Prints element data to \a *out.
|
---|
80 | * \param *out outstream
|
---|
81 | */
|
---|
82 | bool element::Output(ostream * const out) const
|
---|
83 | {
|
---|
84 | if (out != NULL) {
|
---|
85 | *out << name << "\t" << symbol << "\t" << period << "\t" << group << "\t" << block << "\t" << Z << "\t" << mass << "\t" << CovalentRadius << "\t" << VanDerWaalsRadius << endl;
|
---|
86 | //*out << Z << "\t" << fixed << setprecision(11) << showpoint << mass << "g/mol\t" << name << "\t" << symbol << "\t" << endl;
|
---|
87 | return true;
|
---|
88 | } else
|
---|
89 | return false;
|
---|
90 | };
|
---|
91 |
|
---|
92 | /** Prints element data to \a *out.
|
---|
93 | * \param *out outstream
|
---|
94 | * \param No cardinal number of element
|
---|
95 | * \param NoOfAtoms total number of atom of this element type
|
---|
96 | */
|
---|
97 | bool element::Checkout(ostream * const out, const int Number, const int NoOfAtoms) const
|
---|
98 | {
|
---|
99 | if (out != NULL) {
|
---|
100 | *out << "Ion_Type" << Number << "\t" << NoOfAtoms << "\t" << Z << "\t1.0\t3\t3\t" << fixed << setprecision(11) << showpoint << mass << "\t" << name << "\t" << symbol <<endl;
|
---|
101 | return true;
|
---|
102 | } else
|
---|
103 | return false;
|
---|
104 | };
|
---|
105 |
|
---|
106 | atomicNumber_t element::getNumber() const{
|
---|
107 | return Z;
|
---|
108 | }
|
---|
109 |
|
---|
110 | double element::getMass() const
|
---|
111 | {
|
---|
112 | return mass;
|
---|
113 | }
|
---|
114 |
|
---|
115 | double element::getCovalentRadius() const
|
---|
116 | {
|
---|
117 | return CovalentRadius;
|
---|
118 | }
|
---|
119 |
|
---|
120 | double element::getElectronegativity() const
|
---|
121 | {
|
---|
122 | return Electronegativity;
|
---|
123 | }
|
---|
124 |
|
---|
125 | double element::getVanDerWaalsRadius() const
|
---|
126 | {
|
---|
127 | return VanDerWaalsRadius;
|
---|
128 | }
|
---|
129 |
|
---|
130 | int element::getAtomicNumber() const
|
---|
131 | {
|
---|
132 | return Z;
|
---|
133 | }
|
---|
134 |
|
---|
135 | double element::getValence() const
|
---|
136 | {
|
---|
137 | return Valence;
|
---|
138 | }
|
---|
139 |
|
---|
140 | int element::getNoValenceOrbitals() const
|
---|
141 | {
|
---|
142 | return NoValenceOrbitals;
|
---|
143 | }
|
---|
144 |
|
---|
145 | double element::getHBondDistance(const int i) const
|
---|
146 | {
|
---|
147 | ASSERT((i>=0) && (i<3), "Access to element::HBondDistance out of bounds.");
|
---|
148 | return HBondDistance[i];
|
---|
149 | }
|
---|
150 |
|
---|
151 | double element::getHBondAngle(const int i) const
|
---|
152 | {
|
---|
153 | ASSERT((i>=0) && (i<3), "Access to element::HBondAngle out of bounds.");
|
---|
154 | return HBondAngle[i];
|
---|
155 | }
|
---|
156 |
|
---|
157 | string &element::getSymbol(){
|
---|
158 | return symbol;
|
---|
159 | }
|
---|
160 |
|
---|
161 | const string &element::getSymbol() const{
|
---|
162 | return symbol;
|
---|
163 | }
|
---|
164 |
|
---|
165 | void element::setSymbol(const std::string &temp)
|
---|
166 | {
|
---|
167 | symbol = temp;
|
---|
168 | }
|
---|
169 |
|
---|
170 | std::string &element::getName(){
|
---|
171 | return name;
|
---|
172 | }
|
---|
173 |
|
---|
174 | const std::string &element::getName() const{
|
---|
175 | return name;
|
---|
176 | }
|
---|
177 |
|
---|
178 | void element::setName(const std::string &temp)
|
---|
179 | {
|
---|
180 | name = temp;
|
---|
181 | }
|
---|
182 |
|
---|
183 | std::ostream &operator<<(std::ostream &ost,const element &elem){
|
---|
184 | ost << elem.getName() << "(" << elem.getNumber() << ")";
|
---|
185 | return ost;
|
---|
186 | }
|
---|