diff options
author | 2019-08-26 09:59:33 -0700 | |
---|---|---|
committer | 2019-08-26 09:59:33 -0700 | |
commit | 6fcf8b97a6bb432b08c456a275d551d53b998be6 (patch) | |
tree | 574b35c269f77c396fad367672d4fabd00265ac0 /Source/Evolve/WarpXEvolveEM.cpp | |
parent | 0db327debd4a4f704ac5059a9038ef183e22eee4 (diff) | |
download | WarpX-6fcf8b97a6bb432b08c456a275d551d53b998be6.tar.gz WarpX-6fcf8b97a6bb432b08c456a275d551d53b998be6.tar.zst WarpX-6fcf8b97a6bb432b08c456a275d551d53b998be6.zip |
Clean up of RZ mutlimode, mostly removing imaginary part of mode 0
Diffstat (limited to 'Source/Evolve/WarpXEvolveEM.cpp')
-rw-r--r-- | Source/Evolve/WarpXEvolveEM.cpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/Source/Evolve/WarpXEvolveEM.cpp b/Source/Evolve/WarpXEvolveEM.cpp index 7573b572e..5313fb418 100644 --- a/Source/Evolve/WarpXEvolveEM.cpp +++ b/Source/Evolve/WarpXEvolveEM.cpp @@ -356,7 +356,7 @@ WarpX::OneStep_sub1 (Real curtime) RestrictRhoFromFineToCoarsePatch(fine_lev); ApplyFilterandSumBoundaryJ(fine_lev, PatchType::fine); NodalSyncJ(fine_lev, PatchType::fine); - ApplyFilterandSumBoundaryRho(fine_lev, PatchType::fine, 0, 2*n_rz_azimuthal_modes); + ApplyFilterandSumBoundaryRho(fine_lev, PatchType::fine, 0, rho_fp[fine_lev]->nComp()); NodalSyncRho(fine_lev, PatchType::fine, 0, 2); EvolveB(fine_lev, PatchType::fine, 0.5*dt[fine_lev]); @@ -383,7 +383,7 @@ WarpX::OneStep_sub1 (Real curtime) PushParticlesandDepose(coarse_lev, curtime); StoreCurrent(coarse_lev); AddCurrentFromFineLevelandSumBoundary(coarse_lev); - AddRhoFromFineLevelandSumBoundary(coarse_lev, 0, n_rz_azimuthal_modes); + AddRhoFromFineLevelandSumBoundary(coarse_lev, 0, rho_cp[coarse_lev]->nComp()/2); EvolveB(fine_lev, PatchType::coarse, dt[fine_lev]); EvolveF(fine_lev, PatchType::coarse, dt[fine_lev], DtType::FirstHalf); @@ -410,7 +410,7 @@ WarpX::OneStep_sub1 (Real curtime) RestrictRhoFromFineToCoarsePatch(fine_lev); ApplyFilterandSumBoundaryJ(fine_lev, PatchType::fine); NodalSyncJ(fine_lev, PatchType::fine); - ApplyFilterandSumBoundaryRho(fine_lev, PatchType::fine, 0, 2*n_rz_azimuthal_modes); + ApplyFilterandSumBoundaryRho(fine_lev, PatchType::fine, 0, rho_fp[fine_lev]->nComp()/2); NodalSyncRho(fine_lev, PatchType::fine, 0, 2); EvolveB(fine_lev, PatchType::fine, 0.5*dt[fine_lev]); @@ -436,7 +436,7 @@ WarpX::OneStep_sub1 (Real curtime) // by only half a coarse step (second half) RestoreCurrent(coarse_lev); AddCurrentFromFineLevelandSumBoundary(coarse_lev); - AddRhoFromFineLevelandSumBoundary(coarse_lev, n_rz_azimuthal_modes, n_rz_azimuthal_modes); + AddRhoFromFineLevelandSumBoundary(coarse_lev, rho_cp[coarse_lev]->nComp()/2, rho_cp[coarse_lev]->nComp()/2); EvolveE(fine_lev, PatchType::coarse, dt[fine_lev]); FillBoundaryE(fine_lev, PatchType::coarse); |