Boundary Condition : PEC (#1767)

* Read boundary and set periodicity, enumerate BC types, added support for periodic

* separate particle and field boudnary structs

* PEC Etangential Bnormal and setting default pml=0

* eol

* add PEC and FieldBoundary routines

* eol

* remove duplicate struct

* update guard cells for PEC

* host device functions

* remove duplicate function

* fix compilation error by adding ;

* temporary logic to set PEC boundary

* eol

* add documentation for PEC and abort if PEC is used for RZ

* documentation for PEC functions

* remove print statements

* Adding field CI test for PEC for 3D with and without MR

* fix eol

* Apply suggestions from code review

* Update Docs/source/usage/parameters.rst

* Update Source/BoundaryConditions/WarpX_PEC.cpp

* particle PEC CI test

* rename analysis scripts for PEC

* use iside == 0 as condition

* ParallelFor over guard cells. Set field values on boundary and guard in the kernel

* ijk_guard not needed anymore

* remove unnecessary nbody generated by yt

* pi and clight are built-in

* cleanup

* eol

* fix bug in computing ig and ijk_valid

* bug fix in guard-cell update in PEC changes benchmark for particles

* fix bug in the ig compute

* reset particle PEC benchmark

* fix logic to set PEC with current interface

* adding a none type for RZ

* reset benchmarks for 2dLaserInjection, 2dgalilean_hybrid and 2dcomoving CI tests due to differences in PEC reflection

* ensure noz shape factor is used for RZ and k=0 for RZ

* fix XZ/RZ logic for Ey,Etheta,By,Btheta. Add ncomponent for RZ, generically.

* fix the field update for corner cases and mixed boundaries. rename ijk_valid to ijk_mirror since ijk_mirror need not be in valid region

* set GuardCell and Boundary cell flags

* fix compilation bug

* suppress unused variable warning

* eol whitespace

* fix doc

* Apply suggestions from code review

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

* RZ abort message

* reset benchmark after fixing guard-cell updates

* fix typo

* RZ doc

* Update Source/Utils/WarpXUtil.cpp

* Update Examples/Tests/PEC/inputs_particle_PEC_3d

* changes suggested from code-review

* set Enormal and Btangential in guard cell across PEC boundary

* add more info about guard cell update in the docs

* fix typo

* ijk mirror only if ig>0

* update benchmarks

* ignore unused k

* temporarily revert guardcell damping for Cartesian

* fix typo

* dont set Etangential and Bnormal to zero if field is not nodal on boundary

* set bindary int as boolean

* fix high to hi

* set default nox to 1

* particle shape for PEC test

* update benchmarks for 2D comoving and galilean

* reset benchmark for RZ test-cases

* Add in benchmarks generated by @EZoni

* Changing transverse boundaries to periodic for psatd galilean cases

* resetting benchmark for comoving_hybrid_2d and galilean_hybrid_2d

* at-least one guard cell must be filled for PEC

* abort if PEC is used for PSATD

* use new boundary interface to set none for rmin and rmax, since pec does not work for psatd

* Update Source/WarpX.cpp

* add doc mentioning PEC does not work for PSATD

* missing semicolon

* Update Source/BoundaryConditions/WarpX_PEC.H

Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>

* Update Source/BoundaryConditions/WarpX_PEC.H

* lengthy comment for normal and tangential components at domain boundaries

* eol fix

* grammar

* field_hi

* reset benchmarks

Co-authored-by: Andy Nonaka <AJNonaka@lbl.gov>
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>

1 files changed, 1 insertions, 0 deletions

diff --git a/Source/Utils/WarpXAlgorithmSelection.cpp b/Source/Utils/WarpXAlgorithmSelection.cpp
index a0b90efec..cdf896eb2 100644
--- a/Source/Utils/WarpXAlgorithmSelection.cpp
+++ b/Source/Utils/WarpXAlgorithmSelection.cpp
@@ -81,6 +81,7 @@ const std::map<std::string, int> FieldBCType_algo_to_int = {
     {"pec",      FieldBoundaryType::PEC},
     {"pmc",      FieldBoundaryType::PMC},
     {"damped",   FieldBoundaryType::Damped},
+    {"none",     FieldBoundaryType::None},
     {"default",  FieldBoundaryType::PML}
 };