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
128
129
130
131
132
133
134
|
/////////////////////////////////////////////////////////////
// //
// Copyright (c) 2003-2011 by The University of Queensland //
// Earth Systems Science Computational Centre (ESSCC) //
// http://www.uq.edu.au/esscc //
// //
// Primary Business: Brisbane, Queensland, Australia //
// Licensed under the Open Software License version 3.0 //
// http://www.opensource.org/licenses/osl-3.0.php //
// //
/////////////////////////////////////////////////////////////
#include "Geometry/SimpleBlock.h"
#include <stdexcept>
#include <fstream>
#include <sstream>
#include <iomanip>
#include <boost/limits.hpp>
namespace esys
{
namespace lsm
{
template <typename TmplParticle>
SimpleBlockGenerator<TmplParticle>::SimpleBlockGenerator(
unsigned int numX,
unsigned int numY,
unsigned int numZ,
double radius
)
: m_radius(radius),
m_dimCounts(numX, numY, numZ)
{
}
template <typename TmplParticle>
Vec3 SimpleBlockGenerator<TmplParticle>::getPos(const Vec3L &idx)
{
return
Vec3(
idx[0]*2.0*getRadius() + getRadius(),
idx[1]*2.0*getRadius() + getRadius(),
idx[2]*2.0*getRadius() + ((m_dimCounts[2] > 1) ? getRadius() : 0.0)
);
}
template <typename TmplParticle>
int SimpleBlockGenerator<TmplParticle>::getId(const Vec3L &idx)
{
return
idx[0]
+
idx[1]*m_dimCounts[0]
+
idx[2]*m_dimCounts[0]*m_dimCounts[1];
}
template <typename TmplParticle>
template <typename TmplParticleCollection>
void
SimpleBlockGenerator<TmplParticle>::createParticles(
TmplParticleCollection &particleCollection
)
{
Vec3L idx(0, 0, 0);
for (idx[2]=0; idx[2] < m_dimCounts[2]; (idx[2])++)
{
for (idx[1]=0; idx[1] < m_dimCounts[1]; (idx[1])++)
{
for (idx[0]=0; idx[0] < m_dimCounts[0]; (idx[0])++)
{
particleCollection.createParticle(
TmplParticle(
getPos(idx),
getRadius(),
getId(idx),
0
)
);
}
}
}
}
template <typename TmplParticle>
SimpleBlockGenerator<TmplParticle>::~SimpleBlockGenerator()
{
}
template <typename TmplParticle>
double SimpleBlockGenerator<TmplParticle>::getRadius() const
{
return m_radius;
}
template <typename TmplParticle>
SimpleBlock<TmplParticle>::SimpleBlock(
unsigned int numX,
unsigned int numY,
unsigned int numZ,
double radius
)
: ParticleCollection<TmplParticle>(),
m_generator(numX, numY, numZ, radius)
{
createParticles();
}
template <typename TmplParticle>
SimpleBlock<TmplParticle>::~SimpleBlock()
{
}
template <typename TmplParticle>
void SimpleBlock<TmplParticle>::createParticles()
{
m_generator.createParticles(*this);
}
template <typename TmplParticle>
double SimpleBlock<TmplParticle>::getRadius() const
{
return m_generator.getRadius();
}
}
}
|