| 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 | } | 
|---|