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
|
#include "../benchmark/ArgParser.h"
#include "ShapeTester.h"
#include "VecGeom/volumes/PlacedVolume.h"
#include "VecGeom/volumes/SExtru.h"
template <typename ImplT>
int runTester(ImplT const *shape, int npoints, bool debug, bool stat);
int main(int argc, char *argv[])
{
OPTION_INT(npoints, 10000);
OPTION_BOOL(debug, false);
OPTION_BOOL(stat, false);
int N = 20;
vecgeom::Precision x[N], y[N];
for (size_t i = 0; i < (size_t)N; ++i) {
x[i] = 4 * std::sin(i * (2. * M_PI) / N);
y[i] = 4 * std::cos(i * (2. * M_PI) / N);
}
auto volume = new vecgeom::SimpleSExtru("test_VecGeomSExtru", N, x, y, -5., 5.);
volume->Print();
ShapeTester<vecgeom::VPlacedVolume> tester;
tester.setDebug(debug);
tester.setStat(stat);
tester.SetMaxPoints(npoints);
int errCode = tester.Run(volume);
std::cout << "Final Error count for Shape *** " << volume->GetName() << "*** = " << errCode << "\n";
std::cout << "=========================================================\n";
if (volume) delete volume;
return 0;
}
|
