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
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
|
#include "VecGeom/volumes/UnplacedBox.h"
#include "VecGeom/management/VolumeFactory.h"
#include "VecGeom/volumes/SpecializedBox.h"
#include "VecGeom/base/RNG.h"
#include "VecGeom/base/Utils3D.h"
#include <stdio.h>
namespace vecgeom {
inline namespace VECGEOM_IMPL_NAMESPACE {
#ifndef VECCORE_CUDA
SolidMesh *UnplacedBox::CreateMesh3D(Transformation3D const &trans, const size_t nFaces) const
{
SolidMesh *sm = new SolidMesh();
sm->ResetMesh(8, 6);
Vector3D<Precision> const &box = fBox.fDimensions;
const Utils3D::Vec_t vertices[] = {{-box[0], -box[1], -box[2]}, {-box[0], box[1], -box[2]}, {box[0], box[1], -box[2]},
{box[0], -box[1], -box[2]}, {-box[0], -box[1], box[2]}, {-box[0], box[1], box[2]},
{box[0], box[1], box[2]}, {box[0], -box[1], box[2]}};
sm->SetVertices(vertices, 8);
sm->TransformVertices(trans);
sm->AddPolygon(4, {0, 1, 2, 3}, true);
sm->AddPolygon(4, {4, 7, 6, 5}, true);
sm->AddPolygon(4, {0, 4, 5, 1}, true);
sm->AddPolygon(4, {1, 5, 6, 2}, true);
sm->AddPolygon(4, {2, 6, 7, 3}, true);
sm->AddPolygon(4, {3, 7, 4, 0}, true);
return sm;
}
#endif
VECCORE_ATT_HOST_DEVICE
void UnplacedBox::Print() const
{
printf("UnplacedBox {%.2f, %.2f, %.2f}", x(), y(), z());
}
void UnplacedBox::Print(std::ostream &os) const
{
os << "UnplacedBox {" << x() << ", " << y() << ", " << z() << "}";
}
Vector3D<Precision> UnplacedBox::SamplePointOnSurface() const
{
Vector3D<Precision> p(dimensions());
Precision S[3] = {p[1] * p[2], p[0] * p[2], p[0] * p[1]};
Precision rand = (S[0] + S[1] + S[2]) * RNG::Instance().uniform(-1.0, 1.0);
int axis = 0, direction = rand < 0.0 ? -1 : 1;
rand = std::abs(rand);
// rand is guaranteed (by the contract of RNG::Instance().uniform(-1.0, 1.0);)
// to be less than one of the S[axis] since it starts
// at a random number between 0 and (S[0] + S[1] + S[2])
// and each iteration exits or substracts, respectively, S[0] and then S[1]
// so on the 3rd iteration we have axis==2 and rand <= S[2]
// Note that an automated tools (clang-tidy for example) will not be
// able to detect this guarantee and complains about a possible out of
// bound access (or garbage value).
while (rand > S[axis])
rand -= S[axis], axis++;
p[0] = (axis == 0) ? direction * x() : p[0] * RNG::Instance().uniform(-1.0, 1.0);
p[1] = (axis == 1) ? direction * y() : p[1] * RNG::Instance().uniform(-1.0, 1.0);
p[2] = (axis == 2) ? direction * z() : p[2] * RNG::Instance().uniform(-1.0, 1.0);
return p;
}
#ifndef VECCORE_CUDA
template <TranslationCode trans_code, RotationCode rot_code>
VPlacedVolume *UnplacedBox::Create(LogicalVolume const *const logical_volume,
Transformation3D const *const transformation, VPlacedVolume *const placement)
{
if (placement) {
new (placement) SpecializedBox<trans_code, rot_code>(logical_volume, transformation);
return placement;
}
return new SpecializedBox<trans_code, rot_code>(logical_volume, transformation);
}
VPlacedVolume *UnplacedBox::SpecializedVolume(LogicalVolume const *const volume,
Transformation3D const *const transformation,
const TranslationCode trans_code, const RotationCode rot_code,
VPlacedVolume *const placement) const
{
return VolumeFactory::CreateByTransformation<UnplacedBox>(volume, transformation, trans_code, rot_code, placement);
}
#else
template <TranslationCode trans_code, RotationCode rot_code>
VECCORE_ATT_DEVICE
VPlacedVolume *UnplacedBox::Create(LogicalVolume const *const logical_volume,
Transformation3D const *const transformation, const int id, const int copy_no,
const int child_id, VPlacedVolume *const placement)
{
if (placement) {
new (placement) SpecializedBox<trans_code, rot_code>(logical_volume, transformation, id, copy_no, child_id);
return placement;
}
return new SpecializedBox<trans_code, rot_code>(logical_volume, transformation, id, copy_no, child_id);
}
VECCORE_ATT_DEVICE
VPlacedVolume *UnplacedBox::SpecializedVolume(LogicalVolume const *const volume,
Transformation3D const *const transformation,
const TranslationCode trans_code, const RotationCode rot_code,
const int id, const int copy_no, const int child_id,
VPlacedVolume *const placement) const
{
return VolumeFactory::CreateByTransformation<UnplacedBox>(volume, transformation, trans_code, rot_code, id, copy_no,
child_id, placement);
}
#endif
#ifdef VECGEOM_CUDA_INTERFACE
DevicePtr<cuda::VUnplacedVolume> UnplacedBox::CopyToGpu(DevicePtr<cuda::VUnplacedVolume> const in_gpu_ptr) const
{
return CopyToGpuImpl<UnplacedBox>(in_gpu_ptr, x(), y(), z());
}
DevicePtr<cuda::VUnplacedVolume> UnplacedBox::CopyToGpu() const
{
return CopyToGpuImpl<UnplacedBox>();
}
#endif // VECGEOM_CUDA_INTERFACE
} // namespace VECGEOM_IMPL_NAMESPACE
#ifdef VECCORE_CUDA
namespace cxx {
template size_t DevicePtr<cuda::UnplacedBox>::SizeOf();
template void DevicePtr<cuda::UnplacedBox>::Construct(const Precision x, const Precision y, const Precision z) const;
} // namespace cxx
#endif
} // namespace vecgeom
|
