Merge pull request #609 from RevathiJambunathan/ParticlesEBParser

Externally applied E/B fields on the particles using parser

1 files changed, 21 insertions, 0 deletions

diff --git a/Source/Particles/PhysicalParticleContainer.H b/Source/Particles/PhysicalParticleContainer.H
index 0192ffb37..74d1a0f62 100644
--- a/Source/Particles/PhysicalParticleContainer.H
+++ b/Source/Particles/PhysicalParticleContainer.H
@@ -47,6 +47,26 @@ public:
                            amrex::Real t, amrex::Real dt) override;
 #endif // WARPX_DO_ELECTROSTATIC
 
+    /**
+     * \brief Apply external E and B fields on the particles. The E and B
+     * fields could be defined as a constant or using a parser for reading
+     * in a mathematical expression. The default value for the E- and B-fields
+     * is (0.0,0.0,0.0).
+     *
+     * \param[in,out] Exp-Bzp pointer to fields on particles modified based
+     *                  on external E and B
+     * \param[in] xp,yp,zp arrays of particle positions required to compute
+     *                  mathematical expression for the external fields
+     *                  using parser.
+     */
+    void AssignExternalFieldOnParticles ( WarpXParIter& pti,
+                                         RealVector& Exp, RealVector& Eyp,
+                                         RealVector& Ezp, RealVector& Bxp,
+                                         RealVector& Byp, RealVector& Bzp,
+                  amrex::Gpu::ManagedDeviceVector<amrex::ParticleReal> xp,
+                  amrex::Gpu::ManagedDeviceVector<amrex::ParticleReal> yp,
+                  amrex::Gpu::ManagedDeviceVector<amrex::ParticleReal> zp, int lev);
+
     virtual void FieldGather (int lev,
                               const amrex::MultiFab& Ex,
                               const amrex::MultiFab& Ey,
@@ -318,6 +338,7 @@ protected:
     //radiation reaction
     bool do_classical_radiation_reaction = false;
 
+
 #ifdef WARPX_QED
     // A flag to enable quantum_synchrotron process for leptons
     bool m_do_qed_quantum_sync = false;