AMReX_MLCellABecLap_2D_K.H

Go to the documentation of this file.

#ifndef AMREX_MLCELLABECLAP_2D_K_H_
#define AMREX_MLCELLABECLAP_2D_K_H_
#include <AMReX_Config.H>
 
namespace amrex {
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
int coarsen_overset_mask (Box const& bx, Array4<int> const& cmsk, Array4<int const> const& fmsk)
    noexcept
{
    int nerrors = 0;
    const auto lo = amrex::lbound(bx);
    const auto hi = amrex::ubound(bx);
    for (int j = lo.y; j <= hi.y; ++j) {
    for (int i = lo.x; i <= hi.x; ++i) {
        int ii = 2*i;
        int jj = 2*j;
        cmsk(i,j,0) = fmsk(ii,jj,0) + fmsk(ii+1,jj,0) + fmsk(ii,jj+1,0) + fmsk(ii+1,jj+1,0);
        if (cmsk(i,j,0) == 4) {
            cmsk(i,j,0) = 1;
        } else if (cmsk(i,j,0) != 0) {
            ++nerrors;
        }
    }}
    return nerrors;
}
 
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
void coarsen_overset_mask (int i, int j, int, Array4<int> const& cmsk,
                           Array4<int const> const& fmsk) noexcept
{
    int ii = 2*i;
    int jj = 2*j;
    cmsk(i,j,0) = fmsk(ii,jj,0) + fmsk(ii+1,jj,0) + fmsk(ii,jj+1,0) + fmsk(ii+1,jj+1,0);
    if (cmsk(i,j,0) == 4) {
        cmsk(i,j,0) = 1;
    }
}
 
}
 
#endif