AMReX_Mask.H

Go to the documentation of this file.

 
#ifndef AMREX_MASK_H_
#define AMREX_MASK_H_
#include <AMReX_Config.H>
 
#include <AMReX_BaseFab.H>
#include <AMReX_FArrayBox.H>
 
namespace amrex {
 
class Mask final
    :
    public BaseFab<int>
{
public:
 
    Mask () noexcept = default;
 
    explicit Mask (Arena* ar) noexcept;
 
    Mask (const Box& bx, int nc, Arena* ar);
 
    explicit Mask (const Box& bx,
                   int        nc = 1,
                   bool       alloc=true,
                   bool       shared=false,
                   Arena*     ar = nullptr);
 
    explicit Mask (std::istream& is);
 
    explicit Mask (Array4<int> const& a) noexcept : BaseFab<int>(a) {}
 
    explicit Mask (Array4<int> const& a, IndexType t) noexcept : BaseFab<int>(a,t) {}
 
    explicit Mask (Array4<int const> const& a) noexcept : BaseFab<int>(a) {}
 
    explicit Mask (Array4<int const> const& a, IndexType t) noexcept : BaseFab<int>(a,t) {}
 
    ~Mask () noexcept override = default;
 
    Mask (Mask&& rhs) noexcept = default;
 
    Mask (Mask const& rhs, MakeType make_type, int scomp, int ncomp);
 
    Mask (const Mask&) = delete;
    Mask& operator= (const Mask&) = delete;
    Mask& operator= (Mask&&) = delete;
 
 
    friend std::istream& operator>> (std::istream&, Mask&);
 
    void readFrom (std::istream&);
 
    friend std::ostream& operator<< (std::ostream&, const Mask&);
 
    void writeOn (std::ostream&) const;
    //
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& operator&= (const Mask& src) noexcept { return And<run_on>(src); }
 
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& And (const Mask& src) noexcept;
 
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& And (const Mask& src,
               int         srccomp,
               int         destcomp,
               int         numcomp = 1) noexcept;
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& And (const Mask& src,
               const Box&  subbox,
               int         srccomp,
               int         destcomp,
               int         numcomp = 1) noexcept;
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& And (const Mask& src,
               const Box&  srcbox,
               const Box&  destbox,
               int         srccomp,
               int         destcomp,
               int         numcomp = 1) noexcept;
 
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& operator|= (const Mask& src) noexcept { return Or<run_on>(src); }
 
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& Or (const Mask& src) noexcept;
 
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& Or (const Mask& src,
              int         srccomp,
              int         destcomp,
              int         numcomp = 1) noexcept;
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& Or (const Mask& src,
              const Box&  subbox,
              int         srccomp,
              int         destcomp,
              int         numcomp = 1) noexcept;
    template <RunOn run_on AMREX_DEFAULT_RUNON>
    Mask& Or (const Mask& src,
              const Box&  srcbox,
              const Box&  destbox,
              int         srccomp,
              int         destcomp,
              int         numcomp = 1) noexcept;
};
 
template <RunOn run_on>
Mask&
Mask::And (const Mask& src) noexcept
{
    return this->And<run_on>(src,domain,domain,0,0,nvar);
}
 
template <RunOn run_on>
Mask&
Mask::And (const Mask& src,
           int         srccomp,
           int         destcomp,
           int         numcomp) noexcept
{
    return this->And<run_on>(src,domain,domain,srccomp,destcomp,numcomp);
}
 
template <RunOn run_on>
Mask&
Mask::And (const Mask& src,
           const Box&  subbox,
           int         srccomp,
           int         destcomp,
           int         numcomp) noexcept
{
    return this->And<run_on>(src,subbox,subbox,srccomp,destcomp,numcomp);
}
 
template <RunOn run_on>
Mask&
Mask::And (const Mask& src,
           const Box&  srcbox,
           const Box&  destbox,
           int         srccomp,
           int         destcomp,
           int         numcomp) noexcept
{
    auto const& d = this->array();
    auto const& s = src.const_array();
    const auto dlo = amrex::lbound(destbox);
    const auto slo = amrex::lbound(srcbox);
    const Dim3 offset{slo.x-dlo.x,slo.y-dlo.y,slo.z-dlo.z};
    AMREX_HOST_DEVICE_PARALLEL_FOR_4D_FLAG(run_on, destbox, numcomp, i, j, k, n,
    {
        d(i,j,k,n+destcomp) = d(i,j,k,n+destcomp) ? s(i+offset.x,j+offset.y,k+offset.z,n+srccomp) : 0;
    });
    return *this;
}
 
template <RunOn run_on>
Mask&
Mask::Or (const Mask& src) noexcept
{
    return this->Or<run_on>(src,domain,domain,0,0,nvar);
}
 
template <RunOn run_on>
Mask&
Mask::Or (const Mask& src,
          int         srccomp,
          int         destcomp,
          int         numcomp) noexcept
{
    return this->Or<run_on>(src,domain,domain,srccomp,destcomp,numcomp);
}
 
template <RunOn run_on>
Mask&
Mask::Or (const Mask& src,
          const Box&  subbox,
          int         srccomp,
          int         destcomp,
          int         numcomp) noexcept
{
    return this->Or<run_on>(src,subbox,subbox,srccomp,destcomp,numcomp);
}
 
template <RunOn run_on>
Mask&
Mask::Or (const Mask& src,
          const Box&  srcbox,
          const Box&  destbox,
          int         srccomp,
          int         destcomp,
          int         numcomp) noexcept
{
    auto const& d = this->array();
    auto const& s = src.const_array();
    const auto dlo = amrex::lbound(destbox);
    const auto slo = amrex::lbound(srcbox);
    const Dim3 offset{slo.x-dlo.x,slo.y-dlo.y,slo.z-dlo.z};
    AMREX_HOST_DEVICE_PARALLEL_FOR_4D_FLAG(run_on, destbox, numcomp, i, j, k, n,
    {
        d(i,j,k,n+destcomp) = d(i,j,k,n+destcomp) ? 1: s(i+offset.x,j+offset.y,k+offset.z,n+srccomp);
    });
    return *this;
}
 
}
 
#endif /*_MASK_H_*/