cluster/tree03.hh

Go to the documentation of this file.

/* Cluster tree 03
 */

#ifndef clusterTree03_hh
#define clusterTree03_hh

#include "toolbox/hashMap.hh"
#include "cluster/tree.hh"

namespace cluster {

  // ********************************************************** TreeTraits03 **
 
  template<class Node>
  class TreeTraits03 : public TreeTraits<Node> {
  public:
    typedef typename TreeTraits<Node>::F F;

    static Node* newNode(const concepts::Element<F>& elm, const BBall<F>& ball,
       Node* lnk) {return new Node(elm, ball, lnk);}
    static Node* newNode(uint idx, Node& chld, const concepts::Real3d& c,
       concepts::Real r, uint nlf = 0) {
      return new Node(idx, chld, c, r, nlf);
    }

    static Node* child(const Node* nd) {return nd->child();}
    static Node* link(const Node* nd) {return nd->link();}
    static Node*& link(Node* nd) {return nd->link();}
  };
  
  // ************************************************************ CebyCenter **

  template<class Node>
  class CebyCenter {
  public:
    typedef TreeTraits03<Node> Traits;
    typedef typename Traits::F F;

    inline CebyCenter(uint nvtx);
    inline ~CebyCenter() {delete[] vtx_;}
  
    void operator()(const Node* lfset, concepts::Real3d& c, concepts::Real& r,
        const concepts::Real3d* (&e)[4], uint& ne);

  private:
    static concepts::Real eps_;
    uint nvtx_;
    const concepts::Real3d** vtx_;
    __gnu_cxx::hash_map<uint, uint> vtxtag_;
    uint tag_;
  };

  template<class Node>
  inline CebyCenter<Node>::CebyCenter(uint nvtx)
    : nvtx_(nvtx), vtxtag_(), tag_(0) {
    vtx_ = new const concepts::Real3d*[nvtx_];
  }

  // **************************************************************** Tree03 **

  template<class CNode>
  class Tree03 : public Tree<CNode> {
  public:
    typedef TreeTraits03<CNode> Traits;
    typedef typename Traits::F  F;

    Tree03(const concepts::Space<F>& spc, const BBall<F>& ball,
     uint nlfmin = 1);
    inline ~Tree03() {destructor_(root_);  delete root_;}

    inline uint nclst() const {return nclst_;}
    inline uint nleaf() const {return nlf_;}
    inline const CNode* root() const {return root_;}

    void sketch(std::ostream& os, const CNode* lnk = 0, uint l = 0);

  private:
    uint nlfmin_;
    uint nlf_;
    uint nclst_;
    CNode* root_;

    CNode* constructor_(CNode* lfset, CebyCenter<CNode>& cby);
    void destructor_(CNode* clst);
    CNode* concat_(CNode* a, CNode* b);
  };

} // namespace cluster

#endif // clusterTree03_hh