Source/ParticleIO.cpp


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


#include <ParticleContainer.H>
#include <WarpX.H>

using namespace amrex;

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::Checkpoint (const std::string& dir, 
				    const std::string& name,
				    bool is_checkpoint,
                                    const Vector<std::string>& varnames) const
{
    for (unsigned i = 0, n = allcontainers.size(); i < n; ++i) {
	std::string namei = name + std::to_string(i);
	allcontainers[i]->Checkpoint(dir, namei, is_checkpoint, varnames);
    }
}

void
MultiParticleContainer::Restart (const std::string& dir, const std::string& name)
{
    for (unsigned i = 0, n = allcontainers.size(); i < n; ++i) {
	std::string namei = name + std::to_string(i);
	allcontainers[i]->Restart(dir, namei);
    }
}

void
MultiParticleContainer::ReadHeader (std::istream& is) 
{
    for (auto& pc : allcontainers) {
	pc->ReadHeader(is);
    }
}

void
MultiParticleContainer::WriteHeader (std::ostream& os) const
{
    for (const auto& pc : allcontainers) {
	pc->WriteHeader(os);
    }
}