FTMAtomicUF.h

Go to the documentation of this file.

 
 
#pragma once
 
#include <memory>
#include <vector>
 
#include "FTMDataTypes.h"
#include "FTMStructures.h"
 
namespace ttk {
  namespace ftm {
 
    class AtomicUF {
    private:
      unsigned rank_{};
      AtomicUF *parent_{};
      SharedData data_;
 
    public:
      inline explicit AtomicUF(SimplexId extrema = nullVertex)
        : data_(extrema) {
      }
 
      // heavy recursif
      inline AtomicUF *find() {
        if(parent_ == nullptr)
          return this;
        else {
          decltype(parent_) tmp = parent_->find();
#ifdef TTK_ENABLE_OPENMP
#pragma omp atomic write
#endif
          parent_ = tmp;
 
          return parent_;
        }
      }
 
      // Shared data get/set
 
      inline SimplexId getExtrema() const {
        return data_.extrema;
      }
 
      inline CurrentState *getFirstState() {
#ifndef TTK_ENABLE_KAMIKAZE
        if(data_.states.size() != 1) {
          std::cout << "AtomicUF :: getFirstState : nb state 1 !=  "
                    << data_.states.size() << std::endl;
        }
#endif
        return data_.states[0];
      }
 
      inline CurrentState *getState(const std::size_t id) {
#ifndef TTK_ENABLE_KAMIKAZE
        if(id >= data_.states.size()) {
          std::cout << "AtomicUF :: getState : " << id
                    << " >= " << data_.states.size() << std::endl;
        }
#endif
        return data_.states[id];
      }
 
      inline FTMAtomicVector<CurrentState *> &getStates() {
        return data_.states;
      }
 
      inline size_t getNbStates() const {
        return data_.states.size();
      }
 
      inline FTMAtomicVector<idSuperArc> &getOpenedArcs() {
        return data_.openedArcs;
      }
 
      inline void clearOpenedArcs() {
        data_.openedArcs.reset();
      }
 
      inline void setExtrema(SimplexId v) {
        data_.extrema = v;
      }
 
      inline void reserveData(const size_t &s) {
        data_.reserve(s);
      }
 
      inline void addState(CurrentState *s) {
        data_.addState(s);
      }
 
      inline void addArcToClose(idSuperArc a) {
        data_.addArc(a);
      }
 
      inline CurrentState *mergeStates() {
        CurrentState *s = data_.states[0];
        const auto &nbState = data_.states.size();
 
        for(valence i = 1; i < (valence)nbState; ++i) {
          s->merge(*data_.states[i]);
        }
 
        data_.states.reset(1);
        return s;
      }
 
      // UF get / set
 
      inline int getRank() const {
        return rank_;
      }
 
      inline void setRank(const int &rank) {
#ifdef TTK_ENABLE_OPENMP
#pragma omp atomic write
#endif
        rank_ = rank;
      }
 
      inline void setParent(AtomicUF *parent) {
#ifdef TTK_ENABLE_OPENMP
#pragma omp atomic write
#endif
        parent_ = parent;
      }
 
      static inline AtomicUF *makeUnion(AtomicUF *uf0, AtomicUF *uf1) {
        uf0 = uf0->find();
        uf1 = uf1->find();
 
        if(uf0 == uf1) {
          return uf0;
        } else if(uf0->getRank() > uf1->getRank()) {
          uf1->setParent(uf0);
          uf0->data_.merge(uf1->data_);
          return uf0;
        } else if(uf0->getRank() < uf1->getRank()) {
          uf0->setParent(uf1);
          uf1->data_.merge(uf0->data_);
          return uf1;
        } else {
          uf1->setParent(uf0);
          uf0->setRank(uf0->getRank() + 1);
          uf0->data_.merge(uf1->data_);
          return uf0;
        }
      }
 
      static inline AtomicUF *makeUnion(std::vector<AtomicUF *> &sets) {
        AtomicUF *n = nullptr;
 
        if(!sets.size())
          return nullptr;
 
        if(sets.size() == 1)
          return sets[0];
 
        for(int i = 0; i < (int)sets.size() - 1; i++)
          n = makeUnion(sets[i], sets[i + 1]);
 
        return n;
      }
 
      inline bool operator<(const AtomicUF &other) const {
        return rank_ < other.rank_;
      }
 
      inline bool operator>(const AtomicUF &other) const {
        return rank_ > other.rank_;
      }
    };
 
  } // namespace ftm
} // namespace ttk