From 8e2ecfe27b9d464dcd40031ae5ebc4e3781963e2 Mon Sep 17 00:00:00 2001
From: Yin-YinjianZhao <yinjianzhao@lbl.gov>
Date: Mon, 14 Oct 2019 10:28:20 -0700
Subject: Add Higuera-Cary pusher

---
 Source/Particles/RigidInjectedParticleContainer.cpp | 8 ++++++++
 1 file changed, 8 insertions(+)

(limited to 'Source/Particles/RigidInjectedParticleContainer.cpp')

diff --git a/Source/Particles/RigidInjectedParticleContainer.cpp b/Source/Particles/RigidInjectedParticleContainer.cpp
index 891ade76d..aeaac7194 100644
--- a/Source/Particles/RigidInjectedParticleContainer.cpp
+++ b/Source/Particles/RigidInjectedParticleContainer.cpp
@@ -13,6 +13,7 @@
 #include <WarpXAlgorithmSelection.H>
 #include <UpdateMomentumBoris.H>
 #include <UpdateMomentumVay.H>
+#include <UpdateMomentumHigueraCary.H>
 
 using namespace amrex;
 
@@ -443,6 +444,13 @@ RigidInjectedParticleContainer::PushP (int lev, Real dt,
                               Expp[i], Eypp[i], Ezpp[i], Bxpp[i], Bypp[i], Bzpp[i], q, m, dt);
                     }
                 );
+            } else if (WarpX::particle_pusher_algo == ParticlePusherAlgo::HigueraCary) {
+                amrex::ParallelFor( pti.numParticles(),
+                    [=] AMREX_GPU_DEVICE (long i) {
+                        UpdateMomentumHigueraCary( uxpp[i], uypp[i], uzpp[i],
+                              Expp[i], Eypp[i], Ezpp[i], Bxpp[i], Bypp[i], Bzpp[i], q, m, dt);
+                    }
+                );
             } else {
               amrex::Abort("Unknown particle pusher");
             };
-- 
cgit v1.2.3