| [bcf653] | 1 | /* | 
|---|
|  | 2 | * Project: MoleCuilder | 
|---|
|  | 3 | * Description: creates and alters molecular systems | 
|---|
| [0aa122] | 4 | * Copyright (C)  2010-2012 University of Bonn. All rights reserved. | 
|---|
| [94d5ac6] | 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/>. | 
|---|
| [bcf653] | 21 | */ | 
|---|
|  | 22 |  | 
|---|
| [d74077] | 23 | /* | 
|---|
|  | 24 | * BoundaryPointSet.cpp | 
|---|
|  | 25 | * | 
|---|
|  | 26 | *  Created on: Jul 29, 2010 | 
|---|
|  | 27 | *      Author: heber | 
|---|
|  | 28 | */ | 
|---|
|  | 29 |  | 
|---|
| [bf3817] | 30 | // include config.h | 
|---|
|  | 31 | #ifdef HAVE_CONFIG_H | 
|---|
|  | 32 | #include <config.h> | 
|---|
|  | 33 | #endif | 
|---|
|  | 34 |  | 
|---|
| [9eb71b3] | 35 | //#include "CodePatterns/MemDebug.hpp" | 
|---|
| [bbbad5] | 36 |  | 
|---|
| [d74077] | 37 | #include "BoundaryPointSet.hpp" | 
|---|
|  | 38 |  | 
|---|
|  | 39 | #include <iostream> | 
|---|
|  | 40 |  | 
|---|
|  | 41 | #include "BoundaryLineSet.hpp" | 
|---|
|  | 42 | #include "BoundaryTriangleSet.hpp" | 
|---|
| [6f0841] | 43 | #include "Atom/TesselPoint.hpp" | 
|---|
| [d74077] | 44 |  | 
|---|
| [ad011c] | 45 | #include "CodePatterns/Assert.hpp" | 
|---|
|  | 46 | #include "CodePatterns/Info.hpp" | 
|---|
|  | 47 | #include "CodePatterns/Log.hpp" | 
|---|
| [8f4df1] | 48 | #include "LinearAlgebra/Vector.hpp" | 
|---|
| [ad011c] | 49 | #include "CodePatterns/Verbose.hpp" | 
|---|
| [d74077] | 50 |  | 
|---|
|  | 51 | using namespace std; | 
|---|
|  | 52 |  | 
|---|
|  | 53 | // ======================================== Points on Boundary ================================= | 
|---|
|  | 54 |  | 
|---|
|  | 55 | /** Constructor of BoundaryPointSet. | 
|---|
|  | 56 | */ | 
|---|
|  | 57 | BoundaryPointSet::BoundaryPointSet() : | 
|---|
| [97b825] | 58 | LinesCount(0), | 
|---|
|  | 59 | value(0.), | 
|---|
|  | 60 | Nr(-1) | 
|---|
| [d74077] | 61 | { | 
|---|
| [ce7bfd] | 62 | //Info FunctionInfo(__func__); | 
|---|
| [47d041] | 63 | LOG(1, "Adding noname."); | 
|---|
| [d74077] | 64 | } | 
|---|
|  | 65 | ; | 
|---|
|  | 66 |  | 
|---|
|  | 67 | /** Constructor of BoundaryPointSet with Tesselpoint. | 
|---|
|  | 68 | * \param *Walker TesselPoint this boundary point represents | 
|---|
|  | 69 | */ | 
|---|
|  | 70 | BoundaryPointSet::BoundaryPointSet(TesselPoint * const Walker) : | 
|---|
| [97b825] | 71 | LinesCount(0), | 
|---|
|  | 72 | node(Walker), | 
|---|
|  | 73 | value(0.), | 
|---|
| [735b1c] | 74 | Nr(Walker->getNr()) | 
|---|
| [d74077] | 75 | { | 
|---|
| [ce7bfd] | 76 | //Info FunctionInfo(__func__); | 
|---|
|  | 77 | LOG(5, "DEBUG: Adding Node " << *Walker); | 
|---|
| [d74077] | 78 | } | 
|---|
|  | 79 | ; | 
|---|
|  | 80 |  | 
|---|
|  | 81 | /** Destructor of BoundaryPointSet. | 
|---|
|  | 82 | * Sets node to NULL to avoid removing the original, represented TesselPoint. | 
|---|
|  | 83 | * \note When removing point from a class Tesselation, use RemoveTesselationPoint() | 
|---|
|  | 84 | */ | 
|---|
|  | 85 | BoundaryPointSet::~BoundaryPointSet() | 
|---|
|  | 86 | { | 
|---|
| [ce7bfd] | 87 | //Info FunctionInfo(__func__); | 
|---|
| [47d041] | 88 | //LOG(0, "Erasing point Nr. " << Nr << "."); | 
|---|
| [d74077] | 89 | if (!lines.empty()) | 
|---|
| [47d041] | 90 | ELOG(2, "Memory Leak! I " << *this << " am still connected to some lines."); | 
|---|
| [d74077] | 91 | node = NULL; | 
|---|
|  | 92 | } | 
|---|
|  | 93 | ; | 
|---|
|  | 94 |  | 
|---|
|  | 95 | /** Add a line to the LineMap of this point. | 
|---|
|  | 96 | * \param *line line to add | 
|---|
|  | 97 | */ | 
|---|
|  | 98 | void BoundaryPointSet::AddLine(BoundaryLineSet * const line) | 
|---|
|  | 99 | { | 
|---|
| [ce7bfd] | 100 | //Info FunctionInfo(__func__); | 
|---|
|  | 101 | LOG(5, "DEBUG: Adding " << *this << " to line " << *line << "."); | 
|---|
| [d74077] | 102 | if (line->endpoints[0] == this) { | 
|---|
|  | 103 | lines.insert(LinePair(line->endpoints[1]->Nr, line)); | 
|---|
|  | 104 | } else { | 
|---|
|  | 105 | lines.insert(LinePair(line->endpoints[0]->Nr, line)); | 
|---|
|  | 106 | } | 
|---|
|  | 107 | LinesCount++; | 
|---|
|  | 108 | } | 
|---|
|  | 109 |  | 
|---|
| [ba5215] | 110 | const std::string& BoundaryPointSet::getName() const | 
|---|
|  | 111 | { | 
|---|
|  | 112 | ASSERT(node != NULL, "BoundaryPointSet::getName() - internal node is NULL."); | 
|---|
|  | 113 | return node->getName(); | 
|---|
|  | 114 | } | 
|---|
|  | 115 |  | 
|---|
|  | 116 | const Vector& BoundaryPointSet::getPosition() const | 
|---|
|  | 117 | { | 
|---|
|  | 118 | ASSERT(node != NULL, "BoundaryPointSet::getPosition() - internal node is NULL."); | 
|---|
|  | 119 | return node->getPosition(); | 
|---|
|  | 120 | } | 
|---|
|  | 121 |  | 
|---|
|  | 122 | const int& BoundaryPointSet::getNr() const | 
|---|
|  | 123 | { | 
|---|
|  | 124 | ASSERT(node != NULL, "BoundaryPointSet::getPosition() - internal node is NULL."); | 
|---|
| [735b1c] | 125 | return node->getNr(); | 
|---|
| [ba5215] | 126 | } | 
|---|
|  | 127 |  | 
|---|
|  | 128 | TesselPoint *BoundaryPointSet::getTesselPoint() | 
|---|
|  | 129 | { | 
|---|
|  | 130 | return node; | 
|---|
|  | 131 | } | 
|---|
|  | 132 |  | 
|---|
| [d74077] | 133 | ostream & operator <<(ostream &ost, const BoundaryPointSet &a) | 
|---|
|  | 134 | { | 
|---|
| [b8d215] | 135 | ost << "[" << a.Nr << "|" << a.getName() << "]"; // " at " << a.getPosition() << "]"; | 
|---|
| [d74077] | 136 | return ost; | 
|---|
|  | 137 | } | 
|---|
| [ba5215] | 138 |  | 
|---|