Variation of macroscopic properties for E-update (#1016)

* Adding macro-E Push and new file

* Add macroEvolveE, call it, and include algo selection in utils

* fix eol

* Fixing bug in macroE for sigma method 1

* changing MacroEvolveE to MacroscopicEvolveE

* add class for macroscopicproperties and an object in WarpX class

* fix eol

* adding templated ckc call with comment that EvolveE is same for yee and ckc

* add header file pointing to ckc algorithm

* adding obejct m_macroscopic_properties to access sigma,eps,mu

* some cleaning

* Adding comments

* adding documentation

* spelling wandroff to wendroff

* fixing eol

* eol

* const in the right place. Thanks bot!

* profiler for macroscopic evolveE

* MultiFab macroproperties with constant init, templated macroEvolveE,

* call macroparameter init

* eol fix

* add parser for macroscopic properties

* fix eol

* adding input file

* __device__ lambda cannot be private in class

* fix grown tilebox declaration for init data

* [skip ci] some more merge conflicts

* some comments

* removing redundant calls to Macroscopic EvolveE

* [skip ci] fix growntilebox for initializing macro mf

* clean and fix BackwardEuler call for ckc

* commenting out old alpha and beta implementations

* temporarily commiting local changes with calls in Interp.

* fixing a typo

* clean and add documentation

* remove the test input file

* fix typo

* eol fix

* Update Docs/source/running_cpp/parameters.rst

Co-authored-by: Andy Nonaka <AJNonaka@lbl.gov>

* PR suggestions

* Update Source/FieldSolver/WarpXPushFieldsEM.cpp

Co-authored-by: Andy Nonaka <AJNonaka@lbl.gov>

* removing unnecessary includes

* adding 2D initialization for stag arrays

* removing init_style input parameter for material properties

* eol fix

* Adding dE/dt eq curl of (B/mu)

* PhysConst ep0 and mu0

* Add functor for field access for macro and vacuum (B/mu)

* fix eol

* Update Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/FieldAccessorFunctors.H

Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>

* Apply suggestions from code review

Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>

* Update Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CartesianCKCAlgorithm.H

Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>

* Apply suggestions from code review

Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>

* fixing compilation errors and removing Gpu ManagedVector

Co-authored-by: Andy Nonaka <AJNonaka@lbl.gov>
Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>

1 files changed, 7 insertions, 3 deletions

diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CartesianYeeAlgorithm.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CartesianYeeAlgorithm.H
index 305fb3507..2ec6a3d22 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CartesianYeeAlgorithm.H
+++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CartesianYeeAlgorithm.H
@@ -9,6 +9,7 @@
 #define WARPX_FINITE_DIFFERENCE_ALGORITHM_CARTESIAN_YEE_H_
 
 #include "Utils/WarpXConst.H"
+#include "FieldAccessorFunctors.H"
 
 #include <AMReX.H>
 #include <AMReX_REAL.H>
@@ -68,9 +69,10 @@ struct CartesianYeeAlgorithm {
 
     /**
      * Perform derivative along x on a nodal grid, from a cell-centered field `F`*/
+    template< typename T_Field>
     AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
     static amrex::Real DownwardDx (
-        amrex::Array4<amrex::Real> const& F,
+        T_Field const& F,
         amrex::Real const * const coefs_x, int const /*n_coefs_x*/,
         int const i, int const j, int const k, int const ncomp=0 ) {
 
@@ -103,9 +105,10 @@ struct CartesianYeeAlgorithm {
 
     /**
      * Perform derivative along y on a nodal grid, from a cell-centered field `F`*/
+    template< typename T_Field>
     AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
     static amrex::Real DownwardDy (
-        amrex::Array4<amrex::Real> const& F,
+        T_Field const& F,
         amrex::Real const * const coefs_y, int const n_coefs_y,
         int const i, int const j, int const k, int const ncomp=0 ) {
 
@@ -143,9 +146,10 @@ struct CartesianYeeAlgorithm {
 
     /**
      * Perform derivative along z on a nodal grid, from a cell-centered field `F`*/
+    template< typename T_Field>
     AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
     static amrex::Real DownwardDz (
-        amrex::Array4<amrex::Real> const& F,
+        T_Field const& F,
         amrex::Real const * const coefs_z, int const /*n_coefs_z*/,
         int const i, int const j, int const k, int const ncomp=0 ) {