Staircased embedded boundaries in the YEE solver (#1881)

* Added staircased embedded boundaris to the YEE solver

* adding spherical resonating cavity test

* adding functions for fields initialization

* style adjustments

* fixing tabs

* fixed name of analysis script

* fixed name of analysis script

* fixed a few wrong preprocessor directives

* workaround for missing boost

* Revert "workaround for missing boost"

This reverts commit 601f9eb2ec6f8c2100304379b2bea1c6cf9d1851.

* another workaround for missing boost

* getting rid of boost by depending on c++17

* Removed a few unused variables

* adding USE_EB to addToCompileString for EB testing

* removed tabs

* fixing the inputs name for EB sphere test

* shortened the test

* zero padding the names of the images

* adjusted two for loops

* removed some unused variables

* improving the fields initialization

* removed the sphere test and implemented the cube test

* fixed edges lengths computation and added comments

* Fixed the case of all_regular geometries

* fixing a bug that was breaking some tests

* adding test folder

* fixed the default values for the EB cube test

* simplified the analysis script

* fixed cubic resonator default results

* inputting the plot file name from command line

* fixing the diag name

* Fixed a bug in edges initialization

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

* Adding comments to the staircased yee solver (thanks Remi)

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

* fixed the cube resonator test

* removed an unused import

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

1 files changed, 4 insertions, 0 deletions

diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
index c5041a53f..f434ad3d0 100644
--- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
+++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H
@@ -38,11 +38,13 @@ class FiniteDifferenceSolver
 
         void EvolveB ( std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Bfield,
                        std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Efield,
+                       std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& face_areas,
                        int lev, amrex::Real const dt );
 
         void EvolveE ( std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Efield,
                        std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Bfield,
                        std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Jfield,
+                       std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& edge_lengths,
                        std::unique_ptr<amrex::MultiFab> const& Ffield,
                        int lev, amrex::Real const dt );
 
@@ -150,6 +152,7 @@ class FiniteDifferenceSolver
         void EvolveBCartesian (
             std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Bfield,
             std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Efield,
+            std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& face_areas,
             int lev, amrex::Real const dt );
 
         template< typename T_Algo >
@@ -157,6 +160,7 @@ class FiniteDifferenceSolver
             std::array< std::unique_ptr<amrex::MultiFab>, 3 >& Efield,
             std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Bfield,
             std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& Jfield,
+            std::array< std::unique_ptr<amrex::MultiFab>, 3 > const& edge_lengths,
             std::unique_ptr<amrex::MultiFab> const& Ffield,
             int lev, amrex::Real const dt );