Skip to content

Commit

Permalink
remove implementation for smoothing
Browse files Browse the repository at this point in the history
  • Loading branch information
@RevathiJambunathan
RevathiJambunathan committed May 24, 2024
1 parent 19566fa commit 0f33db5
Showing 1 changed file with 70 additions and 70 deletions.
140 changes: 70 additions & 70 deletions Source/WarpX.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(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<amrex::Real>(kk-k)-ngbuffer*tanh_midpoint)
// // // / ((1.-tanh_midpoint)*(ngbuffer/3.));
// // // wtmsk(i,j,k) = std::tanh(arg)*0.5+0.5;
// // // return;
// // // }
// // //}
// // }
// //});
//}
}
}
}
Expand Down

0 comments on commit 0f33db5

Please sign in to comment.