Line data Source code
1 : #include "crpropa/module/OutputShell.h"
2 : #include "crpropa/Units.h"
3 :
4 : #include <iomanip>
5 :
6 : namespace crpropa {
7 :
8 0 : void ShellOutput::process(Candidate* c) const {
9 0 : #pragma omp critical(ShellOutput)
10 : {
11 : std::cout << std::fixed << std::showpoint << std::setprecision(3)
12 : << std::setw(6);
13 0 : std::cout << c->getTrajectoryLength() / Mpc << " Mpc, ";
14 0 : std::cout << c->getRedshift() << ", ";
15 0 : std::cout << c->current.getId() << ", ";
16 0 : std::cout << c->current.getEnergy() / EeV << " EeV, ";
17 0 : std::cout << c->current.getPosition() / Mpc << " Mpc, ";
18 0 : std::cout << c->current.getDirection();
19 : std::cout << std::endl;
20 : }
21 0 : }
22 :
23 0 : std::string ShellOutput::getDescription() const {
24 0 : return "Shell output";
25 : }
26 :
27 0 : void ShellOutput1D::process(Candidate* c) const {
28 0 : #pragma omp critical(ShellOutput)
29 : {
30 : std::cout << std::fixed << std::showpoint << std::setprecision(3)
31 : << std::setw(6);
32 0 : std::cout << c->current.getPosition().x / Mpc << " Mpc, ";
33 0 : std::cout << c->getRedshift() << ", ";
34 0 : std::cout << c->current.getId() << ", ";
35 0 : std::cout << c->current.getEnergy() / EeV << " EeV";
36 : std::cout << std::endl;
37 : }
38 0 : }
39 :
40 0 : std::string ShellOutput1D::getDescription() const {
41 0 : return "Shell output for 1D";
42 : }
43 :
44 0 : void ShellPropertyOutput::process(Candidate* c) const {
45 : Candidate::PropertyMap::const_iterator i = c->properties.begin();
46 0 : #pragma omp critical(ShellOutput)
47 : {
48 0 : for ( ; i != c->properties.end(); i++) {
49 0 : std::cout << " " << i->first << ", " << i->second << std::endl;
50 : }
51 : }
52 0 : }
53 :
54 0 : std::string ShellPropertyOutput::getDescription() const {
55 0 : return "Shell property output";
56 : }
57 :
58 : } // namespace crpropa
|