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
|
#ifndef WARPX_FielIO_H_
#define WARPX_FielIO_H_
#include <WarpX.H>
#ifdef WARPX_USE_OPENPMD
#include <openPMD/openPMD.hpp>
#endif
using namespace amrex;
void
PackPlotDataPtrs (Vector<const MultiFab*>& pmf,
const std::array<std::unique_ptr<MultiFab>,3>& data);
void
AverageAndPackVectorField( MultiFab& mf_avg,
const std::array< std::unique_ptr<MultiFab>, 3 >& vector_field,
const int dcomp, const int ngrow );
void
AverageAndPackScalarField( MultiFab& mf_avg,
const MultiFab & scalar_field,
const int dcomp, const int ngrow );
void
WriteRawField( const MultiFab& F, const DistributionMapping& dm,
const std::string& filename,
const std::string& level_prefix,
const std::string& field_name,
const int lev, const bool plot_guards );
void
WriteZeroRawField( const MultiFab& F, const DistributionMapping& dm,
const std::string& filename,
const std::string& level_prefix,
const std::string& field_name,
const int lev, const int ng );
void
WriteCoarseScalar( const std::string field_name,
const std::unique_ptr<MultiFab>& F_cp,
const std::unique_ptr<MultiFab>& F_fp,
const DistributionMapping& dm,
const std::string& filename,
const std::string& level_prefix,
const int lev, const bool plot_guards,
const int icomp=0 );
void
WriteCoarseVector( const std::string field_name,
const std::unique_ptr<MultiFab>& Fx_cp,
const std::unique_ptr<MultiFab>& Fy_cp,
const std::unique_ptr<MultiFab>& Fz_cp,
const std::unique_ptr<MultiFab>& Fx_fp,
const std::unique_ptr<MultiFab>& Fy_fp,
const std::unique_ptr<MultiFab>& Fz_fp,
const DistributionMapping& dm,
const std::string& filename,
const std::string& level_prefix,
const int lev, const bool plot_guards );
std::unique_ptr<MultiFab>
getInterpolatedScalar(
const MultiFab& F_cp, const MultiFab& F_fp,
const DistributionMapping& dm, const int r_ratio,
const Real* dx, const int ngrow );
std::array<std::unique_ptr<MultiFab>, 3>
getInterpolatedVector(
const std::unique_ptr<MultiFab>& Fx_cp,
const std::unique_ptr<MultiFab>& Fy_cp,
const std::unique_ptr<MultiFab>& Fz_cp,
const std::unique_ptr<MultiFab>& Fx_fp,
const std::unique_ptr<MultiFab>& Fy_fp,
const std::unique_ptr<MultiFab>& Fz_fp,
const DistributionMapping& dm,
const int r_ratio, const Real* dx,
const int ngrow );
void
coarsenCellCenteredFields(
Vector<MultiFab>& coarse_mf, Vector<Geometry>& coarse_geom,
const Vector<MultiFab>& source_mf, const Vector<Geometry>& source_geom,
int coarse_ratio, int finest_level );
#ifdef WARPX_USE_OPENPMD
void
setOpenPMDUnit( openPMD::Mesh mesh, const std::string field_name );
std::vector<std::uint64_t>
getReversedVec( const IntVect& v );
std::vector<double>
getReversedVec( const Real* v );
void
WriteOpenPMDFields( const std::string& filename,
const std::vector<std::string>& varnames,
const MultiFab& mf, const Geometry& geom,
const int iteration, const double time );
#endif // WARPX_USE_OPENPMD
#endif // WARPX_FielIO_H_
|
