From 0e14e871a29027ded4258b079abd63ada94bc874 Mon Sep 17 00:00:00 2001
From: Remi Lehe <remi.lehe@normalesup.org>
Date: Mon, 27 Jan 2020 16:08:21 -0800
Subject: Implement nodal solver

---
 Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

(limited to 'Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp')

diff --git a/Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp b/Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp
index 72d3a1135..7a3f5aa9b 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp
+++ b/Source/FieldSolver/FiniteDifferenceSolver/EvolveB.cpp
@@ -5,6 +5,7 @@
 #else
     #include "FiniteDifferenceAlgorithms/YeeAlgorithm.H"
     #include "FiniteDifferenceAlgorithms/CKCAlgorithm.H"
+    #include "FiniteDifferenceAlgorithms/NodalAlgorithm.H"
 #endif
 #include <AMReX_Gpu.H>
 
@@ -20,7 +21,9 @@ void FiniteDifferenceSolver::EvolveB ( VectorField& Bfield,
     if (m_fdtd_algo == MaxwellSolverAlgo::Yee){
         EvolveBCylindrical <CylindricalYeeAlgorithm> ( Bfield, Efield, dt );
 #else
-    if (m_fdtd_algo == MaxwellSolverAlgo::Yee){
+    if (m_do_nodal) {
+        EvolveBCartesian <NodalAlgorithm> ( Bfield, Efield, dt );
+    } else if (m_fdtd_algo == MaxwellSolverAlgo::Yee) {
         EvolveBCartesian <YeeAlgorithm> ( Bfield, Efield, dt );
     } else if (m_fdtd_algo == MaxwellSolverAlgo::CKC) {
         EvolveBCartesian <CKCAlgorithm> ( Bfield, Efield, dt );
-- 
cgit v1.2.3