1#ifndef AMREX_ML_CELL_LINOP_H_
2#define AMREX_ML_CELL_LINOP_H_
3#include <AMReX_Config.H>
24 using FAB =
typename MF::fab_type;
25 using RT =
typename MF::value_type;
46 void setLevelBC (
int amrlev,
const MF* levelbcdata,
47 const MF* robinbc_a =
nullptr,
48 const MF* robinbc_b =
nullptr,
49 const MF* robinbc_f =
nullptr) final;
63 void updateSolBC (
int amrlev,
const MF& crse_bcdata)
const;
64 void updateCorBC (
int amrlev,
const MF& crse_bcdata)
const;
67 const MLMGBndryT<MF>* bndry=
nullptr,
bool skip_fillboundary=
false)
const;
78 IntVect const& nghost)
const override;
81 const MF& fine_sol,
const MF& fine_rhs)
override;
83 void apply (
int amrlev,
int mglev, MF& out, MF& in,
BCMode bc_mode,
85 void smooth (
int amrlev,
int mglev, MF& sol,
const MF& rhs,
86 bool skip_fillboundary=
false) const final;
89 const MF* crse_bcdata=
nullptr) override;
91 void prepareForFluxes (
int amrlev, const MF* crse_bcdata =
nullptr) override;
94 BCMode bc_mode, const MF* crse_bcdata=
nullptr) final;
97 void reflux (
int crse_amrlev,
98 MF& res, const MF& crse_sol, const MF&,
99 MF&, MF& fine_sol, const MF&) const final;
100 void compFlux (
int amrlev, const
Array<MF*,AMREX_SPACEDIM>& fluxes,
101 MF& sol,
Location loc) const override;
102 void compGrad (
int amrlev, const
Array<MF*,AMREX_SPACEDIM>& grad,
103 MF& sol,
Location loc) const override;
108 MF const& rhs) const override;
114 RT xdoty (
int amrlev,
int mglev, const MF&
x, const MF& y,
bool local) const final;
120 virtual
void Fapply (
int amrlev,
int mglev, MF& out, const MF& in) const = 0;
121 virtual
void Fsmooth (
int amrlev,
int mglev, MF& sol, const MF& rhs,
int redblack) const = 0;
123 const
Array<
FAB*,AMREX_SPACEDIM>& flux,
124 const
FAB& sol,
Location loc,
int face_only=0) const = 0;
127 const
Array<MF*,AMREX_SPACEDIM>& ,
131 RT normInf (
int amrlev, MF
const& mf,
bool local)
const override;
135 void avgDownResAmr (
int clev, MF& cres, MF
const& fres)
const override;
177 LinOpBCType crse_fine_bc_type);
186 return bcond[mfi][icomp];
189 return bcloc[mfi][icomp];
271 Array4<int const> mlo;
272 Array4<int const> mhi;
283 Box
const& box() const noexcept {
return bx; }
287template <
typename MF>
294 bctl_dv(bctl.local_size()*ncomp),
299 bcond[mfi].resize(ncomp);
300 bcloc[mfi].resize(ncomp);
306template <
typename MF>
315 LinOpBCType crse_fine_bc_type)
324 const Box& bx = mfi.validbox();
325 for (
int icomp = 0; icomp < m_ncomp; ++icomp) {
327 BCTuple & bctag = bcond[mfi][icomp];
329 lobc[icomp], hibc[icomp],
330 dx, ratio, interior_bloc,
331 domain_bloc_lo, domain_bloc_hi,
341 for (
int icomp = 0; icomp < m_ncomp; ++icomp) {
343 for (
int m = 0; m < 2*AMREX_SPACEDIM; ++m) {
344 tmp[m].type = bcond[mfi][icomp][m];
345 tmp[m].location = bcloc[mfi][icomp][m];
354template <
typename MF>
360template <
typename MF>
373template <
typename MF>
392 this->
m_dmap[amrlev][mglev],
408 this->
m_dmap[amrlev][mglev],
409 this->
m_geom[amrlev][mglev],
410 face, 0, ngrow, extent, 1,
true);
415 for (
int amrlev = 0; amrlev < this->m_num_amr_levels-1; ++amrlev)
420 this->
m_dmap[amrlev+1][0],
422 this->
m_geom[amrlev+1][0],
424 ratio, amrlev+1, ncomp);
431#if (AMREX_SPACEDIM != 3)
436template <
typename MF>
462 const int in_rad = 0;
463 const int out_rad = 1;
464 const int extent_rad = 2;
469 (cba, this->
m_dmap[amrlev][0], in_rad, out_rad, extent_rad, ncomp);
474 const int in_rad = 0;
475 const int out_rad = 1;
476 const int extent_rad = 2;
481 (cba, this->
m_dmap[amrlev][0], in_rad, out_rad, extent_rad, ncomp);
488 m_bndry_cor[amrlev] = std::make_unique<MLMGBndryT<MF>>
490 MF bc_data(this->m_grids[amrlev][0], this->m_dmap[amrlev][0], ncomp, 1);
501 BCType::Dirichlet)}});
511 m_bcondloc[amrlev][mglev] = std::make_unique<BndryCondLoc>(this->
m_grids[amrlev][mglev],
512 this->
m_dmap[amrlev][mglev],
518template <
typename MF>
521 const MF* robinbc_b,
const MF* robinbc_f)
532 if (a_levelbcdata ==
nullptr) {
538 const MF& bcdata = (a_levelbcdata ==
nullptr) ?
zero : *a_levelbcdata;
554 const int in_rad = 0;
555 const int out_rad = 1;
556 const int extent_rad = 2;
557 const IntVect crse_ratio = br_ref_ratio;
561 (cba, this->
m_dmap[amrlev][0], in_rad, out_rad, extent_rad, ncomp);
572 bcdata, 0, 0, ncomp, br_ref_ratio,
579 m_bndry_sol[amrlev]->setPhysBndryValues(bcdata,0,0,ncomp);
585 m_bndry_sol[amrlev]->setPhysBndryValues(bcdata,0,0,ncomp);
593 const Real* dx = this->
m_geom[amrlev][0].CellSize();
604 AMREX_ASSERT(robinbc_a !=
nullptr && robinbc_b !=
nullptr && robinbc_f !=
nullptr);
608 const Box& domain = this->
m_geom[amrlev][0].Domain();
612#pragma omp parallel if (Gpu::notInLaunchRegion())
615 Box const& vbx = mfi.validbox();
619 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
622 bool outside_domain_lo = !(domain.
contains(blo));
623 bool outside_domain_hi = !(domain.
contains(bhi));
624 if ((!outside_domain_lo) && (!outside_domain_hi)) {
continue; }
625 for (
int icomp = 0; icomp < ncomp; ++icomp) {
627 if (this->
m_lobc_orig[icomp][idim] == LinOpBCType::Robin && outside_domain_lo)
631 rbc(i,j,k,0) = ra(i,j,k,icomp);
632 rbc(i,j,k,1) = rb(i,j,k,icomp);
633 rbc(i,j,k,2) = rf(i,j,k,icomp);
636 if (this->
m_hibc_orig[icomp][idim] == LinOpBCType::Robin && outside_domain_hi)
640 rbc(i,j,k,0) = ra(i,j,k,icomp);
641 rbc(i,j,k,1) = rb(i,j,k,icomp);
642 rbc(i,j,k,2) = rf(i,j,k,icomp);
651template <
typename MF>
658template <
typename MF>
667 this->
m_geom[amrlev-1][0].periodicity());
674template <
typename MF>
682 this->
m_geom[amrlev-1][0].periodicity());
689template <
typename MF>
696 BL_ASSERT(mglev == 0 || bc_mode == BCMode::Homogeneous);
697 BL_ASSERT(bndry !=
nullptr || bc_mode == BCMode::Homogeneous);
702 if (!skip_fillboundary) {
703 in.FillBoundary(0, ncomp, this->
m_geom[amrlev][mglev].periodicity(), cross);
706 int flagbc = bc_mode == BCMode::Inhomogeneous;
707 const int imaxorder = this->
maxorder;
709 const Real* dxinv = this->
m_geom[amrlev][mglev].InvCellSize();
710 const RT dxi =
static_cast<RT>(dxinv[0]);
711 const RT dyi = (AMREX_SPACEDIM >= 2) ?
static_cast<RT>(dxinv[1]) :
RT(1.0);
712 const RT dzi = (AMREX_SPACEDIM == 3) ?
static_cast<RT>(dxinv[2]) :
RT(1.0);
714 const auto& maskvals =
m_maskvals[amrlev][mglev];
715 const auto& bcondloc = *
m_bcondloc[amrlev][mglev];
718 const auto& foo = foofab.const_array();
724 "non-cross stencil not support for gpu");
732 tags.reserve(in.local_size()*AMREX_SPACEDIM*ncomp);
735 const Box& vbx = mfi.validbox();
736 const auto& iofab = in.array(mfi);
737 const auto & bdlv = bcondloc.bndryLocs(mfi);
738 const auto & bdcv = bcondloc.bndryConds(mfi);
740 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
741 if (idim != hidden_direction) {
744 const auto& bvlo = (bndry !=
nullptr) ?
746 const auto& bvhi = (bndry !=
nullptr) ?
748 for (
int icomp = 0; icomp < ncomp; ++icomp) {
750 maskvals[olo].const_array(mfi),
751 maskvals[ohi].const_array(mfi),
752 bdlv[icomp][olo], bdlv[icomp][ohi],
754 bdcv[icomp][olo], bdcv[icomp][ohi],
755 vbx.
length(idim), icomp, idim});
766 mllinop_apply_bc_x(0, i, j, k, tag.blen, tag.fab,
767 tag.mask_lo, tag.bctype_lo, tag.bcloc_lo, tag.bcval_lo,
768 imaxorder, dxi, flagbc, tag.comp);
769 mllinop_apply_bc_x(1, i+tag.blen+1, j, k, tag.blen, tag.fab,
770 tag.mask_hi, tag.bctype_hi, tag.bcloc_hi, tag.bcval_hi,
771 imaxorder, dxi, flagbc, tag.comp);
773#
if (AMREX_SPACEDIM > 1)
775#
if (AMREX_SPACEDIM > 2)
779 mllinop_apply_bc_y(0, i, j, k, tag.blen, tag.fab,
780 tag.mask_lo, tag.bctype_lo, tag.bcloc_lo, tag.bcval_lo,
781 imaxorder, dyi, flagbc, tag.comp);
782 mllinop_apply_bc_y(1, i, j+tag.blen+1, k, tag.blen, tag.fab,
783 tag.mask_hi, tag.bctype_hi, tag.bcloc_hi, tag.bcval_hi,
784 imaxorder, dyi, flagbc, tag.comp);
786#
if (AMREX_SPACEDIM > 2)
790 imaxorder, dzi, flagbc, tag.
comp);
793 imaxorder, dzi, flagbc, tag.
comp);
800 if (cross || tensorop)
803#pragma omp parallel if (Gpu::notInLaunchRegion())
807 const Box& vbx = mfi.validbox();
808 const auto& iofab = in.array(mfi);
810 const auto & bdlv = bcondloc.bndryLocs(mfi);
811 const auto & bdcv = bcondloc.bndryConds(mfi);
813 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim)
815 if (hidden_direction == idim) {
continue; }
820 const int blen = vbx.
length(idim);
821 const auto& mlo = maskvals[olo].array(mfi);
822 const auto& mhi = maskvals[ohi].array(mfi);
823 const auto& bvlo = (bndry !=
nullptr) ? bndry->
bndryValues(olo).const_array(mfi) : foo;
824 const auto& bvhi = (bndry !=
nullptr) ? bndry->
bndryValues(ohi).const_array(mfi) : foo;
825 for (
int icomp = 0; icomp < ncomp; ++icomp) {
826 const BoundCond bctlo = bdcv[icomp][olo];
827 const BoundCond bcthi = bdcv[icomp][ohi];
828 const RT bcllo = bdlv[icomp][olo];
829 const RT bclhi = bdlv[icomp][ohi];
833 imaxorder, dxi, flagbc, icomp);
836 imaxorder, dxi, flagbc, icomp);
837 }
else if (idim == 1) {
840 imaxorder, dyi, flagbc, icomp);
843 imaxorder, dyi, flagbc, icomp);
847 imaxorder, dzi, flagbc, icomp);
850 imaxorder, dzi, flagbc, icomp);
859 amrex::Abort(
"amrex_mllinop_apply_bc not available when BL_NO_FORT=TRUE");
861 if constexpr (std::is_same_v<Real,RT>) {
867 const Box& vbx = mfi.validbox();
869 const auto & bdlv = bcondloc.bndryLocs(mfi);
870 const auto & bdcv = bcondloc.bndryConds(mfi);
883 const auto& fsfab = (bndry !=
nullptr) ? bndry->
bndryValues(ori)[mfi] : foofab;
885 const Mask& m = maskvals[ori][mfi];
888 BL_TO_FORTRAN_ANYD(in[mfi]),
889 BL_TO_FORTRAN_ANYD(m),
891 BL_TO_FORTRAN_ANYD(fsfab),
892 imaxorder, dxinv, flagbc, ncomp, cross);
902template <
typename MF>
906 const Box& dombx = this->
m_geom[0].back().Domain();
909 const int N = old_ba.
size();
912 for (
int i = 0; i < N; ++i)
920 IntVect big =
b.smallEnd() + grid_size - 1;
923#if (AMREX_SPACEDIM == 3)
924 for (
int kk = 0; kk < nblks[2]; ++kk) {
926#if (AMREX_SPACEDIM >= 2)
927 for (
int jj = 0; jj < nblks[1]; ++jj) {
929 for (
int ii = 0; ii < nblks[0]; ++ii)
936#if (AMREX_SPACEDIM >= 2)
939#if (AMREX_SPACEDIM == 3)
944 std::sort(bv.begin(), bv.end());
945 bv.erase(std::unique(bv.begin(), bv.end()), bv.end());
952template <
typename MF>
961template <
typename MF>
967 Dim3 ratio3 = {2,2,2};
970 ratio3.
y = ratio[1];,
971 ratio3.
z = ratio[2];);
975 auto const& finema =
fine.arrays();
976 auto const& crsema =
crse.const_arrays();
983 finema[box_no](i,j,k,n) += crsema[box_no](ic,jc,kc,n);
990#pragma omp parallel if (Gpu::notInLaunchRegion())
994 const Box& bx = mfi.tilebox();
1002 ffab(i,j,k,n) += cfab(ic,jc,kc,n);
1008template <
typename MF>
1012 const int ncomp = this->
getNComp();
1014 const Geometry& crse_geom = this->
Geom(amrlev,fmglev+1);
1030 cfine.
define(cba,
fine.DistributionMap(), ncomp, ng);
1031 cfine.setVal(
RT(0.0));
1036 bool isEB =
fine.hasEBFabFactory();
1047#pragma omp parallel if (Gpu::notInLaunchRegion())
1051 const Box& bx = mfi.tilebox();
1052 const auto& ff =
fine.array(mfi);
1053 const auto& cc = cmf->array(mfi);
1058 const auto& flag = (*flags)[mfi];
1065 mlmg_eb_cc_interp_r<2>(tbx, ff, cc, flg, ncomp);
1076 const bool call_lincc =
true;
1080#if (AMREX_SPACEDIM == 3)
1101template <
typename MF>
1106 const int ncomp = this->
getNComp();
1107 const int refratio = this->
AMRRefRatio(famrlev-1);
1117#pragma omp parallel if (Gpu::notInLaunchRegion())
1121 const Box& bx = mfi.tilebox();
1122 auto const& ff =
fine.array(mfi);
1123 auto const& cc =
crse.const_array(mfi);
1128 const auto& flag = (*flags)[mfi];
1138 mlmg_eb_cc_interp_r<2>(tbx, ff, cc, flg, ncomp);
1146 mlmg_eb_cc_interp_r<4>(tbx, ff, cc, flg, ncomp);
1151 amrex::Abort(
"mlmg_eb_cc_interp: only refratio 2 and 4 are supported");
1162 const bool call_lincc =
true;
1184 amrex::Abort(
"mlmg_lin_cc_interp: only refratio 2 and 4 are supported");
1190template <
typename MF>
1193 const MF& fine_sol,
const MF& fine_rhs)
1196 const int ncomp = this->
getNComp();
1201template <
typename MF>
1207 applyBC(amrlev, mglev, in, bc_mode, s_mode, bndry);
1208 Fapply(amrlev, mglev, out, in);
1211template <
typename MF>
1214 bool skip_fillboundary)
const
1217 for (
int redblack = 0; redblack < 2; ++redblack)
1219 applyBC(amrlev, mglev, sol, BCMode::Homogeneous, StateMode::Solution,
1220 nullptr, skip_fillboundary);
1221 Fsmooth(amrlev, mglev, sol, rhs, redblack);
1222 skip_fillboundary =
false;
1226template <
typename MF>
1229 const MF* crse_bcdata)
1231 BL_PROFILE(
"MLCellLinOp::solutionResidual()");
1232 const int ncomp = this->
getNComp();
1233 if (crse_bcdata !=
nullptr) {
1236 const int mglev = 0;
1237 apply(amrlev, mglev, resid,
x, BCMode::Inhomogeneous, StateMode::Solution,
1241 MF::Xpay(resid,
RT(-1.0),
b, 0, 0, ncomp,
IntVect(0));
1244template <
typename MF>
1248 if (crse_bcdata !=
nullptr) {
1253template <
typename MF>
1256 BCMode bc_mode,
const MF* crse_bcdata)
1258 BL_PROFILE(
"MLCellLinOp::correctionResidual()");
1259 const int ncomp = this->
getNComp();
1260 if (bc_mode == BCMode::Inhomogeneous)
1267 apply(amrlev, mglev, resid,
x, BCMode::Inhomogeneous, StateMode::Correction,
1273 apply(amrlev, mglev, resid,
x, BCMode::Homogeneous, StateMode::Correction,
nullptr);
1276 MF::Xpay(resid, Real(-1.0),
b, 0, 0, ncomp,
IntVect(0));
1279template <
typename MF>
1282 MF&, MF& fine_sol,
const MF&)
const
1289 const int ncomp = this->
getNComp();
1291 const int fine_amrlev = crse_amrlev+1;
1293 Real dt = Real(1.0);
1294 const Real* crse_dx = this->
m_geom[crse_amrlev][0].CellSize();
1295 const Real* fine_dx = this->
m_geom[fine_amrlev][0].CellSize();
1297 const int mglev = 0;
1298 applyBC(fine_amrlev, mglev, fine_sol, BCMode::Inhomogeneous, StateMode::Solution,
1305#pragma omp parallel if (Gpu::notInLaunchRegion())
1314 if (fluxreg.CrseHasWork(mfi))
1316 const Box& tbx = mfi.tilebox();
1321 Elixir elify = flux[1].elixir();,
1322 Elixir elifz = flux[2].elixir(););
1323 FFlux(crse_amrlev, mfi, pflux, crse_sol[mfi], Location::FaceCentroid);
1324 fluxreg.CrseAdd(mfi, cpflux, crse_dx, dt,
RunOn::Gpu);
1334 if (fluxreg.FineHasWork(mfi))
1336 const Box& tbx = mfi.tilebox();
1337 const int face_only =
true;
1342 Elixir elify = flux[1].elixir();,
1343 Elixir elifz = flux[2].elixir(););
1344 FFlux(fine_amrlev, mfi, pflux, fine_sol[mfi], Location::FaceCentroid, face_only);
1345 fluxreg.FineAdd(mfi, cpflux, fine_dx, dt,
RunOn::Gpu);
1350 fluxreg.Reflux(res);
1353template <
typename MF>
1360 const int mglev = 0;
1361 const int ncomp = this->
getNComp();
1362 applyBC(amrlev, mglev, sol, BCMode::Inhomogeneous, StateMode::Solution,
1369#pragma omp parallel if (Gpu::notInLaunchRegion())
1376 const Box& tbx = mfi.tilebox();
1381 Elixir elify = flux[1].elixir();,
1382 Elixir elifz = flux[2].elixir(););
1383 FFlux(amrlev, mfi, pflux, sol[mfi], loc);
1384 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
1385 const Box& nbx = mfi.nodaltilebox(idim);
1386 auto const& dst = fluxes[idim]->array(mfi);
1387 auto const& src = pflux[idim]->const_array();
1390 dst(i,j,k,n) = src(i,j,k,n);
1397template <
typename MF>
1404 if (sol.nComp() > 1) {
1405 amrex::Abort(
"MLCellLinOp::compGrad called, but only works for single-component solves");
1408 const int mglev = 0;
1409 applyBC(amrlev, mglev, sol, BCMode::Inhomogeneous, StateMode::Solution,
1412 const int ncomp = this->
getNComp();
1415 const RT dyi =
static_cast<RT>(this->
m_geom[amrlev][mglev].InvCellSize(1));,
1416 const RT dzi =
static_cast<RT>(this->
m_geom[amrlev][mglev].InvCellSize(2)););
1418#pragma omp parallel if (Gpu::notInLaunchRegion())
1423 const Box& ybx = mfi.nodaltilebox(1);,
1424 const Box& zbx = mfi.nodaltilebox(2););
1425 const auto& s = sol.array(mfi);
1427 const auto& gy = grad[1]->array(mfi);,
1428 const auto& gz = grad[2]->array(mfi););
1432 gx(i,j,k,n) = dxi*(s(i,j,k,n) - s(i-1,j,k,n));
1434#if (AMREX_SPACEDIM >= 2)
1437 gy(i,j,k,n) = dyi*(s(i,j,k,n) - s(i,j-1,k,n));
1440#if (AMREX_SPACEDIM == 3)
1443 gz(i,j,k,n) = dzi*(s(i,j,k,n) - s(i,j,k-1,n));
1451template <
typename MF>
1456#if (AMREX_SPACEDIM != 3)
1459 const int ncomp = rhs.nComp();
1461 bool cc = rhs.ixType().cellCentered(0);
1465 const RT probxlo =
static_cast<RT>(geom.
ProbLo(0));
1468#pragma omp parallel if (Gpu::notInLaunchRegion())
1472 const Box& tbx = mfi.tilebox();
1473 auto const& rhsarr = rhs.array(mfi);
1474#if (AMREX_SPACEDIM == 1)
1478 RT rc = probxlo + (i+
RT(0.5))*dx;
1479 rhsarr(i,j,k,n) *= rc*rc;
1484 RT re = probxlo + i*dx;
1485 rhsarr(i,j,k,n) *= re*re;
1488#elif (AMREX_SPACEDIM == 2)
1492 RT rc = probxlo + (i+
RT(0.5))*dx;
1493 rhsarr(i,j,k,n) *= rc;
1498 RT re = probxlo + i*dx;
1499 rhsarr(i,j,k,n) *= re;
1507template <
typename MF>
1512#if (AMREX_SPACEDIM != 3)
1515 const int ncomp = rhs.nComp();
1517 bool cc = rhs.ixType().cellCentered(0);
1521 const RT probxlo =
static_cast<RT>(geom.
ProbLo(0));
1524#pragma omp parallel if (Gpu::notInLaunchRegion())
1528 const Box& tbx = mfi.tilebox();
1529 auto const& rhsarr = rhs.array(mfi);
1530#if (AMREX_SPACEDIM == 1)
1534 RT rcinv =
RT(1.0)/(probxlo + (i+
RT(0.5))*dx);
1535 rhsarr(i,j,k,n) *= rcinv*rcinv;
1540 RT re = probxlo + i*dx;
1541 RT reinv = (re==
RT(0.0)) ?
RT(0.0) :
RT(1.)/re;
1542 rhsarr(i,j,k,n) *= reinv*reinv;
1545#elif (AMREX_SPACEDIM == 2)
1549 RT rcinv =
RT(1.0)/(probxlo + (i+
RT(0.5))*dx);
1550 rhsarr(i,j,k,n) *= rcinv;
1555 RT re = probxlo + i*dx;
1556 RT reinv = (re==
RT(0.0)) ?
RT(0.0) :
RT(1.)/re;
1557 rhsarr(i,j,k,n) *= reinv;
1565template <
typename MF>
1572 const int ncomp = this->
getNComp();
1577 if (factory && !factory->isAllRegular())
1579 if constexpr (std::is_same<MF,MultiFab>()) {
1580 const MultiFab& vfrac = factory->getVolFrac();
1581 for (
int c = 0; c < ncomp; ++c) {
1586 amrex::Abort(
"TODO: MLMG with EB only works with MultiFab");
1592 for (
int c = 0; c < ncomp; ++c) {
1602template <
typename MF>
1607 const int ncomp = this->
getNComp();
1608 for (
int c = 0; c < ncomp; ++c) {
1609 rhs.plus(-
offset[c], c, 1);
1612 if (!rhs.isAllRegular()) {
1613 if constexpr (std::is_same<MF,MultiFab>()) {
1616 amrex::Abort(
"amrex::EB_set_covered only works with MultiFab");
1622template <
typename MF>
1626 BL_PROFILE(
"MLCellLinOp::prepareForSolve()");
1628 const int imaxorder = this->
maxorder;
1629 const int ncomp = this->
getNComp();
1635 const auto& bcondloc = *
m_bcondloc[amrlev][mglev];
1636 const auto& maskvals =
m_maskvals[amrlev][mglev];
1638 const RT dxi =
static_cast<RT>(this->
m_geom[amrlev][mglev].InvCellSize(0));
1639 const RT dyi =
static_cast<RT>((AMREX_SPACEDIM >= 2) ? this->
m_geom[amrlev][mglev].InvCellSize(1) : Real(1.0));
1640 const RT dzi =
static_cast<RT>((AMREX_SPACEDIM == 3) ? this->
m_geom[amrlev][mglev].InvCellSize(2) : Real(1.0));
1642 auto& undrrelxr = this->
m_undrrelxr[amrlev][mglev];
1643 MF foo(this->
m_grids[amrlev][mglev], this->
m_dmap[amrlev][mglev], ncomp, 0,
MFInfo().SetAlloc(
false));
1648 (factory) ? &(factory->getMultiEBCellFlagFab()) :
nullptr;
1649 auto area = (factory) ? factory->getAreaFrac()
1657 if (factory && !factory->isAllRegular()) {
1658#if defined(AMREX_USE_CUDA) && defined(_WIN32)
1659 if (!std::is_same<MF,MultiFab>()) {
1661 if constexpr (!std::is_same<MF,MultiFab>()) {
1663 amrex::Abort(
"MLCellLinOp with EB only works with MultiFab");
1666 tags.reserve(foo.local_size()*AMREX_SPACEDIM*ncomp);
1670 const Box& vbx = mfi.validbox();
1672 const auto & bdlv = bcondloc.bndryLocs(mfi);
1673 const auto & bdcv = bcondloc.bndryConds(mfi);
1677 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim)
1684 for (
int icomp = 0; icomp < ncomp; ++icomp) {
1685 tags.emplace_back(MLMGPSEBTag<RT>{undrrelxr[olo].array(mfi),
1686 undrrelxr[ohi].array(mfi),
1688 maskvals[olo].const_array(mfi),
1689 maskvals[ohi].const_array(mfi),
1690 bdlv[icomp][olo], bdlv[icomp][ohi],
1692 bdcv[icomp][olo], bdcv[icomp][ohi],
1693 vbx.
length(idim), icomp, idim});
1700 [=]
AMREX_GPU_DEVICE (
int i,
int j,
int k, MLMGPSEBTag<RT>
const& tag)
noexcept
1705 mllinop_comp_interp_coef0_x_eb
1706 (0, i , j, k, tag.blen, tag.flo, tag.mlo, tag.ap,
1707 tag.bctlo, tag.bcllo, imaxorder, dxi, tag.comp);
1708 mllinop_comp_interp_coef0_x_eb
1709 (1, i+tag.blen+1, j, k, tag.blen, tag.fhi, tag.mhi, tag.ap,
1710 tag.bcthi, tag.bclhi, imaxorder, dxi, tag.comp);
1712#
if (AMREX_SPACEDIM > 1)
1714#
if (AMREX_SPACEDIM > 2)
1718 mllinop_comp_interp_coef0_y_eb
1719 (0, i, j , k, tag.blen, tag.flo, tag.mlo, tag.ap,
1720 tag.bctlo, tag.bcllo, imaxorder, dyi, tag.comp);
1721 mllinop_comp_interp_coef0_y_eb
1722 (1, i, j+tag.blen+1, k, tag.blen, tag.fhi, tag.mhi, tag.ap,
1723 tag.bcthi, tag.bclhi, imaxorder, dyi, tag.comp);
1725#
if (AMREX_SPACEDIM > 2)
1727 mllinop_comp_interp_coef0_z_eb
1728 (0, i, j, k , tag.blen, tag.flo, tag.mlo, tag.ap,
1729 tag.bctlo, tag.bcllo, imaxorder, dzi, tag.comp);
1730 mllinop_comp_interp_coef0_z_eb
1731 (1, i, j, k+tag.blen+1, tag.blen, tag.fhi, tag.mhi, tag.ap,
1732 tag.bcthi, tag.bclhi, imaxorder, dzi, tag.comp);
1740 (0, i , j, k, tag.blen, tag.flo, tag.mlo,
1741 tag.bctlo, tag.bcllo, imaxorder, dxi, tag.comp);
1743 (1, i+tag.blen+1, j, k, tag.blen, tag.fhi, tag.mhi,
1744 tag.bcthi, tag.bclhi, imaxorder, dxi, tag.comp);
1746#if (AMREX_SPACEDIM > 1)
1748#if (AMREX_SPACEDIM > 2)
1753 (0, i, j , k, tag.blen, tag.flo, tag.mlo,
1754 tag.bctlo, tag.bcllo, imaxorder, dyi, tag.comp);
1756 (1, i, j+tag.blen+1, k, tag.blen, tag.fhi, tag.mhi,
1757 tag.bcthi, tag.bclhi, imaxorder, dyi, tag.comp);
1759#if (AMREX_SPACEDIM > 2)
1762 (0, i, j, k , tag.blen, tag.flo, tag.mlo,
1763 tag.bctlo, tag.bcllo, imaxorder, dzi, tag.comp);
1765 (1, i, j, k+tag.blen+1, tag.blen, tag.fhi, tag.mhi,
1766 tag.bcthi, tag.bclhi, imaxorder, dzi, tag.comp);
1777 tags.reserve(foo.local_size()*AMREX_SPACEDIM*ncomp);
1781 const Box& vbx = mfi.validbox();
1783 const auto & bdlv = bcondloc.bndryLocs(mfi);
1784 const auto & bdcv = bcondloc.bndryConds(mfi);
1786 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim)
1788 if (idim != hidden_direction) {
1791 for (
int icomp = 0; icomp < ncomp; ++icomp) {
1793 undrrelxr[ohi].array(mfi),
1794 maskvals[olo].const_array(mfi),
1795 maskvals[ohi].const_array(mfi),
1796 bdlv[icomp][olo], bdlv[icomp][ohi],
1798 bdcv[icomp][olo], bdcv[icomp][ohi],
1799 vbx.
length(idim), icomp, idim});
1810 mllinop_comp_interp_coef0_x
1811 (0, i , j, k, tag.blen, tag.flo, tag.mlo,
1812 tag.bctlo, tag.bcllo, imaxorder, dxi, tag.comp);
1813 mllinop_comp_interp_coef0_x
1814 (1, i+tag.blen+1, j, k, tag.blen, tag.fhi, tag.mhi,
1815 tag.bcthi, tag.bclhi, imaxorder, dxi, tag.comp);
1817#
if (AMREX_SPACEDIM > 1)
1819#
if (AMREX_SPACEDIM > 2)
1823 mllinop_comp_interp_coef0_y
1824 (0, i, j , k, tag.blen, tag.flo, tag.mlo,
1825 tag.bctlo, tag.bcllo, imaxorder, dyi, tag.comp);
1826 mllinop_comp_interp_coef0_y
1827 (1, i, j+tag.blen+1, k, tag.blen, tag.fhi, tag.mhi,
1828 tag.bcthi, tag.bclhi, imaxorder, dyi, tag.comp);
1830#
if (AMREX_SPACEDIM > 2)
1851 const Box& vbx = mfi.validbox();
1853 const auto & bdlv = bcondloc.bndryLocs(mfi);
1854 const auto & bdcv = bcondloc.bndryConds(mfi);
1859 for (
int idim = 0; idim < AMREX_SPACEDIM; ++idim)
1861 if (idim == hidden_direction) {
continue; }
1866 const int blen = vbx.
length(idim);
1867 const auto& mlo = maskvals[olo].array(mfi);
1868 const auto& mhi = maskvals[ohi].array(mfi);
1869 const auto& flo = undrrelxr[olo].array(mfi);
1870 const auto& fhi = undrrelxr[ohi].array(mfi);
1871 for (
int icomp = 0; icomp < ncomp; ++icomp) {
1872 const BoundCond bctlo = bdcv[icomp][olo];
1873 const BoundCond bcthi = bdcv[icomp][ohi];
1874 const auto bcllo = bdlv[icomp][olo];
1875 const auto bclhi = bdlv[icomp][ohi];
1878 if constexpr (!std::is_same<MF,MultiFab>()) {
1879 amrex::Abort(
"MLCellLinOp with EB only works with MultiFab");
1881 auto const& ap = area[idim]->const_array(mfi);
1883 mllinop_comp_interp_coef0_x_eb
1884 (0, blo, blen, flo, mlo, ap, bctlo, bcllo,
1885 imaxorder, dxi, icomp);
1886 mllinop_comp_interp_coef0_x_eb
1887 (1, bhi, blen, fhi, mhi, ap, bcthi, bclhi,
1888 imaxorder, dxi, icomp);
1889 }
else if (idim == 1) {
1890 mllinop_comp_interp_coef0_y_eb
1891 (0, blo, blen, flo, mlo, ap, bctlo, bcllo,
1892 imaxorder, dyi, icomp);
1893 mllinop_comp_interp_coef0_y_eb
1894 (1, bhi, blen, fhi, mhi, ap, bcthi, bclhi,
1895 imaxorder, dyi, icomp);
1897 mllinop_comp_interp_coef0_z_eb
1898 (0, blo, blen, flo, mlo, ap, bctlo, bcllo,
1899 imaxorder, dzi, icomp);
1900 mllinop_comp_interp_coef0_z_eb
1901 (1, bhi, blen, fhi, mhi, ap, bcthi, bclhi,
1902 imaxorder, dzi, icomp);
1910 (0, blo, blen, flo, mlo, bctlo, bcllo,
1911 imaxorder, dxi, icomp);
1913 (1, bhi, blen, fhi, mhi, bcthi, bclhi,
1914 imaxorder, dxi, icomp);
1915 }
else if (idim == 1) {
1917 (0, blo, blen, flo, mlo, bctlo, bcllo,
1918 imaxorder, dyi, icomp);
1920 (1, bhi, blen, fhi, mhi, bcthi, bclhi,
1921 imaxorder, dyi, icomp);
1924 (0, blo, blen, flo, mlo, bctlo, bcllo,
1925 imaxorder, dzi, icomp);
1927 (1, bhi, blen, fhi, mhi, bcthi, bclhi,
1928 imaxorder, dzi, icomp);
1939template <
typename MF>
1944 const int ncomp = this->
getNComp();
1953template <
typename MF>
1957 const int ncomp = this->
getNComp();
1960 for (
int ilev = 0; ilev < this->
NAMRLevels()-1; ++ilev) {
1964 *y[this->
NAMRLevels()-1], 0, ncomp, nghost,
true);
1969template <
typename MF>
1973 const int ncomp = this->
getNComp();
1976 for (
int ilev = 0; ilev < this->
NAMRLevels()-1; ++ilev) {
1981 return std::sqrt(result);
1984template <
typename MF>
1997 auto f = [&] (
int amrlev,
int mglev) {
2000 if (factory && !factory->isAllRegular())
2002 if constexpr (std::is_same<MF,MultiFab>()) {
2006 amrex::Abort(
"MLCellLinOp with EB only works with MultiFab");
2029 if (factory && !factory->isAllRegular())
2046template <
typename MF>
2050 const int ncomp = this->
getNComp();
2051 const int finest_level = this->
NAMRLevels() - 1;
2055 if (factory && !factory->isAllRegular()) {
2056#if defined(AMREX_USE_CUDA) && defined(_WIN32)
2057 if (!std::is_same<MF,MultiFab>()) {
2059 if constexpr (!std::is_same<MF,MultiFab>()) {
2061 amrex::Abort(
"MLCellLinOpT with EB only works with MultiFab");
2063 const MultiFab& vfrac = factory->getVolFrac();
2064 if (amrlev == finest_level) {
2074 return std::abs(ma[box_no](i,j,k,n)
2075 *vfrac_ma[box_no](i,j,k));
2081#pragma omp parallel reduction(max:norm)
2084 Box const& bx = mfi.tilebox();
2085 auto const& fab = mf.const_array(mfi);
2089 norm = std::max(
norm, std::abs(fab(i,j,k,n)*v(i,j,k)));
2096 auto const& ma = mf.const_arrays();
2104 if (mask_ma[box_no](i,j,k)) {
2105 return std::abs(ma[box_no](i,j,k,n)
2106 *vfrac_ma[box_no](i,j,k));
2115#pragma omp parallel reduction(max:norm)
2118 Box const& bx = mfi.tilebox();
2119 auto const& fab = mf.const_array(mfi);
2125 norm = std::max(
norm, std::abs(fab(i,j,k,n)*v(i,j,k)));
2135 if (amrlev == finest_level) {
2146template <
typename MF>
2151 for (
int falev = this->
NAMRLevels()-1; falev > 0; --falev)
2154 if (!sol[falev].isAllRegular()) {
2155 if constexpr (std::is_same<MF,MultiFab>()) {
2158 amrex::Abort(
"EB_average_down only works with MultiFab");
2168template <
typename MF>
2173 if (!fres.isAllRegular()) {
2174 if constexpr (std::is_same<MF,MultiFab>()) {
2178 amrex::Abort(
"EB_average_down only works with MultiFab");
2188template <
typename MF>
2196 const int ncomp = this->
getNComp();
2202 this->
m_geom[amrlev][0]);
2216template <
typename MF>
#define BL_PROFILE(a)
Definition AMReX_BLProfiler.H:551
#define AMREX_ALWAYS_ASSERT_WITH_MESSAGE(EX, MSG)
Definition AMReX_BLassert.H:49
#define BL_ASSERT(EX)
Definition AMReX_BLassert.H:39
#define AMREX_ASSERT(EX)
Definition AMReX_BLassert.H:38
#define AMREX_ALWAYS_ASSERT(EX)
Definition AMReX_BLassert.H:50
#define AMREX_FORCE_INLINE
Definition AMReX_Extension.H:119
#define AMREX_HOST_DEVICE_PARALLEL_FOR_3D(...)
Definition AMReX_GpuLaunchMacrosC.nolint.H:110
#define AMREX_GPU_LAUNCH_HOST_DEVICE_LAMBDA_RANGE(TN, TI, block)
Definition AMReX_GpuLaunchMacrosC.nolint.H:4
#define AMREX_HOST_DEVICE_PARALLEL_FOR_4D(...)
Definition AMReX_GpuLaunchMacrosC.nolint.H:111
#define AMREX_GPU_DEVICE
Definition AMReX_GpuQualifiers.H:18
#define AMREX_GPU_HOST_DEVICE
Definition AMReX_GpuQualifiers.H:20
Array4< int const > offset
Definition AMReX_HypreMLABecLap.cpp:1089
Box cbx
Definition AMReX_HypreMLABecLap.cpp:1091
Array4< Real > fine
Definition AMReX_InterpFaceRegister.cpp:90
Array4< int const > mask
Definition AMReX_InterpFaceRegister.cpp:93
Array4< Real const > crse
Definition AMReX_InterpFaceRegister.cpp:92
#define AMREX_LOOP_4D(bx, ncomp, i, j, k, n, block)
Definition AMReX_Loop.nolint.H:16
void amrex_mllinop_apply_bc(const int *lo, const int *hi, amrex_real *phi, const int *philo, const int *phihi, const int *mask, const int *mlo, const int *mhi, int cdir, int bct, amrex_real bcl, const amrex_real *bcval, const int *blo, const int *bhi, int maxorder, const amrex_real *dxinv, int inhomog, int nc, int cross)
#define AMREX_D_TERM(a, b, c)
Definition AMReX_SPACE.H:129
if(!(yy_init))
Definition amrex_iparser.lex.nolint.H:935
const FabSetT< MF > & bndryValues(Orientation face) const noexcept
Return FabSet on given face.
Definition AMReX_BndryData.H:77
Maintain an identifier for boundary condition types.
Definition AMReX_BoundCond.H:20
A collection of Boxes stored in an Array.
Definition AMReX_BoxArray.H:550
void define(const Box &bx)
Initialize the BoxArray from a single box. It is an error if the BoxArray has already been initialize...
BoxArray & refine(int refinement_ratio)
Refine each Box in the BoxArray to the specified ratio.
BoxArray & coarsen(int refinement_ratio)
Coarsen each Box in the BoxArray to the specified ratio.
Long size() const noexcept
Return the number of boxes in the BoxArray.
Definition AMReX_BoxArray.H:597
A class for managing a List of Boxes that share a common IndexType. This class implements operations ...
Definition AMReX_BoxList.H:52
static AMREX_GPU_HOST_DEVICE BoxND TheUnitBox() noexcept
This static member function returns a constant reference to an object of type BoxND representing the ...
Definition AMReX_Box.H:737
AMREX_GPU_HOST_DEVICE double d_numPts() const noexcept
Returns the number of points contained in the BoxND. This is intended for use only in diagnostic mess...
Definition AMReX_Box.H:366
AMREX_GPU_HOST_DEVICE IntVectND< dim > length() const noexcept
Return the length of the BoxND.
Definition AMReX_Box.H:146
AMREX_GPU_HOST_DEVICE bool contains(const IntVectND< dim > &p) const noexcept
Returns true if argument is contained within BoxND.
Definition AMReX_Box.H:204
const Real * CellSize() const noexcept
Returns the cellsize for each coordinate direction.
Definition AMReX_CoordSys.H:71
Calculates the distribution of FABs to MPI processes.
Definition AMReX_DistributionMapping.H:41
Definition AMReX_EBFabFactory.H:24
const MultiFab & getVolFrac() const noexcept
Definition AMReX_EBFabFactory.H:55
An Array of FortranArrayBox(FAB)-like Objects.
Definition AMReX_FabArray.H:344
Array4< typename FabArray< FAB >::value_type const > const_array(const MFIter &mfi) const noexcept
Definition AMReX_FabArray.H:584
MultiArray4< typename FabArray< FAB >::value_type const > const_arrays() const noexcept
Definition AMReX_FabArray.H:646
Definition AMReX_FabFactory.H:50
Rectangular problem domain geometry.
Definition AMReX_Geometry.H:73
const Box & Domain() const noexcept
Returns our rectangular domain.
Definition AMReX_Geometry.H:210
Periodicity periodicity() const noexcept
Definition AMReX_Geometry.H:355
GpuArray< int, AMREX_SPACEDIM > isPeriodicArray() const noexcept
Definition AMReX_Geometry.H:347
const Real * ProbLo() const noexcept
Returns the lo end of the problem domain in each dimension.
Definition AMReX_Geometry.H:178
Definition AMReX_Tuple.H:93
Definition AMReX_GpuElixir.H:13
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE bool allGT(const IntVectND< dim > &rhs) const noexcept
Returns true if this is greater than argument for all components. NOTE: This is NOT a strict weak ord...
Definition AMReX_IntVect.H:418
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE constexpr std::size_t size() noexcept
Definition AMReX_IntVect.H:725
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE IntVectND & setVal(int i, int val) noexcept
Set i'th coordinate of IntVectND to val.
Definition AMReX_IntVect.H:274
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE constexpr IntVectND< dim > TheCellVector() noexcept
This static member function returns a reference to a constant IntVectND object, all of whose dim argu...
Definition AMReX_IntVect.H:711
An InterpBndryData object adds to a BndryData object the ability to manipulate and set the data store...
Definition AMReX_InterpBndryData.H:40
a one-thingy-per-box distributed object
Definition AMReX_LayoutData.H:13
Definition AMReX_MFIter.H:57
bool isValid() const noexcept
Is the iterator valid i.e. is it associated with a FAB?
Definition AMReX_MFIter.H:141
Definition AMReX_MLCellLinOp.H:167
const Vector< BCTuple > & bndryConds(const MFIter &mfi) const noexcept
Definition AMReX_MLCellLinOp.H:179
GpuArray< BCTL, 2 *AMREX_SPACEDIM > const * getBCTLPtr(const MFIter &mfi) const noexcept
Definition AMReX_MLCellLinOp.H:191
const BCTuple & bndryConds(const MFIter &mfi, int icomp) const noexcept
Definition AMReX_MLCellLinOp.H:185
Gpu::DeviceVector< GpuArray< BCTL, 2 *AMREX_SPACEDIM > > bctl_dv
Definition AMReX_MLCellLinOp.H:198
BndryCondLoc(const BoxArray &ba, const DistributionMapping &dm, int ncomp)
Definition AMReX_MLCellLinOp.H:288
LayoutData< GpuArray< BCTL, 2 *AMREX_SPACEDIM > * > bctl
Definition AMReX_MLCellLinOp.H:197
const Vector< RealTuple > & bndryLocs(const MFIter &mfi) const noexcept
Definition AMReX_MLCellLinOp.H:182
int m_ncomp
Definition AMReX_MLCellLinOp.H:199
void setLOBndryConds(const Geometry &geom, const Real *dx, const Vector< Array< BCType, AMREX_SPACEDIM > > &lobc, const Vector< Array< BCType, AMREX_SPACEDIM > > &hibc, IntVect const &ratio, const RealVect &interior_bloc, const Array< Real, AMREX_SPACEDIM > &domain_bloc_lo, const Array< Real, AMREX_SPACEDIM > &domain_bloc_hi, LinOpBCType crse_fine_bc_type)
Definition AMReX_MLCellLinOp.H:309
LayoutData< Vector< RealTuple > > bcloc
Definition AMReX_MLCellLinOp.H:196
LayoutData< Vector< BCTuple > > bcond
Definition AMReX_MLCellLinOp.H:195
const RealTuple & bndryLocs(const MFIter &mfi, int icomp) const noexcept
Definition AMReX_MLCellLinOp.H:188
Definition AMReX_MLCellLinOp.H:21
void smooth(int amrlev, int mglev, MF &sol, const MF &rhs, bool skip_fillboundary=false) const final
Smooth.
Definition AMReX_MLCellLinOp.H:1213
virtual void Fsmooth(int amrlev, int mglev, MF &sol, const MF &rhs, int redblack) const =0
Vector< RT > getSolvabilityOffset(int amrlev, int mglev, MF const &rhs) const override
get offset for fixing solvability
Definition AMReX_MLCellLinOp.H:1567
typename MF::fab_type FAB
Definition AMReX_MLCellLinOp.H:24
void averageDownSolutionRHS(int camrlev, MF &crse_sol, MF &crse_rhs, const MF &fine_sol, const MF &fine_rhs) override
Average-down data from fine AMR level to coarse AMR level.
Definition AMReX_MLCellLinOp.H:1192
void averageDownAndSync(Vector< MF > &sol) const override
Definition AMReX_MLCellLinOp.H:2148
BoxArray makeNGrids(int grid_size) const
Definition AMReX_MLCellLinOp.H:904
Vector< YAFluxRegisterT< MF > > m_fluxreg
Definition AMReX_MLCellLinOp.H:211
virtual void applyBC(int amrlev, int mglev, MF &in, BCMode bc_mode, StateMode s_mode, const MLMGBndryT< MF > *bndry=nullptr, bool skip_fillboundary=false) const
Definition AMReX_MLCellLinOp.H:691
void updateSolBC(int amrlev, const MF &crse_bcdata) const
Definition AMReX_MLCellLinOp.H:660
Vector< std::unique_ptr< MLMGBndryT< MF > > > m_bndry_sol
Definition AMReX_MLCellLinOp.H:153
MLCellLinOpT< MF > & operator=(const MLCellLinOpT< MF > &)=delete
void compGrad(int amrlev, const Array< MF *, AMREX_SPACEDIM > &grad, MF &sol, Location loc) const override
Compute gradients of the solution.
Definition AMReX_MLCellLinOp.H:1399
void avgDownResAmr(int clev, MF &cres, MF const &fres) const override
Definition AMReX_MLCellLinOp.H:2170
void reflux(int crse_amrlev, MF &res, const MF &crse_sol, const MF &, MF &, MF &fine_sol, const MF &) const final
Reflux at AMR coarse/fine boundary.
Definition AMReX_MLCellLinOp.H:1281
RT dotProductPrecond(Vector< MF const * > const &x, Vector< MF const * > const &y) const final
Definition AMReX_MLCellLinOp.H:1955
virtual void Fapply(int amrlev, int mglev, MF &out, const MF &in) const =0
void defineBC()
Definition AMReX_MLCellLinOp.H:438
void endPrecondBC() override
Definition AMReX_MLCellLinOp.H:2218
void update() override
Update for reuse.
Definition AMReX_MLCellLinOp.H:653
typename MF::value_type RT
Definition AMReX_MLCellLinOp.H:25
LinOpBCType BCType
Definition AMReX_MLCellLinOp.H:27
void correctionResidual(int amrlev, int mglev, MF &resid, MF &x, const MF &b, BCMode bc_mode, const MF *crse_bcdata=nullptr) final
Compute residual for the residual-correction form, resid = b - L(x)
Definition AMReX_MLCellLinOp.H:1255
MLCellLinOpT(const MLCellLinOpT< MF > &)=delete
typename MLLinOpT< MF >::BCMode BCMode
Definition AMReX_MLCellLinOp.H:28
RT normInf(int amrlev, MF const &mf, bool local) const override
Definition AMReX_MLCellLinOp.H:2048
Array< RT, 2 *BL_SPACEDIM > RealTuple
Definition AMReX_MLCellLinOp.H:164
virtual bool isCrossStencil() const
Definition AMReX_MLCellLinOp.H:60
void updateCorBC(int amrlev, const MF &crse_bcdata) const
Definition AMReX_MLCellLinOp.H:676
Array< BoundCond, 2 *BL_SPACEDIM > BCTuple
Definition AMReX_MLCellLinOp.H:165
void interpolation(int amrlev, int fmglev, MF &fine, const MF &crse) const override
Add interpolated coarse MG level data to fine MG level data.
Definition AMReX_MLCellLinOp.H:963
virtual void addInhomogNeumannFlux(int, const Array< MF *, AMREX_SPACEDIM > &, MF const &, bool) const
Definition AMReX_MLCellLinOp.H:126
RT xdoty(int amrlev, int mglev, const MF &x, const MF &y, bool local) const final
x dot y, used by the bottom solver
Definition AMReX_MLCellLinOp.H:1941
void prepareForSolve() override
Definition AMReX_MLCellLinOp.H:1624
void unapplyMetricTerm(int amrlev, int mglev, MF &rhs) const final
unapply metric terms if there are any
Definition AMReX_MLCellLinOp.H:1509
void apply(int amrlev, int mglev, MF &out, MF &in, BCMode bc_mode, StateMode s_mode, const MLMGBndryT< MF > *bndry=nullptr) const override
Apply the linear operator, out = L(in)
Definition AMReX_MLCellLinOp.H:1203
void applyMetricTerm(int amrlev, int mglev, MF &rhs) const final
apply metric terms if there are any
Definition AMReX_MLCellLinOp.H:1453
Vector< std::unique_ptr< BndryRegisterT< MF > > > m_crse_cor_br
Definition AMReX_MLCellLinOp.H:157
RT norm2Precond(Vector< MF const * > const &x) const final
Definition AMReX_MLCellLinOp.H:1971
void beginPrecondBC() override
Definition AMReX_MLCellLinOp.H:2190
Vector< Vector< std::unique_ptr< BndryCondLoc > > > m_bcondloc
Definition AMReX_MLCellLinOp.H:201
void setGaussSeidel(bool flag) noexcept
Definition AMReX_MLCellLinOp.H:58
void defineAuxData()
Definition AMReX_MLCellLinOp.H:375
Vector< std::unique_ptr< MF > > m_robin_bcval
Definition AMReX_MLCellLinOp.H:145
Vector< Vector< Array< MultiMask, 2 *AMREX_SPACEDIM > > > m_maskvals
Definition AMReX_MLCellLinOp.H:207
void define(const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info=LPInfo(), const Vector< FabFactory< FAB > const * > &a_factory={})
Definition AMReX_MLCellLinOp.H:362
~MLCellLinOpT() override=default
virtual bool isTensorOp() const
Definition AMReX_MLCellLinOp.H:61
void setInterpBndryHalfWidth(int w)
Definition AMReX_MLCellLinOp.H:147
void restriction(int, int, MF &crse, MF &fine) const override
Restriction onto coarse MG level.
Definition AMReX_MLCellLinOp.H:954
Vector< Vector< BndryRegisterT< MF > > > m_undrrelxr
Definition AMReX_MLCellLinOp.H:204
int m_interpbndry_halfwidth
Definition AMReX_MLCellLinOp.H:223
MLCellLinOpT()
Definition AMReX_MLCellLinOp.H:355
MLCellLinOpT(MLCellLinOpT< MF > &&)=delete
void interpAssign(int amrlev, int fmglev, MF &fine, MF &crse) const override
Overwrite fine MG level data with interpolated coarse data.
Definition AMReX_MLCellLinOp.H:1010
void fixSolvabilityByOffset(int amrlev, int mglev, MF &rhs, Vector< RT > const &offset) const override
Definition AMReX_MLCellLinOp.H:1604
Vector< std::unique_ptr< MLMGBndryT< MF > > > m_bndry_cor
Definition AMReX_MLCellLinOp.H:156
typename MLLinOpT< MF >::Location Location
Definition AMReX_MLCellLinOp.H:30
void prepareForFluxes(int amrlev, const MF *crse_bcdata=nullptr) override
Definition AMReX_MLCellLinOp.H:1246
void solutionResidual(int amrlev, MF &resid, MF &x, const MF &b, const MF *crse_bcdata=nullptr) override
Compute residual for solution.
Definition AMReX_MLCellLinOp.H:1228
Vector< Vector< RT > > m_volinv
Definition AMReX_MLCellLinOp.H:221
bool m_has_metric_term
Definition AMReX_MLCellLinOp.H:151
bool m_use_gauss_seidel
Definition AMReX_MLCellLinOp.H:213
Vector< std::unique_ptr< MLMGBndryT< MF > > > m_bndry_sol_zero
Definition AMReX_MLCellLinOp.H:159
Vector< std::unique_ptr< BndryRegisterT< MF > > > m_crse_sol_br
Definition AMReX_MLCellLinOp.H:154
void interpolationAmr(int famrlev, MF &fine, const MF &crse, IntVect const &nghost) const override
Interpolation between AMR levels.
Definition AMReX_MLCellLinOp.H:1103
void compFlux(int amrlev, const Array< MF *, AMREX_SPACEDIM > &fluxes, MF &sol, Location loc) const override
Compute fluxes.
Definition AMReX_MLCellLinOp.H:1355
void computeVolInv() const
Definition AMReX_MLCellLinOp.H:1986
Vector< std::unique_ptr< iMultiFab > > m_norm_fine_mask
Definition AMReX_MLCellLinOp.H:209
virtual void FFlux(int amrlev, const MFIter &mfi, const Array< FAB *, AMREX_SPACEDIM > &flux, const FAB &sol, Location loc, int face_only=0) const =0
typename MLLinOpT< MF >::StateMode StateMode
Definition AMReX_MLCellLinOp.H:29
void setLevelBC(int amrlev, const MF *levelbcdata, const MF *robinbc_a=nullptr, const MF *robinbc_b=nullptr, const MF *robinbc_f=nullptr) final
Set boundary conditions for given level. For cell-centered solves only.
Definition AMReX_MLCellLinOp.H:520
bool needsUpdate() const override
Does it need update if it's reused?
Definition AMReX_MLCellLinOp.H:53
Definition AMReX_MLLinOp.H:98
const MF * m_coarse_data_for_bc
Definition AMReX_MLLinOp.H:637
Vector< Vector< std::unique_ptr< FabFactory< FAB > > > > m_factory
Definition AMReX_MLLinOp.H:610
int NAMRLevels() const noexcept
Return the number of AMR levels.
Definition AMReX_MLLinOp.H:563
IntVect m_ixtype
Definition AMReX_MLLinOp.H:598
Array< Real, AMREX_SPACEDIM > m_domain_bloc_hi
Definition AMReX_MLLinOp.H:631
RealVect m_coarse_bc_loc
Definition AMReX_MLLinOp.H:636
virtual bool needsUpdate() const
Does it need update if it's reused?
Definition AMReX_MLLinOp.H:266
FabFactory< FAB > const * Factory(int amr_lev, int mglev=0) const noexcept
Definition AMReX_MLLinOp.H:648
Vector< Array< BCType, AMREX_SPACEDIM > > m_hibc_orig
Definition AMReX_MLLinOp.H:576
bool hasRobinBC() const noexcept
Definition AMReX_MLLinOp.H:1276
Array< Real, AMREX_SPACEDIM > m_domain_bloc_lo
Definition AMReX_MLLinOp.H:630
Vector< Vector< BoxArray > > m_grids
Definition AMReX_MLLinOp.H:608
virtual int getNComp() const
Return number of components.
Definition AMReX_MLLinOp.H:261
Vector< int > m_amr_ref_ratio
Definition AMReX_MLLinOp.H:593
Vector< int > m_num_mg_levels
Definition AMReX_MLLinOp.H:595
Vector< Vector< DistributionMapping > > m_dmap
Definition AMReX_MLLinOp.H:609
IntVect m_coarse_data_crse_ratio
Definition AMReX_MLLinOp.H:635
const Vector< int > & AMRRefRatio() const noexcept
Return AMR refinement ratios.
Definition AMReX_MLLinOp.H:643
Vector< Array< BCType, AMREX_SPACEDIM > > m_lobc_orig
Definition AMReX_MLLinOp.H:575
bool needsCoarseDataForBC() const noexcept
Needs coarse data for bc?
Definition AMReX_MLLinOp.H:177
virtual void update()
Update for reuse.
Definition AMReX_MLLinOp.H:268
bool m_precond_mode
Definition AMReX_MLLinOp.H:640
bool hasHiddenDimension() const noexcept
Definition AMReX_MLLinOp.H:706
const Geometry & Geom(int amr_lev, int mglev=0) const noexcept
Definition AMReX_MLLinOp.H:568
Vector< Array< BCType, AMREX_SPACEDIM > > m_lobc
Definition AMReX_MLLinOp.H:571
int hiddenDirection() const noexcept
Definition AMReX_MLLinOp.H:707
Vector< Vector< Geometry > > m_geom
first Vector is for amr level and second is mg level
Definition AMReX_MLLinOp.H:607
Box compactify(Box const &b) const noexcept
Definition AMReX_MLLinOp.H:1283
int maxorder
Definition AMReX_MLLinOp.H:588
void define(const Vector< Geometry > &a_geom, const Vector< BoxArray > &a_grids, const Vector< DistributionMapping > &a_dmap, const LPInfo &a_info, const Vector< FabFactory< FAB > const * > &a_factory, bool eb_limit_coarsening=true)
Definition AMReX_MLLinOp.H:768
Vector< IntVect > mg_coarsen_ratio_vec
Definition AMReX_MLLinOp.H:604
LPInfo info
Definition AMReX_MLLinOp.H:584
Vector< Array< BCType, AMREX_SPACEDIM > > m_hibc
Definition AMReX_MLLinOp.H:572
int NMGLevels(int amrlev) const noexcept
Return the number of MG levels at given AMR level.
Definition AMReX_MLLinOp.H:566
LinOpBCType m_coarse_fine_bc_type
Definition AMReX_MLLinOp.H:634
int m_num_amr_levels
Definition AMReX_MLLinOp.H:592
Definition AMReX_MLMGBndry.H:12
static void setBoxBC(RealTuple &bloc, BCTuple &bctag, const Box &bx, const Box &domain, const Array< LinOpBCType, AMREX_SPACEDIM > &lo, const Array< LinOpBCType, AMREX_SPACEDIM > &hi, const Real *dx, IntVect const &ratio, const RealVect &interior_bloc, const Array< Real, AMREX_SPACEDIM > &domain_bloc_lo, const Array< Real, AMREX_SPACEDIM > &domain_bloc_hi, const GpuArray< int, AMREX_SPACEDIM > &is_periodic, LinOpBCType a_crse_fine_bc_type)
Definition AMReX_MLMGBndry.H:107
Definition AMReX_Mask.H:28
A collection (stored as an array) of FArrayBox objects.
Definition AMReX_MultiFab.H:38
Real sum(int comp=0, bool local=false) const
Returns the sum of component "comp" over the MultiFab – no ghost cells are included.
Definition AMReX_MultiFab.cpp:1266
An Iterator over the Orientation of Faces of a Box.
Definition AMReX_Orientation.H:135
Encapsulation of the Orientation of the Faces of a Box.
Definition AMReX_Orientation.H:29
@ low
Definition AMReX_Orientation.H:34
@ high
Definition AMReX_Orientation.H:34
Definition AMReX_PODVector.H:262
void reserve(size_type a_capacity)
Definition AMReX_PODVector.H:663
iterator begin() noexcept
Definition AMReX_PODVector.H:617
iterator end() noexcept
Definition AMReX_PODVector.H:621
void push_back(const T &a_value)
Definition AMReX_PODVector.H:572
A Real vector in SpaceDim-dimensional space.
Definition AMReX_RealVect.H:32
This class is a thin wrapper around std::vector. Unlike vector, Vector::operator[] provides bound che...
Definition AMReX_Vector.H:27
void copyAsync(HostToDevice, InIter begin, InIter end, OutIter result) noexcept
A host-to-device copy routine. Note this is just a wrapper around memcpy, so it assumes contiguous st...
Definition AMReX_GpuContainers.H:233
static constexpr HostToDevice hostToDevice
Definition AMReX_GpuContainers.H:98
void streamSynchronize() noexcept
Definition AMReX_GpuDevice.H:237
bool inLaunchRegion() noexcept
Definition AMReX_GpuControl.H:86
bool notInLaunchRegion() noexcept
Definition AMReX_GpuControl.H:87
void Sum(T &v, MPI_Comm comm)
Definition AMReX_ParallelReduce.H:204
void Max(KeyValuePair< K, V > &vi, MPI_Comm comm)
Definition AMReX_ParallelReduce.H:126
MPI_Comm CommunicatorSub() noexcept
sub-communicator for current frame
Definition AMReX_ParallelContext.H:70
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mlmg_lin_cc_interp_r2(Box const &bx, Array4< T > const &ff, Array4< T const > const &cc, int nc) noexcept
Definition AMReX_MLMG_2D_K.H:9
Definition AMReX_Amr.cpp:49
std::enable_if_t< std::is_integral_v< T > > ParallelFor(TypeList< CTOs... > ctos, std::array< int, sizeof...(CTOs)> const &runtime_options, T N, F &&f)
Definition AMReX_CTOParallelForImpl.H:191
void average_down(const MultiFab &S_fine, MultiFab &S_crse, const Geometry &fgeom, const Geometry &cgeom, int scomp, int ncomp, int rr)
Definition AMReX_MultiFabUtil.cpp:314
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mlmg_lin_cc_interp_r2(Box const &bx, Array4< T > const &ff, Array4< T const > const &cc, int nc) noexcept
Definition AMReX_MLMG_1D_K.H:9
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > surroundingNodes(const BoxND< dim > &b, int dir) noexcept
Returns a BoxND with NODE based coordinates in direction dir that encloses BoxND b....
Definition AMReX_Box.H:1399
ReduceData< Ts... >::Type ParReduce(TypeList< Ops... > operation_list, TypeList< Ts... > type_list, FabArray< FAB > const &fa, IntVect const &nghost, F &&f)
Parallel reduce for MultiFab/FabArray.
Definition AMReX_ParReduce.H:47
void EB_set_covered(MultiFab &mf, Real val)
Definition AMReX_EBMultiFabUtil.cpp:21
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mlmg_lin_cc_interp_r4(Box const &bx, Array4< T > const &ff, Array4< T const > const &cc, int nc) noexcept
Definition AMReX_MLMG_1D_K.H:27
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > adjCellLo(const BoxND< dim > &b, int dir, int len=1) noexcept
Returns the cell centered BoxND of length len adjacent to b on the low end along the coordinate direc...
Definition AMReX_Box.H:1591
bool isMFIterSafe(const FabArrayBase &x, const FabArrayBase &y)
Definition AMReX_MFIter.H:219
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > coarsen(const BoxND< dim > &b, int ref_ratio) noexcept
Coarsen BoxND by given (positive) refinement ratio. NOTE: if type(dir) = CELL centered: lo <- lo/rati...
Definition AMReX_Box.H:1304
void EB_average_down(const MultiFab &S_fine, MultiFab &S_crse, const MultiFab &vol_fine, const MultiFab &vfrac_fine, int scomp, int ncomp, const IntVect &ratio)
Definition AMReX_EBMultiFabUtil.cpp:336
IntVectND< AMREX_SPACEDIM > IntVect
Definition AMReX_BaseFwd.H:30
AMREX_GPU_HOST_DEVICE constexpr GpuTupleElement< I, GpuTuple< Ts... > >::type & get(GpuTuple< Ts... > &tup) noexcept
Definition AMReX_Tuple.H:179
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > grow(const BoxND< dim > &b, int i) noexcept
Grow BoxND in all directions by given amount.
Definition AMReX_Box.H:1211
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > adjCell(const BoxND< dim > &b, Orientation face, int len=1) noexcept
Similar to adjCellLo and adjCellHi; operates on given face.
Definition AMReX_Box.H:1634
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void ignore_unused(const Ts &...)
This shuts up the compiler about unused variables.
Definition AMReX.H:111
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE T norm(const GpuComplex< T > &a_z) noexcept
Return the norm (magnitude squared) of a complex number.
Definition AMReX_GpuComplex.H:344
bool TilingIfNotGPU() noexcept
Definition AMReX_MFIter.H:12
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_apply_bc_x(int side, Box const &box, int blen, Array4< T > const &phi, Array4< int const > const &mask, BoundCond bct, T bcl, Array4< T const > const &bcval, int maxorder, T dxinv, int inhomog, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:14
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_comp_interp_coef0_x(int side, Box const &box, int blen, Array4< T > const &f, Array4< int const > const &mask, BoundCond bct, T bcl, int maxorder, T dxinv, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:329
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_apply_bc_z(int side, Box const &box, int blen, Array4< T > const &phi, Array4< int const > const &mask, BoundCond bct, T bcl, Array4< T const > const &bcval, int maxorder, T dzinv, int inhomog, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:224
void Abort(const std::string &msg)
Print out message to cerr and exit via abort().
Definition AMReX.cpp:225
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_comp_interp_coef0_y(int side, Box const &box, int blen, Array4< T > const &f, Array4< int const > const &mask, BoundCond bct, T bcl, int maxorder, T dyinv, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:410
FAB::value_type Dot(FabArray< FAB > const &x, int xcomp, FabArray< FAB > const &y, int ycomp, int ncomp, IntVect const &nghost, bool local=false)
Compute dot products of two FabArrays.
Definition AMReX_FabArrayUtility.H:1554
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > adjCellHi(const BoxND< dim > &b, int dir, int len=1) noexcept
Similar to adjCellLo but builds an adjacent BoxND on the high end.
Definition AMReX_Box.H:1612
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > shift(const BoxND< dim > &b, int dir, int nzones) noexcept
Return a BoxND with indices shifted by nzones in dir direction.
Definition AMReX_Box.H:1372
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_comp_interp_coef0_z(int side, Box const &box, int blen, Array4< T > const &f, Array4< int const > const &mask, BoundCond bct, T bcl, int maxorder, T dzinv, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:491
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mllinop_apply_bc_y(int side, Box const &box, int blen, Array4< T > const &phi, Array4< int const > const &mask, BoundCond bct, T bcl, Array4< T const > const &bcval, int maxorder, T dyinv, int inhomog, int icomp) noexcept
Definition AMReX_MLLinOp_K.H:119
iMultiFab makeFineMask(const BoxArray &cba, const DistributionMapping &cdm, const BoxArray &fba, const IntVect &ratio, int crse_value, int fine_value)
Definition AMReX_MultiFabUtil.cpp:607
std::array< T, N > Array
Definition AMReX_Array.H:24
Definition AMReX_Array4.H:61
Definition AMReX_Dim3.H:12
int x
Definition AMReX_Dim3.H:12
int z
Definition AMReX_Dim3.H:12
int y
Definition AMReX_Dim3.H:12
Definition AMReX_Array.H:34
Definition AMReX_MLLinOp.H:35
bool has_metric_term
Definition AMReX_MLLinOp.H:43
StateMode
Definition AMReX_MLLinOp.H:86
BCMode
Definition AMReX_MLLinOp.H:85
Location
Definition AMReX_MLLinOp.H:87
FabArray memory allocation information.
Definition AMReX_FabArray.H:66
Definition AMReX_MFIter.H:20
MFItInfo & SetDynamic(bool f) noexcept
Definition AMReX_MFIter.H:34
MFItInfo & EnableTiling(const IntVect &ts=FabArrayBase::mfiter_tile_size) noexcept
Definition AMReX_MFIter.H:29
Definition AMReX_MLCellLinOp.H:140
BoundCond type
Definition AMReX_MLCellLinOp.H:141
RT location
Definition AMReX_MLCellLinOp.H:142
Definition AMReX_MLCellLinOp.H:227
T bcloc_hi
Definition AMReX_MLCellLinOp.H:234
Box bx
Definition AMReX_MLCellLinOp.H:235
Array4< T const > bcval_hi
Definition AMReX_MLCellLinOp.H:230
Array4< int const > mask_hi
Definition AMReX_MLCellLinOp.H:232
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Box const & box() const noexcept
Definition AMReX_MLCellLinOp.H:243
int comp
Definition AMReX_MLCellLinOp.H:239
int dir
Definition AMReX_MLCellLinOp.H:240
Array4< T const > bcval_lo
Definition AMReX_MLCellLinOp.H:229
int blen
Definition AMReX_MLCellLinOp.H:238
Array4< T > fab
Definition AMReX_MLCellLinOp.H:228
BoundCond bctype_lo
Definition AMReX_MLCellLinOp.H:236
T bcloc_lo
Definition AMReX_MLCellLinOp.H:233
Array4< int const > mask_lo
Definition AMReX_MLCellLinOp.H:231
BoundCond bctype_hi
Definition AMReX_MLCellLinOp.H:237
Definition AMReX_MLCellLinOp.H:247
int blen
Definition AMReX_MLCellLinOp.H:257
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Box const & box() const noexcept
Definition AMReX_MLCellLinOp.H:262
BoundCond bcthi
Definition AMReX_MLCellLinOp.H:256
int comp
Definition AMReX_MLCellLinOp.H:258
Array4< T > flo
Definition AMReX_MLCellLinOp.H:248
Array4< int const > mlo
Definition AMReX_MLCellLinOp.H:250
int dir
Definition AMReX_MLCellLinOp.H:259
Array4< T > fhi
Definition AMReX_MLCellLinOp.H:249
T bcllo
Definition AMReX_MLCellLinOp.H:252
Box bx
Definition AMReX_MLCellLinOp.H:254
BoundCond bctlo
Definition AMReX_MLCellLinOp.H:255
Array4< int const > mhi
Definition AMReX_MLCellLinOp.H:251
T bclhi
Definition AMReX_MLCellLinOp.H:253
Struct for holding types.
Definition AMReX_TypeList.H:12