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
|
/* Copyright 2019-2020 Andrew Myers, Axel Huebl,
* Maxence Thevenet
*
* This file is part of WarpX.
*
* License: BSD-3-Clause-LBNL
*/
#include "Particles/ElementaryProcess/QEDPairGeneration.H"
#include "WarpX.H"
#include <AMReX_Box.H>
#include <AMReX_FArrayBox.H>
#include <AMReX_IntVect.H>
#include <algorithm>
#include <array>
PairGenerationTransformFunc::
PairGenerationTransformFunc (BreitWheelerGeneratePairs const generate_functor,
const WarpXParIter& a_pti, int lev, amrex::IntVect ngEB,
amrex::FArrayBox const& exfab,
amrex::FArrayBox const& eyfab,
amrex::FArrayBox const& ezfab,
amrex::FArrayBox const& bxfab,
amrex::FArrayBox const& byfab,
amrex::FArrayBox const& bzfab,
amrex::Vector<amrex::Real> v_galilean,
int a_offset)
: m_generate_functor(generate_functor)
{
m_get_position = GetParticlePosition(a_pti, a_offset);
m_get_externalEB = GetExternalEBField(a_pti, a_offset);
m_ex_arr = exfab.array();
m_ey_arr = eyfab.array();
m_ez_arr = ezfab.array();
m_bx_arr = bxfab.array();
m_by_arr = byfab.array();
m_bz_arr = bzfab.array();
m_ex_type = exfab.box().ixType();
m_ey_type = eyfab.box().ixType();
m_ez_type = ezfab.box().ixType();
m_bx_type = bxfab.box().ixType();
m_by_type = byfab.box().ixType();
m_bz_type = bzfab.box().ixType();
amrex::Box box = a_pti.tilebox();
box.grow(ngEB);
const std::array<amrex::Real,3>& dx = WarpX::CellSize(std::max(lev, 0));
m_dx_arr = {dx[0], dx[1], dx[2]};
// Lower corner of tile box physical domain (take into account Galilean shift)
amrex::Real cur_time = WarpX::GetInstance().gett_new(lev);
const auto& time_of_last_gal_shift = WarpX::GetInstance().time_of_last_gal_shift;
amrex::Real time_shift = (cur_time - time_of_last_gal_shift);
amrex::Array<amrex::Real,3> galilean_shift = { v_galilean[0]*time_shift, v_galilean[1]*time_shift, v_galilean[2]*time_shift };
const std::array<amrex::Real, 3>& xyzmin = WarpX::LowerCorner(box, galilean_shift, lev);
m_xyzmin_arr = {xyzmin[0], xyzmin[1], xyzmin[2]};
m_galerkin_interpolation = WarpX::galerkin_interpolation;
m_nox = WarpX::nox;
m_n_rz_azimuthal_modes = WarpX::n_rz_azimuthal_modes;
m_lo = amrex::lbound(box);
}
|
