| 1 | /*
 | 
|---|
| 2 |  * Project: MoleCuilder
 | 
|---|
| 3 |  * Description: creates and alters molecular systems
 | 
|---|
| 4 |  * Copyright (C)  2014 Frederik Heber. All rights reserved.
 | 
|---|
| 5 |  * 
 | 
|---|
| 6 |  *
 | 
|---|
| 7 |  *   This file is part of MoleCuilder.
 | 
|---|
| 8 |  *
 | 
|---|
| 9 |  *    MoleCuilder is free software: you can redistribute it and/or modify
 | 
|---|
| 10 |  *    it under the terms of the GNU General Public License as published by
 | 
|---|
| 11 |  *    the Free Software Foundation, either version 2 of the License, or
 | 
|---|
| 12 |  *    (at your option) any later version.
 | 
|---|
| 13 |  *
 | 
|---|
| 14 |  *    MoleCuilder is distributed in the hope that it will be useful,
 | 
|---|
| 15 |  *    but WITHOUT ANY WARRANTY; without even the implied warranty of
 | 
|---|
| 16 |  *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 | 
|---|
| 17 |  *    GNU General Public License for more details.
 | 
|---|
| 18 |  *
 | 
|---|
| 19 |  *    You should have received a copy of the GNU General Public License
 | 
|---|
| 20 |  *    along with MoleCuilder.  If not, see <http://www.gnu.org/licenses/>.
 | 
|---|
| 21 |  */
 | 
|---|
| 22 | 
 | 
|---|
| 23 | /** \file ion.cpp
 | 
|---|
| 24 |  * 
 | 
|---|
| 25 |  * Function implementations for the class ion.
 | 
|---|
| 26 |  * 
 | 
|---|
| 27 |  */
 | 
|---|
| 28 | 
 | 
|---|
| 29 | // include config.h
 | 
|---|
| 30 | #ifdef HAVE_CONFIG_H
 | 
|---|
| 31 | #include <config.h>
 | 
|---|
| 32 | #endif
 | 
|---|
| 33 | 
 | 
|---|
| 34 | //#include "CodePatterns/MemDebug.hpp"
 | 
|---|
| 35 | 
 | 
|---|
| 36 | #include "CodePatterns/Assert.hpp"
 | 
|---|
| 37 | #include "CodePatterns/Log.hpp"
 | 
|---|
| 38 | #include "ion.hpp"
 | 
|---|
| 39 | 
 | 
|---|
| 40 | /************************************* Functions for class ion **********************************/
 | 
|---|
| 41 | ion::ion(const element &_element, const int _ionization) :
 | 
|---|
| 42 |     element(_element),
 | 
|---|
| 43 |     ionization(_ionization)
 | 
|---|
| 44 | {
 | 
|---|
| 45 |   // only outer shell can be ionized
 | 
|---|
| 46 |   ASSERT( abs(ionization) <= NoValenceOrbitals,
 | 
|---|
| 47 |       "ion::ion() - ionization ("+toString(ionization)+
 | 
|---|
| 48 |       ") must not exceed valence ("+toString(NoValenceOrbitals)+").");
 | 
|---|
| 49 |   ASSERT (ionization != 0,
 | 
|---|
| 50 |      "ion::ion() - ionization must not be zero.");
 | 
|---|
| 51 | }
 | 
|---|
| 52 | 
 | 
|---|
| 53 | ion::ion(const element* const _element, const int _ionization) :
 | 
|---|
| 54 |     element(*_element),
 | 
|---|
| 55 |     ionization(_ionization)
 | 
|---|
| 56 | {
 | 
|---|
| 57 |   // only outer shell can be ionized
 | 
|---|
| 58 |   ASSERT( abs(ionization) <= NoValenceOrbitals,
 | 
|---|
| 59 |       "ion::ion() - ionization ("+toString(ionization)+
 | 
|---|
| 60 |       ") must not exceed valence ("+toString(NoValenceOrbitals)+").");
 | 
|---|
| 61 |   ASSERT (ionization != 0,
 | 
|---|
| 62 |      "ion::ion() - ionization must not be zero.");
 | 
|---|
| 63 | }
 | 
|---|
| 64 | 
 | 
|---|
| 65 | ion::ion() :
 | 
|---|
| 66 |     element(),
 | 
|---|
| 67 |     ionization(0)
 | 
|---|
| 68 | {
 | 
|---|
| 69 | 
 | 
|---|
| 70 | }
 | 
|---|
| 71 | 
 | 
|---|
| 72 | ion::~ion()
 | 
|---|
| 73 | {}
 | 
|---|
| 74 | 
 | 
|---|
| 75 | ion& ion::operator=(const element &other)
 | 
|---|
| 76 | {
 | 
|---|
| 77 |   if (this != &other) { // check for self-assignment
 | 
|---|
| 78 |     dynamic_cast<element &>(*this) = dynamic_cast<const element &>(other);
 | 
|---|
| 79 |   }
 | 
|---|
| 80 |   return *this;
 | 
|---|
| 81 | }
 | 
|---|
| 82 | 
 | 
|---|
| 83 | ion& ion::operator=(const ion &other)
 | 
|---|
| 84 | {
 | 
|---|
| 85 |   if (this != &other) { // check for self-assignment
 | 
|---|
| 86 |     dynamic_cast<element &>(*this) = dynamic_cast<const element &>(other);
 | 
|---|
| 87 |     const_cast<int &>(ionization) = other.ionization;
 | 
|---|
| 88 |   }
 | 
|---|
| 89 |   return *this;
 | 
|---|
| 90 | }
 | 
|---|
| 91 | 
 | 
|---|
| 92 | atomicNumber_t ion::getAtomicNumber() const
 | 
|---|
| 93 | {
 | 
|---|
| 94 |   return Z - ionization;
 | 
|---|
| 95 | }
 | 
|---|
| 96 | 
 | 
|---|
| 97 | double ion::getCharge() const
 | 
|---|
| 98 | {
 | 
|---|
| 99 |   return (double)ionization;
 | 
|---|
| 100 | }
 | 
|---|
| 101 | 
 | 
|---|
| 102 | double ion::getValence() const
 | 
|---|
| 103 | {
 | 
|---|
| 104 |   return Valence - ionization;
 | 
|---|
| 105 | }
 | 
|---|
| 106 | 
 | 
|---|
| 107 | int ion::getNoValenceOrbitals() const
 | 
|---|
| 108 | {
 | 
|---|
| 109 |   return NoValenceOrbitals - ionization;
 | 
|---|
| 110 | }
 | 
|---|
| 111 | 
 | 
|---|
| 112 | bool ion::operator==(const element &other) const {
 | 
|---|
| 113 |   return dynamic_cast<const element &>(*this) == other;
 | 
|---|
| 114 | }
 | 
|---|
| 115 | bool ion::operator!=(const element &other) const {
 | 
|---|
| 116 |   return !(*this == other);
 | 
|---|
| 117 | }
 | 
|---|
| 118 | bool ion::operator==(const ion &other) const {
 | 
|---|
| 119 |   bool status = true;
 | 
|---|
| 120 |   status &= dynamic_cast<const element &>(*this) == other;
 | 
|---|
| 121 |   status &= ionization == other.ionization;
 | 
|---|
| 122 |   return status;
 | 
|---|
| 123 | }
 | 
|---|