diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2021-10-21 10:01:39 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-10-21 10:01:39 -0700 |
| commit | 7d011e8f733c9bcb1de8d44c27bf912655d36068 (patch) | |
| tree | ac70d3e05ef087f36c3fd56918bd182cd14c6295 /Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp | |
| parent | 5500567e17314557c44c84bf38b135aeb75d9a99 (diff) | |
| download | WarpX-7d011e8f733c9bcb1de8d44c27bf912655d36068.tar.gz WarpX-7d011e8f733c9bcb1de8d44c27bf912655d36068.tar.zst WarpX-7d011e8f733c9bcb1de8d44c27bf912655d36068.zip | |
Div(B) Cleaning: Change Normalization of G (#2429)
* PSATD with div(B) Cleaning: Change Normalization of G
* Use New Normalization also in RZ
* Use New Normalization also in FDTD
* Reset Benchmark of Test divb_cleaning_3d
Diffstat (limited to 'Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp')
| -rw-r--r-- | Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp index c497ab618..60a30e6f6 100644 --- a/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp +++ b/Source/FieldSolver/SpectralSolver/SpectralAlgorithms/PsatdAlgorithmRZ.cpp @@ -268,11 +268,11 @@ PsatdAlgorithmRZ::pushSpectralFields(SpectralFieldDataRZ & f) { const Complex k_dot_B = -I * (kr * (Bp_old - Bm_old) + I * kz * Bz_old); - fields(i,j,k,Bp_m) += -c2 * kr * 0.5_rt * S_ck * G_old; - fields(i,j,k,Bm_m) += c2 * kr * 0.5_rt * S_ck * G_old; - fields(i,j,k,Bz_m) += I * c2 * kz * S_ck * G_old; + fields(i,j,k,Bp_m) += -kr * 0.5_rt * S_ck * G_old; + fields(i,j,k,Bm_m) += kr * 0.5_rt * S_ck * G_old; + fields(i,j,k,Bz_m) += I * kz * S_ck * G_old; - fields(i,j,k,G_m) = C * G_old + I * S_ck * k_dot_B; + fields(i,j,k,G_m) = C * G_old + I * c2 * S_ck * k_dot_B; } if (time_averaging) @@ -315,9 +315,9 @@ PsatdAlgorithmRZ::pushSpectralFields(SpectralFieldDataRZ & f) if (divb_cleaning) { - fields(i,j,k,Bp_avg_m) += -c2 * kr * 0.5_rt * ep0 * X1 * G_old; - fields(i,j,k,Bm_avg_m) += c2 * kr * 0.5_rt * ep0 * X1 * G_old; - fields(i,j,k,Bz_avg_m) += I * c2 * ep0 * X1 * G_old * kz; + fields(i,j,k,Bp_avg_m) += -kr * 0.5_rt * ep0 * X1 * G_old; + fields(i,j,k,Bm_avg_m) += kr * 0.5_rt * ep0 * X1 * G_old; + fields(i,j,k,Bz_avg_m) += I * ep0 * X1 * G_old * kz; } } } |