| 1 | /* | 
|---|
| 2 | * Project: MoleCuilder | 
|---|
| 3 | * Description: creates and alters molecular systems | 
|---|
| 4 | * Copyright (C)  2010-2012 University of Bonn. 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 | /* | 
|---|
| 24 | * BoundaryPointSet.cpp | 
|---|
| 25 | * | 
|---|
| 26 | *  Created on: Jul 29, 2010 | 
|---|
| 27 | *      Author: heber | 
|---|
| 28 | */ | 
|---|
| 29 |  | 
|---|
| 30 | // include config.h | 
|---|
| 31 | #ifdef HAVE_CONFIG_H | 
|---|
| 32 | #include <config.h> | 
|---|
| 33 | #endif | 
|---|
| 34 |  | 
|---|
| 35 | #include "CodePatterns/MemDebug.hpp" | 
|---|
| 36 |  | 
|---|
| 37 | #include "BoundaryPointSet.hpp" | 
|---|
| 38 |  | 
|---|
| 39 | #include <iostream> | 
|---|
| 40 |  | 
|---|
| 41 | #include "BoundaryLineSet.hpp" | 
|---|
| 42 | #include "BoundaryTriangleSet.hpp" | 
|---|
| 43 | #include "Atom/TesselPoint.hpp" | 
|---|
| 44 |  | 
|---|
| 45 | #include "CodePatterns/Assert.hpp" | 
|---|
| 46 | #include "CodePatterns/Info.hpp" | 
|---|
| 47 | #include "CodePatterns/Log.hpp" | 
|---|
| 48 | #include "LinearAlgebra/Vector.hpp" | 
|---|
| 49 | #include "CodePatterns/Verbose.hpp" | 
|---|
| 50 |  | 
|---|
| 51 | using namespace std; | 
|---|
| 52 |  | 
|---|
| 53 | // ======================================== Points on Boundary ================================= | 
|---|
| 54 |  | 
|---|
| 55 | /** Constructor of BoundaryPointSet. | 
|---|
| 56 | */ | 
|---|
| 57 | BoundaryPointSet::BoundaryPointSet() : | 
|---|
| 58 | LinesCount(0), | 
|---|
| 59 | value(0.), | 
|---|
| 60 | Nr(-1) | 
|---|
| 61 | { | 
|---|
| 62 | //Info FunctionInfo(__func__); | 
|---|
| 63 | LOG(1, "Adding noname."); | 
|---|
| 64 | } | 
|---|
| 65 | ; | 
|---|
| 66 |  | 
|---|
| 67 | /** Constructor of BoundaryPointSet with Tesselpoint. | 
|---|
| 68 | * \param *Walker TesselPoint this boundary point represents | 
|---|
| 69 | */ | 
|---|
| 70 | BoundaryPointSet::BoundaryPointSet(TesselPoint * const Walker) : | 
|---|
| 71 | LinesCount(0), | 
|---|
| 72 | node(Walker), | 
|---|
| 73 | value(0.), | 
|---|
| 74 | Nr(Walker->getNr()) | 
|---|
| 75 | { | 
|---|
| 76 | //Info FunctionInfo(__func__); | 
|---|
| 77 | LOG(5, "DEBUG: Adding Node " << *Walker); | 
|---|
| 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 | { | 
|---|
| 87 | //Info FunctionInfo(__func__); | 
|---|
| 88 | //LOG(0, "Erasing point Nr. " << Nr << "."); | 
|---|
| 89 | if (!lines.empty()) | 
|---|
| 90 | ELOG(2, "Memory Leak! I " << *this << " am still connected to some lines."); | 
|---|
| 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 | { | 
|---|
| 100 | //Info FunctionInfo(__func__); | 
|---|
| 101 | LOG(5, "DEBUG: Adding " << *this << " to line " << *line << "."); | 
|---|
| 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 |  | 
|---|
| 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."); | 
|---|
| 125 | return node->getNr(); | 
|---|
| 126 | } | 
|---|
| 127 |  | 
|---|
| 128 | TesselPoint *BoundaryPointSet::getTesselPoint() | 
|---|
| 129 | { | 
|---|
| 130 | return node; | 
|---|
| 131 | } | 
|---|
| 132 |  | 
|---|
| 133 | ostream & operator <<(ostream &ost, const BoundaryPointSet &a) | 
|---|
| 134 | { | 
|---|
| 135 | ost << "[" << a.Nr << "|" << a.getName() << " at " << a.getPosition() << "]"; | 
|---|
| 136 | return ost; | 
|---|
| 137 | } | 
|---|
| 138 |  | 
|---|