1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
/////////////////////////////////////////////////////////////
// //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing //
// http://earth.uq.edu.au/centre-geoscience-computing //
// //
// Primary Business: Brisbane, Queensland, Australia //
// Licensed under the Open Software License version 3.0 //
// http://www.apache.org/licenses/LICENSE-2.0 //
// //
/////////////////////////////////////////////////////////////
#include "MeshReader.h"
#include "MeshReaderImpl.h"
// -- STL includes --
#include <stdexcept>
using std::runtime_error;
namespace esys
{
namespace lsm //!Lattice Solid Model namespace.
{
// ==== NodeReader functions ==========
/*!
construct NodeReader
\param instream the stream to read from
*/
NodeReader::NodeReader(std::istream& instream)
: IterativeReader<IStreamIterator<MeshNodeData> >(instream)
{}
/*!
initialize node reader
*/
void NodeReader::initialise()
{
string token;
while ((!getIStream().eof()) && (token != "3D-Nodes")) {
getIStream() >> token;
}
if (token != "3D-Nodes") {
throw runtime_error("Could not find '3D-Nodes' token in stream.");
}
int nr_nodes = 0;
getIStream() >> nr_nodes;
setNumElements(nr_nodes);
IterativeReader<IStreamIterator<MeshNodeData> >::initialise();
}
// ==== TriReader Functions ===========
/*!
construct TriReader
\param instream the stream to read from
*/
TriReader::TriReader(std::istream& instream)
: IterativeReader<IStreamIterator<MeshTriData> >(instream)
{}
/*!
initialize triangle reader
*/
void TriReader::initialise()
{
string token;
while ((!getIStream().eof()) && (token != "Tri3")) {
getIStream() >> token;
}
if (token != "Tri3") {
throw runtime_error("Could not find 'Tri3' token in stream.");
}
int nr_nodes = 0;
getIStream() >> nr_nodes;
setNumElements(nr_nodes);
IterativeReader<IStreamIterator<MeshTriData> >::initialise();
}
// ==== MeshReader functions ==========
/*!
construct a mesh reader to read from file
\param filename the name of the file to read from
*/
MeshReader::MeshReader(const string& filename)
: m_impl_ptr(new Impl(filename))
{
m_impl_ptr->initialise();
}
/*!
destroy a mesh reader. Closes all accociated files.
*/
MeshReader::~MeshReader()
{
delete m_impl_ptr;
}
/*!
return iterator for nodes in file/stream
*/
MeshReader::NodeIterator &MeshReader::getNodeIterator()
{
return m_impl_ptr->m_node_reader_ptr->getIterator();
}
/*!
return iterator for triangles in file/stream
*/
MeshReader::TriIterator &MeshReader::getTriIterator()
{
return m_impl_ptr->m_tri_reader_ptr->getIterator();
}
} // end namespace lsm
} // end namespace esys
|
