From 091f7efd85e76b00510ea9dd909cbbebe47095e6 Mon Sep 17 00:00:00 2001 From: David Grote Date: Thu, 7 Jan 2021 08:33:23 -0800 Subject: Implemented update without rho in RZ spectral solver (#1569) * Implemented update without rho in RZ spectral solver * Updated documentation for update_with_rho for RZ * Tiny fix in GalileanPsatdAlgorithmRZ reordering declarations * In Langmuir_multi_rz_psatd, set update_with_rho = 1 --- .../SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H') diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H index 3c8d3d60e..18a7f2dcd 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/GalileanPsatdAlgorithmRZ.H @@ -21,7 +21,8 @@ class GalileanPsatdAlgorithmRZ : public SpectralBaseAlgorithmRZ int const n_rz_azimuthal_modes, int const norder_z, bool const nodal, const amrex::Array& v_galilean, - amrex::Real const dt_step); + amrex::Real const dt_step, + bool const update_with_rho); // Redefine functions from base class virtual void pushSpectralFields (SpectralFieldDataRZ & f) override final; virtual int getRequiredNumberOfFields () const override final { @@ -63,9 +64,10 @@ class GalileanPsatdAlgorithmRZ : public SpectralBaseAlgorithmRZ // Note that dt and v_galilean are saved to use in InitializeSpectralCoefficients amrex::Real const m_dt; amrex::Array m_v_galilean; + bool m_update_with_rho; SpectralRealCoefficients C_coef, S_ck_coef; - SpectralComplexCoefficients Theta2_coef, X1_coef, X2_coef, X3_coef, X4_coef; + SpectralComplexCoefficients Theta2_coef, T_rho_coef, X1_coef, X2_coef, X3_coef, X4_coef; }; -- cgit v1.2.3