diff options
| author |  Remi Lehe <remi.lehe@normalesup.org>
| 2021-08-20 17:11:27 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-08-20 17:11:27 -0700 |
| commit | 9f1e291bf86221e41c8495658c140d249d35da65 (patch) | |
| tree | 5b41aa596dfa25abd49cd735bb70876122cf6653 /Source/Utils/WarpXUtil.cpp | |
| parent | 45d4826633451a0f7b5003994e654f3e7950b78c (diff) | |
| download | WarpX-9f1e291bf86221e41c8495658c140d249d35da65.tar.gz WarpX-9f1e291bf86221e41c8495658c140d249d35da65.tar.zst WarpX-9f1e291bf86221e41c8495658c140d249d35da65.zip | |
Allow Silver-Mueller boundary conditions to only be applied in certain directions (#2220)
* Apply Silver-Mueller boundary only in requested direction
* Include a test with independent Silver-Mueller boundary conditions
* Correct typo in implementation
Diffstat (limited to 'Source/Utils/WarpXUtil.cpp')
| -rw-r--r-- | Source/Utils/WarpXUtil.cpp | 32 |
1 files changed, 0 insertions, 32 deletions
diff --git a/Source/Utils/WarpXUtil.cpp b/Source/Utils/WarpXUtil.cpp index 6124e8afa..6a1b5259a 100644 --- a/Source/Utils/WarpXUtil.cpp +++ b/Source/Utils/WarpXUtil.cpp @@ -655,38 +655,6 @@ void ReadBCParams () } } } - // temporarily check : If silver mueller is selected for one boundary, it should be - // selected at all valid boundaries. - for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) { - if (WarpX::field_boundary_lo[idim] == FieldBoundaryType::Absorbing_SilverMueller || - WarpX::field_boundary_hi[idim] == FieldBoundaryType::Absorbing_SilverMueller){ -#if (AMREX_SPACEDIM == 3) - AMREX_ALWAYS_ASSERT_WITH_MESSAGE( - (WarpX::field_boundary_lo[0] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[0] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_lo[2] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[2] == FieldBoundaryType::Absorbing_SilverMueller) - , " The current implementation requires silver-mueller boundary condition to be applied at all boundaries!"); -#else -#ifndef WARPX_DIM_RZ - AMREX_ALWAYS_ASSERT_WITH_MESSAGE( - (WarpX::field_boundary_lo[0] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[0] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller) - , " The current implementation requires silver-mueller boundary condition to be applied at all boundaries!"); -#else - AMREX_ALWAYS_ASSERT_WITH_MESSAGE( - (WarpX::field_boundary_hi[0] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_lo[1] == FieldBoundaryType::Absorbing_SilverMueller)&& - (WarpX::field_boundary_hi[1] == FieldBoundaryType::Absorbing_SilverMueller) - , " The current implementation requires silver-mueller boundary condition to be applied at all boundaries!"); -#endif -#endif - } - } #ifdef WARPX_DIM_RZ // Ensure code aborts if PEC is specified at r=0 for RZ AMREX_ALWAYS_ASSERT_WITH_MESSAGE( WarpX::field_boundary_lo[0] == FieldBoundaryType::None, |