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
|
/* Copyright 2019 Andrew Myers, Axel Huebl, David Grote
* Luca Fedeli, Maxence Thevenet, Revathi Jambunathan
* Weiqun Zhang, levinem, Yinjian Zhao
*
* This file is part of WarpX.
*
* License: BSD-3-Clause-LBNL
*/
#include "Particles/ParticleIO.H"
#include "Particles/MultiParticleContainer.H"
#include "Particles/ParticleBuffer.H"
#include "Particles/PhysicalParticleContainer.H"
#include "Particles/RigidInjectedParticleContainer.H"
#include "Particles/SpeciesPhysicalProperties.H"
#include "Particles/WarpXParticleContainer.H"
#include "Utils/WarpXConst.H"
#include "Utils/WarpXProfilerWrapper.H"
#include "WarpX.H"
#include <AMReX_BLassert.H>
#include <AMReX_Config.H>
#include <AMReX_Extension.H>
#include <AMReX_GpuControl.H>
#include <AMReX_GpuLaunch.H>
#include <AMReX_GpuQualifiers.H>
#include <AMReX_PODVector.H>
#include <AMReX_ParIter.H>
#include <AMReX_ParticleIO.H>
#include <AMReX_REAL.H>
#include <AMReX_Vector.H>
#include <array>
#include <istream>
#include <memory>
#include <string>
#include <vector>
using namespace amrex;
void
RigidInjectedParticleContainer::ReadHeader (std::istream& is)
{
is >> charge >> mass;
WarpX::GotoNextLine(is);
int nlevs;
is >> nlevs;
WarpX::GotoNextLine(is);
AMREX_ASSERT(zinject_plane_levels.size() == 0);
for (int i = 0; i < nlevs; ++i)
{
int zinject_plane_tmp;
is >> zinject_plane_tmp;
zinject_plane_levels.push_back(zinject_plane_tmp);
WarpX::GotoNextLine(is);
}
}
void
RigidInjectedParticleContainer::WriteHeader (std::ostream& os) const
{
// no need to write species_id
os << charge << " " << mass << "\n";
int nlevs = zinject_plane_levels.size();
os << nlevs << "\n";
for (int i = 0; i < nlevs; ++i)
{
os << zinject_plane_levels[i] << "\n";
}
}
void
WarpXParticleContainer::ReadHeader (std::istream& is)
{
is >> charge >> mass;
WarpX::GotoNextLine(is);
}
void
WarpXParticleContainer::WriteHeader (std::ostream& os) const
{
// no need to write species_id
os << charge << " " << mass << "\n";
}
void
MultiParticleContainer::Restart (const std::string& dir)
{
// note: all containers is sorted like this
// - species_names
// - laser_names
// we don't need to read back the laser particle charge/mass
for (unsigned i = 0, n = species_names.size(); i < n; ++i) {
allcontainers.at(i)->Restart(dir, species_names.at(i));
}
}
void
MultiParticleContainer::ReadHeader (std::istream& is)
{
// note: all containers is sorted like this
// - species_names
// - laser_names
// we don't need to read back the laser particle charge/mass
for (unsigned i = 0, n = species_names.size(); i < n; ++i) {
allcontainers.at(i)->ReadHeader(is);
}
}
void
MultiParticleContainer::WriteHeader (std::ostream& os) const
{
// note: all containers is sorted like this
// - species_names
// - laser_names
// we don't need to read back the laser particle charge/mass
for (unsigned i = 0, n = species_names.size(); i < n; ++i) {
allcontainers.at(i)->WriteHeader(os);
}
}
void
PhysicalParticleContainer::ConvertUnits (ConvertDirection convert_direction)
{
WARPX_PROFILE("PhysicalParticleContainer::ConvertUnits()");
// Account for the special case of photons
const auto t_mass =
this->AmIA<PhysicalSpecies::photon>() ? PhysConst::m_e : this->getMass();
particlesConvertUnits(convert_direction, this, t_mass);
}
|
