Implement nodal solver

1 files changed, 4 insertions, 1 deletions

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 );