| [31bd9c] | 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 | }
 | 
|---|