From 0f33db576fcdb245bf98b824b4355d1b54c46fb6 Mon Sep 17 00:00:00 2001 From: RevathiJambunathan Date: Thu, 23 May 2024 20:37:57 -0400 Subject: [PATCH] remove implementation for smoothing --- Source/WarpX.cpp | 140 +++++++++++++++++++++++------------------------ 1 file changed, 70 insertions(+), 70 deletions(-) diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index 5d11ea5c209..d92d2d53606 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -3142,76 +3142,76 @@ WarpX::BuildBufferMasks () const Box tbx = mfi.growntilebox(); BuildBufferMasksInBox( tbx, (*bmasks)[mfi], tmp[mfi], ngbuffer ); } - if (ipass==0) continue; - amrex::MultiFab* weight_gbuffer = interp_weight_gbuffer[lev].get(); - // Using tmp to also set weights in the interp_weight_gbuffer multifab - for (MFIter mfi(*weight_gbuffer, true); mfi.isValid(); ++mfi) - { - const Box& tbx = mfi.tilebox(IntVect::TheNodeVector(),weight_gbuffer->nGrowVect()); - auto const& gmsk = tmp[mfi].const_array(); - auto const& bmsk = (*bmasks)[mfi].array(); - auto const& wtmsk = (*weight_gbuffer)[mfi].array(); - SetWeightsInGatherBuffer(tbx, wtmsk, gmsk, bmsk, ngbuffer, do_interpolate, tanh_midpoint); - //amrex::ParallelFor(tbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) { - // wtmsk(i,j,k) = 0._rt; - // if (bmsk(i,j,k) == 0 && do_interpolate) { - // if(gmsk(i,j,k)==0) { - // wtmsk(i,j,k) = 0.; - // return; - // } - // for (int ii = i-1; ii >= i-ngbuffer; --ii) { - // if (gmsk(ii,j,k)==0) { - // amrex::Real arg = (static_cast(i-ii)-ngbuffer*tanh_midpoint) - // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // wtmsk(i,j,k) = std::tanh(arg)*0.5 + 0.5; - // amrex::Print() << " i edge wt is " << wtmsk(i,j,k) << "\n"; - // return; - // } - // } - // for (int ii = i+1; ii <= i+ngbuffer; ++ii) { - // if (gmsk(ii,j,k)==0) { - // amrex::Real arg = (static_cast(ii-i)-ngbuffer*tanh_midpoint) - // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; - // amrex::Print() << " wt is " << wtmsk(i,j,k) << "\n"; - // return; - // } - // } - // //for (int jj = j-1; jj >= j-ngbuffer; --jj) { - // // if (gmsk(i,jj,k)==0) { - // // amrex::Real arg = (static_cast(j-jj)-ngbuffer*tanh_midpoint) - // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // // wtmsk(i,j,k) = std::tanh(arg)*0.5 + 0.5; - // // return; - // // } - // //} - // //for (int jj = j+1; jj <= j+ngbuffer; ++jj) { - // // if (gmsk(i,jj,k)==0) { - // // amrex::Real arg = (static_cast(jj - j)-ngbuffer*tanh_midpoint) - // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; - // // return; - // // } - // //} - // //for (int kk = k-1; kk >= k-ngbuffer; --kk) { - // // if (gmsk(i,j,kk)==0) { - // // amrex::Real arg = (static_cast(k-kk)-ngbuffer*tanh_midpoint) - // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; - // // return; - // // } - // //} - // //for (int kk = k+1; kk <= k+ngbuffer; ++kk) { - // // if (gmsk(i,j,kk)==0) { - // // amrex::Real arg = (static_cast(kk-k)-ngbuffer*tanh_midpoint) - // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); - // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; - // // return; - // // } - // //} - // } - //}); - } + //if (ipass==0) continue; + //amrex::MultiFab* weight_gbuffer = interp_weight_gbuffer[lev].get(); + //// Using tmp to also set weights in the interp_weight_gbuffer multifab + //for (MFIter mfi(*weight_gbuffer, true); mfi.isValid(); ++mfi) + //{ + // const Box& tbx = mfi.tilebox(IntVect::TheNodeVector(),weight_gbuffer->nGrowVect()); + // auto const& gmsk = tmp[mfi].const_array(); + // auto const& bmsk = (*bmasks)[mfi].array(); + // auto const& wtmsk = (*weight_gbuffer)[mfi].array(); + // SetWeightsInGatherBuffer(tbx, wtmsk, gmsk, bmsk, ngbuffer, do_interpolate, tanh_midpoint); + // //amrex::ParallelFor(tbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) { + // // wtmsk(i,j,k) = 0._rt; + // // if (bmsk(i,j,k) == 0 && do_interpolate) { + // // if(gmsk(i,j,k)==0) { + // // wtmsk(i,j,k) = 0.; + // // return; + // // } + // // for (int ii = i-1; ii >= i-ngbuffer; --ii) { + // // if (gmsk(ii,j,k)==0) { + // // amrex::Real arg = (static_cast(i-ii)-ngbuffer*tanh_midpoint) + // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // wtmsk(i,j,k) = std::tanh(arg)*0.5 + 0.5; + // // amrex::Print() << " i edge wt is " << wtmsk(i,j,k) << "\n"; + // // return; + // // } + // // } + // // for (int ii = i+1; ii <= i+ngbuffer; ++ii) { + // // if (gmsk(ii,j,k)==0) { + // // amrex::Real arg = (static_cast(ii-i)-ngbuffer*tanh_midpoint) + // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; + // // amrex::Print() << " wt is " << wtmsk(i,j,k) << "\n"; + // // return; + // // } + // // } + // // //for (int jj = j-1; jj >= j-ngbuffer; --jj) { + // // // if (gmsk(i,jj,k)==0) { + // // // amrex::Real arg = (static_cast(j-jj)-ngbuffer*tanh_midpoint) + // // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // // wtmsk(i,j,k) = std::tanh(arg)*0.5 + 0.5; + // // // return; + // // // } + // // //} + // // //for (int jj = j+1; jj <= j+ngbuffer; ++jj) { + // // // if (gmsk(i,jj,k)==0) { + // // // amrex::Real arg = (static_cast(jj - j)-ngbuffer*tanh_midpoint) + // // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; + // // // return; + // // // } + // // //} + // // //for (int kk = k-1; kk >= k-ngbuffer; --kk) { + // // // if (gmsk(i,j,kk)==0) { + // // // amrex::Real arg = (static_cast(k-kk)-ngbuffer*tanh_midpoint) + // // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; + // // // return; + // // // } + // // //} + // // //for (int kk = k+1; kk <= k+ngbuffer; ++kk) { + // // // if (gmsk(i,j,kk)==0) { + // // // amrex::Real arg = (static_cast(kk-k)-ngbuffer*tanh_midpoint) + // // // / ((1.-tanh_midpoint)*(ngbuffer/3.)); + // // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5; + // // // return; + // // // } + // // //} + // // } + // //}); + //} } } }