570 int nlev =
static_cast<int>(vars_new.size());
571 int nplane =
static_cast<int>(
m_bnd_rbx.size());
573 bool interpolate =
true;
576 for (
int iplane(0); iplane<nplane; ++iplane) {
577 int dir =
m_dir[iplane];
578 amrex::RealBox bnd_rbx =
m_bnd_rbx[iplane];
582 for (
int ilev(nlev-1); ilev>=0; --ilev) {
585 amrex::MultiFab mf_cc_vel;
586 auto ba = vars_new[ilev][
Vars::cons].boxArray();
587 auto dm = vars_new[ilev][
Vars::cons].DistributionMap();
588 mf_cc_vel.define(ba, dm, AMREX_SPACEDIM, amrex::IntVect(1,1,1));
589 average_face_to_cellcenter(mf_cc_vel,0,
590 amrex::Array<const amrex::MultiFab*,3>{&vars_new[ilev][
Vars::xvel],
595 amrex::MultiFab mf_cc_data;
596 mf_cc_data.define(ba, dm, ncomp, 1);
598 #pragma omp parallel if (amrex::Gpu::notInLaunchRegion())
600 for (amrex::MFIter mfi(mf_cc_data, amrex::TilingIfNotGPU()); mfi.isValid(); ++mfi) {
601 const amrex::Box& tbx = mfi.tilebox();
602 auto const& dfab = mf_cc_data.array(mfi);
603 auto const& tfab = vars_new[ilev][
Vars::cons].array(mfi);
604 auto const& wfab = mf_cc_vel.array(mfi);
605 amrex::ParallelFor(tbx, [=] AMREX_GPU_DEVICE(
int i,
int j,
int k) noexcept
607 dfab(i,j,k,0) = tfab(i,j,k,1)/tfab(i,j,k,0);
608 dfab(i,j,k,1) = std::sqrt(wfab(i,j,k,0)*wfab(i,j,k,0)
609 + wfab(i,j,k,1)*wfab(i,j,k,1)
610 + wfab(i,j,k,2)*wfab(i,j,k,2)) ;
615 m_lev[iplane] = ilev;
616 m_ps_mf[iplane] = get_slice_data(dir, point, mf_cc_data, geom[ilev],
617 0, ncomp, interpolate, bnd_rbx);
620 auto min_bnd_bx =
m_ps_mf[iplane]->boxArray().minimalBox();
621 amrex::Box bnd_bx =
getIndexBox(bnd_rbx, geom[ilev]);
622 if (bnd_bx == min_bnd_bx) {
break; }
amrex::Real Real
Definition: ERF_ShocInterface.H:19
@ xvel
Definition: ERF_IndexDefines.H:141
@ cons
Definition: ERF_IndexDefines.H:140
@ zvel
Definition: ERF_IndexDefines.H:143
@ yvel
Definition: ERF_IndexDefines.H:142
amrex::Box getIndexBox(const amrex::RealBox &real_box, const amrex::Geometry &geom)
Definition: ERF_SampleData.H:551