Skip to content

Commit

Permalink
EOL
Browse files Browse the repository at this point in the history
  • Loading branch information
RevathiJambunathan committed Feb 21, 2024
1 parent 6351ddc commit 4c6d9c6
Show file tree
Hide file tree
Showing 5 changed files with 4 additions and 6 deletions.
1 change: 0 additions & 1 deletion Tests/Particles/Mapped/InitHelix.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,6 @@ InitHelix (amrex::MultiFab& a_xyz_loc, amrex::Geometry& geom)
Real xi = (static_cast<int>(i)) * d_xi;
Real eta = (static_cast<int>(j)) * d_eta;
Real zeta = (static_cast<int>(k)) * d_zeta;

loc_arr(i,j,k,0) = cx + (0.1 * xi + 0.1) * cos(tpi * eta);
loc_arr(i,j,k,1) = cy + (0.1 * xi + 0.1) * sin(tpi * eta);
//loc_arr(i,j,k,2) = cz + ( xi * 0.1 - 0.1 ) + zeta * 0.1 * (2. + 1. * xi); 3D torus
Expand Down
1 change: 0 additions & 1 deletion Tests/Particles/Mapped/InitTorus.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,6 @@ InitTorus (amrex::MultiFab& a_xyz_loc, amrex::Geometry& geom)
Real xi = (static_cast<int>(i)) * d_xi;
Real eta = (static_cast<int>(j)) * d_eta;
Real zeta = (static_cast<int>(k)) * d_zeta;

loc_arr(i,j,k,0) = cx + (0.1 * xi + 0.1) * cos(tpi * eta);
loc_arr(i,j,k,1) = cy + (0.1 * xi + 0.1) * sin(tpi * eta);
loc_arr(i,j,k,2) = cz + ( xi * 0.1 - 0.1 ) + zeta * 0.1 * (2. + 1. * xi); //3D torus
Expand Down
2 changes: 1 addition & 1 deletion Tests/Particles/Mapped/MappedPC.H
Original file line number Diff line number Diff line change
Expand Up @@ -265,7 +265,7 @@ public:
}
}}}
}
if (!inside) {
if (!inside) {
// // looping over neighboring cells
// for (int lo_k = -2; lo_k <=2; ++lo_k){
// for (int lo_j = -2; lo_j <=2; ++lo_j){
Expand Down
2 changes: 1 addition & 1 deletion Tests/Particles/Mapped/MappedPC.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -274,7 +274,7 @@ MappedPC::AdvectWithUND (MultiFab& vel_nd, int lev, Real dt, const MultiFab& a_x
amrex::Real rc = sqrt((p.pos(0)-0.5)*(p.pos(0)-0.5) + (p.pos(1)-0.5)*(p.pos(1)-0.5));
amrex::Real theta = std::atan2( (p.pos(2) - 0.5) , (p.pos(1) - 0.5));
amrex::Print() << "TO " << p.pos(0) << " " << p.pos(1) << " " << p.pos(AMREX_SPACEDIM-1) << " " << rc << " " << theta <<
" WITH VEL " << v[0] << " " << v[1] << " " << v[2] << std::endl;
" WITH VEL " << v[0] << " " << v[1] << " " << v[2] << std::endl;
#endif
update_mapped_idata(p,plo,dxi,loc_arr);
// amrex::AllPrintToFile("partpos.txt") << p.pos(0) << " " << p.pos(1) << " " << p.pos(2)<< " \n";
Expand Down
4 changes: 2 additions & 2 deletions Tests/Particles/Mapped/main.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -291,7 +291,7 @@ WriteGenericPlotfileHeaderWithMapped (std::ostream &HeaderFile,
amrex::Vector<IntVect> ref_ratio(nlevels);
for (int i = 0 ; i <nlevels; ++i) {
ref_ratio[i] = amrex::IntVect::TheUnitVector();
}
}
BL_ASSERT(nlevels <= bArray.size());
BL_ASSERT(nlevels <= ref_ratio.size()+1);
BL_ASSERT(nlevels <= level_steps.size());
Expand Down Expand Up @@ -682,7 +682,7 @@ void Test()
GetVecOfConstPtrs(mfarr),
GetVecOfConstPtrs(mf_nd),
varname,
nt, level_steps, versionName, levelPrefix, mfPrefix, extra_dirs, geomarr);
nt, level_steps, versionName, levelPrefix, mfPrefix, extra_dirs, geomarr);
mapped_pc.WritePlotFile(plotfilename,"particles");

//WriteSingleLevelPlotfile("plt_grid",a_xyz_loc,{"gridmap",AMREX_D_DECL("x1","y1","z1")},geom,0.0,0);
Expand Down

0 comments on commit 4c6d9c6

Please sign in to comment.