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 | /*
|
---|
9 | * BoundaryPointSet.cpp
|
---|
10 | *
|
---|
11 | * Created on: Jul 29, 2010
|
---|
12 | * Author: heber
|
---|
13 | */
|
---|
14 |
|
---|
15 | // include config.h
|
---|
16 | #ifdef HAVE_CONFIG_H
|
---|
17 | #include <config.h>
|
---|
18 | #endif
|
---|
19 |
|
---|
20 | #include "CodePatterns/MemDebug.hpp"
|
---|
21 |
|
---|
22 | #include "BoundaryPointSet.hpp"
|
---|
23 |
|
---|
24 | #include <iostream>
|
---|
25 |
|
---|
26 | #include "BoundaryLineSet.hpp"
|
---|
27 | #include "BoundaryTriangleSet.hpp"
|
---|
28 | #include "TesselPoint.hpp"
|
---|
29 |
|
---|
30 | #include "CodePatterns/Assert.hpp"
|
---|
31 | #include "CodePatterns/Info.hpp"
|
---|
32 | #include "CodePatterns/Log.hpp"
|
---|
33 | #include "LinearAlgebra/Vector.hpp"
|
---|
34 | #include "CodePatterns/Verbose.hpp"
|
---|
35 |
|
---|
36 | using namespace std;
|
---|
37 |
|
---|
38 | // ======================================== Points on Boundary =================================
|
---|
39 |
|
---|
40 | /** Constructor of BoundaryPointSet.
|
---|
41 | */
|
---|
42 | BoundaryPointSet::BoundaryPointSet() :
|
---|
43 | LinesCount(0),
|
---|
44 | value(0.),
|
---|
45 | Nr(-1)
|
---|
46 | {
|
---|
47 | Info FunctionInfo(__func__);
|
---|
48 | DoLog(1) && (Log() << Verbose(1) << "Adding noname." << endl);
|
---|
49 | }
|
---|
50 | ;
|
---|
51 |
|
---|
52 | /** Constructor of BoundaryPointSet with Tesselpoint.
|
---|
53 | * \param *Walker TesselPoint this boundary point represents
|
---|
54 | */
|
---|
55 | BoundaryPointSet::BoundaryPointSet(TesselPoint * const Walker) :
|
---|
56 | LinesCount(0),
|
---|
57 | node(Walker),
|
---|
58 | value(0.),
|
---|
59 | Nr(Walker->ParticleInfo_nr)
|
---|
60 | {
|
---|
61 | Info FunctionInfo(__func__);
|
---|
62 | DoLog(1) && (Log() << Verbose(1) << "Adding Node " << *Walker << endl);
|
---|
63 | }
|
---|
64 | ;
|
---|
65 |
|
---|
66 | /** Destructor of BoundaryPointSet.
|
---|
67 | * Sets node to NULL to avoid removing the original, represented TesselPoint.
|
---|
68 | * \note When removing point from a class Tesselation, use RemoveTesselationPoint()
|
---|
69 | */
|
---|
70 | BoundaryPointSet::~BoundaryPointSet()
|
---|
71 | {
|
---|
72 | Info FunctionInfo(__func__);
|
---|
73 | //Log() << Verbose(0) << "Erasing point ParticleInfo_nr. " << Nr << "." << endl;
|
---|
74 | if (!lines.empty())
|
---|
75 | DoeLog(2) && (eLog() << Verbose(2) << "Memory Leak! I " << *this << " am still connected to some lines." << endl);
|
---|
76 | node = NULL;
|
---|
77 | }
|
---|
78 | ;
|
---|
79 |
|
---|
80 | /** Add a line to the LineMap of this point.
|
---|
81 | * \param *line line to add
|
---|
82 | */
|
---|
83 | void BoundaryPointSet::AddLine(BoundaryLineSet * const line)
|
---|
84 | {
|
---|
85 | Info FunctionInfo(__func__);
|
---|
86 | DoLog(1) && (Log() << Verbose(1) << "Adding " << *this << " to line " << *line << "." << endl);
|
---|
87 | if (line->endpoints[0] == this) {
|
---|
88 | lines.insert(LinePair(line->endpoints[1]->Nr, line));
|
---|
89 | } else {
|
---|
90 | lines.insert(LinePair(line->endpoints[0]->Nr, line));
|
---|
91 | }
|
---|
92 | LinesCount++;
|
---|
93 | }
|
---|
94 |
|
---|
95 | const std::string& BoundaryPointSet::getName() const
|
---|
96 | {
|
---|
97 | ASSERT(node != NULL, "BoundaryPointSet::getName() - internal node is NULL.");
|
---|
98 | return node->getName();
|
---|
99 | }
|
---|
100 |
|
---|
101 | const Vector& BoundaryPointSet::getPosition() const
|
---|
102 | {
|
---|
103 | ASSERT(node != NULL, "BoundaryPointSet::getPosition() - internal node is NULL.");
|
---|
104 | return node->getPosition();
|
---|
105 | }
|
---|
106 |
|
---|
107 | const int& BoundaryPointSet::getNr() const
|
---|
108 | {
|
---|
109 | ASSERT(node != NULL, "BoundaryPointSet::getPosition() - internal node is NULL.");
|
---|
110 | return node->ParticleInfo_nr;
|
---|
111 | }
|
---|
112 |
|
---|
113 | TesselPoint *BoundaryPointSet::getTesselPoint()
|
---|
114 | {
|
---|
115 | return node;
|
---|
116 | }
|
---|
117 |
|
---|
118 | ostream & operator <<(ostream &ost, const BoundaryPointSet &a)
|
---|
119 | {
|
---|
120 | ost << "[" << a.Nr << "|" << a.getName() << " at " << a.getPosition() << "]";
|
---|
121 | return ost;
|
---|
122 | }
|
---|
123 |
|
---|