AMReX_Cluster.H

Go to the documentation of this file.

 
#ifndef AMREX_Cluster_H_
#define AMREX_Cluster_H_
#include <AMReX_Config.H>
 
#include <AMReX_BoxList.H>
#include <AMReX_REAL.H>
 
#include <list>
 
namespace amrex {
 
class BoxDomain;
class BoxArray;
class ClusterList;
 
 
class Cluster
{
public:
 
    Cluster () noexcept = default;
 
    Cluster (IntVect* a, Long len) noexcept;
 
    Cluster (Cluster& c, const Box& b);
 
    ~Cluster () = default;
    Cluster (const Cluster&) = delete;
    Cluster (Cluster&&) = delete;
    Cluster& operator= (const Cluster&) = delete;
    Cluster& operator= (Cluster&&) = delete;
 
    [[nodiscard]] const Box& box () const noexcept { return m_bx; }
 
    [[nodiscard]] bool ok () const noexcept { return m_bx.d_numPts() > 0. && m_ar != nullptr && m_len > 0; }
 
    [[nodiscard]] Long numTag () const noexcept { return m_len; }
 
    [[nodiscard]] Long numTag (const Box& b) const noexcept;
 
    Cluster* chop ();
 
 
    Cluster* new_chop ();
 
    void distribute (ClusterList&     clst,
                     const BoxDomain& bd);
 
    [[nodiscard]] Real eff () const noexcept {
        BL_ASSERT(ok());
        return static_cast<Real>(double(numTag()) / m_bx.d_numPts());
    }
 
private:
 
    void minBox () noexcept;
 
    Box      m_bx;
    IntVect* m_ar = nullptr;
    Long     m_len = 0;
};
 
 
class ClusterList
{
public:
 
    ClusterList () = default;
 
    ClusterList (IntVect* pts, Long len);
 
    ~ClusterList ();
 
    ClusterList (const ClusterList&) = delete;
    ClusterList (ClusterList&&) = delete;
    ClusterList& operator= (const ClusterList&) = delete;
    ClusterList& operator= (ClusterList&&) = delete;
 
    [[nodiscard]] int length () const { return static_cast<int>(lst.size()); }
 
    void append (Cluster* c) { lst.push_back(c); }
 
    [[nodiscard]] BoxArray boxArray () const;
 
    void boxArray (BoxArray& ba) const;
 
    [[nodiscard]] BoxList boxList() const;
 
    void boxList (BoxList& blst) const;
 
    void chop (Real eff);
 
    void new_chop (Real eff);
 
    void intersect (BoxArray& ba);
 
private:
 
    std::list<Cluster*> lst;
};
 
}
 
#endif /*_Cluster_H_*/