From 9035ee165054e25aedf98d97f16786d8d5f2965a Mon Sep 17 00:00:00 2001
From: Remi Lehe <remi.lehe@normalesup.org>
Date: Fri, 10 Jan 2020 10:52:31 -0800
Subject: Started implementing finite difference solver

Added Yee algorithm
---
 Source/FieldSolver/WarpXPushFieldsEM.cpp | 7 +++++++
 1 file changed, 7 insertions(+)

(limited to 'Source/FieldSolver/WarpXPushFieldsEM.cpp')

diff --git a/Source/FieldSolver/WarpXPushFieldsEM.cpp b/Source/FieldSolver/WarpXPushFieldsEM.cpp
index 4848b051e..1137da741 100644
--- a/Source/FieldSolver/WarpXPushFieldsEM.cpp
+++ b/Source/FieldSolver/WarpXPushFieldsEM.cpp
@@ -111,6 +111,13 @@ WarpX::EvolveB (int lev, amrex::Real a_dt)
 void
 WarpX::EvolveB (int lev, PatchType patch_type, amrex::Real a_dt)
 {
+    if (patch_type == PatchType::fine) {
+        fdtd_solver_fp->EvolveB( Bfield_fp[lev], Efield_fp[lev], a_dt );
+    } else {
+        fdtd_solver_cp->EvolveB( Bfield_cp[lev], Efield_cp[lev], a_dt );
+    }
+
+    // Goes into initializer
     const int patch_level = (patch_type == PatchType::fine) ? lev : lev-1;
     const std::array<Real,3>& dx = WarpX::CellSize(patch_level);
     const Real dtsdx = a_dt/dx[0], dtsdy = a_dt/dx[1], dtsdz = a_dt/dx[2];
-- 
cgit v1.2.3