aboutsummaryrefslogtreecommitdiff
path: root/Source/Parallelization/WarpXComm.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'Source/Parallelization/WarpXComm.cpp')
-rw-r--r--Source/Parallelization/WarpXComm.cpp78
1 files changed, 60 insertions, 18 deletions
diff --git a/Source/Parallelization/WarpXComm.cpp b/Source/Parallelization/WarpXComm.cpp
index 31f651f00..cb87abe4a 100644
--- a/Source/Parallelization/WarpXComm.cpp
+++ b/Source/Parallelization/WarpXComm.cpp
@@ -980,14 +980,29 @@ WarpX::ApplyFilterandSumBoundaryJ (int lev, PatchType patch_type)
const auto& period = Geom(glev).periodicity();
auto& j = (patch_type == PatchType::fine) ? current_fp[lev] : current_cp[lev];
for (int idim = 0; idim < 3; ++idim) {
+ IntVect ng = j[idim]->nGrowVect();
+ IntVect ng_depos_J = get_ng_depos_J();
+ if (WarpX::do_current_centering)
+ {
+#if (AMREX_SPACEDIM == 2)
+ ng_depos_J[0] += WarpX::current_centering_nox / 2;
+ ng_depos_J[1] += WarpX::current_centering_noz / 2;
+#elif (AMREX_SPACEDIM == 3)
+ ng_depos_J[0] += WarpX::current_centering_nox / 2;
+ ng_depos_J[1] += WarpX::current_centering_noy / 2;
+ ng_depos_J[2] += WarpX::current_centering_noz / 2;
+#endif
+ }
if (use_filter) {
- IntVect ng = j[idim]->nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J.min(ng);
MultiFab jf(j[idim]->boxArray(), j[idim]->DistributionMap(), j[idim]->nComp(), ng);
bilinear_filter.ApplyStencil(jf, *j[idim], lev);
- WarpXSumGuardCells(*(j[idim]), jf, period, 0, (j[idim])->nComp());
+ WarpXSumGuardCells(*(j[idim]), jf, period, ng_depos_J, 0, (j[idim])->nComp());
} else {
- WarpXSumGuardCells(*(j[idim]), period, 0, (j[idim])->nComp());
+ ng_depos_J.min(ng);
+ WarpXSumGuardCells(*(j[idim]), period, ng_depos_J, 0, (j[idim])->nComp());
}
}
}
@@ -1019,11 +1034,25 @@ WarpX::AddCurrentFromFineLevelandSumBoundary (int lev)
MultiFab mf(current_fp[lev][idim]->boxArray(),
current_fp[lev][idim]->DistributionMap(), current_fp[lev][idim]->nComp(), 0);
mf.setVal(0.0);
+ IntVect ng = current_cp[lev+1][idim]->nGrowVect();
+ IntVect ng_depos_J = get_ng_depos_J();
+ if (WarpX::do_current_centering)
+ {
+#if (AMREX_SPACEDIM == 2)
+ ng_depos_J[0] += WarpX::current_centering_nox / 2;
+ ng_depos_J[1] += WarpX::current_centering_noz / 2;
+#elif (AMREX_SPACEDIM == 3)
+ ng_depos_J[0] += WarpX::current_centering_nox / 2;
+ ng_depos_J[1] += WarpX::current_centering_noy / 2;
+ ng_depos_J[2] += WarpX::current_centering_noz / 2;
+#endif
+ }
if (use_filter && current_buf[lev+1][idim])
{
// coarse patch of fine level
- IntVect ng = current_cp[lev+1][idim]->nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J.min(ng);
MultiFab jfc(current_cp[lev+1][idim]->boxArray(),
current_cp[lev+1][idim]->DistributionMap(), current_cp[lev+1][idim]->nComp(), ng);
bilinear_filter.ApplyStencil(jfc, *current_cp[lev+1][idim], lev);
@@ -1036,35 +1065,38 @@ WarpX::AddCurrentFromFineLevelandSumBoundary (int lev)
MultiFab::Add(jfb, jfc, 0, 0, current_buf[lev+1][idim]->nComp(), ng);
mf.ParallelAdd(jfb, 0, 0, current_buf[lev+1][idim]->nComp(), ng, IntVect::TheZeroVector(), period);
- WarpXSumGuardCells(*current_cp[lev+1][idim], jfc, period, 0, current_cp[lev+1][idim]->nComp());
+ WarpXSumGuardCells(*current_cp[lev+1][idim], jfc, period, ng_depos_J, 0, current_cp[lev+1][idim]->nComp());
}
else if (use_filter) // but no buffer
{
// coarse patch of fine level
- IntVect ng = current_cp[lev+1][idim]->nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_J.min(ng);
MultiFab jf(current_cp[lev+1][idim]->boxArray(),
current_cp[lev+1][idim]->DistributionMap(), current_cp[lev+1][idim]->nComp(), ng);
bilinear_filter.ApplyStencil(jf, *current_cp[lev+1][idim], lev);
mf.ParallelAdd(jf, 0, 0, current_cp[lev+1][idim]->nComp(), ng, IntVect::TheZeroVector(), period);
- WarpXSumGuardCells(*current_cp[lev+1][idim], jf, period, 0, current_cp[lev+1][idim]->nComp());
+ WarpXSumGuardCells(*current_cp[lev+1][idim], jf, period, ng_depos_J, 0, current_cp[lev+1][idim]->nComp());
}
else if (current_buf[lev+1][idim]) // but no filter
{
+ ng_depos_J.min(ng);
MultiFab::Add(*current_buf[lev+1][idim],
*current_cp [lev+1][idim], 0, 0, current_buf[lev+1][idim]->nComp(),
current_cp[lev+1][idim]->nGrowVect());
mf.ParallelAdd(*current_buf[lev+1][idim], 0, 0, current_buf[lev+1][idim]->nComp(),
current_buf[lev+1][idim]->nGrowVect(), IntVect::TheZeroVector(),
period);
- WarpXSumGuardCells(*(current_cp[lev+1][idim]), period, 0, current_cp[lev+1][idim]->nComp());
+ WarpXSumGuardCells(*(current_cp[lev+1][idim]), period, ng_depos_J, 0, current_cp[lev+1][idim]->nComp());
}
else // no filter, no buffer
{
+ ng_depos_J.min(ng);
mf.ParallelAdd(*current_cp[lev+1][idim], 0, 0, current_cp[lev+1][idim]->nComp(),
current_cp[lev+1][idim]->nGrowVect(), IntVect::TheZeroVector(),
period);
- WarpXSumGuardCells(*(current_cp[lev+1][idim]), period, 0, current_cp[lev+1][idim]->nComp());
+ WarpXSumGuardCells(*(current_cp[lev+1][idim]), period, ng_depos_J, 0, current_cp[lev+1][idim]->nComp());
}
MultiFab::Add(*current_fp[lev][idim], mf, 0, 0, current_fp[lev+1][idim]->nComp(), 0);
}
@@ -1096,14 +1128,18 @@ void
WarpX::ApplyFilterandSumBoundaryRho (int /*lev*/, int glev, amrex::MultiFab& rho, int icomp, int ncomp)
{
const auto& period = Geom(glev).periodicity();
+ IntVect ng = rho.nGrowVect();
+ IntVect ng_depos_rho = get_ng_depos_rho();
if (use_filter) {
- IntVect ng = rho.nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho.min(ng);
MultiFab rf(rho.boxArray(), rho.DistributionMap(), ncomp, ng);
bilinear_filter.ApplyStencil(rf, rho, glev, icomp, 0, ncomp);
- WarpXSumGuardCells(rho, rf, period, icomp, ncomp );
+ WarpXSumGuardCells(rho, rf, period, ng_depos_rho, icomp, ncomp );
} else {
- WarpXSumGuardCells(rho, period, icomp, ncomp);
+ ng_depos_rho.min(ng);
+ WarpXSumGuardCells(rho, period, ng_depos_rho, icomp, ncomp);
}
}
@@ -1134,11 +1170,14 @@ WarpX::AddRhoFromFineLevelandSumBoundary(int lev, int icomp, int ncomp)
rho_fp[lev]->DistributionMap(),
ncomp, 0);
mf.setVal(0.0);
+ IntVect ng = rho_cp[lev+1]->nGrowVect();
+ IntVect ng_depos_rho = get_ng_depos_rho();
if (use_filter && charge_buf[lev+1])
{
// coarse patch of fine level
- IntVect ng = rho_cp[lev+1]->nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho.min(ng);
MultiFab rhofc(rho_cp[lev+1]->boxArray(),
rho_cp[lev+1]->DistributionMap(), ncomp, ng);
bilinear_filter.ApplyStencil(rhofc, *rho_cp[lev+1], lev, icomp, 0, ncomp);
@@ -1150,19 +1189,21 @@ WarpX::AddRhoFromFineLevelandSumBoundary(int lev, int icomp, int ncomp)
MultiFab::Add(rhofb, rhofc, 0, 0, ncomp, ng);
mf.ParallelAdd(rhofb, 0, 0, ncomp, ng, IntVect::TheZeroVector(), period);
- WarpXSumGuardCells( *rho_cp[lev+1], rhofc, period, icomp, ncomp );
+ WarpXSumGuardCells( *rho_cp[lev+1], rhofc, period, ng_depos_rho, icomp, ncomp );
}
else if (use_filter) // but no buffer
{
- IntVect ng = rho_cp[lev+1]->nGrowVect();
ng += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho += bilinear_filter.stencil_length_each_dir-1;
+ ng_depos_rho.min(ng);
MultiFab rf(rho_cp[lev+1]->boxArray(), rho_cp[lev+1]->DistributionMap(), ncomp, ng);
bilinear_filter.ApplyStencil(rf, *rho_cp[lev+1], lev, icomp, 0, ncomp);
mf.ParallelAdd(rf, 0, 0, ncomp, ng, IntVect::TheZeroVector(), period);
- WarpXSumGuardCells( *rho_cp[lev+1], rf, period, icomp, ncomp );
+ WarpXSumGuardCells( *rho_cp[lev+1], rf, period, ng_depos_rho, icomp, ncomp );
}
else if (charge_buf[lev+1]) // but no filter
{
+ ng_depos_rho.min(ng);
MultiFab::Add(*charge_buf[lev+1],
*rho_cp[lev+1], icomp, icomp, ncomp,
rho_cp[lev+1]->nGrowVect());
@@ -1170,14 +1211,15 @@ WarpX::AddRhoFromFineLevelandSumBoundary(int lev, int icomp, int ncomp)
ncomp,
charge_buf[lev+1]->nGrowVect(), IntVect::TheZeroVector(),
period);
- WarpXSumGuardCells(*(rho_cp[lev+1]), period, icomp, ncomp);
+ WarpXSumGuardCells(*(rho_cp[lev+1]), period, ng_depos_rho, icomp, ncomp);
}
else // no filter, no buffer
{
+ ng_depos_rho.min(ng);
mf.ParallelAdd(*rho_cp[lev+1], icomp, 0, ncomp,
rho_cp[lev+1]->nGrowVect(), IntVect::TheZeroVector(),
period);
- WarpXSumGuardCells(*(rho_cp[lev+1]), period, icomp, ncomp);
+ WarpXSumGuardCells(*(rho_cp[lev+1]), period, ng_depos_rho, icomp, ncomp);
}
MultiFab::Add(*rho_fp[lev], mf, 0, icomp, ncomp, 0);
NodalSyncRho(lev+1, PatchType::coarse, icomp, ncomp);
generated by cgit v1.2.3 (git 2.39.1) at 2025-08-09 11:22:18 +0000