diff options
author | 2021-04-19 08:27:19 -0700 | |
---|---|---|
committer | 2021-04-19 08:27:19 -0700 | |
commit | 0beb477484381a012d9e80f326c4cb51926b5e5c (patch) | |
tree | be3b5e0d9374e2cb718582b4fdd64c13a1aa55ff /Source/Parallelization | |
parent | e0d1e2494f1190f1b42088273409887a459987bd (diff) | |
download | WarpX-0beb477484381a012d9e80f326c4cb51926b5e5c.tar.gz WarpX-0beb477484381a012d9e80f326c4cb51926b5e5c.tar.zst WarpX-0beb477484381a012d9e80f326c4cb51926b5e5c.zip |
Divergence cleaning for PSATD in PML (#1600)
* Add div(E)/div(B) cleaning options for PSATD in PMLs
* Pass missing flags to spectral solver in PML
* Duplicate MPI exchange and communication functions for G
* Use separate parameters for div cleaning in PMLs
* Add asserts for features that are not implemented
* Do not need to duplicate MPI exchange functions for G
* Add short documentation for new input parameters
* Set new parameters true by default with PSATD solver
* Add CI test for PML div cleaning with PSATD
* Use new syntax <diag_name>.intervals in new input file
* Reset benchmark of new CI test
* Always synchronize nodal points of G MultiFab
* Fix few warnings in 2D build
* Update Benchmark of pml_psatd_dive_divb_cleaning
* Improve Documentation of warpx.do_pml_dive_cleaning
Co-authored-by: Neïl Zaim <49716072+NeilZaim@users.noreply.github.com>
* Improve Documentation of warpx.do_pml_divb_cleaning
Co-authored-by: Neïl Zaim <49716072+NeilZaim@users.noreply.github.com>
* Improve Abort Message
* Clean Up
Co-authored-by: Neïl Zaim <49716072+NeilZaim@users.noreply.github.com>
Diffstat (limited to 'Source/Parallelization')
-rw-r--r-- | Source/Parallelization/WarpXComm.cpp | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/Source/Parallelization/WarpXComm.cpp b/Source/Parallelization/WarpXComm.cpp index 992583862..d92cffb11 100644 --- a/Source/Parallelization/WarpXComm.cpp +++ b/Source/Parallelization/WarpXComm.cpp @@ -1142,6 +1142,7 @@ void WarpX::NodalSyncPML (int lev, PatchType patch_type) const auto& pml_E = (patch_type == PatchType::fine) ? pml[lev]->GetE_fp() : pml[lev]->GetE_cp(); const auto& pml_B = (patch_type == PatchType::fine) ? pml[lev]->GetB_fp() : pml[lev]->GetB_cp(); const auto& pml_F = (patch_type == PatchType::fine) ? pml[lev]->GetF_fp() : pml[lev]->GetF_cp(); + const auto& pml_G = (patch_type == PatchType::fine) ? pml[lev]->GetG_fp() : pml[lev]->GetG_cp(); // Always synchronize nodal points const auto& period = Geom(lev).periodicity(); @@ -1154,6 +1155,9 @@ void WarpX::NodalSyncPML (int lev, PatchType patch_type) if (pml_F) { pml_F->OverrideSync(period); } + if (pml_G) { + pml_G->OverrideSync(period); + } } } |