region.hpp

#ifndef ABLATELIBRARY_REGION_HPP
#define ABLATELIBRARY_REGION_HPP
 
#include <petsc.h>
#include <memory>
#include <string>
#include <utilities/petscUtilities.hpp>
#include <vector>
namespace ablate::domain {
 
class Region {
   public:
    inline const static std::shared_ptr<Region> ENTIREDOMAIN = {};
 
   private:
    const std::string name;
    const PetscInt value;
    std::size_t id;
 
   public:
    explicit Region(const std::string& name = {}, int value = 1);
 
    [[nodiscard]] inline const std::size_t& GetId() const { return id; }
 
    [[nodiscard]] inline const std::string& GetName() const { return name; }
 
    [[nodiscard]] inline const PetscInt& GetValue() const { return value; }
 
    [[nodiscard]] inline std::string ToString() const { return name + ":" + std::to_string(value); };
 
    void CreateLabel(DM dm, DMLabel& regionLabel, PetscInt& regionValue) const;
 
    static void GetLabel(const std::shared_ptr<Region>& region, DM dm, DMLabel& regionLabel, PetscInt& regionValue);
 
    void GetLabel(DM dm, DMLabel& regionLabel, PetscInt& regionValue);
 
    static bool InRegion(const std::shared_ptr<Region>& region, DM dm, PetscInt point);
 
    inline bool InRegion(DM dm, PetscInt point) const {
        PetscInt ptValue;
        DMGetLabelValue(dm, name.c_str(), point, &ptValue) >> utilities::PetscUtilities::checkError;
        return ptValue == value;
    }
 
    void CheckForLabel(DM dm) const;
 
    void CheckForLabel(DM dm, MPI_Comm comm) const;
 
    inline bool operator==(const Region& otherRegion) const { return id == otherRegion.id; }
 
    inline bool operator<(const Region& otherRegion) const { return id < otherRegion.id; }
};
 
std::ostream& operator<<(std::ostream& os, const Region& region);
 
std::ostream& operator<<(std::ostream& os, const std::shared_ptr<Region>& region);
 
}  // namespace ablate::domain
#endif  // ABLATELIBRARY_REGION_HPP