Method boundaryDterm

Source/radiation/HypreExtMultiABec.cpp:1020–1071

void HypreExtMultiABec::boundaryDterm(int level,
                                      MultiFab* Dterm,
                                      MultiFab& Soln,
                                      int icomp)
{
  BL_PROFILE("HypreExtMultiABec::boundaryDterm");

  const Box& domain = bd[level]->getDomain();
#ifdef _OPENMP
#pragma omp parallel
#endif
  for (MFIter mfi(Soln); mfi.isValid(); ++mfi) {
    int i = mfi.index();
    const Box &reg = grids[level][i];
    for (OrientationIter oitr; oitr; oitr++) {
      int cdir(oitr());
      int idim = oitr().coordDir();
      const RadBoundCond &bct = bd[level]->bndryConds(oitr())[i];
      const Real      &bcl = bd[level]->bndryLocs(oitr())[i];
      const FArrayBox       &bcv  = bd[level]->bndryValues(oitr())[mfi];
      const Mask      &msk = bd[level]->bndryMasks(oitr(), i);

      if (reg[oitr()] == domain[oitr()]) {
        Array4<int const> tf_arr;
        int bctype = bct;
        if (bd[level]->mixedBndry(oitr())) {
          const BaseFab<int> &tf = *(bd[level]->bndryTypes(oitr())[i]);
          tf_arr = tf.array();
          bctype = -1;
        }
        HABEC::hdterm3(Dterm[idim][mfi].array(),
                       Soln[mfi].array(icomp),
                       reg,
                       cdir, bctype, tf_arr, bcl,
                       bcv.array(bdcomp),
                       msk.array(),
                       (*d2coefs[level])[idim][mfi].array(),
                       geom[level].CellSize());
      }
      else {
          HABEC::hdterm(Dterm[idim][mfi].array(),
                        Soln[mfi].array(icomp),
                        reg,
                        cdir, bct, bcl,
                        bcv.array(bdcomp),
                        msk.array(),
                        (*d2coefs[level])[idim][mfi].array(),
                        geom[level].CellSize());
      }
    }
  }
}