From 0564feb5041728173716ef251f1d66d37c314770 Mon Sep 17 00:00:00 2001 From: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> Date: Thu, 6 Aug 2020 06:30:14 -0700 Subject: Galilean PSATD: current correction and rho-free formulation (#1151) * Introduce option to update E with/without rho * Clean up * Implement current correction for Galilean PSATD (needs bug fix) * Include equations in docs * Fix EOL whitespaces error * Small clean-up * Implement Galilean PSATD update without rho * Clean up * Fix bug in current correction * Fix EOL whitespaces * Clean up * Fix unused import * Remove unused variable * [skip CI] Improve docs * Clean up style * Fix EOL whitespaces * Fix EOL whitespaces * Clean up style * Revert analysis script to old status * [skip CI] Clean up style * Make equations more human-readable and improve comments * 2D test with current correction works * Temporary build fix as in #1197 * 3D test with current correction works * Rename th and th_star as theta and theta_star * Fix a couple of wrong comments * Add vertical spaces to improve readability * Improve documentation * Function CurrentCorrection is now pure * 2D benchmark fields data are now correct * Add limits of coefficients for nu=-1 * Change default of update_with_rho for Galilean PSATD --- .../SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp | 8 ++++++++ 1 file changed, 8 insertions(+) (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp') diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp index d2f087706..4bc147cd8 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PMLPsatdAlgorithm.cpp @@ -155,6 +155,14 @@ void PMLPsatdAlgorithm::InitializeSpectralCoefficients ( } }; +void +PMLPsatdAlgorithm::CurrentCorrection (SpectralFieldData& field_data, + std::array,3>& current, + const std::unique_ptr& rho) +{ + amrex::Abort("Current correction not implemented for PML PSATD"); +} + void PMLPsatdAlgorithm::VayDeposition (SpectralFieldData& field_data, std::array,3>& current) -- cgit v1.2.3