diff options
| author |  David Grote <grote1@llnl.gov>
| 2020-07-27 09:19:47 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2020-07-27 09:19:47 -0700 |
| commit | df0ec5715c081e71f7f4bcd13a1a895fcd942436 (patch) | |
| tree | 70e203bd6a8e06a0eda68695ca0eb1d7687c5f6b /Source/WarpX.cpp | |
| parent | f68983f03ff24d6207c76a6a0ca30ed87ce1adad (diff) | |
| download | WarpX-df0ec5715c081e71f7f4bcd13a1a895fcd942436.tar.gz WarpX-df0ec5715c081e71f7f4bcd13a1a895fcd942436.tar.zst WarpX-df0ec5715c081e71f7f4bcd13a1a895fcd942436.zip | |
Added k-space filter for RZ spectral solver (#1006)
* Added k-space filter for RZ spectral solver
* Added SpectralBinomialFilter files for RZ spectral solver
* Added RZspectral binomial filter to CMakeLists.txt
* Update Docs/source/running_cpp/parameters.rst
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* Update Docs/source/running_cpp/parameters.rst
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* Update Source/FieldSolver/SpectralSolver/SpectralBinomialFilter.H
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* Update Source/FieldSolver/SpectralSolver/SpectralBinomialFilter.cpp
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* Update Source/FieldSolver/SpectralSolver/SpectralFieldDataRZ.H
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* Fixed literals in SpectralBinomialFilter.cpp
* For RZ spectral, apply filter to rho old and new
* Added SpectralBinomialFilter::InitFilterArray
* For SpectralBinomialFilter, combine R and Z into one routine
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* For SpectralBinomialFilter, combine R and Z into one routine, part 2
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
* For SpectralBinomialFilter, combine R and Z into one routine, part 3
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
Co-authored-by: Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
Diffstat (limited to 'Source/WarpX.cpp')
| -rw-r--r-- | Source/WarpX.cpp | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index e4ac79cec..967f69fdd 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -85,6 +85,9 @@ bool WarpX::use_fdtd_nci_corr = false; int WarpX::l_lower_order_in_v = true; bool WarpX::use_filter = false; +bool WarpX::use_kspace_filter = false; +bool WarpX::use_filter_compensation = false; + bool WarpX::serialize_ics = false; bool WarpX::refine_plasma = false; @@ -443,6 +446,8 @@ WarpX::ReadParameters () // Read filter and fill IntVect filter_npass_each_dir with // proper size for AMREX_SPACEDIM pp.query("use_filter", use_filter); + pp.query("use_kspace_filter", use_kspace_filter); + pp.query("use_filter_compensation", use_filter_compensation); Vector<int> parse_filter_npass_each_dir(AMREX_SPACEDIM,1); pp.queryarr("filter_npass_each_dir", parse_filter_npass_each_dir); filter_npass_each_dir[0] = parse_filter_npass_each_dir[0]; @@ -948,6 +953,9 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm realspace_ba.grow(1, ngE[1]); // add guard cells only in z spectral_solver_fp[lev].reset( new SpectralSolverRZ( realspace_ba, dm, n_rz_azimuthal_modes, noz_fft, do_nodal, dx_vect, dt[lev], lev ) ); + if (use_kspace_filter) { + spectral_solver_fp[lev]->InitFilter(filter_npass_each_dir, use_filter_compensation); + } # else if ( fft_periodic_single_box == false ) { realspace_ba.grow(ngE); // add guard cells @@ -1066,6 +1074,9 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm realspace_ba.grow(1, ngE[1]); // add guard cells only in z spectral_solver_cp[lev].reset( new SpectralSolverRZ( realspace_ba, dm, n_rz_azimuthal_modes, noz_fft, do_nodal, cdx_vect, dt[lev], lev ) ); + if (use_kspace_filter) { + spectral_solver_cp[lev]->InitFilter(filter_npass_each_dir, use_filter_compensation); + } # else realspace_ba.grow(ngE); // add guard cells spectral_solver_cp[lev].reset( new SpectralSolver( realspace_ba, dm, |